restart=1;
weightcost=0;
data=train2;
batchify8;
rbmlearn;

testdata=test2;
getscores;
scoresmodel2data2=testscores;

testdata=test3;
getscores;
scoresmodel2data3=testscores;


restart=1;
data=train3;
batchify8;
rbmlearn;

testdata=test2;
getscores;
scoresmodel3data2=testscores;

testdata=test3;
getscores;
scoresmodel3data3=testscores;

plot23test;

posscores=scoresmodel2data2-scoresmodel3data2;
negscores=scoresmodel2data3-scoresmodel3data3;
%% if you save posscores and negscores for this model pair, you can
%% combine models by just adding them to posscores and negscores 
%% produced by  another model pair.
%% It might be better to average rahter than add so that the numbers
%% do not get too big, because the code for getoverlap is fragile.


overlap=getoverlap(posscores,negscores);

fprintf(1, 'total negative log prob on test data = %4.2f bits \n', overlap);