search

HCBravoLab / metagenomeSeq

Statistical analysis for sparse high-throughput sequencing
64 stars 20 forks source link

ref1 in function cumNormStatFast() - Means or Medians? #63

Open andrebolerbarros opened 5 years ago

andrebolerbarros commented 5 years ago

I've been reviewing the function cumNormStatFast, in order to get the CSS normalization. When going through, I found that ref1 is determined as rowMeans; however, in the original formula, ref1 refers to q-l, defined as median of _l_th quantile across samples.

Is that supposed to be like this? Or am I just understanding it wrongly?

Thanks!

jnpaulson commented 5 years ago

Thanks! The function is correct.

Best,

On Tue, Jan 15, 2019 at 7:52 PM André Barros notifications@github.com wrote:

I've been reviewing the function cumNormStatFast, in order to get the CSS normalization. When going through, I found that ref1 is determined as rowMeans; however, in the original formula, ref1 refers to q-l, defined as median of _l_th quantile across samples.

Is that supposed to be like this? Or am I just understanding it wrongly?

Thanks!

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/HCBravoLab/metagenomeSeq/issues/63, or mute the thread https://github.com/notifications/unsubscribe-auth/AC5hgwdbtsnLC3xv1nvMJGnZG77j9WzNks5vDiNQgaJpZM4aBle3 .