在R中快速与eCDF集成


10

我有一个形式为的积分方程, 其中是经验cdf,是一个函数。我有一个压缩映射,所以我尝试使用Banach不动点定理序列来求解积分方程。˚F Ñ

Ť1个X=0XGŤ1个ÿ dF^ñÿ
F^ñG

但是,这在R中运行非常缓慢,我想这是因为我一次又一次地对使用sum()函数进行积分。XF^ñ

有没有一种更快的方法可以将经验分布与诸如Integrated()之类的函数结合使用?


6
尽管这实际上是一个R问题而不是一个统计问题(因此可能属于stackoverflow)...您可以发布代码吗?在R中,通常有很多机会可以提高运行时性能,而又看不到代码,很难确定哪一个适用。
jbowman

Answers:


14

限定经验分布函数 它遵循 因此,您无需使用即可解决此问题。这种代码∫ -

F^ñŤ=1个ñ一世=1个ñ一世[X一世Ť
-GŤdF^ñŤ=1个ñ一世=1个ñGX一世
integrate()R
x <- rnorm(10^6)
g <- function(t) exp(t) # say
mean(g(x))

应该是超快速的,因为它是矢量化的。


请注意,我添加了一个相关的问题,即为什么函数相对于经验分布的积分是在观察点处评估的函数的平均值。math.stackexchange.com/questions/2340290/...
的TexMex
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.