I am new with Rstudio so maybe this is a stupid question.
I have to conduct a one-way ANOVA (with appropriate constrain) on the chick weight with respect to the six feed supplements, without using the aov or lm functions.
First of all I am trying to create a incindence matrix with ones and zeros. I used the function model.matrix but now I have one column to less. I used the code: model.matrix(~feed) and I obtained the next image.
In this image you see the first 10 rows of the matrix. But I am missing the casein column and also the order is changed.
Could please some help me out
Thanks in advance