Rの二項比の区間推定のカバレッジ確率を計算するための一般関数を記述しようとしています。 Wald、Clopper-Pearson、様々なプライオリティのHPD間隔。R関数の引数としての関数パラメータ
理想的には、間隔を計算するために使用するメソッドを引数として取ることができる関数が1つあります。次に私の質問:別の関数に関数を引数として含めるにはどうしたらいいですか?正確な(X、n)は、適切な間隔を計算するだけの関数である
# 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)
}
:例として
は、正確なClopperピアソン間隔でIは、以下の機能を有します。私はカバレッジ確率ではなく間隔計算の各方法のための別々のカバレッジ関数を計算するための機能を有するように
Coverage <- function(n, FUN, ...)
...
# interval
int <- FUN(...)
したいと思います。これを行う標準的な方法はありますか?私は説明を見つけることができませんでした。
おかげで、Rは ジェームズ
https://github.com/hadle y/devtools/wiki/Functional-programming – hadley
役に立つリンクをありがとう。 – jatotterdell