R语言中执行diag(Gl)=diag(Gl)+1是否等价于直接修改原矩阵Gl

在R语言环境下，diag(Gl) = diag(Gl) + 1 操作等价于直接修改原矩阵Gl本身，不会生成新的独立副本。

原理说明 #

• R中的diag()属于支持赋值操作的「替换函数」，你写的赋值语句本质是调用了底层的 `diag<-`() 函数，这类函数的设计目标就是直接修改传入对象的对应部分，而非返回新对象。
• 你可以通过下方的测试代码验证行为：

# 初始化3x3全零矩阵
Gl <- matrix(0, nrow = 3, ncol = 3)
print("修改前的Gl：")
print(Gl)

# 执行对角线加1操作
diag(Gl) = diag(Gl) + 1
print("修改后的Gl：")
print(Gl)

运行后可以看到，原矩阵Gl的对角线元素全部变为1，其余元素保持不变，原对象已被直接修改。

注意事项 #

上述行为仅在Gl为matrix类型时生效，如果Gl是数据框data.frame类型，执行diag()赋值会直接报错，操作前可以用class(Gl)确认变量类型。

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