-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcorrMatMeasures.m
69 lines (54 loc) · 1.82 KB
/
corrMatMeasures.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
% Measurements on the correlation matrices
clear all;
% load('files\newGeneMat_5RPKM.mat');
load('files\asdGeneIDs_5RPKM.mat');
load('files\szGeneIDs_5RPKM.mat');
load('files\ndGeneIDs_5RPKM.mat');
load('files\hkGeneIDs_5RPKM.mat');
% for i = 1 : 4
%
% % ASDs
% load(['results\diffNetworks\abs_abs\asdDiffNet_' num2str(i+1) num2str(i) '.mat']);
% temp1 = triu(asdDiffMat,1);
% nonZ = find(temp1 ~= 0);
% temp2 = temp1(nonZ);
% N = isnan(temp2);
% NN = find(N == 1);
% temp2(NN) = [];
% asdM(i) = mean(temp2);
% clear asdDiffMat; clear temp1; clear temp2; clear nonZ;
%
% % SZ
% load(['results\diffNetworks\abs_abs\szDiffNet_' num2str(i+1) num2str(i) '.mat']);
% temp1 = triu(szDiffMat,1);
% nonZ = find(temp1 ~= 0);
% temp2 = temp1(nonZ);
% N = isnan(temp2);
% NN = find(N == 1);
% temp2(NN) = [];
% szM(i) = mean(temp2);
% clear szDiffMat; clear temp1; clear temp2; clear nonZ;
%
% % ND
% load(['results\diffNetworks\abs_abs\ndDiffNet_' num2str(i+1) num2str(i) '.mat']);
% temp1 = triu(ndDiffMat,1);
% nonZ = find(temp1 ~= 0);
% temp2 = temp1(nonZ);
% N = isnan(temp2);
% NN = find(N == 1);
% temp2(NN) = [];
% ndM(i) = mean(temp2);
% clear ndDiffMat; clear temp1; clear temp2; clear nonZ;
%
% end
%%%------------------------------------------------------------------------
%%%------------------------------------------------------------------------
%%% Hirarchial Clustering
%%%------------------------------------------------------------------------
% load(['results\geneCorr\Avg(log2)\ageGroup1_corrMat.mat']);
%
% Z = linkage(corrMat);
% save('results\geneCorr\Avg(log2)\ageGroup1_Z.mat', 'Z');
load('results\geneCorr\Avg(log2)\ageGroup1_Z.mat');
% set(0,'RecursionLimit',1000);
dendrogram(zNEW)