[R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Luciane Maria Pilotto
lutipilotto em yahoo.com.br
Quinta Dezembro 18 22:52:56 BRST 2014
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 <jalvesaq em gmail.com> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Para: r-br em 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 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.
Mais detalhes sobre a lista de discussão R-br