Heterozygosity filter does not remove samples below -3SD

[mvochteloo]

In the filter_het.R script samples for which the proportion of non-homozygous sites below -3 or above +3 standard deviations are correctly identified als failed (het_fail). However, for the het_pass samples below -3 standard deviations are not filtered out.

het_fail <- subset(het, (het$HET_RATE < mean(het$HET_RATE) - 3*sd(het$HET_RATE)) | (het$HET_RATE > mean(het$HET_RATE) +3 * sd(het$HET_RATE)));
het_pass <- subset(het, (het$HET_RATE < mean(het$HET_RATE) - 3*sd(het$HET_RATE)) | (het$HET_RATE < mean(het$HET_RATE) +3 * sd(het$HET_RATE)));