function analyze_predictions(sbgrp1,sbgrp2)

pairstring = [num2str(sbgrp1) '.' num2str(sbgrp2)];

cd('D:\Bio Terrorism\Data\Hourly\Extended data\');
totvol = load('volume.All.agg','Ascii');
sbgrp12data = load(['count.' pairstring '.agg'],'Ascii');
if isempty(sbgrp12data)
    fprintf('The pair information could not be loaded');
end;
sbgrp1data = load(['volume.' num2str(sbgrp1) '.agg'],'Ascii');
if isempty(sbgrp1data)
    fprintf('The information for the 1st subgroup could not be loaded');
end;
sbgrp2data = load(['volume.' num2str(sbgrp2) '.agg'],'Ascii');
if isempty(sbgrp2data)
    fprintf('The information for the 2nd subgroup could not be loaded');
end;

cd('D:\Bio Terrorism\Matlab\');
[dates12,counts12,missindex12] = fillingapwithperiodicities(sbgrp12data(:,1),sbgrp12data(:,2)./totvol(:,2),[num2str(sbgrp1) '.' num2str(sbgrp2)]);
[dates1,counts1,missindex1] = fillingapwithperiodicities(sbgrp1data(:,1),sbgrp1data(:,2)./totvol(:,2),num2str(sbgrp1));
[dates2,counts2,missindex2] = fillingapwithperiodicities(sbgrp2data(:,1),sbgrp2data(:,2)./totvol(:,2),num2str(sbgrp2));

predstd(counts12,1,missindex12,pairstring);
predstd(counts1,1,missindex1,num2str(sbgrp1));
predstd(counts2,1,missindex2,num2str(sbgrp2));