5 betaLassoMLE(:,:,r)=inv(rhoLassoMLE(:,:,r,L)) * phiLassoMLE(:,:,r,L);
8 dessinRho = zeros(max(K),size(rhoLassoMLE,1),max(ind));
13 k=size(find(piLassoMLE(:,L)~=0),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')
25 set(gca, 'FontSize', 30)
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')
35 set(gca, 'FontSize', 30)
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')
45 set(gca, 'FontSize', 30)
48 for z=1:size(rhoLassoMLE,1)
49 dessinRho(r,z,L) = inv(rhoLassoMLE(z,z,r,L)^2);
57 % colormap(cmapRouge(end:-1:1,:));
59 % colormap(cmapBleu(end:-1:1,:));
60 % else colormap(cmapVert(end:-1:1,:));
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])
75 % colormap(cmapRouge(end:-1:1,:));
77 % colormap(cmapBleu(end:-1:1,:));
78 % else colormap(cmapVert(end:-1:1,:));
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])
92 colormap(gr(end:-1:1,:));
93 axes('Position', [0.14 0.1 0.8 0.8], 'Visible', 'off');
95 % set(gca, 'FontSize', 40)
96 % caxis([0 max(max(dessinRho(:,:,L)))])