r - 如何生成对称随机矩阵？

Question

我想生成一个应该是对称的随机矩阵。

我试过这个:

matrix(sample(0:1, 25, TRUE), 5, 5)

但它不一定是对称的。

我怎样才能做到这一点？

Answer 1

另一个非常有趣的机会是基于以下数学事实:如果 A 是某个矩阵，则 A 乘以其转置始终是对称的。

> A <- matrix(runif(25), 5, 5)
> A %*% t(A)
         [,1]      [,2]      [,3]      [,4]      [,5]
[1,] 1.727769 1.0337816 1.2195505 1.4661507 1.1041355
[2,] 1.033782 1.0037048 0.7368944 0.9073632 0.7643080
[3,] 1.219551 0.7368944 1.8383986 1.3309980 0.9867812
[4,] 1.466151 0.9073632 1.3309980 1.3845322 1.0034140
[5,] 1.104135 0.7643080 0.9867812 1.0034140 0.9376534