function out=handMetaReport(data) varList=[6 7 8 9 10 11 12 14 16 17 18]; out=cell(length(varList),2); for i=1:length(varList) [C,~,IC]=unique(data(:,varList(i))); numTypes=size(C,1); tally=zeros(numTypes,1); for j=1:numTypes tally(j)=sum(IC==j); end if isequal(class(C{1,1}),'double') C1=table2array(C); C2=C1; C2(isnan(C1))=[]; C2=[C2;NaN]; tally(isnan(C1))=[]; tally=[tally;sum(isnan(C1))]; report=array2table([C2 tally]); report.Properties.VariableNames(1)=C.Properties.VariableNames(1); report.Properties.VariableNames(2)={'tally'}; else report=[C array2table(tally)]; end if report{end,2}==0 report(end,:)=[]; end report=flipud(sortrows(report,2)); display(report); out{i,1}=data.Properties.VariableNames(varList(i)); out{i,2}=report; end num_flies=size(data,1) total_turns=nansum(data.numTurns) mean_handedness=nanmean(data.handedness) sem_handedness=nanstd(data.handedness)/sqrt(num_flies)