RUV with K more than 1.
1
0
Entering edit mode
star ▴ 10
@star-15676
Last seen 1 day ago
Netherlands

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)
RUVSeq • 81 views
ADD COMMENT
0
Entering edit mode
ATpoint ★ 2.1k
@atpoint-13662
Last seen 16 hours ago
Germany

If you want to use k=10 then you need the first 10 RUV factors. Basically you would combine colData with that matrix that RUVseq returns and then include the columns (=factors) into the design. See for example this analysis from the DESeq2 developer using k=5:

https://github.com/mikelove/preNivolumabOnNivolumab/blob/main/preNivolumabOnNivolumab.Rmd#L112-L144

ADD COMMENT

Login before adding your answer.

Traffic: 289 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6