我为一些微分方程(deSolve)创建了一个数组,并希望将变量(n = 1000)的几个主题(n = 250)的值放入数组中,但是现在我使用两个for循环和它需要大量的运行时间。我不知道如何改善我的代码。有什么建议吗?
array.fruit <- array(NA, dim=c(X+1, length(yini)+1, 250, 1000),
dimnames = list(NULL, c("time","Glycogen", "FatMass", "LeanMass", "ECF", "BW", "ATh", "BMI"))) #here I created an empty array
for (o in 1:250){
for (i in 1:1000){
array.fruit[,,o,i] <- ode(func=model.fruit, y=yini, parms=pars, times=times1)
}
}