[R-br] Material Regressão logística múltipla

[Leonardo Ferreira Fontenelle]

Para incorporar os pesos amostrais é necessário usar o pacote survey. Se
não me engano, o peso do polr é de precisão analítica.

-- 
Leonardo Ferreira Fontenelle
http://lattes.cnpq.br/9234772336296638

Em Qui 15 jan. 2015, às 10:10, Luciane Maria Pilotto escreveu:
> Olá, 
> 
> estou tentando usar a função "polr" mas não está rodando com os pesos
> amostrais. Será que podes ajudar?
> 
> Link do banco reduzido:
> https://www.dropbox.com/s/9qrdf4mhd6tzypi/id3.rda?dl=0
> 
> id$q13 <- as.factor(id$q13)
> mp <- polr(q13 ~ q11 , method="logistic", weights = bwgr_et, data=id,
> Hess=TRUE)
> 
> #Erro que está dando:
> Erro em optim(s0, fmin, gmin, method = "BFGS", ...) : 
>   valor inicial em vmmin não é finito
> Além disso: Mensagens de aviso perdidas:
> 1: In eval(expr, envir, enclos) :
>   non-integer #successes in a binomial glm!
> 2: glm.fit: fitted probabilities numerically 0 or 1 occurred 
> 
> summary(mp)
> ordinal.or.display(mp)
> 
> __________________________________________________
> Luciane Maria Pilotto
> 
> --------------------------------------------
> Em qua, 14/1/15, Felipe <felipe.e.barletta em gmail.com> escreveu:
> 
>  Assunto: Re: [R-br] Material Regressão logística múltipla
>  Para: r-br em listas.c3sl.ufpr.br, "ELIAS MENDES COSTA MENDES"
>  <eliasmccosta em yahoo.com.br>
>  Data: Quarta-feira, 14 de Janeiro de 2015, 17:20
>  
>  
>      Boa tarde, não sei se já lhe responderam.
>  
>      Talvez o script abaixo lhe ajude:
>  
>      
>  
>      require(MASS)
>  
>     
>   melhora<-rep(c("ac","alg","nenh"),4)
>  
>       sexo<-c(1,1,1,1,1,1,0,0,0,0,0,0)
>  
>       trat<-c(1,1,1,0,0,0,1,1,1,0,0,0)
>  
>       Freq<-c(16,5,6,6,7,19,5,2,7,1,0,10)
>  
>       artrite<-cbind(melhora,sexo,trat)
>  
>       artrite<-as.data.frame(artrite)
>  
>       attach(artrite)
>  
>       options(contrasts = c("contr.treatment",
>  "contr.poly"))
>  
>       ajust1 <- polr(melhora ~ sexo + trat + sexo*trat,
>  weights =
>      Freq, data=artrite)
>  
>       ajust1
>  
>       summary(ajust1)
>  
>       ajust2 <- polr(melhora ~ sexo + trat, weights=
>  Freq, data =
>      artrite)
>  
>       ajust2
>  
>       summary(ajust2)
>  
>       ajust2$fitted.values
>  
>      
>  
>      
>  
>      Uma boa forma de procurar algo que deseja referente ao R
>  é através
>      do comando:
>  
>      >RSiteSearch('Assunto')
>  
>      
>  
>      
>  
>      
>  
>      On 08-01-2015 17:17,
>  ELIAS MENDES COSTA
>        MENDES wrote:
>  
>      
>      
>        
>          Boa tarde
>  pessoal,
>              Algum de vocês tem material sobre regressão
>  logística
>              múltipla? (De preferencia usando o
>  R) 
>          
>  
>            
>          Desde
>              de já agradeço 
>           
>          
>  
>          
>          Att,
>          Elias Mendes Costa
>  
>              
>  
>            
>        
>        
>  
>        
>        
>  
>        _______________________________________________
>  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.