我有以下内容:7次治疗,6次治疗,4次重复,并且产量。我使用rep命令生成处理,位置和复制的模式。我现在需要生成yield数据。对此的一个重要部分是该图可以由4,8或12个植物组成(图解)。该计划首先初始化将要采样的1,000,000株植物的载体。我使用for循环填充yield变量。
for (i in 1:reps) y[i]<-mean(sample(vector,plotsize))
yield<-append(yield,y1)
for (i in 1:reps) y[i]<-mean(5+sample(vector,plotsize))
yield<-append(yield,y1)
for (i in 1:reps) y[i]<-mean(10+sample(vector,plotsize))
yield<-append(yield,y1)
第一对语句用于控件。每次治疗比最后一次治疗大5个单位。继续这个过程给出了一个有效的答案。我可以将代码重写为:
for (j in 1:7) {
for (i in 1:reps) y[i]<-mean((j*5)+sample(vector,plotsize))
yield<-append(yield,y1)
}
是否有办法获得更快执行的相同结果?