<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content="text/html; charset=windows-1252" http-equiv=Content-Type>
<META name=GENERATOR content="MSHTML 8.00.6001.19258">
<STYLE></STYLE>
</HEAD>
<BODY bgColor=#ffffff>
<DIV><FONT color=#0000ff face=Arial>Tentei rodar o script, mas a partir daí não
roda mais:</FONT></DIV>
<DIV><FONT size=2 face=Arial></FONT> </DIV>
<DIV><FONT size=2 face=Arial>banco$F6rec=as.factor(banco$F6rec) <BR>Erro em
is.factor(x) : objeto 'banco' não encontrado<BR></DIV></FONT>
<BLOCKQUOTE
style="BORDER-LEFT: #000000 2px solid; PADDING-LEFT: 5px; PADDING-RIGHT: 0px; MARGIN-LEFT: 5px; MARGIN-RIGHT: 0px">
<DIV><FONT size=2 face=Arial></FONT><BR> </DIV>
<DIV>Será que o <SPAN
style="BACKGROUND-COLOR: rgb(255,255,255); FONT-FAMILY: 'trebuchet ms',sans-serif; FONT-SIZE: 13px">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
style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex"
class=gmail_quote>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>
<P>
<HR>
<P></P>_______________________________________________<BR>R-br mailing
list<BR>R-br@listas.c3sl.ufpr.br<BR>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<BR>Leia
o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo
reproduzível.</BLOCKQUOTE></BODY></HTML>