function out=handMetaMoments(in) in=double(in); in(isnan(in))=[]; numReps=1000; numFlies=size(in,1); numMoments=20; out=zeros(numReps,numMoments,5); h=waitbar(0,'progress'); for i=1:numReps whichTemp=ceil(rand(numFlies,1)*numFlies); X=in(whichTemp); mu=mean(X); for j=3:numMoments % out(i,:)=[mean(dataTemp) var(dataTemp) skewness(dataTemp) kurtosis(dataTemp)]; out(i,j,1)=mean((X-mu).^j)/mean((X-mu).^2).^(j/2); end out(i,1,1)=mean(X); out(i,2,1)=var(X); waitbar(i/numReps,h); end for n=[4 3 2 1] for i=1:numReps whichTemp=ceil(rand(10^n,1)*numFlies); X=in(whichTemp); mu=mean(X); for j=3:numMoments out(i,j,6-n)=mean((X-mu).^j)/mean((X-mu).^2).^(j/2); end out(i,1,6-n)=mean(X); out(i,2,6-n)=var(X); waitbar(i/numReps,h); end end close(h);