Question: Intepretation of Z scores for gene set enrichment analysis in MAST
gravatar for siajunren
16 months ago by
siajunren0 wrote:

I performed zlmfit on a matrix with MAST for 2 populations of cells, A and B, with population A as intercept. Then I run gseaAfterBoot on the zlfit object after performing bootstrapping with bootVcov1. If I get a positive value Z score for a particular gene set, does it mean that genes in this geneset are upregulated in population A compared to B? Similarly, a negative Z score mean that genes in this geneset are downregulated in A compared to B. Or is it the other way round? Not sure of the directionality here, because I am not sure what the model coefficients that were used to calculate the Z score mean.

ADD COMMENTlink modified 15 months ago by Andrew_McDavid150 • written 16 months ago by siajunren0
gravatar for Andrew_McDavid
15 months ago by
Andrew_McDavid150 wrote:

The model coefficients that are being tested for enrichment are the coefficients from the zlmfit objects. With A as a ref, then positive coefficients in the ZlmFit object means that the mean of B > mean of A. The GSEA is calculating the average ZlmFit coefficient inside the set. So a positive Z-score in gseaAfterBoot means that B > A, on average, in the gene set.


## "1.3.0"
vb1 = subset(vbetaFA, ncells==1)
vb1 = vb1[,freq(vb1)>.1][1:15,]
zf = zlm(~Stim.Condition, vb1)
coef(zf, 'D') #discrete components

boots = bootVcov1(zf, 5)
sets = list(A = 1:2, #positive coefficients 
B = 3:6) #negative
gsea = gseaAfterBoot(zf, boots, sets, CoefficientHypothesis('Stim.ConditionUnstim'))
calcZ(gsea)[,,'Z'] #zscore is positive for A, negative for B
summary(gsea) #pretty-printer
ADD COMMENTlink written 15 months ago by Andrew_McDavid150
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 370 users visited in the last hour