function out=handMetaEyeColorAnalysis(in)

eyeColors={'++','+m','mm','+-','m-','--'};

numReps=1000;
out=zeros(length(eyeColors),3,numReps);

for n=1:numReps
    
    for i=1:length(eyeColors)
        which=find(strcmp(in.eyeColor,eyeColors{i}));
        N=length(which);
        whichFlies=ceil(rand(N,1)*N);
        mu=nanmean(in.handedness(which(whichFlies)));
        std=nanstd(in.handedness(which(whichFlies)));
        out(i,1,n)=std/mu;
        mu=nanmean(in.numTurns(which(whichFlies)));
        std=nanstd(in.numTurns(which(whichFlies)));
        out(i,2,n)=std/mu;
        mu=nanmean(in.switchiness(which(whichFlies)));
        std=nanstd(in.switchiness(which(whichFlies)));
        out(i,3,n)=std/mu;
    end
    
end