Entering edit mode
4 weeks ago
IM
•
0
Dear all,
Could SVA be executed in the following way?
mod <- model.matrix(~1, data=as.data.frame(t(annot_)))
sva_fit <- sva(data_, mod, method="two-step")
I found this approach in a paper (https://doi.org/10.1038/s41591-019-0457-8), but I'm not entirely sure how SVA works with the "intercept-only" model. Any insights would be appreciated!