[R-br] Razão de prevalência

Terça Abril 19 22:00:15 BRT 2016

Prezados...

Executei o comando sugerido pelo Mauricio

*mod1<-glm(PoliFar~ClassEco+Educ,family="poisson"(link="log"),data=dmha2)*

 e estou encontrando a seguinte mensagem de erro

*Error in if (any(y < 0)) stop("negative values not allowed for the
'Poisson' family") : valor ausente onde TRUE/FALSE necessárioAlém disso:
Warning message:In Ops.factor(y, 0) : ‘<’ not meaningful for factors*

Meus dados estão no seguinte formato: (*data frame dmha2*):

*ClassEco     Educ     PoliFar*

*A                 0-2         S*

*D-E              3-5         NB                 12+         S*

*A                 6-8          S*

*C                 9-11         N*

*D-E              0-2          S*
*...                 ...          ...*

Vocês sabem o que está errado?

Att.,

*Emerson*

Em 19 de abril de 2016 19:03, Emerson Cotta Bodevan <bodevan.ec em gmail.com>
escreveu:

> Obrigado Mauricio.
>
> Vou executar aqui.
>
> Att.,
>
>
>
> *Emerson*
>
> Em 19 de abril de 2016 19:01, Mauricio Cardeal <mcardeal2010 em gmail.com>
> escreveu:
>
>> Emerson, você pode tentar assim:
>>
>> modelo=glm(variavel dependente ~ variaveis
>> independentes,family="poisson"(link="log"))
>> exp(modelo$coefficients)
>> summary(modelo)
>>
>> Para os intervalos de confiança use a função confint:
>>
>> exp(confint(modelo)) # 95% CI for exponentiated coefficients
>>
>> Mauricio Cardeal
>> UFBA
>>
>>
>> Em 19/04/2016 18:37, Emerson Cotta Bodevan escreveu:
>>
>> Olá Marco!
>>
>> Obrigado pela dica. Depois de ler as referências e sugestões, vi que
>> (apesar de minha resposta ser dicotômica) , como meu estudo é transversal,
>> o melhor é usar a Regressão de Poisson com variância Robusta.
>>
>> A dica do Maurício Cardeal. Agora... como construo o IC da Razão de
>> Prevalência?
>>
>> Agradeço a todos.
>>
>>
>> *Emerson*
>>
>> Em 19 de abril de 2016 17:22, Marco Nunes <nunes.ma em outlook.com>
>> escreveu:
>>
>>> Olá Emerson
>>>
>>> O uso da razão de prevalência na regressão logística já foi reapondido.
>>> Para o cálculo de razão de prevalência (assim como risco relatico e
>>> razão de odds também) utilizo o pacote "epiR".
>>> Encaminho um tutorial abaixo.
>>>
>>>
>>> # Exemplo de utilização do epi.2by2 para calculo da razão de prevalencia
>>>
>>> library(epiR)
>>>
>>> dat <- matrix(c(13,2163,5,3349), nrow = 2, byrow = TRUE)
>>> rownames(dat) <- c("DF+", "DF-")
>>> colnames(dat) <- c("FUS+", "FUS-")
>>> dat
>>> epi.2by2(dat = as.table(dat), method = "cross.sectional", conf.level =
>>> 0.95, units = 100,  homogeneity = "breslow.day", outcome = "as.columns")
>>>
>>>
>>> Prof. Dr. Marco Antonio Prado Nunes
>>> Departamento de Medicina/UFS
>>>
>>> ------------------------------
>>> Date: Sat, 9 Apr 2016 15:49:02 -0300
>>> From: bodevan.ec em gmail.com
>>> To: r-br em listas.c3sl.ufpr.br
>>> Subject: Re: [R-br] Razão de prevalência
>>>
>>>
>>> Obrigado César.
>>>
>>> Vou verificar.
>>>
>>> Abraço,
>>>
>>>
>>> *Emerson*
>>>
>>> Em 8 de abril de 2016 14:30, Cesar Rabak <cesar.rabak em gmail.com>
>>> escreveu:
>>>
>>> Olá Emerson,
>>>
>>> Embora a referência apresentada pelo Leonardo seja recente e
>>> interessante, pode acontecer que devido a exigências do veículo que você
>>> publicará, orientação ou chefia, o cálculo tenha que seguir com a RL
>>> regular para os IC da RP.
>>>
>>> Nesse caso, recomendo a leitura de um artigo mais "clássico" e que tem
>>> boas referências e código e R:
>>> <http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf>
>>> http://www.scielo.br/pdf/csp/v31n3/0102-311X-csp-31-03-00487.pdf
>>>
>>> HTH
>>> --
>>> Cesar Rabak
>>>
>>> 2016-04-07 15:21 GMT-03:00 Emerson Cotta Bodevan <
>>> <bodevan.ec em gmail.com>bodevan.ec em gmail.com>:
>>>
>>> Prezados, boa tarde.
>>>
>>> Estou usando regressão logística, tendo como resposta a ocorrência ou
>>> não de uma doença comum (alta prevalência).
>>>
>>> Como calcular a *razão de prevalência* e seus intervalos de confiança?
>>> Estudo transversal.
>>>
>>> Algum pacote adequado?
>>>
>>> Agradeço qualquer ajuda.
>>>
>>> *Emerson*
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem ( <http://www.leg.ufpr.br/r-br-guia>
>>> http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>>>
>>>
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem ( <http://www.leg.ufpr.br/r-br-guia>
>>> http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
>>>
>>>
>>>
>>> _______________________________________________ R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de
>>> postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo
>>> reproduz�vel.
>>>
>>> _______________________________________________
>>> R-br mailing list
>>> R-br em listas.c3sl.ufpr.br
>>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>>> código mínimo reproduzível.
>>>
>>
>>
>>
>> _______________________________________________
>> R-br mailing listR-br em listas.c3sl.ufpr.brhttps://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forne�a c�digo m�nimo reproduz�vel.
>>
>>
>>
>> _______________________________________________
>> R-br mailing list
>> R-br em listas.c3sl.ufpr.br
>> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
>> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
>> código mínimo reproduzível.
>>
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160419/12201b42/attachment.html>