R design.matrix проблема - упал столбец в матрице дизайна?

У меня возникла странная проблема при попытке настроить матрицу проектирования для последующего анализа парных дифференциальных выражений на данных RNAseq. Для матрицы дизайна у меня есть и информация о доноре, и каждое условие:

group<-factor(y$samples$group) #44 samples, 6 different conditions
sample<-factor(y$samples$samples) #44 samples, 11 different donors.
design<- model.matrix(~0+sample+group)
head(design)

      Donor11.CD8    Donor12.CD8    Donor14.CD8   Donor15.CD8   Donor16.CD8
1             1             0             0             0             0
2             1             0             0             0             0
3             1             0             0             0             0
4             1             0             0             0             0
5             1             0             0             0             0
6             1             0             0             0             0
       Donor17.CD8    Donor18.CD8   Donor19.CD8    Donor20.CD8   Donor3.CD8
1             0             0             0             0            0
2             0             0             0             0            0
3             0             0             0             0            0
4             0             0             0             0            0
5             0             0             0             0            0
6             0             0             0             0            0
        Donor4.CD8        Treatment2         Treatment3  Treatment4 Treatment5
1            0                  0                  0           0             0
2            0                  0                  0           0             1
3            0                  0                  0           1             0
4            0                  0                  0           0             0
5            0                  0                  1           0             0
6            0                  1                  0           0             0
        Treatment6
1             1
2             0
3             0
4             0
5             0
6             0
> 

Проблема в том, что я, кажется, теряю условие (лечение 1), когда формирую матрицу дизайна, и я не уверен почему.

Заранее большое спасибо за помощь!

1 ответ

Решение

Это не проблема. Обработка 1 обозначена всеми 0 для столбцов в матрице дизайна. Посмотрите на строку 4 - ноль для процедур со 2 по 6. Это означает, что это лечение 1. Это называется "контрастом лечения", потому что коэффициенты в модели противопоставляют названное лечение "базовому" уровню, в данном случае базовому уровню такое лечение1.

Другие вопросы по тегам