FAScaleRat=55.9686/100;
%str_Data255.meanRBC2TIS
%str_Data255.meanTiss2Gas/FAScaleRat
%str_Data255.meanRBC2Gas/FAScaleRat
%str_Data255.meanDiss2Gas/FAScaleRat

%             255,  259, 266, 268, 281, 284, 285, 293, 316, 379, 384, 391
%              0,    0,   0,   0,   0,   0,   0,   0,   0,   1,   0,   0
%              W,    W,   W,   W,   W,   W,   W,   W,   W,   W,   A,   W
 
%             290,  313, 438, 446, 388, 444, 306, 291, 303, 469, 470, 437
%              1,    0,   1,   1,   1,   0,   0,   0,   1,   1,   1,   0
%              W,    W,   W,   W,   AA,  W,   W,   AA,  


%            538,  282, 292, 371, 476, 499, 520, 522
%             0,    0,   0,   1,   1,   0,   0,   0
%             W,    W,   W,   W,   W,   W,   W,   W


%            305,  271, 451, 442, 321, 314, 439, 399
%             1,    1,   0,   0,   0,   0,   0,   0
%            API,  W,   W,   W,   W,   W,   W,   W

%[p_FE1_Old,h_FE1_Old] = ranksum(Age_Heal_You,Age_Heal_Old)
%[p_R2T_YO,h_R2T_YO] = ranksum(R2T_Heal_You,R2T_Heal_Old)

%

%            255,  259, 266, 268, 281, 284, 285, 293, 316, 379, 384, 391,
%            379,  384,  387,  391
%Age=23+/-5 yrs
Age_Heal_You=[18,  21,  21,  21,  21,  32,  25,  21,  21,  32,  20,   20];
R2T_Heal_You=[0.27,0.35,0.28,0.27,0.33,0.28,0.32,0.32,0.34,0.38,0.34,0.28];
T2G_Heal_You=[1.34,1.13,1.29,1.25,1.16,1.13,1.01,1.06,1.09,1.21,1.02,1.33];
R2G_Heal_You=[0.36,0.39,0.36,0.34,0.39,0.32,0.33,0.34,0.37,0.46,0.35,0.37];
D2G_Heal_You=[1.55,1.40,1.49,1.44,1.35,1.36,1.29,1.34,1.23,1.49,1.24,1.51];
%Ven_Heal_You=[10, 10,   10,   8,  10,  10,  10,  10,  10,  12, 10.06, 9.9];
Ven_Heal_You=[4,   2,     9,   8,   4,   3,  7,   9,   12,  12, 10.06, 9.9];
FE1_Heal_You=[119, 105, 108,  91,  95,   94,  96,  92,  96,  104,  116, 85];
FVC_Heal_You=[122, 104, 105,  96,  100,  95,  110, 104, 101, 104,  117, 95];
F1V_Heal_You=[101, 100, 102,  94,  94,   99,  87,  88,  94,  102,  99,  88];
F25_Heal_You=[114, 103, 107,  70,  76,   87,  69,  67,  77,  103,  91,  57];
BMI_Heal_You=[24,  24,  23,   21,  18,   25,  22,  22,  18,  22,   23,  30];


%             na,376_in,275?,  , 316,  na, 293,    ,  281

%            290,  313, 438, 446, 388, 444, 306, 291, 303, 469, 470, 437
%Age=24+/-5 yrs
Age_Asth_You=[26,  28,  18,  19,  21,  21,  22,  22,  26,  21,  21,  36];
R2T_Asth_You=[0.51,0.35,0.49,0.41,0.48,0.26,0.33,0.32,0.35,0.57,0.43,0.30];
T2G_Asth_You=[0.61,0.85,1.04,1.42,0.78,1.25,0.83,1.25,0.81,0.90,0.97,1.12];
R2G_Asth_You=[0.31,0.30,0.51,0.58,0.37,0.32,0.27,0.40,0.28,0.52,0.42,0.34];
D2G_Asth_You=[0.79,1.03,1.48,1.82,1.03,1.46,0.96,1.55,0.92,1.27,1.13,1.40];
Ven_Asth_You=[25,  17,  12,  15,  14,  13,  9,   9,   24,  10,  10,   13];
Sev_Asth_You=[2,   2,   3,   4,   3,   4,   3,   4,   4,   4,   3,   3];
DLC_Asth_You=[101, 111, 105, 91,  99,  96,  108, 82,  138, 70,  98,  79];
TLC_Asth_You=[125, 127, 100, 104, 105, 108, 100, 96,  130, 91,  107, 84];
RV_Asth_You= [165, 136, 101, 84,  97,  126, 115, 94,  271, 158, 155, 88];
RTC_Asth_You=[133, 107, 101, 81,  93,  117, 115, 98,  208, 175, 146, 104];
FE1_Asth_You=[80,  100, 96,  105, 71,  75,  71,  67,  77,  59,  96,  74];
FVC_Asth_You=[112, 128, 111, 116, 94,  102, 84,  80,  107, 69,  116, 90];
F1V_Asth_You=[70,  78,  89,  92,  81,  73,  84,  84,  72,  87,  83,  81];
F25_Asth_You=[41,  56,  72,  86,  48,  63,  46,  42,  39,  43,  59,  64];
BMI_Asth_You=[18,  27,  22,  26,  31,  19,  23,  24,  27,  20,  29,  24];


