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)
Comments
0 B
|👍
/👎
0 B
|👍
/👎