How to integrate a distribution function for a specific interval?
信息
此问题已关闭。 请重新打开它进行编辑或回答。
显示 更早的评论
I computed two distribution function:
[distribution_1,x]=ecdf(X); Distribution_2= f(distribution_1);
Now I want to integrate Distribution_2 over the interval [a,b] and tried to use q = integral(fun,xmin,xmax) but I have trouble creating the function handle "fun" in my specific case (since I cannot create the connection to x?)
Thankful for any advice!
采纳的回答
Robert
2017-11-6
To make your empirical CDF into a function that integral can use, you should interpolate distribution_2. Because of the way the CDF is defined, I recommend you interpolate based on the previous value. The first value of x is a duplicate of the second, but for use with interp1(..., 'previous') it should be a (finite) number less than any you plan to use. -realmax should do the trick.
x(1) = -realmax;
f_dist2 = @(new_x) interp1(x, distribution_2, new_x, 'previous', 'extrap');
integral(f_dist2, -1, 1)
3 个评论
Robert
2017-11-6
Here is the code I used for testing:
X = randn(100, 1); % sample data
[distribution_1, x] = ecdf(X);
g=0.1;
h=0.2;
distribution_2 = 1-(1-distribution_1.^(1/(1+g))).^(1+h);
cla
stairs(x, distribution_1)
hold all
stairs(x, distribution_2)
hold off
x(1) = -realmax;
f_dist2 = @(new_x) interp1(x, distribution_2, new_x, 'previous', 'extrap');
f_dist2(0)
f_dist2([-10, 10])
integral(f_dist2, -1, 1)
integral(f_dist2, -10, 10)
Nina
2017-11-7
Nina
2017-11-7
更多回答（1 个）
David Goodmanson
2017-11-6
编辑：David Goodmanson
2017-11-6
Hi Nina,
If this is not a misinterpretation of the problem, then it seems to work. I am assuming that the function ecdf is already on the Matlab path.
g = .1; h = .2;
fun = @(x) 1-(1-ecdf(x).^(1/(1+g))).^(1+h)
integral(@(x) fun(x), 2,3)
function a = ecdf(x) % for demo purposes
a = sin(x);
end
I'm curious about the idea of integrating a cdf, could you perhaps comment on that?
1 个评论
Nina
2017-11-7
此问题已关闭。
Translated by 
