Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
Não consegui usar o comando svyolr, está dando os erros abaixo.
# Utilizando as variáveis criadas q11 e q13 que foram transformadas em fator:
> id$q11 <- as.factor(id$q11)
> id$q13 <- as.factor(id$q13)
> dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL, fpc=NULL,
+ data = NULL, nest = FALSE, weights=~peso,pps=FALSE)
> dstrat
Independent Sampling design (with replacement)
svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL,
fpc = NULL, data = NULL, nest = FALSE, weights = ~peso, pps = FALSE)
> svyolr(q13~q11,design=dstrat, method="logistic")
Erro em svyolr.survey.design2(q13 ~ q11, design = dstrat, method = "logistic") :
response must be a factor
# Utilizando as variáveis quest_11 e quest_13 que são fatores no banco:
> svyolr(quest_13~quest_11,design=dstrat, method="logistic")
Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { :
arguments must have same length
==========================================================
Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot = FALSE)).
Dicas com outros comandos que utilizam os pesos amostrais são bem vindas!!!
__________________________________________________
Luciane Maria Pilotto
Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS
--------------------------------------------
Em ter, 25/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
Para: r-br@listas.c3sl.ufpr.br
Data: Terça-feira, 25 de Novembro de 2014, 17:09
Não serve a http://r-survey.r-forge.r-project.org/survey/html/svyolr.html?
Leonardo Ferreira
Fontenelle
Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br
escreveu:
Olá,
estou tentando rodar regressão logística ordinal e
não estou conseguindo usar os pesos amostrais, está dando
o seguinte erro:
> m11 <- polr(q13 ~ q11 ,
method="logistic", weights = peso, data=id,
Hess=TRUE)
Aviso: glm.fit: fitted probabilities numerically 0 or 1
occurred
Erro em optim(s0, fmin, gmin, method =
"BFGS", ...) :
valor inicial em vmmin não é finito
*A variável q13 possui 5 ctegorias (1,2,3,4,5) e a
variável q11 possui 4 categorias (1,2,3,4).
Tem como solicitar o qui-quadrado, o número de
observações e o pseudo R2 no mesmo comando utilizado para
reg. log. ord.?
Como utilizar o peso amostral nas tabelas de
contingência? Estou usando os seguintes comandos:
table(q11,q13)
rowPercents(table(q11,q13))
ou
CrossTable(q11, q13, digits=2,
expected=TRUE, prop.r=TRUE, prop.c=TRUE,
prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE,
format=c("SPSS"), dnn = c("local
atend","satisfação"))
Atenciosamente,
Luciane M. Pilotto
Doutoranda em Saúde Bucal Coletiva/UFRGS
Enviado do Email do Windows
_______________________________________________
R-br mailing list
R-br@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.
-----Anexo incorporado-----
_______________________________________________
R-br mailing list
R-br@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.
No comando abaixo:
dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL,
fpc=NULL, data = NULL, nest = FALSE, weights=~peso,pps=FALSE)
Você provavelmente quer passar o argumento "data = id", em vez de "data
= NULL".
--
Leonardo Ferreira Fontenelle
http://lattes.cnpq.br/9234772336296638
Em Qui 27 nov. 2014, às 14:58, Luciane Maria Pilotto escreveu:
> Não consegui usar o comando svyolr, está dando os erros abaixo.
>
> # Utilizando as variáveis criadas q11 e q13 que foram transformadas em
> fator:
>
> > id$q11 <- as.factor(id$q11)
> > id$q13 <- as.factor(id$q13)
>
>
> > dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL, fpc=NULL,
> + data = NULL, nest = FALSE, weights=~peso,pps=FALSE)
> > dstrat
> Independent Sampling design (with replacement)
> svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL,
> fpc = NULL, data = NULL, nest = FALSE, weights = ~peso, pps = FALSE)
>
> > svyolr(q13~q11,design=dstrat, method="logistic")
> Erro em svyolr.survey.design2(q13 ~ q11, design = dstrat, method =
> "logistic") :
> response must be a factor
>
>
> # Utilizando as variáveis quest_11 e quest_13 que são fatores no banco:
> > svyolr(quest_13~quest_11,design=dstrat, method="logistic")
> Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { :
> arguments must have same length
>
> ==========================================================
>
> Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no
> entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab
> (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot = FALSE)).
>
> Dicas com outros comandos que utilizam os pesos amostrais são bem
> vindas!!!
>
> __________________________________________________
> Luciane Maria Pilotto
> Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS
>
>
> --------------------------------------------
> Em ter, 25/11/14, Leonardo Ferreira Fontenelle
> <leonardof@leonardof.med.br> escreveu:
>
> Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
> Para: r-br@listas.c3sl.ufpr.br
> Data: Terça-feira, 25 de Novembro de 2014, 17:09
>
>
>
>
>
>
> Não serve a
> http://r-survey.r-forge.r-project.org/survey/html/svyolr.html?
>
>
> Leonardo Ferreira
> Fontenelle
>
>
>
>
> Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br
> escreveu:
>
> Olá,
>
>
> estou tentando rodar regressão logística ordinal e
> não estou conseguindo usar os pesos amostrais, está dando
> o seguinte erro:
>
>
> > m11 <- polr(q13 ~ q11 ,
> method="logistic", weights = peso, data=id,
> Hess=TRUE)
>
> Aviso: glm.fit: fitted probabilities numerically 0 or 1
> occurred
>
> Erro em optim(s0, fmin, gmin, method =
> "BFGS", ...) :
>
> valor inicial em vmmin não é finito
>
>
> *A variável q13 possui 5 ctegorias (1,2,3,4,5) e a
> variável q11 possui 4 categorias (1,2,3,4).
>
>
> Tem como solicitar o qui-quadrado, o número de
> observações e o pseudo R2 no mesmo comando utilizado para
> reg. log. ord.?
>
>
> Como utilizar o peso amostral nas tabelas de
> contingência? Estou usando os seguintes comandos:
>
> table(q11,q13)
>
> rowPercents(table(q11,q13))
>
> ou
>
> CrossTable(q11, q13, digits=2,
> expected=TRUE, prop.r=TRUE, prop.c=TRUE,
>
>
> prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE,
>
>
> format=c("SPSS"), dnn = c("local
> atend","satisfação"))
>
>
>
> Atenciosamente,
> Luciane M. Pilotto
> Doutoranda em Saúde Bucal Coletiva/UFRGS
>
> Enviado do Email do Windows
>
>
>
>
>
> _______________________________________________
>
> R-br mailing list
>
> R-br@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.
>
>
>
>
> -----Anexo incorporado-----
>
> _______________________________________________
> R-br mailing list
> R-br@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@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.
Ah, não Leonardo, já havia alterado e dá o mesmo erro!!
Enviado do Windows
--------------------------------------------
Em qui, 27/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
Para: r-br@listas.c3sl.ufpr.br
Data: Quinta-feira, 27 de Novembro de 2014, 21:11
No comando abaixo:
dstrat<-svydesign(ids=~1,
probs=NULL, strata = NULL, variable = NULL,
fpc=NULL, data = NULL, nest = FALSE,
weights=~peso,pps=FALSE)
Você provavelmente quer passar o argumento
"data = id", em vez de "data
= NULL".
--
Leonardo Ferreira
Fontenelle
http://lattes.cnpq.br/9234772336296638
Em Qui 27 nov. 2014, às
14:58, Luciane Maria Pilotto escreveu:
>
Não consegui usar o comando svyolr, está dando os erros
abaixo.
>
> #
Utilizando as variáveis criadas q11 e q13 que foram
transformadas em
> fator:
>
> > id$q11 <-
as.factor(id$q11)
> > id$q13 <-
as.factor(id$q13)
>
>
> > dstrat<-svydesign(ids=~1,
probs=NULL, strata = NULL, variable = NULL, fpc=NULL,
> + data = NULL, nest = FALSE,
weights=~peso,pps=FALSE)
> > dstrat
> Independent Sampling design (with
replacement)
> svydesign(ids = ~1, probs
= NULL, strata = NULL, variable = NULL,
> fpc = NULL, data = NULL,
nest = FALSE, weights = ~peso, pps = FALSE)
>
> >
svyolr(q13~q11,design=dstrat,
method="logistic")
> Erro em
svyolr.survey.design2(q13 ~ q11, design = dstrat, method
=
> "logistic") :
> response must be a factor
>
>
> # Utilizando as variáveis quest_11 e
quest_13 que são fatores no banco:
>
> svyolr(quest_13~quest_11,design=dstrat,
method="logistic")
> Erro em
tapply(1:NROW(x), list(factor(strata)), function(index) { :
> arguments must have same
length
>
>
==========================================================
>
> Eu tb poderia usar
os pesos amostrais nas tabelas de contingencia, no
> entanto, tb não consigo usar os pesos
nestas. Utilizei o comando crosstab
>
(crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot
= FALSE)).
>
> Dicas
com outros comandos que utilizam os pesos amostrais são
bem
> vindas!!!
>
>
__________________________________________________
> Luciane Maria Pilotto
>
Mestre e Doutoranda em Saúde Bucal Coletiva -
FO/UFRGS
>
>
>
--------------------------------------------
> Em ter, 25/11/14, Leonardo Ferreira
Fontenelle
> <leonardof@leonardof.med.br>
escreveu:
>
>
Assunto: Re: [R-br] Regressão logística ordinal e uso de
pesos amostrais
> Para: r-br@listas.c3sl.ufpr.br
> Data: Terça-feira, 25 de Novembro de
2014, 17:09
>
>
>
>
>
>
> Não serve a
>
http://r-survey.r-forge.r-project.org/survey/html/svyolr.html?
>
>
> Leonardo Ferreira
>
Fontenelle
>
>
>
>
> Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br
> escreveu:
>
> Olá,
>
>
> estou tentando
rodar regressão logística ordinal e
>
não estou conseguindo usar os pesos amostrais, está
dando
> o seguinte erro:
>
>
> > m11 <- polr(q13 ~ q11 ,
> method="logistic", weights =
peso, data=id,
> Hess=TRUE)
>
> Aviso: glm.fit:
fitted probabilities numerically 0 or 1
> occurred
>
> Erro em optim(s0, fmin, gmin, method
=
> "BFGS", ...) :
>
> valor inicial
em vmmin não é finito
>
>
> *A variável
q13 possui 5 ctegorias (1,2,3,4,5) e a
> variável q11 possui 4 categorias
(1,2,3,4).
>
>
> Tem como solicitar o qui-quadrado,
o número de
> observações e o pseudo
R2 no mesmo comando utilizado para
>
reg. log. ord.?
>
>
> Como utilizar
o peso amostral nas tabelas de
>
contingência? Estou usando os seguintes comandos:
>
>
table(q11,q13)
>
> rowPercents(table(q11,q13))
>
> ou
>
> CrossTable(q11,
q13, digits=2,
> expected=TRUE,
prop.r=TRUE, prop.c=TRUE,
>
>
> prop.t=FALSE,prop.chisq=FALSE, chisq =
FALSE,
>
>
>
format=c("SPSS"), dnn = c("local
>
atend","satisfação"))
>
>
>
>
Atenciosamente,
>
Luciane M. Pilotto
> Doutoranda em
Saúde Bucal Coletiva/UFRGS
>
> Enviado do Email do Windows
>
>
>
>
>
>
_______________________________________________
>
> R-br mailing
list
>
> R-br@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.
>
>
>
>
> -----Anexo
incorporado-----
>
>
_______________________________________________
> R-br mailing list
>
R-br@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@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@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.
Se você quer usar pesos de amostragem, o pacote survey é a única forma
que conheço no R. Se você quer usar as variáveis do data frame "id",
precisa passá-lo como argumento para a função svydesign. Por favor,
forneça um código mínimo reprodutível, incluindo acesso a dados e os
dois comandos (svyolr e svydesign com dados).
--
Leonardo Ferreira Fontenelle
http://lattes.cnpq.br/9234772336296638
Em Sex 28 nov. 2014, às 16:27, Luciane Maria Pilotto escreveu:
> Ah, não Leonardo, já havia alterado e dá o mesmo erro!!
>
>
>
>
> Enviado do Windows
>
> --------------------------------------------
> Em qui, 27/11/14, Leonardo Ferreira Fontenelle
> <leonardof@leonardof.med.br> escreveu:
>
> Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
> Para: r-br@listas.c3sl.ufpr.br
> Data: Quinta-feira, 27 de Novembro de 2014, 21:11
>
> No comando abaixo:
>
> dstrat<-svydesign(ids=~1,
> probs=NULL, strata = NULL, variable = NULL,
> fpc=NULL, data = NULL, nest = FALSE,
> weights=~peso,pps=FALSE)
>
> Você provavelmente quer passar o argumento
> "data = id", em vez de "data
> = NULL".
>
> --
> Leonardo Ferreira
> Fontenelle
> http://lattes.cnpq.br/9234772336296638
>
> Em Qui 27 nov. 2014, às
> 14:58, Luciane Maria Pilotto escreveu:
> >
> Não consegui usar o comando svyolr, está dando os erros
> abaixo.
> >
> > #
> Utilizando as variáveis criadas q11 e q13 que foram
> transformadas em
> > fator:
> >
> > > id$q11 <-
> as.factor(id$q11)
> > > id$q13 <-
> as.factor(id$q13)
> >
> >
>
> > > dstrat<-svydesign(ids=~1,
> probs=NULL, strata = NULL, variable = NULL, fpc=NULL,
> > + data = NULL, nest = FALSE,
> weights=~peso,pps=FALSE)
> > > dstrat
> > Independent Sampling design (with
> replacement)
> > svydesign(ids = ~1, probs
> = NULL, strata = NULL, variable = NULL,
> > fpc = NULL, data = NULL,
> nest = FALSE, weights = ~peso, pps = FALSE)
> >
> > >
> svyolr(q13~q11,design=dstrat,
> method="logistic")
> > Erro em
> svyolr.survey.design2(q13 ~ q11, design = dstrat, method
> =
> > "logistic") :
> > response must be a factor
> >
> >
> > # Utilizando as variáveis quest_11 e
> quest_13 que são fatores no banco:
> >
> > svyolr(quest_13~quest_11,design=dstrat,
> method="logistic")
> > Erro em
> tapply(1:NROW(x), list(factor(strata)), function(index) { :
>
> > arguments must have same
> length
> >
> >
> ==========================================================
> >
> > Eu tb poderia usar
> os pesos amostrais nas tabelas de contingencia, no
> > entanto, tb não consigo usar os pesos
> nestas. Utilizei o comando crosstab
> >
> (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot
> = FALSE)).
> >
> > Dicas
> com outros comandos que utilizam os pesos amostrais são
> bem
> > vindas!!!
> >
> >
> __________________________________________________
> > Luciane Maria Pilotto
> >
> Mestre e Doutoranda em Saúde Bucal Coletiva -
> FO/UFRGS
> >
> >
> >
> --------------------------------------------
> > Em ter, 25/11/14, Leonardo Ferreira
> Fontenelle
> > <leonardof@leonardof.med.br>
> escreveu:
> >
> >
> Assunto: Re: [R-br] Regressão logística ordinal e uso de
> pesos amostrais
> > Para: r-br@listas.c3sl.ufpr.br
> > Data: Terça-feira, 25 de Novembro de
> 2014, 17:09
> >
> >
>
> >
> >
> >
> >
> > Não serve a
> >
> http://r-survey.r-forge.r-project.org/survey/html/svyolr.html?
> >
> >
> > Leonardo Ferreira
> >
> Fontenelle
> >
> >
>
> >
> >
> > Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br
> > escreveu:
> >
> > Olá,
> >
> >
> > estou tentando
> rodar regressão logística ordinal e
> >
> não estou conseguindo usar os pesos amostrais, está
> dando
> > o seguinte erro:
> >
> >
> > > m11 <- polr(q13 ~ q11 ,
> > method="logistic", weights =
> peso, data=id,
> > Hess=TRUE)
> >
> > Aviso: glm.fit:
> fitted probabilities numerically 0 or 1
> > occurred
> >
> > Erro em optim(s0, fmin, gmin, method
> =
> > "BFGS", ...) :
> >
> > valor inicial
> em vmmin não é finito
> >
> >
> > *A variável
> q13 possui 5 ctegorias (1,2,3,4,5) e a
> > variável q11 possui 4 categorias
> (1,2,3,4).
> >
> >
>
> > Tem como solicitar o qui-quadrado,
> o número de
> > observações e o pseudo
> R2 no mesmo comando utilizado para
> >
> reg. log. ord.?
> >
> >
> > Como utilizar
> o peso amostral nas tabelas de
> >
> contingência? Estou usando os seguintes comandos:
> >
> >
> table(q11,q13)
> >
> > rowPercents(table(q11,q13))
> >
> > ou
> >
> > CrossTable(q11,
> q13, digits=2,
> > expected=TRUE,
> prop.r=TRUE, prop.c=TRUE,
> >
> >
> > prop.t=FALSE,prop.chisq=FALSE, chisq =
> FALSE,
> >
> >
>
> >
> format=c("SPSS"), dnn = c("local
> >
> atend","satisfação"))
> >
> >
> >
> >
> Atenciosamente,
> >
> Luciane M. Pilotto
> > Doutoranda em
> Saúde Bucal Coletiva/UFRGS
> >
> > Enviado do Email do Windows
> >
> >
> >
> >
> >
> >
> _______________________________________________
> >
> > R-br mailing
> list
> >
> > R-br@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.
> >
> >
> >
> >
> > -----Anexo
> incorporado-----
> >
> >
> _______________________________________________
> > R-br mailing list
> >
> R-br@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@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@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@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.
Boa noite, 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. A função polr tb permite usar os pesos amostrais, porém, não consigo identificar o erro que está acusando. A função svyolr pode ser uma alternativa, mas gostaria de decifrar este este problema para uso em outras situações. Segue abaixo o script e uma parte do banco para testes. Ajudas são bem vindas! ---------------------------------------------------------------------------------------------------------------- ##Banco id3 - banco parcial para testes load("id3.rda") attach(id3) id3$q11 <- as.factor(id3$q11) id3$q13 <- as.factor(id3$q13) table(q11,q13, useNA="always") #Não está informando o pvalor ao usar os pesos amostrais! - variável para peso="bwgr_et" ou "peso" 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) ############################################### #reportando erro em optim(s0, fmin, gmin, method = "BFGS", ...) : valor inicial em vmmin não é finito m1 <- polr(q13 ~ q11 , method="logistic", data=id3, weights=peso, Hess=TRUE) summary(m1) ordinal.or.display(m1) __________________________________________________ Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS -------------------------------------------- Em sáb, 29/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais Para: r-br@listas.c3sl.ufpr.br Data: Sábado, 29 de Novembro de 2014, 22:25 Se você quer usar pesos de amostragem, o pacote survey é a única forma que conheço no R. Se você quer usar as variáveis do data frame "id", precisa passá-lo como argumento para a função svydesign. Por favor, forneça um código mínimo reprodutível, incluindo acesso a dados e os dois comandos (svyolr e svydesign com dados). -- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638 Em Sex 28 nov. 2014, às 16:27, Luciane Maria Pilotto escreveu:
Ah, não Leonardo, já havia alterado e dá o mesmo erro!!
Enviado do Windows
--------------------------------------------
Em qui, 27/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 27 de Novembro de 2014, 21:11
No comando abaixo:
dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL, fpc=NULL, data = NULL, nest = FALSE, weights=~peso,pps=FALSE) Você provavelmente quer passar o argumento
"data = id", em vez de "data
= NULL". -- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638 Em Qui 27 nov. 2014, às 14:58, Luciane Maria Pilotto escreveu: > Não consegui usar o comando svyolr, está dando os erros abaixo. > > # Utilizando as variáveis criadas q11 e q13 que foram transformadas em > fator: >
> > id$q11 <- as.factor(id$q11)
id$q13 <- as.factor(id$q13) > > > > dstrat<-svydesign(ids=~1,
probs=NULL, strata = NULL, variable = NULL, fpc=NULL,
> + data = NULL, nest = FALSE, weights=~peso,pps=FALSE) > > dstrat Independent Sampling design (with replacement) > svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL, fpc = NULL, data = NULL, nest = FALSE, weights = ~peso, pps = FALSE) > >
svyolr(q13~q11,design=dstrat, method="logistic") > Erro em svyolr.survey.design2(q13 ~ q11, design = dstrat, method = > "logistic") : response must be a factor
> # Utilizando as variáveis quest_11 e quest_13 que são fatores no banco: > > svyolr(quest_13~quest_11,design=dstrat, method="logistic") > Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { : arguments must have same length >
========================================================== > > Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no > entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab > (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot = FALSE)). > > Dicas com outros comandos que utilizam os pesos amostrais são
bem
> vindas!!! > >
> Luciane Maria Pilotto > Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS >
> >
> Em ter, 25/11/14, Leonardo Ferreira Fontenelle > <leonardof@leonardof.med.br> escreveu: >
>
Assunto: Re: [R-br] Regressão logística ordinal e uso de
pesos amostrais > Para: r-br@listas.c3sl.ufpr.br > Data: Terça-feira, 25 de Novembro de 2014, 17:09 > > > > >
> Não serve a > http://r-survey.r-forge.r-project.org/survey/html/svyolr.html? > > > Leonardo Ferreira > Fontenelle > > > > > Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br > escreveu: > Olá, > > > estou tentando rodar regressão logística ordinal e > não estou conseguindo usar os pesos amostrais, está dando > o seguinte erro: > > > m11 <- polr(q13 ~ q11 , method="logistic", weights = peso, data=id, Hess=TRUE) > > Aviso: glm.fit: fitted probabilities numerically 0 or 1 > occurred > > Erro em optim(s0, fmin, gmin, method
=
> "BFGS", ...) : > > valor inicial em vmmin não é finito > > *A variável q13 possui 5 ctegorias (1,2,3,4,5) e a > variável q11 possui 4 categorias (1,2,3,4). > > > Tem como solicitar o qui-quadrado, o número de > observações e o pseudo R2 no mesmo comando utilizado para > reg. log. ord.? > > > Como utilizar o peso amostral nas tabelas de > contingência? Estou usando os seguintes comandos: > > table(q11,q13) > > rowPercents(table(q11,q13)) > > ou > CrossTable(q11, q13, digits=2, > expected=TRUE, prop.r=TRUE, prop.c=TRUE, > > > prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE, >
>
> format=c("SPSS"), dnn = c("local > atend","satisfação")) > > > Atenciosamente, > Luciane M. Pilotto > Doutoranda em Saúde Bucal Coletiva/UFRGS > > Enviado do Email do Windows > > > > >
> > R-br mailing list > > R-br@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. > > >
> -----Anexo incorporado----- >
> R-br mailing list > R-br@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@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@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@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@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.
Link do banco: https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0 -------------------------------------------- Em sáb, 13/12/14, Luciane Maria Pilotto <lutipilotto@yahoo.com.br> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@listas.c3sl.ufpr.br Data: Sábado, 13 de Dezembro de 2014, 0:21 Boa noite, 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. A função polr tb permite usar os pesos amostrais, porém, não consigo identificar o erro que está acusando. A função svyolr pode ser uma alternativa, mas gostaria de decifrar este este problema para uso em outras situações. Segue abaixo o script e uma parte do banco para testes. Ajudas são bem vindas! ---------------------------------------------------------------------------------------------------------------- ##Banco id3 - banco parcial para testes load("id3.rda") attach(id3) id3$q11 <- as.factor(id3$q11) id3$q13 <- as.factor(id3$q13) table(q11,q13, useNA="always") #Não está informando o pvalor ao usar os pesos amostrais! - variável para peso="bwgr_et" ou "peso" 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) ############################################### #reportando erro em optim(s0, fmin, gmin, method = "BFGS", ...) : valor inicial em vmmin não é finito m1 <- polr(q13 ~ q11 , method="logistic", data=id3, weights=peso, Hess=TRUE) summary(m1) ordinal.or.display(m1) __________________________________________________ Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS -------------------------------------------- Em sáb, 29/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais Para: r-br@listas.c3sl.ufpr.br Data: Sábado, 29 de Novembro de 2014, 22:25 Se você quer usar pesos de amostragem, o pacote survey é a única forma que conheço no R. Se você quer usar as variáveis do data frame "id", precisa passá-lo como argumento para a função svydesign. Por favor, forneça um código mínimo reprodutível, incluindo acesso a dados e os dois comandos (svyolr e svydesign com dados). -- Leonardo Ferreira Fontenelle http://lattes.cnpq.br/9234772336296638 Em Sex 28 nov. 2014, às 16:27, Luciane Maria Pilotto escreveu:
Ah, não Leonardo, já havia alterado e dá o mesmo erro!!
Enviado
do Windows
--------------------------------------------
Em qui, 27/11/14, Leonardo Ferreira Fontenelle <leonardof@leonardof.med.br> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
Para: r-br@listas.c3sl.ufpr.br Data: Quinta-feira, 27 de Novembro de 2014, 21:11
No comando abaixo:
dstrat<-svydesign(ids=~1,
probs=NULL, strata = NULL, variable =
NULL,
fpc=NULL, data = NULL, nest = FALSE,
weights=~peso,pps=FALSE)
Você provavelmente quer passar o argumento
"data = id", em vez de "data
= NULL".
--
Leonardo Ferreira
Fontenelle http://lattes.cnpq.br/9234772336296638 Em Qui 27 nov. 2014, às
14:58, Luciane Maria Pilotto escreveu:
> Não consegui usar o comando svyolr, está dando os erros abaixo.
> > #
Utilizando as variáveis criadas q11 e q13 que foram
transformadas em
> fator: >
> > id$q11 <-
as.factor(id$q11)
id$q13 <-
as.factor(id$q13)
>
> > >
dstrat<-svydesign(ids=~1,
probs=NULL, strata = NULL, variable = NULL, fpc=NULL, > + data = NULL, nest = FALSE, weights=~peso,pps=FALSE) > > dstrat
Independent Sampling design (with replacement)
> svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL,
fpc = NULL, data = NULL,
nest = FALSE, weights = ~peso, pps = FALSE) >
svyolr(q13~q11,design=dstrat,
method="logistic")
> Erro em
svyolr.survey.design2(q13 ~ q11, design = dstrat, method
= > "logistic") :
response must be a factor
> # Utilizando as variáveis quest_11 e quest_13 que são fatores no banco: > > svyolr(quest_13~quest_11,design=dstrat, method="logistic") > Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { :
arguments must have same
length >
==========================================================
> > Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no > entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab > (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot
= FALSE)). > > Dicas com outros comandos que utilizam os pesos amostrais são
bem
> vindas!!! > >
__________________________________________________
> Luciane Maria Pilotto > Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS
> >
--------------------------------------------
> Em ter, 25/11/14, Leonardo Ferreira Fontenelle > <leonardof@leonardof.med.br> escreveu:
>
Assunto: Re: [R-br] Regressão logística ordinal e uso de
pesos amostrais Para: r-br@listas.c3sl.ufpr.br > Data: Terça-feira, 25 de Novembro de 2014, 17:09 >
>
> > >
> Não serve a
http://r-survey.r-forge.r-project.org/survey/html/svyolr.html?
>
> Leonardo Ferreira > Fontenelle > >
> Em Ter 25 nov. 2014, às 15:36, lutipilotto@yahoo.com.br > escreveu:
Olá, > > > estou tentando
rodar regressão logística ordinal e > não estou conseguindo usar os pesos amostrais, está dando o seguinte erro:
> > > m11 <- polr(q13 ~ q11 ,
method="logistic", weights = peso, data=id,
Hess=TRUE)
> > Aviso: glm.fit: fitted probabilities numerically 0 or 1 occurred > > Erro em optim(s0, fmin, gmin, method
=
> "BFGS", ...) : > > valor inicial em vmmin não é finito > > *A variável
q13 possui 5 ctegorias (1,2,3,4,5) e
a
> variável q11 possui 4 categorias
(1,2,3,4).
>
>
> Tem como solicitar o qui-quadrado, o
número de
> observações e o pseudo R2 no mesmo comando utilizado para > reg. log. ord.? >
>
Como utilizar
o peso amostral nas tabelas de > contingência? Estou usando os seguintes
comandos:
>
>
table(q11,q13)
> > rowPercents(table(q11,q13)) > ou >
CrossTable(q11,
q13, digits=2, > expected=TRUE, prop.r=TRUE, prop.c=TRUE, > >
> prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE,
>
> format=c("SPSS"), dnn = c("local
atend","satisfação"))
> > > Atenciosamente,
>
Luciane M. Pilotto
> Doutoranda em Saúde Bucal Coletiva/UFRGS
> > Enviado do Email do Windows > > > > >
_______________________________________________
>
R-br mailing
list > R-br@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.
>
>
-----Anexo incorporado-----
>
_______________________________________________
> R-br mailing list > R-br@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@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@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@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@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.
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
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@gmail.com> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Para: r-br@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@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.
2014-12-18 19:52 GMT-05:00 Luciane Maria Pilotto <lutipilotto@yahoo.com.br>:
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!
Talvez o p-valor seja diferente nos dois casos porque o código usado para exibi-lo na tela é o mesmo nos dois casos. A sequência é a seguinte: 1. A função crosstab usa a função xtabs() para fazer uma tabela de X & Y considerando os pesos. 2. A tabela é repassada para a função CrossTable(), que cria um objeto com a tabela, valores esperados, teste de qui-quadrado, etc. 3. A função print.CrossTable() formata o objeto para exibição. Ou seja, tanto faz usar crosstab() (com ou sem pesos) como usar CrossTable() diretamente porque quem exibe o p-valor é print.CrossTable(). Eu mantenho o pacote descr e já fiz a correção no github (https://github.com/jalvesaq/descr), mas ainda não enviei para o CRAN. Se você não estiver usando Windows, será fácil compilar o pacote a partir do código fonte. -- Jakson
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@gmail.com> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais |
CROSSTAB
Para: r-br@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
No link https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0 __________________________________________________ Luciane M. Pilotto -------------------------------------------- Em sex, 19/12/14, Mauro Sznelwar <sznelwar@uol.com.br> escreveu: Assunto: [R-br] RES: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@listas.c3sl.ufpr.br Data: Sexta-feira, 19 de Dezembro de 2014, 0:30 Onde tem este arquivo "id3.rda" para rodar?
Obrigado! Vou tentar rodar No link https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0 __________________________________________________ Luciane M. Pilotto Onde tem este arquivo "id3.rda" para rodar? --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com
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@gmail.com> escreveu:
Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais |
CROSSTAB
Para: r-br@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@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@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
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@uol.com.br> escreveu: Assunto: [R-br] RES: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@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@gmail.com> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@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@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@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@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.
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@uol.com.br> escreveu: Assunto: [R-br] RES: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@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@gmail.com> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais | CROSSTAB Para: r-br@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@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@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@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@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
Os comandos estão abaixo. ########################################################### 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@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.
participantes (4)
-
Jakson Alves de Aquino -
Leonardo Ferreira Fontenelle -
Luciane Maria Pilotto -
Mauro Sznelwar