Obrigada pela ajuda! Com estes comandos é possível verificar o pvalue, pois o valor não é zero.
O estranho é que o pvalor só aparece como zero ao usar os pesos amostrais na função crosstab! # com uso dos pesos amostrais > crosstab(q11[getario == "35 a 44 anos"], q13[getario == "35 a 44 anos"], > weight = bwgr_et[getario == "35 a 44 anos"], digits = 1, expected = TRUE, > prop.r=TRUE, prop.c=TRUE, chisq = TRUE, + format = "SPSS", dnn = c("local atend","satisfação"), plot = FALSE,data=id) Conteúdo das células |-------------------------| | Contagem | | Valores esperados | | Percentual por linha | | Percentual por coluna | |-------------------------| ================================================================================ satisfação local atend 1 2 3 4 5 Total -------------------------------------------------------------------------------- 1 1109669 2443218 403767 181185 140021 4277860 1196118.5 2439668.8 380279.6 142113.3 119679.7 25.9 57.1 9.4 4.2 3.3 37.7 35.0 37.8 40.1 48.1 44.1 -------------------------------------------------------------------------------- 2 1573195 3265609 498158 149920 157032 5643914 1578076.5 3218731.1 501714.7 187494.6 157897.2 27.9 57.9 8.8 2.7 2.8 49.8 49.6 50.5 49.4 39.8 49.5 -------------------------------------------------------------------------------- 3 471510 704133 97758 44018 18497 1335916 373531.1 761874.5 118756.0 44380.0 37374.3 35.3 52.7 7.3 3.3 1.4 11.8 14.9 10.9 9.7 11.7 5.8 -------------------------------------------------------------------------------- 4 15438 52362 8089 1489 1611 78989 22085.9 45047.5 7021.7 2624.1 2209.8 19.5 66.3 10.2 1.9 2.0 0.7 0.5 0.8 0.8 0.4 0.5 -------------------------------------------------------------------------------- Total 3169812 6465322 1007772 376612 317161 11336679 28.0 57.0 8.9 3.3 2.8 ================================================================================ Estatísticas para todos os fatores da tabela Pearson's Chi-squared test ------------------------------------------------------------ Qui² = 77492.28 g.l. = 12 p = 0 Frequência esperada mínima: 2209.839 #sem uso do peso amostral > crosstab(q11, q13, digits = 1, expected = TRUE, prop.r=TRUE, prop.c=TRUE, > chisq = TRUE, + format = "SPSS", dnn = c("local atend","satisfação"),plot = FALSE,data=id) Aviso em chisq.test(t, correct = FALSE, ...) : Chi-squared approximation may be incorrect Conteúdo das células |-------------------------| | Contagem | | Valores esperados | | Percentual por linha | | Percentual por coluna | |-------------------------| ============================================================== satisfação local atend 1 2 3 4 5 Total -------------------------------------------------------------- 1 1616 4740 952 315 137 7760 1964.1 4632.3 805.6 242.3 115.7 20.8 61.1 12.3 4.1 1.8 39.3 32.4 40.3 46.5 51.1 46.6 -------------------------------------------------------------- 2 2589 5459 829 235 112 9224 2334.6 5506.3 957.6 288.0 137.5 28.1 59.2 9.0 2.5 1.2 46.8 51.9 46.4 40.5 38.1 38.1 -------------------------------------------------------------- 3 714 1383 227 57 35 2416 611.5 1442.2 250.8 75.4 36.0 29.6 57.2 9.4 2.4 1.4 12.2 14.3 11.7 11.1 9.3 11.9 -------------------------------------------------------------- 4 74 194 40 9 10 327 82.8 195.2 33.9 10.2 4.9 22.6 59.3 12.2 2.8 3.1 1.7 1.5 1.6 2.0 1.5 3.4 -------------------------------------------------------------- Total 4993 11776 2048 616 294 19727 25.3 59.7 10.4 3.1 1.5 ============================================================== Estatísticas para todos os fatores da tabela Pearson's Chi-squared test ------------------------------------------------------------ Qui² = 210.3717 g.l. = 12 p = 2.343732e-38 Frequência esperada mínima: 4.873422 Células com frequências esperada < 5: 1 de 20 (5%) Atenciosamente, Luciane Pilotto -------------------------------------------- Em sáb, 13/12/14, Jakson Alves de Aquino <jalve...@gmail.com> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@listas.c3sl.ufpr.br Data: Sábado, 13 de Dezembro de 2014, 12:51 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 _______________________________________________ R-br mailing list R-br@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. _______________________________________________ R-br mailing list R-br@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.