RUV with K more than 1.
1
0
Entering edit mode
star ▴ 10
@star-15676
Last seen 21 months 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 • 907 views
ADD COMMENT
0
Entering edit mode
ATpoint ★ 4.5k
@atpoint-13662
Last seen 1 day 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: 349 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