Hi,

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