[p_R2T_You,h_R2T_You] = ranksum(R2T_Heal_You,R2T_Asth_You)
[p_T2G_You,h_T2G_You] = ranksum(T2G_Heal_You,T2G_Asth_You)
[p_R2G_You,h_R2G_You] = ranksum(R2G_Heal_You,R2G_Asth_You)
[p_Ven_You,h_Ven_You] = ranksum(Ven_Heal_You,Ven_Asth_You)

% plot out R2T
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [R2T_Heal_You R2T_Asth_You];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(R2T_Heal_You)) ones(size(R2T_Asth_You))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('RBC / Tissue','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age<35','fontsize',22,'fontweight','b')

ylim([0.1 0.6])

hold on;
plot(ones(size(R2T_Heal_You)), R2T_Heal_You,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(R2T_Asth_You))+1, R2T_Asth_You,'b*', 'MarkerSize',10);


% plot out R2G
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [R2G_Heal_You R2G_Asth_You];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(R2G_Heal_You)) ones(size(R2G_Asth_You))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('RBC / Gas','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age<35','fontsize',22,'fontweight','b')

ylim([0.1 0.6])

hold on;
plot(ones(size(R2G_Heal_You)), R2G_Heal_You,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(R2G_Asth_You))+1, R2G_Asth_You,'b*', 'MarkerSize',10);

% plot out T2G
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [T2G_Heal_You T2G_Asth_You];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(T2G_Heal_You)) ones(size(T2G_Asth_You))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('Tissue/Gas','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age<35','fontsize',22,'fontweight','b')

ylim([0.6 1.5])

hold on;
plot(ones(size(T2G_Heal_You)), T2G_Heal_You,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(T2G_Asth_You))+1, T2G_Asth_You,'b*', 'MarkerSize',10);

% plot out ventilation
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [Ven_Heal_You Ven_Asth_You];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(Ven_Heal_You)) ones(size(Ven_Asth_You))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('Ventilation Defect [%]','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age<35','fontsize',22,'fontweight','b')

ylim([0  60])

hold on;
plot(ones(size(Ven_Heal_You)), Ven_Heal_You,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(Ven_Asth_You))+1, Ven_Asth_You,'b*', 'MarkerSize',10);


%{
%            538,  282, 292, 371, 476, 499, 520,  522,  
%Age=59+/-4 yrs
Age_Heal_Old=[60,  55,  59,  59,  60,  58,  52,   67];
R2T_Heal_Old=[0.23,0.22,0.20,0.28,0.28,0.29,0.34,0.27];
T2G_Heal_Old=[1.18,1.43,1.18,1.01,1.02,1.27,1.05,1.19];
R2G_Heal_Old=[0.26,0.32,0.23,0.28,0.29,0.32,0.33,0.33];
D2G_Heal_Old=[1.31,1.61,1.25,1.08,1.20,1.43,1.31,1.40];
Ven_Heal_Old=[38,  17,  18,  17,  24,  24,  20,   17];
%}
%            528,  282, 292, 371, 476, 499, 520,  522,  
%Age=59+/-4 yrs
Age_Heal_Old=[52,  55,  59,  59,  60,  58,  52,   67];
R2T_Heal_Old=[0.28,0.22,0.20,0.28,0.28,0.29,0.34,0.27];
T2G_Heal_Old=[0.92,1.43,1.18,1.01,1.02,1.27,1.05,1.19];
R2G_Heal_Old=[0.26,0.32,0.23,0.28,0.29,0.32,0.33,0.33];
D2G_Heal_Old=[1.10,1.61,1.25,1.08,1.20,1.43,1.31,1.40];
Ven_Heal_Old=[16,  17,  18,  17,  24,  24,  20,   17];
DLC_Heal_Old=[82, 122, 100,  66,  97,  113, 111,  90];
TLC_Heal_Old=[101,112, 104, 104,  113, 107, 115,  110];
RV_Heal_Old= [84, 109, 100, 100,  96,  93,  106,  114];
RTC_Heal_Old=[83, 97,  96,  78,   85,  87,  92,   104];
FE1_Heal_Old=[103,124, 101, 99,   109, 111, 106,  102];
FVC_Heal_Old=[110, 126, 103, 109, 117, 120, 115,  111];
F1V_Heal_Old=[92,  96,  96,  92,  93,  98,  90,   90];
F25_Heal_Old=[78,  110, 89,  75,  87,  111, 76,   73];
BMI_Heal_Old=[22,  29,  25,  26,  28,  27,  21,   23];



