function out=handMetaLeftyAnalysis(in) data=in.handedness; numReps=100; %grand mean analysis h=waitbar(0,'progress'); grandMeanVals=zeros(numReps,1); for i=1:numReps whichFlies=ceil(rand(size(in,1),1)*size(in,1)); grandMeanVals(i)=nanmean(data(whichFlies)); waitbar(i/numReps,h); end close(h); out.grandMeanVals=grandMeanVals; experimenters=unique(in.experimenterID); out.experimenterMeanVals=zeros(size(experimenters,1),2); for i=1:length(experimenters) whichFlies=find(strcmp(experimenters{i},in.experimenterID)); out.experimenterMeanVals(i,1)=nanmean(data(whichFlies)); out.experimenterMeanVals(i,2)=length(whichFlies); end figure; violinPlot(out.experimenterMeanVals(:,1)) scatter(1+0.1*randn(length(experimenters),1),out.experimenterMeanVals(:,1),40,[zeros(length(experimenters),1) out.experimenterMeanVals(:,2) out.experimenterMeanVals(:,2)]/50000); genotypes=unique(in.genotype); out.genotypeMeanVals=zeros(size(genotypes,1),2); for i=1:length(genotypes) whichFlies=find(strcmp(genotypes{i},in.genotype)); out.genotypeMeanVals(i,1)=nanmean(data(whichFlies)); out.genotypeMeanVals(i,2)=length(whichFlies); end figure; violinPlot(out.genotypeMeanVals(:,1)) scatter(1+0.1*randn(length(genotypes),1),out.genotypeMeanVals(:,1),40,[zeros(length(genotypes),1) out.genotypeMeanVals(:,2) out.genotypeMeanVals(:,2)]/800); sexes=unique(in.sex); out.sexMeanVals=zeros(size(sexes,1),2); for i=1:length(sexes) whichFlies=find(strcmp(sexes{i},in.sex)); out.sexMeanVals(i,1)=nanmean(data(whichFlies)); out.sexMeanVals(i,2)=length(whichFlies); end figure; violinPlot(out.sexMeanVals(:,1)) scatter(1+0.1*randn(length(sexes),1),out.sexMeanVals(:,1),40,[zeros(length(sexes),1) out.sexMeanVals(:,2) out.sexMeanVals(:,2)]/100000);