Added a script to generate plots for example4.m without the statistics

toolbox.

examples/example3_plots.m

@@ -1,7 +1,8 @@
-% Note: Only works in MATLAB 2014a or later since "histogram" function
-% was introduced in that version of MATLAB.
-se1=matfile('example3_se1.mat');
-se2=matfile('example3_se2.mat');
+% Run this script after example3.m. Note: This script only works in MATLAB
+% 2014a or later since "histogram" function was introduced in that version
+% of MATLAB.
+se1 = matfile('example3_se1.mat');
+se2 = matfile('example3_se2.mat');
 
 figure(1);
 subplot(2,2,1);

examples/example4_plots.m

@@ -0,0 +1,29 @@
+% Run this script after running example4.m.
+load('example4-theta04-theta2-100trials.mat');
+
+figure(1);
+clf;
+subplot(2,2,1);
+bplot([alpha_n_diff mu_n_diff s_n_diff]);
+set(gca,'YLim',[0 0.002]);
+set(gca,'XTick',1:3,'XTickLabel',{'\alpha','\mu','s'},'FontSize',14);
+ylabel('maximum absolute difference');
+title('Null models');
+
+subplot(2,2,2);
+bplot([alpha_e_diff mu_e_diff s_e_diff])
+set(gca,'YLim',[0 0.002]);
+set(gca,'XTick',1:3,'XTickLabel',{'\alpha','\mu','s'},'FontSize',14);
+ylabel('maximum absolute difference');
+title('Enrichment models');
+
+subplot(2,2,3);
+scatter(log10(bf(:, 1)), log10(bf(:, 2)),'b','*');
+hold on
+plot([-5 30],[-5 30],'-','Color','red');
+hold off
+% ozline = refline([1 0]);
+% ozline.Color = 'r';
+xlabel('Log10 Bayes factors from VARBVS','Interpreter','none');
+ylabel('Log10 Bayes factors from RSS-VARBVSR','Interpreter','none');
+title('');