Asked by Ganesh
on 1 Oct 2012

I have 2 sets of data of the same size (X1, Y1) and (X2,Y2). (X1,Y1) has a linear relationship with zero intercept (Y1=slope*X1). How can I estimate the slope, R^2 value and plot it along with 95% confidence interval lines? Using Corrcoeff, polyfit and polyval gives me a relationship of the form Y1=slope*X1 + c, but I would like to force c=0.

(X2,Y2) has a non-linear relationship of the form Y2=A-Bexp(C*X2). Likewise, I would like to determine A,B & C along with R^2 value and plot it along with 95% confidence interval lines. I know I have to use nlinfit, nlparci but I am new in this and would need your help.

Thanks much.

*No products are associated with this question.*

Answer by Star Strider
on 1 Oct 2012

Edited by Star Strider
on 1 Oct 2012

Accepted Answer

You seem to be acquainted with `nlinfit` and know how to use it, but apparently not with `Anonymous Functions`.

I suggest you use `nlinfit` for both problems. With `B = [B(1)` ... `B(n)]'` a column vector of parameters, your `Y1` function becomes:

Y1 = @(B,x) B.*x;

since with only one parameter in it, you do not need to use subscripts.

Your `Y2` function becomes:

Y2 = @(B,x) B(1)-B(2).*exp(B(3).*x);

These should work as your `modelfun` functions in `nlinfit`.

Ganesh
on 5 Oct 2012

Thank you for providing with the code. It ran smoothly. The only change I made is in the plot command. F2_est is a row array, while FitCI95 is a column array. So, I made a minor change as: plot(X2fit, F2_est'-FitCI95+0.01, ':r', 'LineWidth', 1.5) plot(X2fit, F2_est'+FitCI95, ':r', 'LineWidth', 1.5)

Thanks a ton for your help.

Ganesh
on 5 Oct 2012

Is there a way to find the R^2 value for the fit?

Star Strider
on 5 Oct 2012

** R^2:** Yes, if you use the

It is my pleasure to help!

Log in to comment.

Related Content

MATLAB and Simulink resources for Arduino, LEGO, and Raspberry Pi

Learn moreOpportunities for recent engineering grads.

Apply Today
## 0 Comments

Log in to comment.