Skip to content
Prev 299052 / 398506 Next

design matrix creation in R

On Jul 3, 2012, at 4:10 PM, mmstat at comcast.net wrote:

            
What first part if the next part is ... this ?
That's just generating a factor variable.
Something like this?

model.matrix(~a+b, data= dd, contrasts = list(a="contr.sum",  
b="contr.sum") )

    (Intercept) a1 a2 b1 b2 b3
1            1  1  0  1  0  0
2            1  1  0  0  1  0
3            1  1  0  0  0  1
4            1  1  0 -1 -1 -1
5            1  0  1  1  0  0
6            1  0  1  0  1  0
7            1  0  1  0  0  1
8            1  0  1 -1 -1 -1
9            1 -1 -1  1  0  0
10           1 -1 -1  0  1  0
11           1 -1 -1  0  0  1
12           1 -1 -1 -1 -1 -1
attr(,"assign")
[1] 0 1 1 2 2 2
attr(,"contrasts")
attr(,"contrasts")$a
[1] "contr.sum"

attr(,"contrasts")$b
[1] "contr.sum"