RAS算法生成随机矩阵

Question

试图生成一个随机矩阵，其中包含行和列的约束，对角线上有零。

到目前为止，通过在Ipfp包中使用mipfp函数，我可以生成一个具有行和和约束的随机矩阵，但仍然不能得到等于0的对角线。

InitialM <- array(1,dim=c(5,5))
SumOfRows <- c(40, 60, 20, 10, 10)
SumOfCols <- c(30, 30, 60, 10, 10)
TargetData<- list(SumOfRows, SumOfCols)
Target.list <- list(1,2)
RandomMatrix <- Ipfp(InitialM, Target.list, TargetData, 
   iter=100, tol=1e-10)

有人有什么想法吗？

Answer 1

为了实现您想要做的事情，只需将初始矩阵InitialM的对角线设置为0：

InitialM <- array(1,dim=c(5,5))
for(i in 1:5) {
  InitialM[i,i] <- 0
}

SumOfRows <- c(40, 60, 20, 10, 10)
SumOfCols <- c(30, 30, 60, 10, 10)
TargetData<- list(SumOfRows, SumOfCols)
Target.list <- list(1,2)
RandomMatrix <- Ipfp(InitialM, Target.list, TargetData, 
                     iter=100, tol=1e-10)

因此，最终的RandomMatrix将在对角线上有零，行和上的约束将匹配。