%            305,  271, 451, 442, 321, 314, 439, 399
%Age=58+/-9 yrs
Age_Asth_Old=[45,  53,  69,  59,  49,  57,  60,  72];
R2T_Asth_Old=[0.34,0.20,0.30,0.13,0.42,0.19,0.20,0.23];
T2G_Asth_Old=[1.07,0.97,1.27,1.14,0.80,0.87,1.43,0.93];
R2G_Asth_Old=[0.37,0.19,0.38,0.16,0.33,0.17,0.29,0.22];
D2G_Asth_Old=[1.29,1.06,1.54,1.21,1.03,0.92,1.58,1.01];
Ven_Asth_Old=[18,  39,  18,  17,  45,  60,  30,  55];
Sev_Asth_Old=[1,   2,   2,   4,   4,   4,   4,   4];
DLC_Asth_Old=[77,  80,  104, 82,  110, 76,  98, 82];
TLC_Asth_Old=[97,  132, 108, 85,  108, 124, 96, 135];
RV_Asth_Old= [127, 147, 108, 87,  164, 218, 116, 206];
RTC_Asth_Old=[131, 104, 100, 103, 152, 175, 121, 152];
FE1_Asth_Old=[83,  89,  110, 73,  48,  52,  85,  48];
FVC_Asth_Old=[86,  138, 130, 87,  83,  76,  90,  111];
F1V_Asth_Old=[96,  69,  83,  82,  58,  67,  92,  57];
F25_Asth_Old=[72,  39,  81,  46,  18,  25,  62,  16];
BMI_Asth_Old=[29,  25,  24,  26,  19,  24,  36,  28];


[p_R2T_Old,h_R2T_Old] = ranksum(R2T_Heal_Old,R2T_Asth_Old)
[p_T2G_Old,h_T2G_Old] = ranksum(T2G_Heal_Old,T2G_Asth_Old)
[p_R2G_Old,h_R2G_Old] = ranksum(R2G_Heal_Old,R2G_Asth_Old)
[p_Ven_Old,h_Ven_Old] = ranksum(Ven_Heal_Old,Ven_Asth_Old)


% plot out R2T
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [R2T_Heal_Old R2T_Asth_Old];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(R2T_Heal_Old)) ones(size(R2T_Asth_Old))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('RBC / Tissue','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age>45','fontsize',22,'fontweight','b')

ylim([0.1 0.6])

hold on;
plot(ones(size(R2T_Heal_Old)), R2T_Heal_Old,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(R2T_Asth_Old))+1, R2T_Asth_Old,'b*', 'MarkerSize',10);


% plot out R2G
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [R2G_Heal_Old R2G_Asth_Old];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(R2G_Heal_Old)) ones(size(R2G_Asth_Old))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('RBC / Gas','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age>45','fontsize',22,'fontweight','b')

ylim([0.1 0.6])

hold on;
plot(ones(size(R2G_Heal_Old)), R2G_Heal_Old,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(R2G_Asth_Old))+1, R2G_Asth_Old,'b*', 'MarkerSize',10);

% plot out T2G
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [T2G_Heal_Old T2G_Asth_Old];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(T2G_Heal_Old)) ones(size(T2G_Asth_Old))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('Tissue/Gas','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age>45','fontsize',22,'fontweight','b')

ylim([0.6 1.5])

hold on;
plot(ones(size(T2G_Heal_Old)), T2G_Heal_Old,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(T2G_Asth_Old))+1, T2G_Asth_Old,'b*', 'MarkerSize',10);


% plot out ventilation
hh=figure;
set(hh, 'DefaultTextFontSize', 18);
x = [Ven_Heal_Old Ven_Asth_Old];  % [0.4 0.5 0.1 0.7 0.9];
group = [ones(size(Ven_Heal_Old)) ones(size(Ven_Asth_Old))+1 ]; % [1 1 1 2 2];
bh=boxplot(x, group, 'labels',{'Normal','Asthma'});

