Olá
Estou enviando uma dúvida de um amigo:
Ele está rodando um modelo hierárquico de poisson com variância
robusta
structure(list(AFDIC = c(0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L,
0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L,
1L,
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L,
0L,
0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L), F6rec =
structure(c(1L,
2L, 1L, 5L, 4L, 5L, 5L, 5L, 2L, 5L, 1L, 4L, 5L, 5L, 2L, 2L,
2L,
2L, 1L, 1L, 1L, 1L, 2L, 1L, 4L, 1L, 2L, 1L, 1L, 1L, 1L, 1L,
2L,
1L, 1L, 2L, 4L, 2L, 2L, 2L, 3L, 2L, 3L, 2L, 2L, 5L, 1L, 1L,
2L,
3L, 2L, 2L, 1L, 2L, 4L, 3L, 1L, 3L, 1L, 3L, 1L, 2L, 1L, 1L,
2L,
2L, 3L, 1L, 5L, 1L, 4L, 2L, 1L, 2L, 2L, 1L, 2L, 4L, 2L, 1L,
3L,
5L, 2L, 1L, 3L, 2L, 4L, 1L, 2L, 4L, 4L, 2L, 2L, 1L, 3L, 2L,
4L,
2L, 2L), .Label = c("1", "2", "3", "4", "5"), class = "factor")),
.Names = c("AFDIC",
"F6rec"), row.names = c(NA, -99L), class = "data.frame")
# 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:
banco$F6rec=as.factor(banco$F6rec)
library(epicalc)
library(lmtest)
library(sandwich)
model2=glm(AFDIC~F6rec,family=poisson,data=banco)
coeftest(model2, vcov=sandwich)
z test of coefficients:
Estimate Std. Error
z value Pr(>|z|)
(Intercept) -2.77259 0.68465 -4.0496
5.13e-05 ***
F6rec2 1.49165 0.73551
2.0281 0.04255 *
F6rec3 1.16315 0.93207
1.2479 0.21206
F6rec4 -15.53000 0.74810 -20.7592 <
2.2e-16 ***
F6rec5 0.47000 1.16994
0.4017 0.68788
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’
1
#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.
Vocês têm alguma ideia?
Abraços
--
Sérgio Henrique Almeida da Silva Junior
Doutorando
em Epidemiologia em Saúde Pública
Escola Nacional de Saúde Pública Sérgio
Arouca - ENSP/FIOCRUZ
http://lattes.cnpq.br/1611345552843383
Tel: (21)
94429486/78101651 id:
123*20942