一种在R中创建矩阵的有效方法
一种在R中创建矩阵的有效方法
提问于 2012-06-01 09:38:07
我正在尝试从一个向量中创建一个这样的矩阵：
vec= c(2, 5, 9)
> A
[,1] [,2] [,3] [,4]
[1,] 2 0 0 0
[2,] 5 3 0 0
[3,] 9 7 4 0实际上，第一列总是向量元素，第二列从0开始，然后是(5-2 = 3)，然后第二列的第三个元素是(9-2 = 7)。然后第三列从0开始，然后是0和(9-5 = 4)，最后一列总是零。可以将vec的长度更改为任何数字，例如4，5，... .How我可以编写一个有效的函数或代码来创建此矩阵吗？
回答 2
Stack Overflow用户
回答已采纳
发布于 2012-06-01 10:45:04
我认为这将会做你想要的：
f = function(vec)
{
n = length(vec)
M = matrix(0,n,n+1)
M[,1] = vec
for(i in 1:n) M[,i+1] = c(rep(0,i),vec[-c(1:i)]-vec[i])
return(M)
}
vec = c(2,5,9)
f(vec)
[,1] [,2] [,3] [,4]
[1,] 2 0 0 0
[2,] 5 3 0 0
[3,] 9 7 4 0Stack Overflow用户
发布于 2012-06-01 13:25:01
我不知道效率如何，但这里有两个不使用for循环的解决方案：
n <- length(vec)
A <- replicate(n+1, vec) - cbind(0, t(replicate(n, vec)))
A[upper.tri(A)] <- 0这个比较长，但只创建了一个矩阵
n <- length(vec)
A <- replicate(n, vec)
A <- A - t(A)
A <- cbind(vec, A)
A[upper.tri(A)] <- 0页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/10851976
复制相关文章
相似问题
腾讯云开发者
