我正在运行一个基于个体的模拟,其中每个个体都由一个随机变量mort_rate
控制。 N
的{{1}}数可以根据以下两个参数进行随机化:mort_rate
。在这里,出于说明目的,
mort_distr(N, parm1, parm2)
所以
mort_distr <- function(N, parm1, parm2)
{
v <- rnorm(n=N, mean=parm1)
mort <- v * exp(1./parm2)
data.frame(id=1:N,
mort=mort,
parm1=parm1,
parm2=parm2)
}
我的问题是:为了减少计算时间,我如何在R中预先生成一个3D查找表,其中包含{{1}的每个组合的> mort_distr(N=3, parm1=2., parm2=10.)
id mort parm1 parm2
1 1 1.908670 2 10
2 2 5.351502 2 10
3 3 2.440259 2 10
个N=500
样本}}和mort_rate
,分别从parm1
和parm2
开始变化?换句话说,我希望有一个表格seq(0.,4.,0.2)
,该表格使我能够通过运行例如seq(5.,10.,1.)