[R-br] RES: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Luciane Maria Pilotto
lutipilotto em yahoo.com.br
Segunda Dezembro 22 17:07:08 BRST 2014
Sim, está rodando, mas o valor do p aparece como zero. Utilizando os comandos que o Jakson postou dá certo.
__________________________________________________
Luciane Maria Pilotto
--------------------------------------------
Em seg, 22/12/14, Mauro Sznelwar <sznelwar em uol.com.br> escreveu:
Assunto: [R-br] RES: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Para: r-br em listas.c3sl.ufpr.br
Data: Segunda-feira, 22 de Dezembro de 2014, 2:35
Por que eu não consigo
rodar isto, o seguinte dá certo!
>
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)
Error in chisq.test(t,
correct = FALSE, ...) :
at least one
entry of 'x' must be positive
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.
_______________________________________________
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.
---
Este email foi escaneado pelo Avast
antivírus.
http://www.avast.com
_______________________________________________
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