Question about sva correction: why does IsoformSwitchAnalyzeR tell me I have a lot of SVs but sva actually disagrees?

[bozbezbozzel]

Hi,

I'm following the steps in your vignette (great work by the way) for my data of about 100 NSCLC patients. With importRdata I keep receiving the warning though that sva finds too many sources of variation.

My data is fairly heterogeneous as you'd expect from patient samples from disease such as lung cancer, but I have a few things going for me namely that

• these are all resected tumors that received little or no pre-treatment
• I split them out by histological subtype since that is by far the largest source of variation in the gene level counts
• I added sex and tumor stage as covariates since I know from my other analyses that these do add a little bit of variation
• the quality of my samples is good and my regular RNA expression analyses are going well

Now out of interest I ran sva separately on my transcripts and gene-levels counts, imported with tximport and scaled in the case of transcripts. Using num.sv with be as the method (which is the default as far as I'm aware) gives me 3 SVs for the counts, 1 for the transcripts. That doesn't seem like it would warrant the warnings I keep getting.

I resorted to disabling sva with detectUnwantedEffects = FALSE but now I'm curious where this discrepancy could come from. I did also notice that the guesstimated dtu number is zero-- seems possible but not likely. I'm wondering if I'm doing something wrong and am just not noticing?