Entering edit mode
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)