[R] GLM fitting

BXC (Bendix Carstensen) bxc at steno.dk
Fri Jan 28 11:45:14 CET 2005


You have a perfect separtaion of y by x, i.e.

y == (x>30)

is true for all units.

Bendix Carstensen
----------------------
Bendix Carstensen
Senior Statistician
Steno Diabetes Center
Niels Steensens Vej 2
DK-2820 Gentofte
Denmark
tel: +45 44 43 87 38
mob: +45 30 75 87 38
fax: +45 44 43 07 06
bxc at steno.dk
www.biostat.ku.dk/~bxc
----------------------



> -----Original Message-----
> From: r-help-bounces at stat.math.ethz.ch 
> [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Vito Ricci
> Sent: Friday, January 28, 2005 11:14 AM
> To: r-help at stat.math.ethz.ch
> Subject: [R] GLM fitting
> 
> 
> DeaR R-useRs,
> 
> I'm trying to fit a logist model with these data:
> 
> > dati
>    y  x
> 1  1 37
> 2  1 35
> 3  1 33
> 4  1 40
> 5  1 45
> 6  1 41
> 7  1 42
> 8  0 20
> 9  0 21
> 10 0 25
> 11 0 27
> 12 0 29
> 13 0 18
> 
> I use glm(), having this output:
> 
> > g<-glm(y~x,family=binomial,data=dati)
> Warning messages: 
> 1: Algorithm did not converge in: glm.fit(x = X, y =
> Y, weights = weights, start = start, etastart =
> etastart,  
> 2: fitted probabilities numerically 0 or 1 occurred
> in: glm.fit(x = X, y = Y, weights = weights, start =
> start, etastart = etastart,  
> > g
> 
> Call:  glm(formula = y ~ x, family = binomial, data =
> dati) 
> 
> Coefficients:
> (Intercept)            x  
>     -348.23        11.23  
> 
> Degrees of Freedom: 12 Total (i.e. Null);  11 Residual
> Null Deviance:      17.94 
> Residual Deviance: 7.011e-10    AIC: 4 
> 
> I don't understand the meaning of warning. Can anyone
> help me? Many thanks.
> Cheers,
> Vito
> 
> 
> =====
> Diventare costruttori di soluzioni
> Became solutions' constructors
> 
> "The business of the statistician is to catalyze 
> the scientific learning process."  
> George E. P. Box
> 
> Top 10 reasons to become a Statistician
> 
>      1. Deviation is considered normal
>      2. We feel complete and sufficient
>      3. We are 'mean' lovers
>      4. Statisticians do it discretely and continuously
>      5. We are right 95% of the time
>      6. We can legally comment on someone's posterior distribution
>      7. We may not be normal, but we are transformable
>      8. We never have to say we are certain
>      9. We are honestly significantly different
>     10. No one wants our jobs
> 
> 
> Visitate il portale http://www.modugno.it/
> e in particolare la sezione su Palese  
> http://www.modugno.it/archivio/palese/
> 
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list 
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read 
> the posting guide! http://www.R-project.org/posting-guide.html
>




More information about the R-help mailing list