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@gmail.com>
Para: r-br@listas.c3sl.ufpr.br; Alexandre Santos <alexandresantosbr@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@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@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________
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.

--
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.