[R-br] Dúvida

[Edson Lira]

Boa tarde R-br lista, estou com um pequeno problema.

Tenho uma lista de exames que associo ao laboratório respectivo, vejam um exemplo abaixo:


                                 Descrição        lab Quantidade
1              17 ALFA HIDROXI PROGESTERONA       <NA>          1
2                               ACIDO URICO Bioquímica          1
3                               ACIDO URICO Bioquímica          1
4                               ACIDO URICO Bioquímica          1
5                               ACIDO URICO Bioquímica          1
6                               ACIDO URICO Bioquímica          4
7                               ACIDO URICO Bioquímica          9
8                               ACIDO URICO Bioquímica          1
9                               ACIDO URICO Bioquímica       1480
10                              ACIDO URICO Bioquímica         37
11                          ACIDO VALPROICO       <NA>          1
12  ALANINA AMINOTRANSFERASE (ALT/GPT) IFCC       <NA>          1
13                                 ALBUMINA Bioquímica        205
14                                 ALBUMINA Bioquímica         17
15                       ALFA-FETO-PROTEINA  HORMÔNIOS          3
16                       ALFA-FETO-PROTEINA  HORMÔNIOS          2
17                                  AMILASE Bioquímica         77
18                                  AMILASE Bioquímica          4
19                          ANDROSTENEDIONA       <NA>          1
20                          ANDROSTENEDIONA       <NA>          1


Para gerar os dados de contagem dos exames por laboratório.



 Estou implementando uma rotina para otimizar o trabalho de consolidação dos dados. O problema é que alguns exames que já estão associados ao respectivo laboratório quando não são realizados no período, eles não  têm registros, mas no script estão relacionados, vejam a rotina abaixo:

library(memisc)

acj$lab<-recode(acj$Lab,
"Bioquímica" <-c("ACIDO URICO","ALBUMINA","BILIRRUBINA DIRETA",
"BILIRRUBINA INDIRETA","BILIRRUBINA TOTAL E FRAÇÕES","BILIRRUBINAS TOTAL",
"COLESTEROL", "COLESTEROL HDL", "COLESTEROL LDL", "COLESTEROL VLDL",
"CREATININA","DHL - DESIDROGENASE LÁCTICA", "FERRITINA", "FERRO SÉRICO",
"FERRO SÉRICO LABTEST", "FOSFATASE ALCALINA", 
"FOSFATASE ÁCIDA FRAÇÃO PROSTÁTICA", "FOSFATASE ÁCIDA PROSTÁTICA (RIE)", "FOSFORO",
"GAMA GT- GAMA-GLUTAMIL-TRANSFERASE", "GLICOSE", "GLICOSE COM FLUORETO", "GLICOSE POS-PRANDIAL",
"GLICOSE-6-FOSFATO DEHIDROGENASE (GGFD)", "GLOBULINA", "HEMOGLOBINA GLICOSILADA",
"LIPIDIOS TOTAIS", "LIPIDOGRAMA", "MAGNESIO", "POTÁSSIO", "PROTEINAS TOTAIS",
"SATURAÇÃO DE TRANSFERRINA","SÓDIO","TRANSAMINASE OXALACÉTICA - AST", "TRANSAMINASE PIRÚVICA - ALT",
"TRANSAMINASES", "TRANSAMINASES", "TRANSFERRINA", "ALFA-1 GLICOPROTEINA ÁCIDA",
"ALT-TRANSAMINASE PIRÚVICA", "AMILASE", "AST-TRANSSAMINASE OXALACÉTICA",
"CÁLCIO","CAPACIDADE DE OXIDAÇÃO DO FERRO","COLESTEROL TOTAL","COLESTEROL TOTAL E FRAÇÕES",
 "ESTUDO DO FERRO", "FOSFATASE ÁCIDA TOTAL", "TESTE DE ROG", "TRIGLICERÍDEOS", 
"URÉIA", "LÍTIO"),

"Citoquímica" <-c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM"))

Quando rodo esta rotina o R mostra a mensagem:

Erro em factor(y, levels = 1:max(y), labels = newcodes) : 
  invalid labels; length 13 should be 1 or 12
Além disso: Mensagens de aviso perdidas:
In recode(acj$Lab, "Bioquímica" <- c("ACIDO URICO", "ALBUMINA",  :
  recoding "Citoquímica" <- c("FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM") has no consequences
 
O problema é que não existe nenhum exame associado ao laboratório de Citoquímica (FOSFATASE ALCALINA DE LEUCOCITOS", "IMUNOFENOTIPAGEM). ou seja não foi realizado nenhum exame, portanto, não aparecem no banco, como resolver este problema? 
 


[]'s


Edson Lira
Estatístico
Manaus-Amazonas
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110712/a416be8f/attachment.html>