How to move fitted plot along y axis

Question

Niklas Kurz 2020-6-24

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/554149-how-to-move-fitted-plot-along-y-axis

评论： Niklas Kurz 2020-6-25

Hey, it's me again, this time I want to shift a fit I created with cftool:

[xData, yData] = prepareCurveData( Fqr, Prem );
% Set up fittype and options.
ft = fittype( 'a*atan(2*b*t/3.83^2-t^2)', 'independent', 't', 'dependent', 'y' );
opts = fitoptions( 'Method', 'NonlinearLeastSquares' );
opts.Display = 'Off';
opts.StartPoint = [0.0496544303257421 0.902716109915281];
% Fit model to data.
[fitresult, gof] = fit( xData, yData, ft, opts );
% Plot fit with data.
figure( 'Name', 'untitled fit 1' );
h = plot( fitresult, xData, yData,'x');
set(h,'LineWidth',2,'Markersize',7.5)
legend( h, 'Prem vs. Fqr', 'untitled fit 1', 'Location', 'NorthEast', 'Interpreter', 'none' );
% Label axes
xlabel( 'Fqr', 'Interpreter', 'none' );
ylabel( 'Prem', 'Interpreter', 'none' );
grid on

I manage to relcoate y values just by typing +pi/2 behind yData, but what about the fitted line? That's where I'm struggling.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Dana 2020-6-24

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/554149-how-to-move-fitted-plot-along-y-axis#answer_456643

在 MATLAB Online 中打开

Why not just shift up the raw y data?

[xData, tmp] = prepareCurveData( Fqr, Prem );
yData = tmp+pi/2; % shift up your y data values

Then everything else would be the same as you had.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Dana 2020-6-24

在 MATLAB Online 中打开

Okay, I see now the issue. Because you don't have a constant term in your regression specification, what I suggested before won't work. Before going any further, let me ask you: are you sure you don't want a constant term in your regression? Generally you'd need a compelling reason to not include one. It happens, but it's relatively uncommon.

If you did include one, you'd have something like

ft = fittype( 'c+a*atan(2*b*t/3.83^2-t^2)', 'independent', 't', 'dependent', 'y' );

where now c is the constant term. With that specification, you can then just add pi/2 to your yData as I suggested before, and it won't do anything except shift the data and fitted curve up by a constant pi/2.

If you're sure you don't want a constant term, then a slight modification of what I wrote before would work:

[xData, tmp] = prepareCurveData( Fqr, Prem );
yData = tmp+pi/2;
% Set up fittype and options.
ft = fittype( 'a*atan(2*b*t/3.83^2-t^2)+pi/2', 'independent', 't', 'dependent', 'y' );

So now pi/2 has been added to both the data and the regression specification. Everything after that would be the same as in your original post.

Dana 2020-6-24

When including the estimated constant c in your regression specification, you'll get the exact same shape of fitted curve whether or not you first add pi/2 to yData, the only difference is the curve will be shifted by a constant.

However, you won't generally get the same shape of fitted curve if you're comparing "with constant c, and with pi/2 added to yData" to "without constant c, and without pi/2 added to yData". In that case, only one of your specifications involves estimating a constant, so there's no reason to think you should get the same shape of curve. Hope that clarifies.

Niklas Kurz 2020-6-25