I have used RUVr to remove batch effect in my dataset. Selecting
k=10 looks good in my dataset after checking using PCA plots. Now I want to know which estimated factor (W) should be considered in my design matrix for DE analysis?
r <- RUVr(x, gene, k=10,res) Which W? W-1 or W-10 ? model.matrix(~ Condition + W)