【发布时间】:2014-09-02 01:51:48
【问题描述】:
在进行重要性抽样实验时,我模拟了 Kolmogorov-Smirnov 距离的值
$$ D_n = \max_x |\hat{F}_n(x)-F(x)| $$
其中 $n$ 是原始重要性样本的大小,我想将这些值与 the asymptotic distribution of the Kolmogorov-Smirnov test 或 Kolmogorov distribution 进行比较,即
$$ \sqrt{N} D_n \longrightarrow \sup_{t\in[0,1]}|B(t)| $$
其中 $B$ 是布朗桥。
由于ks.test 依赖于这种渐近分布,它的 cdf 已经存在于 R 中的某个地方,我想知道如何访问它。 R 函数ks.test 包含指令
PVAL <- 1 - if (alternative == "two.sided")
.Call(C_pKolmogorov2x, STATISTIC, n)
但我自己拨打C_pKolmogorov2x 不起作用。
【问题讨论】:
-
在我看来,
kolmim包中的函数pkolm为pKolmogorov2x提供了一个包装器。 -
确实,这解决了我的问题!谢谢。
-
不客气。不能直接调用 C-Function 似乎很奇怪。只是用谷歌搜索了一下,但没有找到关于如何直接访问
pKolmogorov2x的解决方案。 -
确实很奇怪。如果我完全按照包中的方式重新定义
pkolm,它不会识别pKolmogorov2x... -
好的,知道了!尝试以下操作,例如:
.Call(stats:::C_pKolmogorov2x, STATISTIC = 0.04, n = 200),这对我来说产生了 0.1067137 美元。
标签: r