%%Plot PEM vs ALS
load PEMCS_ALSX.txt
PEMCS = PEMCS_ALSX(:,1);
ALSX_CS = PEMCS_ALSX(:,2);
figure
subplot(3,2,1)
plot(PEMCS,ALSX_CS, '.')
title('PEM-CS vs. ALSX (r = 0.4768)')
r1 = corr(PEMCS, ALSX_CS)

load PEMCS_ALSY.txt
%PEMCS = PEMCS_ALSY(:,1);
ALSY_CS = PEMCS_ALSY(:,2);
subplot(3,2,2)
plot(PEMCS,ALSY_CS,'.')
title('PEM-CS vs. ALSY (r=0.7067)')
r2 = corr(PEMCS,ALSY_CS)

load PEMEX_ALSX.txt
PEMEX = PEMEX_ALSX(:,1);
ALSX_EX = PEMEX_ALSX(:,2);
subplot(3,2,3)
plot(PEMEX,ALSX_EX,'.')
title('PEM-EX vs. ALSX (r=0.6026)')
r3 = corr(PEMEX, ALSX_EX)

load PEMEX_ALSY.txt
ALSY_EX = PEMEX_ALSY(:,2);
subplot(3,2,4)
plot(PEMEX,ALSY_EX,'.')
title('PEM-EX vs. ALSY (r=0.7783)')
r4 = corr(PEMEX, ALSY_EX)

load PEMEY_ALSX.txt
PEMEY = PEMEY_ALSX(:,1);
ALSX_EY = PEMEY_ALSX(:,2);
subplot(3,2,5)
plot(PEMEY,ALSX_EY,'.')
title('PEM-EY vs. ALSX (r=0.4978)')
r5 = corr(PEMEY, ALSX_EY)


load PEMEY_ALSY.txt
PEMEY = PEMEY_ALSY(:,1);
ALSY_EY = PEMEY_ALSY(:,2);
subplot(3,2,6)
plot(PEMEY,ALSY_EY,'.')
title('PEM-EY vs. ALSY (r=0.7239)')
r6 = corr(PEMEY, ALSY_EY)

%%PEM vs ISI 
%%Since PEMX and ALSY correlate the most, I'll plot just PEMEX and ISI ETMY

load PEMEX_ETMYX100.txt
PEMEX = PEMEX_ETMYX100(:,1);
ETMYX100 = PEMEX_ETMYX100(:,2);
figure
subplot(4,3,1)
plot(PEMEX, ETMYX100,'.')
title('PEM-EX vs. ISI-ETMY\_X 100-300mHz (r=0.6733)')
r7 = corr(PEMEX, ETMYX100)

load PEMEX_ETMYX300.txt
PEMEX = PEMEX_ETMYX300(:,1);
ETMYX300 = PEMEX_ETMYX300(:,2);
subplot(4,3,2)
plot(PEMEX, ETMYX300,'.')
title('PEM-EX vs. ISI-ETMY\_X 0.3-1 Hz (r=0.7081)')
r8 = corr(PEMEX, ETMYX300)

load PEMEX_ETMYY100.txt
PEMEX = PEMEX_ETMYY100(:,1);
ETMYY100 = PEMEX_ETMYY100(:,2);
subplot(4,3,3)
plot(PEMEX, ETMYY100,'.')
title('PEM-EX vs. ISI-ETMY\_Y 100-300 mHz (r=0.3487)')
r8 = corr(PEMEX, ETMYY100)

load PEMEX_ETMYY300.txt
PEMEX = PEMEX_ETMYY300(:,1);
ETMYY300 = PEMEX_ETMYY300(:,2);
subplot(4,3,4)
plot(PEMEX, ETMYY300,'.')
title('PEM-EX vs. ISI-ETMY\_Y 0.3-1 Hz (r=0.2403)')
r9 = corr(PEMEX, ETMYY300)

load PEMEX_ETMYZ100.txt
PEMEX = PEMEX_ETMYZ100(:,1);
ETMYZ100 = PEMEX_ETMYZ100(:,2);
subplot(4,3,5)
plot(PEMEX, ETMYZ100,'.')
title('PEM-EX vs. ISI-ETMY\_Z 100-300 mHz (r=0.4724 )')
r10 = corr(PEMEX, ETMYZ100)

load PEMEX_ETMYZ300.txt
PEMEX = PEMEX_ETMYZ300(:,1);
ETMYZ300 = PEMEX_ETMYZ300(:,2);
subplot(4,3,6)
plot(PEMEX, ETMYZ300,'.')
title('PEM-EX vs. ISI-ETMY\_Z 0.3-1 Hz (r=0.6829)')
r11 = corr(PEMEX, ETMYZ300)

load PEMEX_ETMYRX100.txt
PEMEX = PEMEX_ETMYRX100(:,1);
ETMYRX100 = PEMEX_ETMYRX100(:,2);
subplot(4,3,7)
plot(PEMEX, ETMYRX100,'.')
title('PEM-EX vs. ISI-ETMY\_RX 100-300 mHz (r= 0.6986)')
r12 = corr(PEMEX, ETMYRX100)

load PEMEX_ETMYRX300.txt
PEMEX = PEMEX_ETMYRX300(:,1);
ETMYRX300 = PEMEX_ETMYRX300(:,2);
subplot(4,3,8)
plot(PEMEX, ETMYRX300,'.')
title('PEM-EX vs. ISI-ETMY\_RX 0.3-1 Hz (r=0.6860)')
r13 = corr(PEMEX, ETMYRX300)

load PEMEX_ETMYRY100.txt
PEMEX = PEMEX_ETMYRY100(:,1);
ETMYRY100 = PEMEX_ETMYRY100(:,2);
subplot(4,3,9)
plot(PEMEX, ETMYRY100,'.')
title('PEM-EX vs. ISI-ETMY\_RY 100-300 mHz (r=0.7469)')
r14 = corr(PEMEX, ETMYRY100)

load PEMEX_ETMYRY300.txt
PEMEX = PEMEX_ETMYRY300(:,1);
ETMYRY300 = PEMEX_ETMYRY300(:,2);
subplot(4,3,10)
plot(PEMEX, ETMYRY300,'.')
title('PEM-EX vs. ISI-ETMY\_RY 0.3-1 Hz (r=0.7075)')
r15 = corr(PEMEX, ETMYRY300)

load PEMEX_ETMYRZ100.txt
PEMEX = PEMEX_ETMYRZ100(:,1);
ETMYRZ100 = PEMEX_ETMYRZ100(:,2);
subplot(4,3,11)
plot(PEMEX, ETMYRZ100,'.')
title('PEM-EX vs. ISI-ETMY\_RZ 100-300 mHz (r=0.6442)')
r16 = corr(PEMEX, ETMYRZ100)

load PEMEX_ETMYRZ300.txt
PEMEX = PEMEX_ETMYRZ300(:,1);
ETMYRZ300 = PEMEX_ETMYRZ300(:,2);
subplot(4,3,12)
plot(PEMEX, ETMYRZ300,'.')
title('PEM-EX vs. ISI-ETMY\_RZ 0.3-1 Hz (r=0.7167)')
r17 = corr(PEMEX, ETMYRZ300)
