Bear with me, I’m new at this.

How do I calculate or extract the p value of coefficients in my GLM in edgeR? I want to gauage empirically the contribution of each of the coefficients to the model

For reference, my design is:

Timepoint <- factor(meta.data_sub$Group3, levels=c("0","1"))
Age <- as.numeric(meta.data_sub$Age)
Sex <- factor(meta.data_sub$Sex2)
Ethnicity <- factor(meta.data_sub$Eth2)
design <- model.matrix(~Age +Sex +Ethnicity +Timepoint)

And I’m using QLFTest

dge <- estimateDisp(dge, design = design)
fit <- glmQLFit(dge, design = design,dispersion=dge$common.dispersion, contrast = contrast)
qlf <- glmQLFTest(fit)