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.