[R] Seeking help with LOGIT model
Christofer Bogaso
bogaso.christofer at gmail.com
Thu Apr 12 13:18:23 CEST 2012
Thanks Ken for your reply. No doubt your english is quite tough!! I
understand something is not normal with the 5th explanatory variable
(se:2872.17069!) However could not understand what you mean by "You
seem to be getting complete separation on X5 "?
Can you please be more elaborate?
Thanks,
On Thu, Apr 12, 2012 at 4:06 PM, ken knoblauch <ken.knoblauch at inserm.fr> wrote:
> Christofer Bogaso <bogaso.christofer <at> gmail.com> writes:
>> Dear all, I am fitting a LOGIT model on this Data...........
> ---- << snip >>---
>> glm(Data[,1] ~ Data[,-1], binomial(link = logit))
>>
>> Call: glm(formula = Data[, 1] ~ Data[, -1], family = binomial(link = logit))
>>
>> Coefficients:
>> (Intercept) Data[, -1]X 1 Data[, -1]X 2 Data[, -1]X 3 Data[,
>> -1]X 4 Data[, -1]X 5
>> 10.99326 0.01943 10.61013 -0.66763
>> 70.98785 17.33126
>>
>> Degrees of Freedom: 43 Total (i.e. Null); 38 Residual
>> Null Deviance: 44.58
>> Residual Deviance: 17.46 AIC: 29.46
>> Warning message:
>> glm.fit: fitted probabilities numerically 0 or 1 occurred
>>
>> However I am getting a warning mesage as "fitted probabilities
>> numerically 0 or 1 occurred". Here my question is,
> have I made any
>> mistakes with my above implementation? I
> s it just because, I have too
>> less number of '0' in my response Variable?
>>
> Look at the output of summary, especially the standard errors.
> You seem to be getting complete
> separation on X5 and X4 doesn,'t look so hot either.
>
> Ken
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list