[R-br] Dúvida em summary em regressão logística

[Alexandre Santos]

Obrigado Vinícius,

       Mas o Benilton já havia chamado a atenção para uma questão crucial: "regressão logística modela a probabilidade de sucesso dado um conjunto de covariaveis", então no modelo original que estava trabalhando existem erros conceituais, uma aproximação do que eu queria que era a obtenção dos coeficientes do intercepto, linear, quadrático e cúbico, ficou:

 pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de mumias
dens<-sort(rep(2^(2:6),10))#### Densidade de pulgoes oferecidos
##
##Prob(Sucesso) = exp(P0+P1*N+P2*N^2+P3*N^3)/(1+exp(P0+P1*N+P2*N^2+P3*N^3))
naoPred<-dens-pred
#
p.model2<-glm(cbind(pred, naoPred)~poly(dens, 3, raw=TRUE), family=binomial)
summary(p.model2)
#

Abs

Alexandre

________________________________
 De: Vinicius Brito Rocha <viniciusbritor em gmail.com>
Para: r-br em listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr em yahoo.com.br> 
Enviadas: Terça-feira, 20 de Novembro de 2012 21:32
Assunto: Re: [R-br] Dúvida em summary em regressão logística
 

Alexandre,

acho que vc chamou o modelo de forma errada.

deveria ser glm(cbind(evento,1-evento)~VA1+VA2+...VAp,data=dados,family=binomial(link = "logit"))

normalmente nos softwares o coeficiente resposta é o do modelo linear, como vc esta usando a funçao logit

log(P(X=evento)/P(X=nao evento))= beta0+beta1*X1+beta2*X2+...+betap*Xp


Abs.

Mas o Walmes provavelmente vai poder te responder melhor.

Abs

Vinicius Brito Rocha


Em 20 de novembro de 2012 16:49, Alexandre Santos <alexandresantosbr em yahoo.com.br> escreveu:

Boa tarde Pessoal,
>       Estou ajustando uma regressão logística e me deparei com a seguinte dúvida:
>      
>###################################################################################################
>#Regressão logística entre a proporção de herbívoros predados e a densidade de herbívoros oferecidos
>#
>pred<-c(1,2,2,3,1,4,2,3,2,3,5,6,5,5,3,7,7,6,2
>,3,15,12,14,12,11,11,11,13,13,13,18,14,27,26
>,17,18,20,22,10,15,29,30,36,40,23,50,30,40,29,52)##Número de herbívoros predados
>dens<-sort(rep(2^(2:6),10))#### Densidade de herbivoros oferecidos
>##
>##Regressão pred/dens = exp(P0+P1*N+P2*N^2+P3*N3)/1+ exp(P0+P1*N+P2*N^2+P3*N3)
>p.model1<-glm(pred/par~dens+I(dens^2)+I(dens^3),family="binomial")
>summary(p.model1)
>#
>
>
>
>
>Minha dúvida é se o coeficiente linear que aparece no summary esta transformado em exp(x)/1+exp(x)
>ou trata-se do valor sem transformação?
>
>
>Obrigado,
>
>-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso - Campus Cáceres
Avenida dos Ramires, s/n
Bairro: Distrito Industrial 
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr em yahoo.com.br alexandre.santos em cas.ifmt.edu.br ======================================================================
>
>
>_______________________________________________
>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.
>


-- 
Vinicius Brito Rocha.
Estatístico e Atuário (IM / UFRJ)
Mestre em Pesquisa Operacional (COPPE / UFRJ)

www.aplicademic.blogspot.com
http://twitter.com/viniciusbritor

"Não se preocupe muito com as suas dificuldades em Matemática, posso assegurar-lhe que as minhas são ainda maiores." - Albert Einstein.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20121120/7a019b89/attachment.html>