require(dplyr)
# Import the data
data("metaData")
# Compute the inter-environment correlation matrix
matCorr <- metaGE.cor(metaData, Threshold = 0.8)
# Fit the Fixed Effect model
FeDF <- metaGE.fit(metaData, matCorr, Method = "Fe")
# Control the FDR (here Benjamini-Hochberg)
Alpha <- 0.05
Signif <- FeDF$PVALUE %>% p.adjust(method = "BH") %>% `<`(Alpha) %>% which
# Draw the z-scores heatmap of the significant markers
heatmap <- metaGE.heatmap(Data = FeDF[Signif,],
Prefix = "Z.")
Run the code above in your browser using DataLab