MATLAB Answers

Fitting a sigmoid curve using curve fitting tool box

390 ビュー (過去 30 日間)
Hari krishnan
Hari krishnan 2020 年 2 月 17 日
編集済み: Srilatha Raghavan 2020 年 7 月 31 日
Hi, I am trying to fit a sigmoid function to the underlying data with the goodness of fit. I was using the curve fitting tool box. Unfortunately, i am not getting an idea to how to go forward with this. I tried to write custom function to fit the data, unfortunately, i am not able to get the required fit. Any help to solve this will be appreciated. X and Y data is shown below.
y= [0.70,1.09,0.96,0.37,0.27,0.35,0.16,0.18,0.18,0.02,0.37,0.67,1.26,1.70,1.85,1.73,1.92,1.72,1.47,1.50,1.46,1.43,2.12,1.6,1.8]
x = [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25]

採用された回答

Alex Sha
Alex Sha 2020 年 2 月 18 日
Hi, Hair, the quailty of your data seems to be bad for sigmoid type of function. If take the function like:
y = y0+a/(1+exp(-(x-x0)/b))^c
Then the fitting result looks like:
  3 件のコメント
Hari krishnan
Hari krishnan 2020 年 2 月 18 日
@ Alex Sha, the data corresponds to population dynamics. So the data captures the trend of a regular logistic growth. I wanted to know what the value of coefficient signifies?

サインインしてコメントする。

その他の回答 (1 件)

John D'Errico
John D'Errico 2020 年 2 月 18 日
編集済み: John D'Errico 2020 年 2 月 18 日
The problem is, there are many functional forms that people think of and call with the generic name "sigmoid". We can find a few of them listed here:
Perhaps the most common one is of the specific form:
f(x) = 1/(1 + exp(-x))
But many people will want to use thing like a transformed cumulative normal, which can be gotten from the erf function. Or possibly an arc tangent curve, or many others. Again, a few of them are shown on the wiki page. The thing is, each of those has a different fundamental shape.
y = [0.70,1.09,0.96,0.37,0.27,0.35,0.16,0.18,0.18,0.02,0.37,0.67,1.26,1.70,1.85,1.73,1.92,1.72,1.47,1.50,1.46,1.43,2.12,1.6,1.8]
x = 1:length(y);
plot(x,y,'o')
Your first problem is that your data is REALLY noisy. The first three points almost look like they are not even from the same curve.
Regardless, lets try to fit a simple sigmoidal model to your data. I'll start with the one I show above, adding some parameters, but keeping the same general form.
ft = fittype('L + (U-L)/(1 + exp(-4*log(3)*(x-xmid)/xscale80))','indep','x')
ft =
General model:
ft(L,U,xmid,xscale80,x) = L + (U-L)/(1 + exp(-4*log(3)*(x-xmid)/xscale80))
mdl = fit(x',y',ft,'start',[.5 ,3, 15,3])
mdl =
General model:
mdl(x) = L + (U-L)/(1 + exp(-4*log(3)*(x-xmid)/xscale80))
Coefficients (with 95% confidence bounds):
L = 0.4242 (0.2428, 0.6057)
U = 1.694 (1.519, 1.869)
xmid = 12.69 (11.91, 13.47)
xscale80 = 1.763 (-1.042, 4.567)
Here, I've named the parameters to be indicative of what they represent in the model. As it turns out, each of the paramters can be simply interpreted here.
That is, L is the lower asymptote, the limit of f(x), as x --> inf. U is the upper asymptote. xmid is the point where the curve goes through 50% of the transition. So xmid is the 50% point.
xscale80 is a little more difficult, but we can express it in the form of the width of the transition. If my arithmetic is correct, we can say the difference between the 10% and the 90% points on the curve will be xscale80. So 80% of the transition occurs between the 10% nd 90% points on that curve. Thus, when xscale80 is a large number, the transition is a slow one. If xscale is small, then the transition approaches a sharp step function.
plot(mdl)
grid on
hold on
plot(x,y,'bo')
xline(mdl.xmid);
xline(mdl.xmid + mdl.xscale80/2);
xline(mdl.xmid - mdl.xscale80/2);
I've drawn vertical lines at the 10% 50% and 90% points on the curve as it goes through the transition.
As you can see, the model does fit quasi-reasonably, although the lower asymptote is a bit off in my eyes. I did say your data was pure crapola, I think? At least those first three data points are dragging the estimate of L off from where it seems it should be. A weighted fit, or a robust regression would help.
mdl = fit(x',y',ft,'Robust','LAR','start',[.5 ,3, 15,3])
Fitting stopped because the number of iterations or function evaluations exceeded the specified maximum.
mdl =
General model:
mdl(x) = L + (U-L)/(1 + exp(-4*log(3)*(x-xmid)/xscale80))
Coefficients (with 95% confidence bounds):
L = 0.2993 (0.1037, 0.4948)
U = 1.72 (1.532, 1.908)
xmid = 12.59 (11.74, 13.44)
xscale80 = 2.374 (-0.7873, 5.536)
plot(mdl)
hold on
grid on
plot(x,y,'bo')
xline(mdl.xmid - mdl.xscale80/2);
xline(mdl.xmid);
xline(mdl.xmid + mdl.xscale80/2);
This fit looks more reasonable. Your data is still pure crapola though. This is hinted at by the width of the confidence intervals on the parameters. As you can see, those bounds on the parameters are all scarily wide.
  7 件のコメント
Srilatha Raghavan
Srilatha Raghavan 2020 年 7 月 31 日
It's the sigmoid function ( the equation) adjustment I am having troubles about. There's always a trade off between the sigmoid curve I got from your equation and my equation adjustment to get my curve to reach a particular height, at my rated x value.

サインインしてコメントする。

Community Treasure Hunt

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

Start Hunting!

Translated by