78c3b91900e3b55d495cf87f73e708675b10e429
[valse.git] / reports / bazar_Emilie / dessinBeta.m
1 niveauColores
2
3 for L=ind
4 for r=1:max(b(:,L))
5 betaLassoMLE(:,:,r)=inv(rhoLassoMLE(:,:,r,L)) * phiLassoMLE(:,:,r,L);
6 end
7 end
8 dessinRho = zeros(max(K),size(rhoLassoMLE,1),max(ind));
9 cpt=1;
10 gr = gray(64);
11
12 for L=ind
13 k=size(find(piLassoMLE(:,L)~=0),1)
14 figure(1)
15 %for r=1:k
16 r=1
17 subplot(1,k+1,r)
18 cpt=cpt+1;
19 colormap(gr(end:-1:1,:));
20 imagesc(abs(betaLassoMLE(:,:,r)))
21 set(gca,'xtick',[],'ytick',[])
22 set(gca, 'FontSize', 30)
23 title(['$\hat{\underline{\beta}_1}$'],'FontSize', 45,'Interpreter','latex')
24 h=colorbar;
25 set(gca, 'FontSize', 30)
26 r=2
27 subplot(1,k+1,r)
28 cpt=cpt+1;
29 colormap(gr(end:-1:1,:));
30 imagesc(abs(betaLassoMLE(:,:,r)))
31 set(gca,'xtick',[],'ytick',[])
32 set(gca, 'FontSize', 30)
33 title(['$\hat{\underline{\beta}_2}$'],'FontSize', 45,'Interpreter','latex')
34 h=colorbar;
35 set(gca, 'FontSize', 30)
36
37 %end
38 subplot(1,k+1,k+1)
39 gr = gray(64);
40 colormap(gr(end:-1:1,:));
41 imagesc(abs(betaLassoMLE(:,:,1)-betaLassoMLE(:,:,2)))
42 set(gca,'xtick',[],'ytick',[])
43 title('$|\hat{\underline{\beta}_1}-\hat{\underline{\beta}}_2|$','FontSize',45,'Interpreter','latex')
44 h=colorbar;
45 set(gca, 'FontSize', 30)
46 figure(L+k+1)
47 for r=1:k
48 for z=1:size(rhoLassoMLE,1)
49 dessinRho(r,z,L) = inv(rhoLassoMLE(z,z,r,L)^2);
50 end
51 end
52
53 %for r=1:k
54 r=1
55 subplot(k,1,r)
56 % if r==1
57 % colormap(cmapRouge(end:-1:1,:));
58 % else if r==2
59 % colormap(cmapBleu(end:-1:1,:));
60 % else colormap(cmapVert(end:-1:1,:));
61 % end
62 % end
63 imagesc(dessinRho(r,:,L))
64 set(gca,'xtick',[],'ytick',[])
65 xlabh = get(gca,'XLabel');
66 set(xlabh,'Position',get(xlabh,'Position') + [0 .15 0])
67 title('$\hat{\Sigma}_1$','FontSize', 55,'Interpreter','latex')
68 xlabel('d_4[1] d_4[2] d_4[3] d_3[1] d_3[2] d_3[3] d_3[4] d_3[5] d_3[6]','FontSize', 30)
69 set(gca, 'FontSize', 45)
70 set(gca,'Position',[0.01 0.1+0.5*(2-r) 0.9 0.2])
71
72 r=2
73 subplot(k,1,r)
74 % if r==1
75 % colormap(cmapRouge(end:-1:1,:));
76 % else if r==2
77 % colormap(cmapBleu(end:-1:1,:));
78 % else colormap(cmapVert(end:-1:1,:));
79 % end
80 % end
81 imagesc(dessinRho(r,:,L))
82 set(gca,'xtick',[],'ytick',[])
83 xlabh = get(gca,'XLabel');
84 set(xlabh,'Position',get(xlabh,'Position') + [0 .15 0])
85 title('$\hat{\Sigma}_2$','FontSize', 55,'Interpreter','latex')
86 xlabel('d_4[1] d_4[2] d_4[3] d_3[1] d_3[2] d_3[3] d_3[4] d_3[5] d_3[6]','FontSize', 30)
87 set(gca, 'FontSize', 45)
88 set(gca,'Position',[0.01 0.1+0.5*(2-r) 0.9 0.2])
89 %end
90 %h=colorbar;
91 gr = gray(64);
92 colormap(gr(end:-1:1,:));
93 axes('Position', [0.14 0.1 0.8 0.8], 'Visible', 'off');
94 % c=colorbar;
95 % set(gca, 'FontSize', 40)
96 % caxis([0 max(max(dessinRho(:,:,L)))])
97 end