Hi all,
I would like to calculate the correlation between two vectors (1*744) for each. but the correlation output gives a matrix (744*744) contains NaNs with no number in the matrix.
Any help will be appreciated.
thanks

 채택된 답변

Birdman
Birdman 2017년 12월 20일

0 개 추천

Try corrcoef
corrcoef(A,B)

댓글 수: 12

Lilya
Lilya 2017년 12월 20일
the obtained result is:
NaN NaN NaN 1
Birdman
Birdman 2017년 12월 20일
Can you share your vectors in a mat file?
Lilya
Lilya 2017년 12월 20일
Find the attachment, please. thank you
If we want to omit NaN values from your data, the two array becomes differently sized. Therefore, a solution could be replacing NaN values with the mean of the array. Only JulTotcurr include NaN values.
JulTotcurr(isnan(JulTotcurr))=mean(JulTotcurr,'omitnan');
corrcoef(JulTotcurr,JulTotwnd)
Lilya
Lilya 2017년 12월 20일
편집: Lilya 2017년 12월 20일
Great! it works But can we obtain it as a vector? to plot the data?
Birdman
Birdman 2017년 12월 20일
What do you want to plot?
Lilya
Lilya 2017년 12월 20일
the values of the correlation
res=corrcoef(JulTotcurr,JulTotwnd);
plot(res(:,1),res(:,2),'o')
Lilya
Lilya 2017년 12월 20일
Totally appreciate it thank you
Jos (10584)
Jos (10584) 2017년 12월 20일
Statistically, it is not recommended to replace the missing values by the mean ("mean imputing"). There are other techniques that are (often) more reliable.
Statistical analyses with missing data is quite common, and there are many pitfalls to be aware of, in order to avoid drawing false conclusions!
Lilya
Lilya 2017년 12월 21일
Thank you for your clarification Jos. Do you recommend one of them?
My pleasure. You can use random sampling, for instance. Replace the NaNs by values sampled from the rest of the data vector. The advantage is that it, in the long run, preserves both the mean and the variance of the original data (whereas mean imputing only preserves the mean). You also might consider bootstrapping this procedure to get a better estimate of the correlation with another vector. If X is your original vector:
Xcopy = X ;
tf = isnan(X) ;
Xcopy(tf) = randsample(X(~tf), sum(tf), true)
You also should analyse whether your missing data is (completely) random or not, as this would influence the inferences you can make.

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

도움말 센터File Exchange에서 NaNs에 대해 자세히 알아보기

질문:

2017년 12월 20일

댓글:

2017년 12월 21일

Community Treasure Hunt

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

Start Hunting!

Translated by