RでeCDFとすばやく統合する


10

Iは、フォームの積分方程式有する F nは経験的累積分布関数であり、gは関数です。私は収縮マッピングを持っているので、バナッハの固定小数点定理シーケンスを使用して積分方程式を解こうとしています。

T1(x)=0xg(T1(y)) dF^n(y)
F^ng

しかし、この非常にゆっくりとRで実行され、私は私がのために合計を使用して()関数を統合していますので、それは考えて何度も繰り返し。xF^n

経験的な分布を使用して、integrate()などの関数と統合するより速い方法はありますか?


6
これは実際には統計の質問ではなくRの質問です(したがって、おそらくstackoverflowに属しています)...コードを投稿できますか?Rでは、多くの場合、実行時のパフォーマンスを大幅に改善する複数の機会があり、コードを確認しないと、どれが適用されるかを判断するのが困難です。
jbowman 2013年

回答:


14

経験分布関数を定義 F NT = 1

F^n(t)=1ni=1nI[xi,)(t),
g(t)dF^n(t)=1ni=1ng(xi).
integrate()R
x <- rnorm(10^6)
g <- function(t) exp(t) # say
mean(g(x))

それはベクトル化されているので超高速でなければなりません。


経験的分布に関する関数の積分が、観測された点で評価された関数の平均である理由に関して、関連する質問を追加しました。math.stackexchange.com/questions/2340290/...
texmex
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.