set(bh(:,:),'linewidth',3);
hh2 = findobj(gca,'Type','text');
set(hh2,'fontweight','b')

%xlabel('Subject','fontsize',22,'fontweight','b')
ylabel('Ventilation Defect [%]','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')

color = ['b', 'r', 'm', 'y'];
h = findobj(gca,'Tag','Box');

title('Age>45','fontsize',22,'fontweight','b')

ylim([0  60])

hold on;
plot(ones(size(Ven_Heal_Old)), Ven_Heal_Old,'k*', 'MarkerSize',10)
hold on;
plot(ones(size(Ven_Asth_Old))+1, Ven_Asth_Old,'b*', 'MarkerSize',10);

%R2T
[R_tlc_r2t, P_tlc_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [TLC_Asth_You,TLC_Asth_Old])
[R_rv_r2t, P_rv_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [RV_Asth_You,RV_Asth_Old])
[R_dlco_r2t, P_dlco_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [DLC_Asth_You,DLC_Asth_Old])
[R_sev_r2t, P_sev_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [Sev_Asth_You,Sev_Asth_Old])
[R_rtc_r2t, P_rtc_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [RTC_Asth_You,RTC_Asth_Old])
[R_fe1_r2t, P_fe1_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [FE1_Asth_You,FE1_Asth_Old])
[R_fvc_r2t, P_fvc_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [FVC_Asth_You,FVC_Asth_Old])
[R_f1v_r2t, P_f1v_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [F1V_Asth_You,F1V_Asth_Old])
[R_f25_r2t, P_f25_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [F25_Asth_You,F25_Asth_Old])
[R_ven_r2t, P_ven_r2t]= corrcoef([R2T_Asth_You,R2T_Asth_Old], [Ven_Asth_You,Ven_Asth_Old])




%T2G
[R_tlc_t2g, P_tlc_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [TLC_Asth_You,TLC_Asth_Old])
[R_rv_t2g, P_rv_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [RV_Asth_You,RV_Asth_Old])
[R_dlco_t2g, P_dlco_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [DLC_Asth_You,DLC_Asth_Old])
[R_sev_t2g, P_sev_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [Sev_Asth_You,Sev_Asth_Old])
[R_rtc_t2g, P_rtc_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [RTC_Asth_You,RTC_Asth_Old])
[R_fe1_t2g, P_fe1_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [FE1_Asth_You,FE1_Asth_Old])
[R_fvc_t2g, P_fvc_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [FVC_Asth_You,FVC_Asth_Old])
[R_f1v_t2g, P_f1v_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [F1V_Asth_You,F1V_Asth_Old])
[R_f25_t2g, P_f25_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [F25_Asth_You,F25_Asth_Old]);
[R_ven_t2g, P_ven_t2g]= corrcoef([T2G_Asth_You,T2G_Asth_Old], [Ven_Asth_You,Ven_Asth_Old])

%R2G
[R_tlc_r2g, P_tlc_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [TLC_Asth_You,TLC_Asth_Old])
[R_rv_r2g, P_rv_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [RV_Asth_You,RV_Asth_Old])
[R_dlco_r2g, P_dlco_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [DLC_Asth_You,DLC_Asth_Old])
[R_sev_r2g, P_sev_r2g] = corrcoef([R2G_Asth_You,R2G_Asth_Old], [Sev_Asth_You,Sev_Asth_Old])
[R_rtc_r2g, P_rtc_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [RTC_Asth_You,RTC_Asth_Old])
[R_fe1_r2g, P_fe1_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [FE1_Asth_You,FE1_Asth_Old])
[R_fvc_r2g, P_fvc_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [FVC_Asth_You,FVC_Asth_Old])
[R_f1v_r2g, P_f1v_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [F1V_Asth_You,F1V_Asth_Old])
[R_f25_r2g, P_f25_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [F25_Asth_You,F25_Asth_Old])
[R_ven_r2g, P_ven_r2g]= corrcoef([R2G_Asth_You,R2G_Asth_Old], [Ven_Asth_You,Ven_Asth_Old])


