我试图从矢量创建这样的矩阵:
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),然后第二列的thirld元素是(9-2 = 7)。然后,第三列以0开头,然后是0和(9-5 = 4),最后一列始终为零。 vec的长度可以更改为任意数字,例如4,5,...。如何编写有效的函数或代码来创建此矩阵?
最佳答案
我认为这将满足您的要求:
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 0
关于r - 在R中创建矩阵的有效方法,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/10851976/