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@gmail.com> escreveu:
Assunto: Re: [R-br] Material Regressão logística múltipla Para: r-br@listas.c3sl.ufpr.br, "ELIAS MENDES COSTA MENDES" <eliasmccosta@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@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@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@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.