Hello everybody,

I aim to use the QSEA R package for the analysis of my MEDIP data. My question concerns the TMM normalisation which is applied by the function qsea::addLibraryFactors. Here, a reference can be either chosen manually or if e.g. ref=1 is set, then the first sample in the sampleTable is automatically chosen.

My data might show some variety, that's why I assume, I have to make a good choice for my reference sample. What are important criteria for the appropriate choice of a reference sample?

I had a look at the function edgeR::calcNormFactors, where there is the default option that "the library whose upper quartile is closest to the mean upper quartile for all the libraries" is chosen as the reference. I would prefer to use this option, as well, but struggle to integrate it into code for qsea::addLibraryFactors / qsea::estimateLibraryFactors.

Could somebody help me with this issue?

explaining the criteria to choose an appropriate reference sample or (even better) tell me how to integrate the code from edgeR::calcNormFactors into qsea:addLibraryFactors ?

Or do I get anything wrong and might don't have concern too much about the choice of the reference sample?

Thank you very much!

The help page for addLibraryFactors says this:

Arguments:

      qs: The qseaSet object

 factors: In case normalization factors have been pre-computed by the
          user, they can be passed with this parameter.  In this case
          QSEA adds this factors to the qseaSet object and does not
          compute normalization factors.

Where the 'factors' argument clearly provides the opportunity to compute the TMM normalization factors yourself and then just add to the qseaSet object. So if you prefer the choice that computeNormFactors makes, it seems easy enough to do?