-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathscc_count_with_cov.m
89 lines (57 loc) · 2.46 KB
/
scc_count_with_cov.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
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
function [ class_id ] = scc_count_with_cov( X,Y,Z_cov,K )
[p1 n] = size(X);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
eucli_dis = norm(X - mean(X,2),'fro')^2;
logit_vec = log ( mean(Y,2)) ;
null_dev_Y = sum(sum( - Y .* logit_vec )) + n * sum( exp( logit_vec) ) + sum(Y .* log(Y+1e-5)) - sum(Y) ;
null_dev_Y = null_dev_Y * 2;
ratio = (eucli_dis) / ( (null_dev_Y) + (eucli_dis ));
fs_weight = [ones(1,p1)* ratio/ (p1) 1- ratio];
%%% Select K and phi in weights
[K_best,phi,w] = select_K_target_gower( [X;Y], fs_weight,@knn_weight_gower_weighted_dense,@knn_weight_gower_concat_2data_weighted_tune,1);
w = w ./ sum(w);
[U_output,Z_output,beta_output] = intclu_l2_poisll_CARP_SOK_with_cov(X,Y,Z_cov,w,K);
%% Rand Index
V_round = round(Z_output,3);
[class_id,iter_cut] = get_cluster_assignment(V_round,w,n,K);
len_V = size(V_round,3);
if length(unique(class_id)) == 1
class_id_mat = zeros(len_V,n);
class_no_vec = zeros(len_V,1);
for i = 1:len_V
[class_no_vec(i), class_id_mat(i,:)] = group_assign_vertice(V_round(:,:,i),w,n);
end
iter_cut = max(find(class_no_vec > 1));
class_id = class_id_mat(iter_cut,:);
end
%%%%%%%%%%%%%%%%% Adaptive
beta = beta_output(:,iter_cut);
theta_output_iter = U_output(end,:,iter_cut);
logit_cov = exp(Z_cov * beta) ;
%%% Adjust for alpha
Y_num = Y;
Y_num(Y_num == 0) = 0.1;
Y_new = Y_num ./ logit_cov';
Y_new(Y_new < 0) = 0;
logit_vec = log ( mean(Y,2)) ;
null_dev_Y = sum(sum( - Y_new .* logit_vec )) + n * sum( exp( logit_vec) ) + sum(Y_new .* log(Y+1e-5)) - sum(Y_new) ;
null_dev_Y = null_dev_Y * 2;
ratio = (eucli_dis) / ( (null_dev_Y) + (eucli_dis ));
fs_weight = [ones(1,p1)* ratio/ (p1) 1- ratio];
[K_best,phi,w] = select_K_target_gower( [X;Y ./ logit_cov' ], fs_weight,@knn_weight_gower_weighted_dense,@knn_weight_gower_concat_2data_weighted_tune,1);
w = w ./ sum(w);
[U_output,Z_output,beta_output] = intclu_l2_poisll_CARP_SOK_with_cov(X,Y,Z_cov,w,K);
%% Rand Index
V_round = round(Z_output,3);
[class_id,iter_cut] = get_cluster_assignment(V_round,w,n,K);
len_V = size(V_round,3);
if length(unique(class_id)) == 1
class_id_mat = zeros(len_V,n);
class_no_vec = zeros(len_V,1);
for i = 1:len_V
[class_no_vec(i), class_id_mat(i,:)] = group_assign_vertice(V_round(:,:,i),w,n);
end
iter_cut = max(find(class_no_vec > 1));
class_id = class_id_mat(iter_cut,:);
end
end