RPKM counts with SVA package & Combat with SVA?
0
0
Entering edit mode
Ahdee ▴ 50
@ahdee-8938
Last seen 28 days ago
United States

Hi everyone.  2 questions concerning the SVA batch correction package.

1. Can I use RPKM values? Calculated with edgeR rpkm() function with log2 transformed.  I know of svaseq() but that seem to only work with counts. 

2. In the manual the author states "It is also possible to use the sva function with the ComBat function to remove both known batch effects and  other  potential  latent  sources  of  variation." So how will this work.  Do I for example run combat first with say batch1, then run the sva with the output matrix?  Thereafter would I do something model.matrix(~0+test+sva_mod$sv+batch1)

 

thanks!

A

batch effect combat sva sva • 672 views
ADD COMMENT

Login before adding your answer.

Traffic: 329 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