Linear model with Constraints
I'm quite new to R and I have a following problem:
I have a simple 2-factor linear model:
Rate~factor1 + factor2 //factor1 has 8 categorical values, factor2 has 6 categories;
model1 <- lm(Rate~factor1+factor2, data=myData)
And want to put const开发者_开发问答raints SUM of factor1 coefficients = 0, the same for factor2.
None of the manuals gives any clue how to do this.
I found a link to similar problem here but it is different and I couldn't figure out how to modify it...
It's described in chapter 6 of MASS (Modern Applied Statistics with S). Use the contrasts
arg of lm
(take a look at ?contr.sum
and ?model.matrix.default
for examples).
精彩评论