[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