How to do curve fitting on data with inflection point?
6 次查看(过去 30 天)
显示 更早的评论
Hello
Dear experts,
I have a scatter plot of data of x and y and I need to fit a curve on them. However, it is obvious that there is no function that can fit the whole range of these data (blue line of the attached picture). So, I need to fit two curves (dashed red line): one is on the data that are less than the inflection point (yellow point), and one for data that are larger than the inflection point. The inflection point can visually be guessed as 0.25, but is there any way that I can calculate the inflection point from the original scatter plot?
Note that I know how to do the curve fitting, and my question is only regarding the inflection point.
Hope this is the appropriate forum to ask this question.
Thank you very much.
1 个评论
Adam Danz
2021-1-14
编辑:Adam Danz
2021-1-14
Why not fit the entire curve to a sigmoid and then extract values less than or greater than x=.25?
Alternatively, but inferior to the suggestion above, would be to separate the curves first and then fit each half to a sigmoid but I can't imagine it would be as good of a fit as fitting the full curve.
采纳的回答
Cris LaPierre
2021-1-14
编辑:Cris LaPierre
2021-1-14
Why do you think there is no equation that can fit this data? This looks like a sigmoid E max model to me. A similar shape is obtained from the Hill equation, too.
Here's a simple example using the Hill equation and trying to match your curve.
mx = 0;
mn = 1;
ec50 = 0.25;
n=5;
x=0:.01:1;
y = mn + (mx-mn)*ec50^n./(ec50^n+x.^n);
plot(x,y)
grid on
6 个评论
Cris LaPierre
2021-1-14
It is possible to numerically integrate it. I'll leave it to others to propose a closed-form solution.
ec50 = 0.25;
n=5;
x=0:.01:1;
y = 1-ec50^n./(ec50^n+x.^n);
A = trapz(x,y)
更多回答(0 个)
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!