r - 使用 for 循环将 Weibull 观测值保存在矩阵中

Question

我编写了代码来使用 for 循环基于 8 个不同的形状参数来模拟 Weibull 值。我想将它保存在 365 x 8 的矩阵中并最终导出到 csv。

我的代码是这样的：

k <- seq(1.4,2.8,0.2)
test <- matrix(rep(NA,365),nrow=365,ncol=8)
for (i in k){
  wind.k <- mapply(rweibull,i,scale.k,MoreArgs=list(n=1))
  test[ ]<-as.matrix(wind.k)
}

结果是 365 x8，但我注意到它只为每列捕获 k=2。任何人都可以帮助修复代码吗？

非常感谢您提前提供的帮助！

Answer 1

Use this:

test <- matrix(rweibull(n=365*8, shape=k, scale=scale.k), nrow=365, ncol=8, byrow=TRUE)

Note that scale.k must be a length 8 vector.