Fitting function to ecdf plot

2020 3 23
1 个回答
回答已采纳
更新时间：2020 3 24
5 次查看（30 天）

0 个投票

Hi
I am looking for a way to fit function g to the ecdf-plot I produce as follows.
T = 1000;
cens = (y>T);
ecdf(y,'Function','survivor', ...
'Censoring',cens,'Bounds','on');
g=fittype('a-b*exp(-c*x)');
I would like to use the bounds determined in ecdf as a weighting factor "w" for the fit.
I tried
h=gcf;
dataObjs = findobj(gcf,'-property','YData');
dmy3=dataObjs.XData;
x=dmy3';
dmy4=dataObjs.YData;
y2=dmy4';
g=fittype('a-b*exp(-c*x)');
f=fit(x,y2,g,'StartPoint',[[ones(size(x)), -exp(-x)]\y2; 1]);
hold on;
Though this is giving me the upper bounds as xData and yData.
Thanks in advance for any suggestion.
Peter

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

 采纳的回答

Peter Mühlenbrock
Peter Mühlenbrock 2020-3-24

0 个投票

Iw ill close this, as I found a solution.
%----access all xy data-----
d = findall(gcf, '-property', 'xdata');
xydatas = arrayfun(@(h) get(h, {'xdata','ydata', 'type'}), d, 'Uniform', 0);
%----adress xy data-----
x1=xydatas{3, 1}{1, 1}' ;
y1=xydatas{3, 1}{1, 2}' ;
w=xydatas{3, 1}{1, 2}'-xydatas{2, 1}{1, 2}';
g=fittype('a*exp(-b*x)+c');
f=fit(x1,y1,g,'StartPoint',[1, 0.02, 0.5],'Weight',w);
hold on;
plot(f)

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 R Language 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by