0

我想进行一个模拟研究,我改变我的输入参数。以下是我的函数的输入:

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] 中的所有值

4

0 回答 0