在R中创建对称矩阵

时间:2013-05-01 04:31:32

标签: r

如何在R中创建此矩阵?实际上,例如,矩阵A[1,1]的fisrt元素等于x^2j的乘法之和,A[1,2]等于x的乘法和yj,...

          x = 1:3
          y = c(2, 4, 6)
          z = c(3, 7, 9)
          j = c(5, 10, 15)
         > A
            [,1] [,2] [,3]
       [1,]  180  360  560
       [2,]  360  720 1120
       [3,]  560 1120 1750

1 个答案:

答案 0 :(得分:5)

我的猜测:

R > mat = as.matrix(cbind(x, y,z))
R > mat
     x y z
[1,] 1 2 3
[2,] 2 4 7
[3,] 3 6 9
R > t(mat)%*%diag(j)%*%mat
    x    y    z
x 180  360  560
y 360  720 1120
z 560 1120 1750