-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathotherdatasets.m
24 lines (24 loc) · 927 Bytes
/
otherdatasets.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
for itr = 2:13
currentdata = alldata{itr};
currentdata = FilterData(currentdata);
[pxx, f, pcascore,explained] = transformdata(currentdata, 2,5);
[idx,P,BIC,gm] = GMMCluster(pcascore(:,1:20),8,0,10);
d2 = mahal(gm,pcascore(:,1:20));
clusters = NewClustering(idx,8,currentdata);
Y = tsne(pxx');
RecreateKmeans(idx,8,pxx,d2,clusters, f, Y, 'otherdatasets/',itr);
close;
end
[pxx, f, pcascore,explained] = transformdata(data, 2,5);
[idx,P,BIC,gm] = GMMCluster(pcascore(:,1:20),8,0,10);
for itr = 1:13
currentdata = alldata{itr};
currentdata = FilterData(currentdata);
[pxx, f, pcascore,explained] = transformdata(currentdata, 2,5);
[idx,P] = GMMPredict(pcascore(:,1:20),gm);
d2 = mahal(gm,pcascore(:,1:20));
clusters = NewClustering(idx,8,currentdata);
Y = tsne(pxx');
RecreateKmeans(idx,8,pxx,d2,clusters, f, Y, 'otherdatasetsinitialmodel/',itr);
close;
end