+++ /dev/null
-niveauColores
-
-for L=ind
- for r=1:max(b(:,L))
- betaLassoMLE(:,:,r)=inv(rhoLassoMLE(:,:,r,L)) * phiLassoMLE(:,:,r,L);
- end
-end
-dessinRho = zeros(max(K),size(rhoLassoMLE,1),max(ind));
-cpt=1;
-gr = gray(64);
-
- for L=ind
- k=size(find(piLassoMLE(:,L)~=0),1)
- figure(1)
- %for r=1:k
- r=1
- subplot(1,k+1,r)
- cpt=cpt+1;
- colormap(gr(end:-1:1,:));
- imagesc(abs(betaLassoMLE(:,:,r)))
- set(gca,'xtick',[],'ytick',[])
- set(gca, 'FontSize', 30)
- title(['$\hat{\underline{\beta}_1}$'],'FontSize', 45,'Interpreter','latex')
- h=colorbar;
- set(gca, 'FontSize', 30)
- r=2
- subplot(1,k+1,r)
- cpt=cpt+1;
- colormap(gr(end:-1:1,:));
- imagesc(abs(betaLassoMLE(:,:,r)))
- set(gca,'xtick',[],'ytick',[])
- set(gca, 'FontSize', 30)
- title(['$\hat{\underline{\beta}_2}$'],'FontSize', 45,'Interpreter','latex')
- h=colorbar;
- set(gca, 'FontSize', 30)
-
- %end
- subplot(1,k+1,k+1)
- gr = gray(64);
- colormap(gr(end:-1:1,:));
- imagesc(abs(betaLassoMLE(:,:,1)-betaLassoMLE(:,:,2)))
- set(gca,'xtick',[],'ytick',[])
- title('$|\hat{\underline{\beta}_1}-\hat{\underline{\beta}}_2|$','FontSize',45,'Interpreter','latex')
- h=colorbar;
- set(gca, 'FontSize', 30)
- figure(L+k+1)
- for r=1:k
- for z=1:size(rhoLassoMLE,1)
- dessinRho(r,z,L) = inv(rhoLassoMLE(z,z,r,L)^2);
- end
- end
-
- %for r=1:k
- r=1
- subplot(k,1,r)
-% if r==1
-% colormap(cmapRouge(end:-1:1,:));
-% else if r==2
-% colormap(cmapBleu(end:-1:1,:));
-% else colormap(cmapVert(end:-1:1,:));
-% end
-% end
- imagesc(dessinRho(r,:,L))
- set(gca,'xtick',[],'ytick',[])
- xlabh = get(gca,'XLabel');
- set(xlabh,'Position',get(xlabh,'Position') + [0 .15 0])
- title('$\hat{\Sigma}_1$','FontSize', 55,'Interpreter','latex')
- 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)
- set(gca, 'FontSize', 45)
- set(gca,'Position',[0.01 0.1+0.5*(2-r) 0.9 0.2])
-
- r=2
- subplot(k,1,r)
-% if r==1
-% colormap(cmapRouge(end:-1:1,:));
-% else if r==2
-% colormap(cmapBleu(end:-1:1,:));
-% else colormap(cmapVert(end:-1:1,:));
-% end
-% end
- imagesc(dessinRho(r,:,L))
- set(gca,'xtick',[],'ytick',[])
- xlabh = get(gca,'XLabel');
- set(xlabh,'Position',get(xlabh,'Position') + [0 .15 0])
- title('$\hat{\Sigma}_2$','FontSize', 55,'Interpreter','latex')
- 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)
- set(gca, 'FontSize', 45)
- set(gca,'Position',[0.01 0.1+0.5*(2-r) 0.9 0.2])
- %end
- %h=colorbar;
- gr = gray(64);
- colormap(gr(end:-1:1,:));
- axes('Position', [0.14 0.1 0.8 0.8], 'Visible', 'off');
-% c=colorbar;
-% set(gca, 'FontSize', 40)
-% caxis([0 max(max(dessinRho(:,:,L)))])
-end
\ No newline at end of file