%SandersvHurlbert.m
% Sanders method based on a misapplication of Theorem 3.5.2
load GEEP
DATA=DATA';
[raren,ESN]=rarefy(DATA,10);
H1=htrans(DATA,1);
SandersH10=H1*10;i=find(SandersH10>=1);li=length(i);
SandersH10(i)=ones(1,li);
SES10=sum(SandersH10')';
H10=htrans(DATA,10);ES10=sum(H10')';
H10(i)
plot(ES10,SES10,'+',[1 8.5]',[1 8.5]','-r');
xlabel('Hurlbert''s E(S_{10}) for GEEP data');
ylabel('Sanders'' ad hoc E(S_{10})')
title('Sanders'' (1968) error in calculating E(S_n)')
figure(gcf)
