我已经生成了一个观察到的矩阵,这里是代码:
obs.matrix <- matrix(c(rep(1,10),rep(2,10)),nrow=10,ncol=2)
现在我想获得3000个置换数据集,在每个数据集中,应该有10个1和10个2,但它们可以位于不同的列中。
我不知道怎么做其余的事。 我试过但失败了。
x =obs.matrix
theta = function(resample){sample(c(1,2),replace = T)}
result <- bootstrap::bootstrap(x,3000,theta)
感谢您的帮助。