[R-br] Regressão logística ordinal e uso de pesos amostrais
Leonardo Ferreira Fontenelle
leonardof em leonardof.med.br
Sábado Novembro 29 22:25:32 BRST 2014
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.
Mais detalhes sobre a lista de discussão R-br