[R-br] Ajuda modelo hierárquico de poisson com variância robusta

[Mauro Sznelwar]

Tentei rodar o script, mas a partir daí não roda mais:

banco$F6rec=as.factor(banco$F6rec) 
Erro em is.factor(x) : objeto 'banco' não encontrado


   
  Será que o  reshape::combine_factor() ajuda?


  Abraços


  2012/6/25 Sérgio Henrique almeida da silva ju <sergio.edfisica em gmail.com>

    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








  -- 
  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





------------------------------------------------------------------------------


  _______________________________________________
  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/20120626/b016fd07/attachment-0001.html>