我对定义单向anova的对比有一个小问题;我尝试了很多方法,但似乎找不到合适的方法。
我有这种设计:
$ Categorie : Factor w/ 15 levels "Culture","culture_loisirs",..: 9 2 12 12 2 3 12 13 9 13 ...
$ RT : num 932 1142 1238 680 983 ...
在我的One-Way Anova中,我想做一个先验对比如下:
cbind(c1,c2,c3)
c1 c2 c3
[1,] -1.0 -1 -1
[2,] -1.0 -1 -1
[3,] -1.0 -1 -1
[4,] -1.0 -1 -1
[5,] -1.0 -1 -1
[6,] 6.5 13 0
[7,] -1.0 -1 -1
[8,] 6.5 0 13
[9,] -1.0 -1 -1
[10,] -1.0 -1 -1
[11,] -1.0 -1 -1
[12,] -1.0 -1 -1
[13,] -1.0 -1 -1
[14,] -1.0 -1 -1
[15,] -1.0 -1 -1
但是当我尝试重新定义数据$ Categorie的对比时,R返回:
contrasts(data$Categorie) <- cbind(c1,c2,c3)
Error in `contrasts<-`(`*tmp*`, value = c(-1, -1, -1, -1, -1, 6.5, -1,
:singular contrast matrix
如果我使用任何两种对比
,则不会出现此问题contrasts(data$Categorie) <- cbind(c1,c2)
contrasts(data$Categorie) <- cbind(c1,c3)
etc.
你对这是为什么有任何线索吗?