为何自定义R语言对比矩阵完全奇异？求原因及影响因素

问题描述 #

有A、B、C、D、E、F共6个条件，尝试在R中构建用于回归模型的对比矩阵，代码如下：

# Compare A, B, C to D, E, F
contrast1 <- c(1/3, 1/3, 1/3,  -1/3, -1/3, -1/3)
# Next compare A and D to C and F.
contrast2 <- c(1/2, 0, -1/2, 1/2, 0, -1/2) 
# Next compare B and E to C and F
contrast3 <- c(0, 1/2, -1/2, 0, 1/2, -1/2) 
# Next compare A and D to B and E
contrast4 <- c(1/2, -1/2, 0, 1/2, -1/2, 0) 
# Create interaction
contrast.interaction <- contrast1 * contrast4

mat.temp <- rbind(constant=1/6, contrast1, contrast2, contrast3, contrast4, contrast.interaction)
mat.temp

矩阵输出结果：

[,1]       [,2]       [,3]       [,4]       [,5]       [,6]
constant             0.1666667  0.1666667  0.1666667  0.1666667  0.1666667  0.1666667
contrast1            0.3333333  0.3333333  0.3333333 -0.3333333 -0.3333333 -0.3333333
contrast2            0.5000000  0.0000000 -0.5000000  0.5000000  0.0000000 -0.5000000
contrast3            0.0000000  0.5000000 -0.5000000  0.0000000  0.5000000 -0.5000000
contrast4            0.5000000 -0.5000000  0.0000000  0.5000000 -0.5000000  0.0000000
contrast.interaction 0.1666667 -0.1666667  0.0000000 -0.1666667  0.1666667  0.0000000

运行solve(mat.temp)时出现报错：

Error in solve.default(mat.temp) : Lapack routine dgesv: system is exactly singular: U[6,6] = 0

已知矩阵奇异，需明确原因及冗余对比项。

奇异原因分析 #

矩阵奇异的核心是行向量线性相关，即某一行可以通过其他行的线性组合得到，导致矩阵的秩小于行数（6）。具体到你的设计：

1. contrast.interaction是contrast1与contrast4的逐元素乘积，但该行并非独立于其他行。通过rankMatrix(mat.temp)可验证，矩阵秩小于6，说明行向量存在冗余。
2. 你的对比设计中，contrast1是前3组与后3组的整体对比，contrast2、contrast3、contrast4是组内/交叉的细分对比，再加上构造的contrast.interaction，这些对比之间并非完全正交，存在线性依赖关系，最终导致矩阵不可逆。

解决方法 #

要让矩阵可逆，必须保证所有行向量线性无关，可通过以下方式调整：

• 移除冗余对比项：替换contrast.interaction为一个独立的、无冗余的对比项。例如添加A与B的直接对比：

  # 替换contrast.interaction为A vs B的独立对比
  contrast5 <- c(1/2, -1/2, 0, 0, 0, 0)
  mat.temp <- rbind(constant=1/6, contrast1, contrast2, contrast3, contrast4, contrast5)
  
  # 验证矩阵秩（需加载Matrix包）
  library(Matrix)
  rankMatrix(mat.temp) # 返回6，说明行向量线性无关
  solve(mat.temp) # 可正常运行
  

• 使用正交对比构造：手动构造两两正交的对比向量（任意两个行向量的点积为0，除自身外），或利用R内置的contrasts()函数生成标准正交对比，确保矩阵满秩。

内容的提问来源于stack exchange，提问作者Another