Hi Dr. Martinez, I am following your tutorial. I was wondering if for the PCA should we use: pca <- prcomp(scale_data, scale=FALSE) #instead of TRUE since each gene was previously scaled across samples? Thank you, Paula