-
Notifications
You must be signed in to change notification settings - Fork 3
/
Copy pathplotROCtehnicalNote.m
executable file
·70 lines (61 loc) · 2.21 KB
/
plotROCtehnicalNote.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
clc; clear all; close all;
load singleROCtechnicalNote.mat
lw = 3; fontsize = 15;
figure;
plot(FPcg,TPcg,'lineWidth',lw);
hold on;
plot(FPsample,TPsample,'lineWidth',lw,'Color','r');
set(gca,'Ytick',[0:0.1:1]);
set(gca,'YTickLabel',[0:0.1:1],'FontSize',fontsize);
set(gca,'Xtick',[0:0.1:1]);
set(gca,'XTickLabel',[0:0.1:1],'FontSize',fontsize);
xlabel('False positive rate','FontSize',20)
ylabel('True positive rate','FontSize',20);
h = legend('$\hat{\mbox{PLV}}_{\mbox{circgauss}}$','$\hat{\mbox{PLV}}_{\mbox{sample}}$');
set(h,'Interpreter','latex','fontsize',13);
axis([0 .5 0.5 1])
load ROCtechnicalNote.mat
% auc : nCoupling x nSigma x nLength
% default values for
% coupling = 0.15, sigma = 0.5, n = 5000
% function of epsilon
cg = squeeze(auc_cg(:,3,3));
sample = squeeze(auc_sample(:,3,3));
aucCompare = [cg sample];
figure; bar(aucCompare,'grouped');
h = legend('$\hat{\mbox{PLV}}_{\mbox{circgauss}}$','$\hat{\mbox{PLV}}_{\mbox{sample}}$');
set(h,'Interpreter','latex','fontsize',13);
axis([0 6 0.5 1.2])
set(gca,'XTickLabel',couplingVec,'FontSize',fontsize);
set(gca,'Ytick',[0:0.5:1]);
set(gca,'YTickLabel',[0:0.5:1],'FontSize',fontsize);
xlabel('\epsilon','FontSize',20)
ylabel('AUC')
% function of sigma
cg = squeeze(auc_cg(3,:,3));
sample = squeeze(auc_sample(3,:,3));
aucCompare = [cg' sample'];
figure; bar(aucCompare,'grouped');
h = legend('$\hat{\mbox{PLV}}_{\mbox{circgauss}}$','$\hat{\mbox{PLV}}_{\mbox{sample}}$');
set(h,'Interpreter','latex','fontsize',13);
axis([0 6 0.5 1.2])
fontsize = 15;
set(gca,'XTickLabel',sigmaVec,'FontSize',fontsize);
set(gca,'Ytick',[0:0.5:1]);
set(gca,'YTickLabel',[0:0.5:1],'FontSize',fontsize);
xlabel('\sigma','FontSize',20)
ylabel('AUC')
% function of sigma
cg = squeeze(auc_cg(3,3,:));
sample = squeeze(auc_sample(3,3,:));
aucCompare = [cg sample];
figure; bar(aucCompare,'grouped');
h = legend('$\hat{\mbox{PLV}}_{\mbox{circgauss}}$','$\hat{\mbox{PLV}}_{\mbox{sample}}$');
set(h,'Interpreter','latex','fontsize',13);
axis([0 6 0.5 1.2])
fontsize = 15;
set(gca,'XTickLabel',nVec,'FontSize',fontsize);
set(gca,'Ytick',[0:0.5:1]);
set(gca,'YTickLabel',[0:0.5:1],'FontSize',fontsize);
xlabel('L','FontSize',20)
ylabel('AUC')