% Ruth's Gel Permeation data for Winter Flounder Lipovitellin.
% A difference plot of Native vs Heat-stable fractions was
% calculated using a cubic spline interpolation function in
% a Matlab script.
clear;
load hslv.prn;
M=hslv;
X = M(:,1);
Y = M(:,2);
Z = M(:,3);
W = M(:,4);
figure(1);
clf;
subplot(3,1,1);
plot(X,W,'y:');
[Nhslv dum]=size(M);
Ns=11;
nhslv = Nhslv-Ns;
Vo=143;
Vt=401;
for i = 1:nhslv,
  yh(i)=W(i+Ns);
  xh(i)=X(i+Ns);
  Kh(i)=(xh(i)-Vo)/(Vt-Vo);
end
subplot(3,1,2);
plot(xh,yh,'y:');
xhs = xh(1):2:xh(nhslv);
[ dum nspli] = size(xhs);
for i = 1:nspli,
  Kspl(i)=(xhs(i)-Vo)/(Vt-Vo);
end
 
yhs = spline(xh,yh,xhs);
subplot(3,1,3);
plot(xhs,yhs,'y:');


load nlv.prn;
M=nlv;
a=0;
b=2.448;
c=0.89;
s=0.8;
X = M(:,1);
Y = M(:,2);
subplot(3,1,1);
hold on;
plot(a+b*(X+s),c*Y,'r-'), hold off;
axis([0 480 -0.05 0.45]);
text(300, 0.35, 'raw data aligned');
[Nnlv dum]=size(M);
Ns=10;
nnlv = Nnlv-Ns;
for i = 1:nnlv,
  yn(i)=c*Y(i+Ns);
  xn(i)=a+b*(X(i+Ns)+s);
  Kn(i)=(xn(i)-Vo)/(Vt-Vo);
end
subplot(3,1,2);
hold on;
plot(xn,yn,'r-'), hold off;
axis([100 480 -0.05 0.45]);
text(300, 0.35, 'truncated data set');

yns = spline(xn,yn,xhs);
subplot(3,1,3);
hold on;
plot(xhs,yns,'r-'), hold off;
axis([100 480 -0.05 0.45]);
text(300, 0.35, 'splined data/common X');
xlabel('Elution Volume');

figure(2);
clf;
subplot(1,1,1);
plot(xhs,yns-yhs,'y-'), hold on;
plot(xn,yn,'r:'), hold off;
figure(3);
clf;
subplot(1,1,1);
plot(Kh,yh,'y-'), hold on;
Vn=0.45;
Vhs=0.533;
yV=[0 0.35];
yV0=[0 0.28];
xVo=[0 0];
xVt=[1 1];
xVn=[Vn Vn];
xVhs=[Vhs Vhs];
plot(xVo,yV,'w:');
plot(xVt,yV0,'w:');
plot(xVn,yV,'w:');
plot(xVhs,yV,'w:');
plot(Kspl,(yns-yhs),'b--');
plot(Kn,yn,'r:');
% Labeling
xkey=[0.8 0.98];
yn=[0.38 0.38];
yhs=[0.355 0.355];
yhl=[0.33 0.33];
plot(xkey,yn,'y-');
plot(xkey,yhs,'r:');
plot(xkey,yhl,'b--'), hold off;
text(0.42, 0.36,'Khl');
text(0.52, 0.36,'Khs');
text(1, 0.38,'n: native');
text(1, 0.355,'hs: heat-stable');
text(1, 0.33,'hl: heat-labile');
text(0.9, 0.305,'hl = n - 1.12 hs');
xlabel('Ke');
ylabel('A280');
axis([-0.1 1.3 -0.01 0.4]);