load Huybers07data.txt
figure
time = Huybers07data(:,1);
Benthic = Huybers07data(:,4);
Planktic = Huybers07data(:,5);

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,Benthic);
subplot(1,2,1);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,P95,'g');
title('Benthic')

[f,P,Pred_norm,Pwhite_norm,P95,P95chi2]=myspectrum(time,Planktic);
subplot(1,2,2);
loglog(f,P,'k'); hold on;
loglog(f,Pred_norm,'r');
loglog(f,Pwhite_norm,'b');
loglog(f,P95,'g');
title('Planktic');


