How do I correlate columns

15 次查看（过去 30 天）
Kris
Kris 2015-4-13
回答： Jonathan Campelli 2015-4-14
I have 2400 x 1015 matrix. I would like to correlate column 1 with column 926. I would then like to correlate column 1 with 927 and so on until I correlate column 1 with column 1015.
I would then like to print the result of each correlation in row 1, in a new file.
Thanks
  1 个评论
Jonathan Campelli
Jonathan Campelli 2015-4-14
I've made a program that's computing each correlation by the formula, which is returning totally different results from my previous code. Just updating you that I'll post the correct code within the next couple hours, in case you have need of it still. "corrcoeff" was yielding different results than those you asked for. Coming soon...
Jonathan Campelli

请先登录，再进行评论。

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

采纳的回答

Jonathan Campelli
Jonathan Campelli 2015-4-14
Hello again, Kris,
This function outputs a 925x90 matrix with individually computed linear correlation coefficients for each of the comparisons you described. It took me a short while to determine precisely what you were seeking to compute and what format you intended:
function [ output_args ] = KrisCorrelation( matrix )
n=size(matrix,1); %Determine number of rows with "size" and"1"
standard=zscore(matrix); %Standardize every column in matrix with "zscore"
x=1:925; %Set first range of columns for correlation
y=926:1015; %Set second range of columns for correlation
%%Use matrix multiplication to multiply column 1 by columns 926:1015, then
%%column 2 by columns 926:1015, and so on until all multiplications are
%%complete. Matrix multiplication sums the values in every column
%%automatically, and this time it leaves one row vector of length 90 for
%%each of the 925 separate loops. This equation completes the equations
%%for finding the correlation coefficient:
coefficients=(standard(:,x)'*standard(:,y))/(n-1);
%Save correlation variable to MAT file:
save('KrisCoefficients', 'coefficients')
%Save correlation variable "correlationResults" to XLSX file:
xlswrite('KrisCoefficients.xlsx', coefficients)
end
Once you hit "Accept Answer" on the solution that meets your needs, we'll know your all set.
Best regards,
Jonathan Campelli

更多回答（2 个）

Chad Greene
Chad Greene 2015-4-13
How's this?
m = rand(2400,1015);
c = corrcoef(m(:,[1 926:1015]));
imagesc(c)
cb = colorbar;
ylabel(cb,'correlation coefficient')
  4 个评论
显示 2更早的评论
Kris
Kris 2015-4-13
I want to know the strength of the relationship between the values in column 1 compared with the values in columns 926, I then want to know the strength of the relationship of the values in column 1 compared with values in column 927 and so on until I reach column 1015. I then want to repeat this operation for columns 1 - 925.
Jonathan Campelli
Jonathan Campelli 2015-4-13
Alright, cooking up some code...

请先登录，再进行评论。

Chad Greene
Chad Greene 2015-4-13
The rand function was only to generate random data to play with.
I don't see how you'll get a 925 x 90 correlation matrix. Here I'll make up some data, but set column 1 to sine wave, then columns 926 to the end will have an increasingly high sine wave to random noise ratio. So correlation increases from columns 926 to 1015, as is evident in the top row of subplot 3:
% Some random data:
m = rand(2400,1015);
% The first column is sine function:
m(:,1) = sind(1:2400);
% Force correlation to come increasingly positive
% from cols 926 to 1015,
mag = 0;
for k = 926:1015
m(:,k) = m(:,k) + mag*sind(1:2400)';
mag = mag+.01;
end
subplot(131)
imagesc(m)
title 'full dataset'
subplot(132)
imagesc(corrcoef(m))
title 'correlation all data'
caxis([-1 1])
subplot(133)
c = corrcoef(m(:,[1 926:1015]));
imagesc(c)
title 'correlation of cols 1, 926:1015'
caxis([-1 1])
rgbmap('red','white','blue')

类别

在 Help Center 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by