我想进行一个模拟研究,我改变我的输入参数。以下是我的函数的输入:
my.function<-function(rho, psi, beta0 = c(50,50), beta1 = c(1,1), sigma=c(1,1)) {
c(rho,psi,beta0,beta1,sigma)
}
我尝试使用 expand.grid 和 mapply 对各种输入值进行模拟
rho= c(0,.5,.99)
psi=c(0,.5,.99)
beta1.1=c(0,.5,1)
beta1.2=c(0,.5,1)
sigma.1= c(0,1)
sigma.2=c(0,1)
my.grid<- expand.grid(rho,psi,beta1.1,beta1.2,sigma.1,sigma.2)
然后我尝试使用mapply:
mapply(my.function, x = my.grid[,1], y=my.grid[,2],
beta1=c(my.grid[,3],my.grid[,4]), sigma= c(my.grid[,5], my.grid[,6]))
这里的问题在于 beta1 项和 sigma 项。而不是治疗
beta1=c(my.grid[,3],my.grid[,4])
sigma= c(my.grid[,5], my.grid[,6]))
作为一对,当我使用 mapply 时,该函数将运行 my.grid[,3] 和 my.grid[,4] 中的所有值