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

[Leonardo Ferreira Fontenelle]

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.