R语言 随机扰动矩阵，保证新矩阵的每行加和 以及 每列加和与原矩阵相同

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有一个0 1矩阵 类似
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    1    1    1
[2,]    1    0    0    1    0
[3,]    0    0    0    1    1
[4,]    0    1    0    0    1

想要随机后  每行及每列的加和不变。用了vegan包中的permatswap函数
> mat
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    1    1    1
[2,]    1    0    0    1    0
[3,]    0    0    0    1    1
[4,]    0    1    0    0    1
> out <- permatswap(mat, times = 9, burnin = 20000, thin = 500, mtype = "prab")
> out$perm[[1]]
     [,1] [,2] [,3] [,4] [,5]
[1,]    1    0    1    1    1
[2,]    1    0    0    1    0
[3,]    0    1    0    0    1
[4,]    0    0    0    1    1
确实能达到我想要的效果，但是我的矩阵有14000x9000这么大，随机起来很慢。
请问大家有没有更高效的办法，拜托了！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：R语言 swap VEGA Mat ERM