Matrix包对称矩阵Dimnames异常：赋值NULL无法移除行/列名？

我碰到过这个版本的Matrix包的坑！这是1.2-18版本里针对对称稀疏矩阵（比如dsCMatrix）的一个特殊行为问题——因为对称矩阵的行、列维度是绑定的，包的旧版本会强制行名和列名保持同步，直接用rownames() <- NULL或者colnames() <- NULL的方式根本不生效，甚至会触发奇怪的互换问题。

先复现一下你遇到的情况，方便对照：

library(Matrix)

# 创建2x2对称稀疏矩阵
m <- Matrix(c(1, 0, 0, 1), nrow = 2, ncol = 2, sparse = TRUE)
m <- as(m, "dsCMatrix")

# 设置列名后，行名自动同步
colnames(m) <- c("a", "b")
print(m)
# 输出里行名也变成了"a"、"b"

# 尝试移除行名，无效
rownames(m) <- NULL
print(m)
# 行名依然存在

# 尝试移除列名，出现异常互换
colnames(m) <- NULL
print(m)
# 列名没消失，Dimnames的行、列元素反而互换了

为什么会这样？ #

在Matrix 1.2-18的旧版本中，对称矩阵的Dimnames属性被设计成强制同步的——当你修改行名或列名时，包会自动把另一个维度的名称也改成一样的。而直接赋值NULL会触发内部的同步逻辑bug，导致出现维度名称互换的异常。

解决办法 #

有两种方案可以选：

方案1：直接修改底层的Dimnames属性

绕过rownames()/colnames()的封装函数，直接操作矩阵的Dimnames属性，这样就能独立设置行、列名：

# 仅保留列名，移除行名
attr(m, "Dimnames") <- list(NULL, c("a", "b"))
print(m)

# 验证结果
rownames(m)  # 返回NULL，行名已移除
colnames(m)  # 返回c("a", "b")，列名正常保留

方案2：升级Matrix包到新版本

这个问题在Matrix包的后续版本（比如1.3.x及以上）已经被修复了，新版本中对称矩阵的行名和列名可以独立修改，赋值NULL也能正常移除对应维度的名称。如果你的环境允许升级，执行以下命令即可：

install.packages("Matrix")

升级之后再测试，就不会再碰到这种奇怪的同步和互换问题了。

内容的提问来源于stack exchange，提问作者Øyvind Langsrud