%ventilation
[R_tlc_vent, P_tlc_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [TLC_Asth_You,TLC_Asth_Old])
[R_rv_vent, P_rv_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [RV_Asth_You,RV_Asth_Old])
[R_dlco_vent, P_dlco_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [DLC_Asth_You,DLC_Asth_Old])
[R_sev_vent, P_sev_vent] = corrcoef([Ven_Asth_You,Ven_Asth_Old], [Sev_Asth_You,Sev_Asth_Old])
[R_t2g_vent, P_t2g_vent] = corrcoef([Ven_Asth_You,Ven_Asth_Old], [T2G_Asth_You,T2G_Asth_Old])
[R_r2t_vent, P_r2t_vent] = corrcoef([Ven_Asth_You,Ven_Asth_Old], [R2T_Asth_You,R2T_Asth_Old])
[R_rtc_vent, P_rtc_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [RTC_Asth_You,RTC_Asth_Old])
[R_fe1_vent, P_fe1_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [FE1_Asth_You,FE1_Asth_Old])
[R_fvc_vent, P_fvc_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [FVC_Asth_You,FVC_Asth_Old])
[R_f1v_vent, P_f1v_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [F1V_Asth_You,F1V_Asth_Old])
[R_f25_vent, P_f25_vent]= corrcoef([Ven_Asth_You,Ven_Asth_Old], [F25_Asth_You,F25_Asth_Old])


    



%correlation
figure
plot([RV_Asth_You], [T2G_Asth_You],'sg', 'Linewidth',3,'MarkerSize',12');
hold on;
plot([RV_Asth_Old],[T2G_Asth_Old], 'db', 'Linewidth',3,'MarkerSize',12');

xlabel('RV %pred','fontsize',22,'fontweight','b')
ylabel('Tissue/Gas Ratio [%]','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')
 %axis([1.5 7e0 0.05 0.6])

p = polyfit([RV_Asth_You,RV_Asth_Old],[T2G_Asth_You,T2G_Asth_Old],1);  
x = [75:300];
f=polyval(p,x) ;
hold on;
plot(x,f,'LineWidth', 2)

lh=legend( 'Younger', 'Older','Fitted', 'Location','NorthWest'   ); 
set( lh, 'fontsize', 18, 'fontweight','b');
M = findobj(lh,'type','line');
set(M,'linewidth',3) % this will also work on a vector of handles 


%correlation
figure
plot([Ven_Asth_You], [R2T_Asth_You],'sg', 'Linewidth',3,'MarkerSize',12');
hold on;
plot([Ven_Asth_Old],[R2T_Asth_Old], 'db', 'Linewidth',3,'MarkerSize',12');

xlabel('Vent Def %','fontsize',22,'fontweight','b')
ylabel('RBC/Tissue Ratio','fontsize',22,'fontweight','b')
%title('Healthy','fontsize',22,'fontweight','b')

set (gca,'fontsize',18,'fontweight','b');
set(gcf,'Color','white')
 %axis([1.5 7e0 0.05 0.6])

p = polyfit([Ven_Asth_You,Ven_Asth_Old],[R2T_Asth_You,R2T_Asth_Old],1);  
x = [0:80];
f=polyval(p,x) ;
hold on;
plot(x,f,'LineWidth', 2)

lh=legend( 'Younger', 'Older','Fitted', 'Location','NorthWest'   ); 
set( lh, 'fontsize', 18, 'fontweight','b');
M = findobj(lh,'type','line');
set(M,'linewidth',3) % this will also work on a vector of handles 

%healthy vs old
[p_FE1_YO,h_FE1_YO] = ranksum(FE1_Asth_You,FE1_Asth_Old)
[p_FVC_YO,h_FVC_YO] = ranksum(FVC_Asth_You,FVC_Asth_Old)
[p_F1V_YO,h_F1V_YO] = ranksum(F1V_Asth_You,F1V_Asth_Old)
[p_F25_YO,h_F25_YO] = ranksum(F25_Asth_You,F25_Asth_Old)
[p_DLC_YO,h_DLC_YO] = ranksum(DLC_Asth_You,DLC_Asth_Old)
[p_TLC_YO,h_TLC_YO] = ranksum(TLC_Asth_You,TLC_Asth_Old)
[p_RV_YO,h_RV_YO] = ranksum(RV_Asth_You,RV_Asth_Old)
[p_RTC_YO,h_RTC_YO] = ranksum(RTC_Asth_You,RTC_Asth_Old)

[p_Ven_YO,h_Ven_YO] = ranksum(Ven_Asth_You,Ven_Asth_Old)
[p_R2T_YO,h_R2T_YO] = ranksum(R2T_Asth_You,R2T_Asth_Old)
[p_T2G_YO,h_T2G_YO] = ranksum(T2G_Asth_You,T2G_Asth_Old)