10

我正在尝试编写一个通用函数来计算 R 中二项式比例的区间估计的覆盖概率。我打算为各种置信区间方法执行此操作,例如 Wald、Clopper-Pearson、HPD 区间用于不同的先验。

理想情况下,我希望有一个函数,可以将用于计算间隔的方法作为参数。那么我的问题是:如何将一个函数作为参数包含在另一个函数中?

例如,对于 Exact Clopper-Pearson 区间,我有以下函数:

# Coverage for Exact interval
ExactCoverage <- function(n) {
p <- seq(0,1,.001)
x <- 0:n

# value of dist
dist <- sapply(p, dbinom, size=n, x=x)

# interval
int <- Exact(x,n)

# indicator function
ind <- sapply(p, function(x) cbind(int[,1] <= x & int[,2] >= x))

list(coverage = apply(ind*dist, 2, sum), p = p)
}

其中 Exact(x,n) 只是一个计算适当间隔的函数。我想拥有

Coverage <- function(n, FUN, ...)
...
# interval
int <- FUN(...)

这样我就有一个函数来计算覆盖概率,而不是为每种间隔计算方法提供一个单独的覆盖函数。有没有标准的方法来做到这一点?我一直无法找到解释。

谢谢,詹姆斯

4

3 回答 3

16

在 R 中,函数可以作为函数参数提供。语法与非函数对象之一匹配。

这是一个示例函数。

myfun <- function(x, FUN) {
  FUN(x)
}

此函数将函数应用于FUN对象x

一些包含从 1 到 10 的数字的向量示例:

vec <- 1:10 

> myfun(vec, mean)
[1] 5.5
> myfun(vec, sum)
[1] 55
> myfun(vec, diff)
[1] 1 1 1 1 1 1 1 1 1

这不仅限于内置函数,还适用于任何函数:

> myfun(vec, function(obj) sum(obj) / length(obj))
[1] 5.5

mymean <- function(obj){
  sum(obj) / length(obj)
}
> myfun(vec, mymean)
[1] 5.5
于 2013-03-22T09:47:16.710 回答
2

您还可以将函数名称存储为字符变量,并使用 do.call() 调用它

> test = c(1:5)
> do.call(mean, list(test))
[1] 3
> 
> func = 'mean'
> do.call(func, list(test))
[1] 3
于 2015-11-22T06:02:29.940 回答
0

哈德利的文字给出了一个很好的(和简单的)例子:

randomise <- function(f) f(runif(1e3))

randomise(mean)
#> [1] 0.5059199
randomise(mean)
#> [1] 0.5029048
randomise(sum)
#> [1] 504.245
于 2020-07-03T10:33:11.763 回答