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

[Jakson Alves de Aquino]

On Fri, Dec 12, 2014 at 06:21:08PM -0800, Luciane Maria Pilotto wrote:
> antes de rodar as regressões preciso fazer as tabelas de
> contingencia utilizando os pesos amostrais com a função
> "crosstab" e não está dando certo. O pvalor está aparecendo como
> 0.

[...]

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