Re: [R] about model.matrix

2015-04-02 Thread Michael Dewey
You cannot have columns which are linearly dependent. Starting a project which uses statistics without having some sort of local statistical backup seems unprofitable. install.packages("fortunes") # if not already done library(fortunes) fortune(122) On 02/04/2015 15:20, Fix Ace wrote: Thank

Re: [R] about model.matrix

2015-04-02 Thread Bert Gunter
This is a big topic. You might try looking for tutorials on "linear models", with "rank" or "rank deficiency" as subtopics. One possible book is: http://www.amazon.com/Linear-Models-Chapman-Statistical-Science/dp/1439887330/ref=sr_1_5?s=books&ie=UTF8&qid=1427987551&sr=1-5&keywords=linear+models+in

Re: [R] about model.matrix

2015-04-02 Thread Fix Ace
Thank you very much for the response. Then what does it mean? I am not a stat person, but have to use it for my project. :( Could you please recommend some readings about it? Thanks a lot! On Wednesday, April 1, 2015 10:58 AM, Michael Dewey wrote: This is really a question about st

Re: [R] about model.matrix

2015-04-01 Thread Michael Dewey
This is really a question about statistics rather than R but see below On 01/04/2015 06:28, Fix Ace wrote: I tried to run the sample code from R: dd <- data.frame(a = gl(3,4), b = gl(4,1,12)) a b 1 1 1 2 1 2 3 1 3 4 1 4 5 2 1 6 2 2 7 2 3 8 2 4 9 3 1 10 3 2 11 3 3 12 3 4 options("contr

[R] about model.matrix

2015-04-01 Thread Fix Ace
I tried to run the sample code from R: dd <- data.frame(a = gl(3,4), b = gl(4,1,12))  a b 1 1 1 2 1 2 3 1 3 4 1 4 5 2 1 6 2 2 7 2 3 8 2 4 9 3 1 10 3 2 11 3 3 12 3 4 options("contrasts") model.matrix(~ a + b, dd)(Intercept) a2 a3 b2 b3 b4 11 0 0 0 0 0 21 0 0