[cid fn] = fnames;

% run CO alg2
co_wk = cell(length(fn),2);
for i=1:length(fn)
    co_wk{i,1} = evco(fn{i},2,60,0);
    co_wk{i,2} = cid(i);
end


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% run CO alg2
co_lili = cell(length(fn),2);
for i=1:length(fn)
    co_lili{i,1} = evco(fn{i},5,60,0);
    co_lili{i,2} = cid(i);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% find k for co_wk
k_wk = zeros(length(co_wk),1);
for i=1:length(co_wk)
    r = co_wk{i,1}(:,3);
    x = co_wk{i,1}(:,4);
    k_wk(i) = calib(r,x,1);
end

% find k for co_lili
k_lili = zeros(length(co_lili),1);
for i=1:length(co_wk)
    r = co_lili{i,1}(:,3);
    x = co_lili{i,1}(:,4);
    k_lili(i) = calib(r,x,1);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% find compliances
C_wk = k_wk;

C_lili = zeros(length(k_lili),1);
for i=1:length(k_lili)
    k = k_lili(i);
    ps = mean(co_lili{i,1}(:,7));
    pd = mean(co_lili{i,1}(:,8));
    C_lili(i) = k/(ps+pd);
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% age vs C
load agegender
z = [];
for i=1:length(co_wk)
    cid = co_wk{i,2};
    ind = find(agegender(:,1)==cid);
    if ind
        z = [z; [cid agegender(ind,3) k_wk(i)]];
    end 
end
z(z(:,2)>89,:) = [];  % get rid of PHI


save tushar k_wk k_lili C_wk C_lili z 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% plotting

figure; plot(z(:,2),z(:,3),'.');
xlabel('age'); ylabel('arterial compliance (windkessel)');

figure; hold on;
plot(C_wk,'.-');
plot(C_lili,'.-r');
xlabel('patients')
ylabel('arterial compliance');
title('wk in blue, lili in red');

figure; hist(k_lili);