% Plot residuals
%
% arPlotResiduals(saveToFile)
%
% saveToFile    [false]

function arPlotResiduals(saveToFile)

global ar

if(isempty(ar))
    error('please initialize by arInit')
end

if(~exist('saveToFile','var'))
    saveToFile = false;
end

% constants
labelfontsize = 12;
labelfonttype = 'TimesNewRoman';

h = myRaiseFigure('Residuals');

tmpres = [];
for m=1:length(ar.model)
    for d=1:length(ar.model(m).data)
        for y=1:size(ar.model(m).data(d).res, 2)
            tmpres1 = ar.model(m).data(d).res(:,y);
            tmpres = [tmpres; tmpres1(~isnan(tmpres1))]; %#ok<AGROW>
        end
    end
end

g = subplot(1,2,1);
mySubplotStyle(g, labelfontsize, labelfonttype);
qqplot(tmpres);
box(g, 'on')
axis(g,'equal')
axis(g,'square')
title(g, 'QQ-plot of sample data versus standard normal');
ylabel(g, 'Quantiles of Input Sample');
xlabel(g, 'Standard Normal Quantiles');

g = subplot(1,2,2);
mySubplotStyle(g, labelfontsize, labelfonttype);
for m=1:length(ar.model)
    for d=1:length(ar.model(m).data)
        for y=1:size(ar.model(m).data(d).res, 2)
            tmpres = ar.model(m).data(d).res(:,y);
            ttmp = ar.model(m).data(d).tExp(~isnan(tmpres));
            tmpres = tmpres(~isnan(tmpres));
			tunique = unique(ttmp); %R2013a compatible
            for jt = 2:length(tunique);
                res1 = tmpres(ttmp==tunique(jt));
                res2 = tmpres(ttmp==tunique(jt-1));
                for jr=1:length(res1)
                    plot(res1(jr), res2, 'b+')
                    hold(g, 'on');
                end
            end
        end
    end
end
plot(xlim*0, ylim, 'r-.');
plot(xlim, ylim*0, 'r-.');
hold(g, 'off');
box(g, 'on')
axis(g,'equal')
axis(g,'square')
title(g, 'Auto-correlation of normalized residuals')
ylabel(g, 'residuals t_i')
xlabel(g, 'residuals t_{i+1}')

if(saveToFile)
    ar.plot.savePath_ResFig = mySaveFigure(h, 'residuals');
end

%% sub-functions



function h = myRaiseFigure(figname)
global ar
openfigs = get(0,'Children');

figcolor = [1 1 1];

ar.plot.time = now;

if(isfield(ar.plot, 'fighandel_res') && ~isempty(ar.plot.fighandel_res) && ...
    ar.plot.fighandel_res ~= 0 && ...
    sum(ar.plot.fighandel_res==openfigs)>0 && ...
    strcmp(get(ar.plot.fighandel_res, 'Name'), figname))

    h = ar.plot.fighandel_res;
    figure(h);
else
    h = figure('Name', figname, 'NumberTitle','off', ...
        'Units', 'normalized', 'Position', ...
        [0.71 0.7 0.25 0.2]);
    set(h,'Color', figcolor);
    ar.plot.fighandel_res = h;
end



function savePath = mySaveFigure(h, name)
savePath = [arSave '/Figures'];

if(~exist(savePath, 'dir'))
    mkdir(savePath)
end

savePath = mypath([savePath '/' name]);

saveas(h, savePath, 'fig');
print('-depsc2', savePath);
eval(['!ps2pdf  -dEPSCrop ' savePath '.eps '  savePath '.pdf']);



function str = mypath(str)
str = strrep(str, ' ', '\ ');
str = strrep(str, '(', '\(');
str = strrep(str, ')', '\)');




function mySubplotStyle(g, labelfontsize, labelfonttype)
set(g, 'FontSize', labelfontsize);
set(g, 'FontName', labelfonttype);