Difference between spearmann correlation and fit command

1 次查看(过去 30 天)
Hello all,
I have some data that I want to analyze (EMG of muscles vs subjects' fatigue) in order to find the relation between two measurs.
At first, I used the 'corr' command to calculate the spearman correlation coefficient between the EMG and the fatigue for each subject. Afterwards, I used the 'fit' command to fit a linear model to my data.
To my surprise, in some cases (about 5% of the cases), the spearman correlation coefficient was positive, but the fitted linear model had a negative slope (or vice versa).
Shouldn't they have the same tendency? Or am I doing something wrong?
I'm attaching my matlab code, and excel file with some example data that shows the spoken discrepancy.
Thank you.

采纳的回答

Peter Perkins
Peter Perkins 2013-10-11
You're (correctly) fitting a linear model with a slope and an intercept. The "model" that a correlation coefficient is equivalent to has no intercept. Look at your figure to see why one is positive and the other negative.
Hope this helps.
  1 个评论
Gena
Gena 2013-10-13
Thank you for your answer, but there is still something that I don't understand. When you are saying that "the model has no intercept" you mean that the model is y=a*x (passes through the origin), if I understand it correctly. But what about other cases (as in the attached file) where the model of the correlation coefficient doesn't pass through the origin (and has a negative value, so as the slope of the fitted line)? What is the difference (after all, in both cases I'm using the same command - "corr")?
Thank you.

请先登录,再进行评论。

更多回答(1 个)

Wayne King
Wayne King 2013-10-11
I would not put any faith in the correlation coefficient, the p-value is 0.9 and the reported rho is close to 0.
  1 个评论
Gena
Gena 2013-10-11
I agree, but that is not the issue. There are some other cases where the correlation is more significant, but the discrepancy still exists.
Thank you for the answer.

请先登录,再进行评论。

类别

Help CenterFile Exchange 中查找有关 Get Started with Curve Fitting Toolbox 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by