[R] FW: Levels and GLM
Kuhn, Max
Max.Kuhn at pfizer.com
Fri Jul 7 22:17:06 CEST 2006
One correction... since you are fitting a logistic model, it is
technically correct to say the "mean value of the linear predictor,"
instead of "mean response".
20 lashes for me.
Max
-----Original Message-----
From: Kuhn, Max
Sent: Friday, July 07, 2006 4:11 PM
To: 'r-help at stat.math.ethz.ch'
Subject: [R] Levels and GLM
jdrapp,
By default, R fits full rank models. If you are coming from SAS, you're
probably used to less than full rank model parameterizations.
>From Section 11.1.1 of "An Introduction to R" at
http://cran.r-project.org/doc/manuals/R-intro.html#Contrasts
there is this:
"What about a k-level factor A? The answer differs for unordered and
ordered factors. For unordered factors k - 1 columns are generated
for the indicators of the second, ..., kth levels of the factor.
(Thus the implicit parameterization is to contrast the response at
each level with that at the first.)"
So level "M" is the "reference cell". Assuming that
data.logistic$Overall is continuous, the intercept is the estimate of
the mean response when maj = "M" and data.logistic$Overall = 0. The
estimate for majN is the difference between the reference cell
(estimated
by the intercept) and the mean response when maj = "N" and
data.logistic$Overall = 0.
You should check out ?model.matrix and ?contrasts.
Max
> I am using the as.factor command to use with glm. When I use the
command
>
> >maj <- as.factor(data.logistic$Majors)
> >maj
>
> I receive the following output:
> [1] M M N M M M M N N M M M N M M M M M M M M M M M N M N N M M N M
> M N M M M M M
> [40] N M N M M N M M M N M N M N M N N N M N M M M M M M N M N M M M
> M M N N M M M
> [79] M M M N N M M N M N M M M M M M M M M M M M M M M N M M M M M N
> M M M M M N M
> [118] M M M N M N N M M M M M M M M N M N M M M M M N M M M M N M M M
> N N M M M N M
> [157] M M M M M M M M M M M M M N M M N N M M N M M M M M M M M M M M
> M M N M N M M
> [196] M N M M M M M M M M N M M M M M M M M N M M M M M M M M M M M M
> M M N M M N N
> [235] M M M M M N M M M M M M N N M M N M M M M M M M M M M M M M M M
> M N M M M M N
> [274] N M M M M M M N M M M M M M M M M M N N M N M M M M M M M M M M
> N M N N M M M
> [313] M M M M M M M N M M M M M N M M M M M M M M M M M M M M M N M M
> M M M M M N M
> [352] M N M N M M N M M M M N M M M M M M M M M M N M M N N
> Levels: M N
>
> When I enter:
>
> > logistic.glm <- glm(data.logistic$X100.Yard.Average ~
data.logistic$Overall + maj, family=binomial)
> > logistic.glm
>
> I receive the following output:
>
> Call: glm(formula = data.logistic$X100.Yard.Average ~
> data.logistic$Overall + maj, family = binomial)
>
> Coefficients:
> (Intercept) data.logistic$Overall majN
> 2.38819 -0.02718 -0.18385
>
> Degrees of Freedom: 377 Total (i.e. Null); 375 Residual
> Null Deviance: 514.5
> Residual Deviance: 410.7 AIC: 416.7
>
> My question: Why is there no output for majM? Any help would be
> greatly appreciated
----------------------------------------------------------------------
LEGAL NOTICE\ Unless expressly stated otherwise, this messag...{{dropped}}
More information about the R-help
mailing list