I have checked the PCA from my sRNAseq data and there is one replicate which goes over 67% away from all my samples (1/18).
I know that the PCA plots the genes with the top variance, and this sample seems to not have crazy values in all the genes (checked random), and with PCA of different biotype class genes (miRNAs and snoRNAs)
Pvars <- rowVars(assay(rld))
topVarGenes <- head(order(rowVars(assay(rld)), decreasing = TRUE),35)
with this I found the genes with the highest variance, but I do not seem to identify the genes that are causing the "crazy" behavior of this sample.
I do not want to throw away the sample because I only have 3 replicates, and the behavior seems to be caused by a few genes. Also the cooks plot shows it as the other samples.
Any suggestion on how to find the genes that are causing the extreme behavior in one replicate?