Acredito que o problema está na criação da variável dummy.<div><br></div><div>Será que o
<span style="font-family:'trebuchet ms',sans-serif;font-size:13px;background-color:rgb(255,255,255)">reshape::combine_factor()</span> ajuda?</div><div><br></div><div>Abraços<br><br><div class="gmail_quote">2012/6/25 Sérgio Henrique almeida da silva ju <span dir="ltr"><<a href="mailto:sergio.edfisica@gmail.com" target="_blank">sergio.edfisica@gmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Olá<div><br></div><div>Estou enviando uma dúvida de um amigo:</div><div><br></div><div>Ele está rodando um modelo hierárquico de poisson com variância robusta</div>
<div><br></div><div>structure(list(AFDIC = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, </div>
<div>0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, </div><div>0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, </div><div>1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, </div>
<div>0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, </div><div>0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, </div><div>0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), F6rec = structure(c(1L, </div>
<div>2L, 1L, 5L, 4L, 5L, 5L, 5L, 2L, 5L, 1L, 4L, 5L, 5L, 2L, 2L, 2L, </div><div>2L, 1L, 1L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 2L, </div><div>1L, 1L, 2L, 4L, 2L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 5L, 1L, 1L, 2L, </div>
<div>3L, 2L, 2L, 1L, 2L, 4L, 3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L, 2L, </div><div>2L, 3L, 1L, 5L, 1L, 4L, 2L, 1L, 2L, 2L, 1L, 2L, 4L, 2L, 1L, 3L, </div><div>5L, 2L, 1L, 3L, 2L, 4L, 1L, 2L, 4L, 4L, 2L, 2L, 1L, 3L, 2L, 4L, </div>
<div>2L, 2L), .Label = c("1", "2", "3", "4", "5"), class = "factor")), .Names = c("AFDIC", </div><div>"F6rec"), row.names = c(NA, -99L), class = "data.frame")</div>
<div><br></div><div># Quando ele roda o modelo o resultado do p-valor é para cada categoria da variável, porém ele gostaria também de obter o valor da variável como um todo:
</div><div><br></div><div>banco$F6rec=as.factor(banco$F6rec) <br clear="all"><div>library(epicalc)</div><div><div>library(lmtest)</div><div>library(sandwich)</div></div><div><br></div><div>model2=glm(AFDIC~F6rec,family=poisson,data=banco)</div>
<div><div>coeftest(model2, vcov=sandwich)</div></div><div><br></div><div><div>z test of coefficients:</div><div><br></div><div> Estimate Std. Error z value Pr(>|z|) </div><div>(Intercept) -2.77259 0.68465 -4.0496 5.13e-05 ***</div>
<div>F6rec2 1.49165 0.73551 2.0281 0.04255 * </div><div>F6rec3 1.16315 0.93207 1.2479 0.21206 </div><div>F6rec4 -15.53000 0.74810 -20.7592 < 2.2e-16 ***</div><div>F6rec5 0.47000 1.16994 0.4017 0.68788 </div>
<div>---</div><div>Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 </div></div><div><br></div><div>#Observe que o resultado é para cada categoria de F6rec (F6rec2, F6rec3, F6rec4, F6rec,5), porém ele gostaria de obter também o p-valor para o F6rec como um todo, quando ele rodar o modelo multivariado ele que decidir se mantem essa variável ou não.</div>
<div><br></div><div>Vocês têm alguma ideia?</div><div><br></div><div>Abraços</div><span class="HOEnZb"><font color="#888888"><div><br></div><div><br></div><div><br></div>-- <br>Sérgio Henrique Almeida da Silva Junior<br>
Doutorando em Epidemiologia em Saúde Pública<br>
Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ<br><a href="http://lattes.cnpq.br/1611345552843383" target="_blank">http://lattes.cnpq.br/1611345552843383</a> <br>Tel: (21) 94429486/78101651 id: 123*20942<br>
<br><br>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>Sérgio Henrique Almeida da Silva Junior<br>Doutorando em Epidemiologia em Saúde Pública<br>Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ<br><a href="http://lattes.cnpq.br/1611345552843383">http://lattes.cnpq.br/1611345552843383</a> <br>
Tel: (21) 94429486/78101651 id: 123*20942<br><br><br>
</div>