edgeR QL with voom weights
1
0
Entering edit mode
@gregorylstone-12225
Last seen 6.2 years ago

I was wondering if it's possible to incorporate the weights calculated by voomWithQualityWeights() in glmQLFit in edgeR. My code is as follows:

v <- voomWithQualityWeights(dge, design, plot=FALSE)
qlfit_weigthed <- glmQLFit(dge$counts, design, weights = v$weights, robust=TRUE, dispersion = dge$trended.dispersion)
lrt_weighted <- glmQLFTest(qlfit_weigthed, contrast=cont)

Adding the weights to glmQLFit gives me several more significant genes, so I want to make sure what I am doing is correct. My hesitation is with which counts to feed into glmQLFit. Should I be entering normalized counts or raw counts? I've read through the reference manual and I couldn't find a specification.

Additionally, is this worth pursuing, or are the weights provided by voomWQW and the weights exprected by glmQLFit not the same type?

Any help would be greatly appreciated. Thanks!

edger limma-voom voom • 1.4k views
ADD COMMENT
3
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 7 hours ago
The city by the bay

No, it doesn't make sense to use the voom weights in QL edgeR. The former describes the precision of the log-CPMs, while the weights in GLM scale the variances of the raw counts. They're totally different things.

You should be supplying the DGEList object to glmQLFit. You can supply the raw counts if you want, but then you need to pass the normalization factors via offsets. edgeR doesn't have a concept of normalized counts.

ADD COMMENT
0
Entering edit mode

Thank you for the clarification!
 

ADD REPLY

Login before adding your answer.

Traffic: 459 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6