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 constraints 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).
If you love us? You can donate to us via Paypal or buy me a coffee so we can maintain and grow! Thank you!
Donate Us With