function out=handMetaSexAndTempVariability(in,varOfInterest)

eval(['data=in.' varOfInterest ';']);

genotypes={'trp';'shi';'r73d06';'r16d01';'np2320';'mb247gal80';'c819';...
    'c465';'c232';'dgrp45';'dgrp535'};

hotIndices=(double(in.expTemp==32) + double(in.expTemp==33))>0;
coolIndices=(double(in.expTemp==22) + double(in.expTemp==23))>0;

out=zeros(length(genotypes),2);

for i=1:length(genotypes)
    tempIndices=strcmp(in.genotype,genotypes{i});

    tempHotIndices=hotIndices.*tempIndices==1;
    tempCoolIndices=coolIndices.*tempIndices==1;

    out(i,1)=nanstd(data(tempCoolIndices))/nanmean(data(tempCoolIndices));
    out(i,2)=nanstd(data(tempHotIndices))/nanmean(data(tempHotIndices));
    
end


figure;
hold on;

for i=1:size(out,1)
    plot([1 2],[out(i,1) out(i,2)],'k.-');
end
xlim([0 3])