I find some protocol will recommend using
pbmc <- ScaleData(pbmc, vars.to.regress = "percent.mt")
to regress out cell cycle influence or other "noise". But I am confused why we can not drop these genes in
feature selection steps. For example, we can drop these genes from high variable gene set, so that protecting our latter PCA and cluster analysis from these influence genes.