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

Segunda Dezembro 22 21:49:37 BRST 2014

Quais os camandos do Jackson?

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.
_______________________________________________
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