复制100次

Question

我想从二项分布中抽取100个大小为187的随机样本作为固定概率值。我在这里附加了我的R代码，在这里我得到了1个随机样本，大小为187。但是我无法将结果存储100次。如何储存？这意味着将有100列，每列将有187个观测值。

mdata=read.csv("c:/germany/st_it.csv",header=TRUE)###susceptibe and infective number on a day t
t<-read.csv("c:/germany/pt_new.csv",header=TRUE)
t1=t$PRO
s=mdata$st  ## susceptible
s1=mdata$it ### infected
q=0.9946289
b=numeric(nrow(mdata))
for (i in 1:(nrow(mdata))){
  b[i]=rbinom(1,s[i],q^(s1[i]))
}
sus=numeric(187)
sus[1]=b[1]
for (j in 2:187) {
  ifelse(t1[j]==t1[j-1], sus[j]<-sus[j-1], sus[j]<-b[j])
}
sus