关于corr函数

2019-07-17 13:58发布


这是小弟编写的一个简单的程序,想得到12153*12153的矩阵,这里是想把原先ppi那个矩阵里面为1的数据利用gene中的数据计算出相关系数替换掉,但是实际运行时不知道为什么没有替换掉,下面是实际的程序
close all;clear all;clc;
load('C:/Users/Administrator/Desktop/ppi_mat_tempo.mat');

load('C:/Users/Administrator/Desktop/gene_microarray_sample.mat');

for c=1:12153
   for d=1:12153
      if ppi_mat_tempo(c,d)==1
       ppi_mat_tempo(c,d)=corr(gene_microarray_sample(c,:)',gene_microarray_sample(d,:)');
      else ppi_mat_tempo(c,d)=0;

      end


   end

end
开始我想是不是编写错误,于是又编了一个小的程序如下:
close all;clear all;clc;

B=[0 1 1;1 0 1;1 1 0];
C=[1 23 3;9 5 6;68 8 9];
for i=1:3
   for j=1:3
      if B(i,j)==1
       B(i,j)=corr(C(i,:)',C(j,:)');
      else
          B(i,j)=0;

      end


   end

end
运行后可以实现替换B矩阵中1的功能,我不知道为什么会这样,是不是因为数据太大的原因(122153*12153)?这个问题又该如何解决呢?求高手指点~~
友情提示: 此问题已得到解决,问题已经关闭,关闭后问题禁止继续编辑,回答。