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

Luciane Maria Pilotto lutipilotto em yahoo.com.br
Sexta Novembro 28 16:27:34 BRST 2014


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.brhttps://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