comparemodels.m

The BNL toolbox is a set of Matlab functions for defining and estimating the parameters of a Bayesi

list={'alarm50','alarm100','alarm200','alarm500','alarm50mis','alarm100mis','alarm200mis','alarm500mis'};
list2={'alarm50test','alarm100test','alarm200test','alarm500test','alarm50testmis','alarm100testmis','alarm200testmis','alarm500testmis'};
loglikunresttrain=zeros(8,1);
loglikunresttest=zeros(8,1);


loglikmaincumtrain=zeros(8,1);
loglikmaincumtest=zeros(8,1);

loglikmaincumtrain50=zeros(8,1);
loglikmaincumtest50=zeros(8,1);

for s=1:length(list)
    set=list{s}
    set2=list2{s}
    eval(['load c:\frank\alarm\twee\unrest_test',set]);
    loglikunresttest(s)=loglikelihood;
    eval(['load c:\frank\alarm\twee\unrest_',set]);
    loglikunresttrain(s)=loglikelihood;
    
 
    eval(['load c:\frank\alarm\twee\main_cum_test',set]);
    loglikmaincumtest(s)=loglikelihood;
    eval(['load c:\frank\alarm\twee\main_cum_',set]);
    loglikmaincumtrain(s)=loglikelihood;
      
    eval(['load c:\frank\alarm\twee\main_cum50_50_test',set]);
    loglikmaincumtest50(s)=loglikelihood;
    eval(['load c:\frank\alarm\twee\main_cum50_50_',set]);
    loglikmaincumtrain50(s)=loglikelihood;

end
plot(1:8,loglikunresttrain./repmat([50 100 200 500 ],1,2)','-b+')
hold on
plot(1:8,loglikmaincumtrain50./repmat([50 100 200 500 ],1,2)','-r+')
plot(1:8,loglikmaincumtrain./repmat([50 100 200 500 ],1,2)','-g+')
plot(1:8,loglikunresttest./repmat([50 100 200 500 ],1,2)','-b')
hold on
plot(1:8,loglikmaincumtest./repmat([50 100 200 500 ],1,2)','-g')
plot(1:8,loglikmaincumtest50./repmat([50 100 200 500 ],1,2)','-r')

subplot(1,2,1),plot(1:4,loglikunresttest(1:4)./repmat([50 100 200 500 ],1,1)','-*')
hold on

subplot(1,2,1),plot(1:4,loglikmaincumtest(1:4)./repmat([50 100 200 500 ],1,1)','-+')

subplot(1,2,1),plot(1:4,loglikmaincumtest50(1:4)./repmat([50 100 200 500 ],1,1)','-^')

subplot(1,2,2),plot(5:8,loglikunresttest(5:8)./repmat([50 100 200 500 ],1,1)','-*')
hold on

subplot(1,2,2),plot(5:8,loglikmaincumtest(5:8)./repmat([50 100 200 500 ],1,1)','-+')

subplot(1,2,2),plot(5:8,loglikmaincumtest50(5:8)./repmat([50 100 200 500 ],1,1)','-^')



subplot(1,2,1),plot(1:4,loglikunresttest(1:4),'-*')
hold on

subplot(1,2,1),plot(1:4,loglikmaincumtest(1:4),'-+')

subplot(1,2,1),plot(1:4,loglikmaincumtest50(1:4),'-^')

subplot(1,2,2),plot(5:8,loglikunresttest(5:8),'-*')
hold on

subplot(1,2,2),plot(5:8,loglikmaincumtest(5:8),'-+')

subplot(1,2,2),plot(5:8,loglikmaincumtest50(5:8),'-^')