求助，使用R生成对称矩阵，sas的方法也可以

希望生成一个5×5的对称矩阵，固定其中一侧数值的最大值、最小值和平均数。是想生成类似下面这样的，例子这个的均值为0.3，范围为0.2-0.4

先生成对角线，再生成上三角或下三角

1. m <- diag(1, 5)
   
2. m[lower.tri(m)] <- runif(length(m[lower.tri(m)]), 0.2, 0.4)
   
3. m[upper.tri(m)] <- m[lower.tri(m)]

复制代码

或者是，如果不能直接生成矩阵的话，有没有办法生成有限制的随机数，限制范围与均值，虽然这样不算是随机数

owenqi 发表于 2021-4-27 23:29
先生成对角线，再生成上三角或下三角

你的代码还要再加一句

1. m[upper.tri(m)] <- t(m)[upper.tri(m)]

复制代码

llb_321 发表于 2021-4-28 09:49
你的代码还要再加一句

感谢两位老师的回复，还想问一下，再限制均值为0.30的且限制两位小数的条件的话，是可行的吗

juicejune 发表于 2021-4-28 10:46
感谢两位老师的回复，还想问一下，再限制均值为0.30的且限制两位小数的条件的话，是可行的吗

把owenqi老师的代码中的这句
m[lower.tri(m)] <- runif(length(m[lower.tri(m)]), 0.2, 0.4)
加一个round()即可限制小数位数，如下
m[lower.tri(m)] <- round(runif(length(m[lower.tri(m)]), 0.2, 0.4),2)