我使用IMIS
包(增量混合重要性抽样)来估计参数。不幸的是,它被写入寻找函数likelihood
,sample.prior
和prior
在它被调用的环境中,所以我不能将它包装在函数中(我的最终目标)。从?IMIS
单因素例如工作得很好,如何在R中填充环境
require(IMIS)
likelihood <- function(theta) exp(-1*sin(3*theta)*sin(theta^2) - 0.1*theta^2)
prior <- function(theta) dnorm(theta, 0, 5)
sample.prior <- function(n) rnorm(n, 0, 5)
result = IMIS(500, 3000, 100, 10)
## also fine using do.call (pertinent below)
result <- do.call(IMIS, args = list(B = 500, B.re = 3000, number_k = 100, D = 10))
,但勿庸置疑,在功能包装并不:
rm(likelihood, prior, sample.prior, result)
imisWrap <- function() {
likelihood <- function(theta) exp(-1*sin(3*theta)*sin(theta^2) - 0.1*theta^2)
prior <- function(theta) dnorm(theta, 0, 5)
sample.prior <- function(n) rnorm(n, 0, 5)
result = IMIS(500, 3000, 100, 10)
return(result)
}
imisWrap() ## can't find sample.prior
我认为的方式解决这个是我的包装,以创造环境(或使用其环境),然后使用do.call
在该环境中运行IMIS
,但我不知道如何创建一个新环境,该环境在i中有likelihood
,prior
,sample.prior
和result
吨。
编辑:使用@ BenBolker的很好的意见我有一个改善,但仍无法正常工作的尝试:
imisWrap2 <- function() {
likelihood <- function(theta) exp(-1*sin(3*theta)*sin(theta^2) - 0.1*theta^2)
prior <- function(theta) dnorm(theta, 0, 5)
sample.prior <- function(n) rnorm(n, 0, 5)
imisEnv <- new.env()
assign("likelihood", likelihood, envir = imisEnv)
assign("sample.prior", sample.prior, envir = imisEnv)
assign("prior", prior, envir = imisEnv)
result = do.call(IMIS,
args = list(B = 500, B.re = 3000, number_k = 100, D = 10),
envir = imisEnv)
return(result)
}
但是,这仍然无法找到的功能。
我想你可以使用'assign'来填充环境。你有没有考虑过让软件包作者改写一些更理智的东西......? –
@BenBolker或者只是使用'trace'实时重写函数。 –