How can I get function expression for Cumulative distribution function

2019 7 29
2 个回答
回答已采纳
更新时间：2019 7 30
5 次查看（30 天）

0 个投票

I have a cumulative distribution function(cdf). I have attached the figure. How can I fit the curve and get the function
expression for it?

请先登录，再回答此问题。

 采纳的回答

Star Strider
Star Strider 2019-7-29

0 个投票

Likely the easiest way to estimate the parameters is with the Statistics and Machine Learning Toolbox fitdist function:
x = randn(1, 25)*2 + 5; % Create Data (Mean = 5, St Dev = 2)
b = fitdist(x(:), 'normal')
producing:
b =
NormalDistribution
Normal distribution
mu = 5.30768 [4.42359, 6.19177]
sigma = 2.1418 [1.67238, 2.97957]
Choose the appropriate distribution function for your data.

6 个评论
显示 4更早的评论 隐藏 4更早的评论

Miley Lee
Miley Lee 2019-7-29
hi, thank you for answer my question.
Actually I already have these mu and sigma value. I would like to derive the function expression for the CDF figure. I have the all the x,y value on the figure.
Star Strider
Star Strider 2019-7-29
If you have the mu and sigma values, use the cdf function.
An alternate option using erf:
p = @(x,mu,sigma) (1 + erf((x-mu)/(sigma*sqrt(2))))/2; % Equivalent to ‘normcdf’
Miley Lee
Miley Lee 2019-7-29
hello,thank you for your answering.
I just type this line in MATLAB and use
plot(p);
to compare with the cdf figure see whether it is fit. But MATLAB can not plot p.
Below is the code I use to plot the above cdf figure. x is a 8760*1 matrix.
[mu,sigma] = normfit(x);
d=cdf('norm',x,mu,sigma);
figure
plot(x,d)
Thank you!
Star Strider
Star Strider 2019-7-29
You probably have to sort ‘x’ first. (Sorting it has no effect on the parameter estimates, and makes a decent-looking plot possible.)
Try this:
x = sort(x);
[mu,sigma] = normfit(x);
d=cdf('norm',x,mu,sigma);
figure
plot(x,d)
Miley Lee
Miley Lee 2019-7-29
hello, I have try the below code:
x=sort(y,'descend');
[mu,sigma] = normfit(x);
d=cdf('norm',x,mu,sigma);
figure
plot(x,d)
hold on;
p = @(x,mu,sigma) (1 + erf((x-mu)/(sigma*sqrt(2))))/2;
plot(p)
However, it still said plot(p) can not achieved. May I ask what is the "erf" mean?
Star Strider
Star Strider 2019-7-29
I provided you with a link to the erf function. It is the error function.
I provided you with the correct code for the plot. Use it!

请先登录，再进行评论。

更多回答（1 个）

Miley Lee
Miley Lee 2019-7-30

0 个投票

hello,sorry for the late reply.
Can I ask where is he correct code for the plot? I haven't find in the answer.
ask.PNG

1 个评论
显示 -1更早的评论 隐藏 -1更早的评论

Star Strider
Star Strider 2019-7-30
It is in this Comment.
Specifically:
x = randn(1, 25)*2 + 5; % Create Data
x = sort(x);
[mu,sigma] = normfit(x);
d=cdf('norm',x,mu,sigma);
figure
plot(x,d)
I did not initially provide the ‘x’ vector, since you have your own. I created it to test my code, and include it here for completeness. (The sort function sorts ascending by default, so there is no need to specify it.)
An alternative using the erf function is:
x = randn(1, 25)*2 + 5; % Create Data
x = sort(x);
[mu,sigma] = normfit(x);
p = @(x,mu,sigma) (1 + erf((x-mu)/(sigma*sqrt(2))))/2; % Equivalent to ‘normcdf’
figure
plot(x,p(x,mu,sigma))
using the approach in Find Cumulative Distribution Function of Normal Distribution.

请先登录，再进行评论。

类别

在 帮助中心 和 File Exchange 中查找有关 Half-Normal Distribution 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by