function out=handMetaGenotypeVariabilityAnalysis(in) data=in.handedness; minN=70; validTemps=in.expTemp==22 | in.expTemp==23; validConds=strcmp('normal',in.expCond); validIndices=validTemps.*validConds; genotypes=unique(in.genotype); out.dgrpVars=[]; out.nondgrpVars=[]; for i=1:length(genotypes) genoIndices=strcmp(genotypes{i},in.genotype); if isequal(genotypes{i}(1:min([4 length(genotypes{i})])),'dgrp') if sum(validIndices.*genoIndices)>=minN out.dgrpVars=[out.dgrpVars;nanstd(data(validIndices.*genoIndices==1))/nanmean(data(validIndices.*genoIndices==1))]; end else if sum(validIndices.*genoIndices)>=minN out.nondgrpVars=[out.nondgrpVars;nanstd(data(validIndices.*genoIndices==1))/nanmean(data(validIndices.*genoIndices==1))]; end end end