% subject A
% import .csv files as tables; rename to acupoint names
     
%LineWidth=0.1;
set(groot,'defaultLineLineWidth',1)

figure
subplot(2,2,1);
plot(st36.TimeD,st36.rST36,'DisplayName','ST36');
hold on;
plot(st36.TimeD,st36.adjacent,'DisplayName','+6mm');
axis([8 22 3 13])
xlabel('Time of day (hours)') 
ylabel('Impedance @5KHz (K\Omega)')
legend('show', 'Location', 'northwest')
title('Subject A. Impedance at right ST-36')
grid

%figure
subplot(2,2,3);
plot(st37.TimeD,st37.rST37,'DisplayName','ST37');
hold on;
plot(st37.TimeD,st37.adjacent,'DisplayName','+6mm');
axis([8 22 3 13])
xlabel('Time of day (hours)') 
ylabel('Impedance @5KHz (K\Omega)')
legend('show', 'Location','northwest')
title('Subject A. Impedance at right ST-37')
grid


%figure
subplot(2,2,4);
plot(gb34.TimeD,gb34.rGB34,'DisplayName','GB34');
hold on;
axis([8 22 3 13])
xlabel('Time of day (hours)') 
ylabel('Impedance @5KHz (K\Omega)')
plot(gb34.TimeD,gb34.adjacent,'DisplayName','+6mm');
legend('show', 'Location','northwest')
title('Subject A. Impedance at right GB-34')
grid

%figure
subplot(2,2,2);
plot(st39.TimeD,st39.rST39,'DisplayName','ST39');
hold on;
axis([8 22 3 13])
xlabel('Time of day (hours)') 
ylabel('Impedance @5KHz (K\Omega)')
plot(st39.TimeD,st39.adjacent,'DisplayName','+6mm');
legend('show', 'Location','northwest')
title('Subject A. Impedance at right ST-39')
grid

