So I am looking for an automated approach to detect outliers in RNA-seq data. I usually looked at a PCA plot and decided visually. Now I would like to automate this. So I have been looking at the
PcaHubert() function in
rrcov package, which then flags suspected outliers as false:
pca <- PcaHubert(t(rna.data)) outliers <- which(pca@flag=='FALSE')
Would this be a good option? Or are there others better suited for RNA-seq data?
Thanks for your input!