update the model selection step. Beginning for the plots

[valse.git] / reports / bazar Emilie / dessinBeta.m

diff --git a/reports/bazar Emilie/dessinBeta.m b/reports/bazar Emilie/dessinBeta.m
new file mode 100644 (file)
index 0000000..78c3b91
--- /dev/null
+++ b/reports/bazar Emilie/dessinBeta.m
@@ -0,0 +1,97 @@
+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