[R-br] Anexando banco: Regressão logística ordinal e uso de pesos amostrais | CROSSTAB
Luciane Maria Pilotto
lutipilotto em yahoo.com.br
Sábado Dezembro 13 00:24:01 BRST 2014
Link do banco: https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0
--------------------------------------------
Em sáb, 13/12/14, Luciane Maria Pilotto <lutipilotto em yahoo.com.br> 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, 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 em leonardof.med.br>
escreveu:
Assunto: Re:
[R-br] Regressão logística ordinal e uso de pesos
amostrais
Para: r-br em 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 em leonardof.med.br>
escreveu:
>
>
Assunto: Re: [R-br]
Regressão logística ordinal e uso de
pesos amostrais
> Para: r-br em 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 em leonardof.med.br>
> escreveu:
>
>
> >
>
Assunto: Re: [R-br]
Regressão logística ordinal e uso
de
> pesos amostrais
>
> Para: r-br em 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 em 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 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.
>
>
>
>
> >
> >
>
> -----Anexo
>
incorporado-----
>
>
> >
>
_______________________________________________
> > 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.
>
_______________________________________________
> 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.
_______________________________________________
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