[R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB

[Luciane Maria Pilotto]

Os comandos estão abaixo.

###########################################################
 Acho que o valor é zero
 mesmo. Olhei o código da função
 print.htest() e descobri que o problema é que
 nela é valor é
 impresso por format.pval()
 enquanto a função print.CrossTable() usa
 cat(). Confira:
 
 ##### Início do código
 
 library(foreign)
 library(descr)
 
 ##Banco id3 - banco parcial para testes
 load("id3.rda")
 
 id3$q11 <- as.factor(id3$q11)
 id3$q13 <- as.factor(id3$q13)
 
 # attach() deve ficar depois
 de finalizada a manipulação dos dados
 attach(id3)
 
 table(q11,q13, useNA="always")
 
 # valor p é informado como
 0
 ct <- crosstab(q11, q13, weight =
 bwgr_et, digits = 1, expected = TRUE,
    
            prop.r=TRUE, prop.c=TRUE, chisq =
 TRUE, format = "SPSS",
        
        dnn = c("local
 atend","satisfação"), plot =
 FALSE,data=id3)
 
 # O teste
 de qui-quadrado está guardado no objeto "ct":
 print(ct$CST$p.value, digits = 16)
 print(ct$CST)
 format.pval(ct$CST$p.value)
 
 ##### FIM do código
 
 Vou corrigir a função print.CrossTable. 
 
 -- 
 Jakson
 Alves de Aquino
 www.lepem.ufc.br/aquino.php
 _______________________________________________
 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.