I performed a surrogate variable analysis on my RNA-seq data using the SVA package in R and want to use the identified surrogate variables as covariates in my DE design formula. There were more than 5 significant SV identified - human samples from different locations and different timepoints, hence high heterogeneity and several confounding factors. The output of sva is anyway limited to the first 5 SV identified and I was unable to successfully increase the number of sv included in the svaseq output file. I am using the code given below:
I couldn't find any helps/options adapting the output in the vignette/R help/forum, hence the question. Thanks in advance for any help!