How to compare the standard errors

1 次查看(过去 30 天)
Lizan
Lizan 2012-11-22
Hi,
I have implemented a code for the slope,
x_ceng = sum(xi)/n; % center of gravity of x
m = sum(xi-x_ceng)*yi/sum(xi-x_ceng)^2; % Slope
di = yi-m*xi-c; % the residuals
D = sum(xi-x_ceng)^2;
% The standard error is then given by;
dm = sqrt(sum(di^2)/D*(n-2))
I have compared m with values of polyfit and it agrees, but I am not certain that the standard error is correct. I am obtaining very low dm that what I would like considering the spread of the data that I obtain.
Is there any function in MATLAB that I can compare dm with?
Many Thanks for any input.
  3 个评论
Lizan
Lizan 2012-11-23
I am curious what is it in eq."dm" that does not look correct for you?
Star Strider
Star Strider 2012-11-23
Initially, there was a right parenthesis missing. (You corrected that in your edit.)
The polyfit documentation outlines a procedure to calculate the covariance matrix from information in the ‘S’ structure, and from that you can calculate the standard errors to check to see how they compare to those you calculated in your code.

请先登录,再进行评论。

回答(1 个)

Star Strider
Star Strider 2012-11-23
If you want to determine the parameter confidence intervals from polyfit, I suggest the File Exchange function polyparci.

产品

Community Treasure Hunt

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

Start Hunting!

Translated by