load aggdb

alg = 1;

rr = cell2mat(tco);
rr = rr(:,3);

% C1 calibrated
for c=1:4
    n = length(F);
    qq = [];
    for i=1:n
        r = tco{i}(:,3);
        x = estco{i,alg}(:,1);
        e = estco{i,alg}(:,2);
        k = calib(r,x,c);
        q = k.*x;
        qq = [qq; q];
    end

    figure; plot(rr,qq,'.')
end