[R] how to understand the interaction of two fixed effects
genomenet at gmail.com
genomenet at gmail.com
Thu Mar 15 01:13:51 CET 2007
Hi There,
I have two questions about how to understand well about the
interaction efffect.
1) Suppose two factors, A and B.
A has n levels, B has m levels.
Why the degree of freedom of interaction effect of A and B ( here I
mean A:B not A*B) is (n-1)*(m-1) not n*m-1?
2) can Lard:Gender be understood as LardGender?
Value=c(709,679,699,657,594,677,592,538,476,508,505,539)
Lard=rep(c("Fresh","Rancid"),each=6)
Gender=rep(c("Male","Male","NONE","NONE","Female","Female"),2)
LardGender=rep(c("FreshMale","FreshNONE","FreshFemale","RancidMale","RancidNONE","RancidFemale"),each=2)
Food=data.frame(Value,Lard, Gender,LardGender)
Why is the following two return outputs of lm() function different?
> lm(Value~Lard+Gender+Lard:Gender,data=Food)
Call:
lm(formula = Value ~ Lard + Gender + Lard:Gender, data = Food)
Coefficients:
(Intercept) LardRancid GenderMale
635.5 -113.5 58.5
GenderNONE LardRancid:GenderMale LardRancid:GenderNONE
42.5 -15.5 -72.5
> lm(Value~Lard+Gender+LardGender,data=Food)
Call:
lm(formula = Value ~ Lard + Gender + LardGender, data = Food)
Coefficients:
(Intercept) LardRancid GenderMale
635.5 -113.5 43.0
GenderNONE LardGenderFreshMale LardGenderFreshNONE
-30.0 15.5 72.5
LardGenderRancidFemale LardGenderRancidMale LardGenderRancidNONE
NA NA NA
Thank you very much!
Fan
More information about the R-help
mailing list