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

Mauro Sznelwar sznelwar em uol.com.br
Sexta Dezembro 19 00:30:53 BRST 2014


Onde tem este arquivo "id3.rda" para rodar?


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
 _______________________________________________
 


---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com



Mais detalhes sobre a lista de discussão R-br