提前感谢您的帮助。我是R的新手，所以请忽略我代码中的任何明显错误。
我试图用我创建的数组替换矩阵中的每个非零值。以下是我到目前为止的情况：
totalsquares<-matrix(data=rbinom(225,1,0.65), nrow=15, ncol=15, byrow=TRUE)
samplesquares <- sum(totalsquares)
totalsquares[totalsquares == 1] <- rpois(samplesquares,30)
urchins<-sum(totalsquares)
for(rowcount in 1:15){
for(colcount in 1:15){
if(totalsquares[rowcount, colcount] != 0){
capprob<-rbeta(1,0.48,0.12)
cmr<-array(data=rbinom((5*totalsquares[rowcount,colcount]), 1, capprob), c(totalsquares[rowcount, colcount],5,samplesquares))
}
}
}
这里的想法是创建一个代表我的采样区域的网格，然后我对它的大约65％进行采样。在我采样的每个方格中，我都有一些基于密度约为30 /平方的海胆。
循环是为捕获标记 - 重新捕获模型创建遇到历史记录。它应该为每个采样的正方形创建一个二项式数据数组，其中维度是采样个体的数量乘以5的平方数。
主要问题是我无法弄清楚如何用代表该样本的数组替换矩阵中的值。我也不认为我的循环工作正常，因为每个数组的每个数组应该有不同的行数。
感谢任何帮助！