Entering edit mode
                    12 months 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!