metaGE.fit | R Documentation |
Quantitative trait loci detection via Fixed or Random effect meta-analysis GWAS procedure.
metaGE.fit(Data, MatCorr, Method, NA.omit = TRUE, DropZScores = FALSE)
Data |
A dataset containing the estimated marker effect and its associated pvalue of each marker (in rows) in each environment (in columns), as obtained from |
MatCorr |
The inter-environments correlation matrix. Can be computed using |
Method |
A string specifying the method to be performed: either " |
NA.omit |
A boolean specifying whether the markers with some |
DropZScores |
A boolean specifying whether the Zscores should be dropped from the dataset or not.( |
Different tests may be performed:
Fixed Effect (Fe), to identify markers with a stable effect across environments.
Random Effect (Re), to identify markers whose effects may be unstable across environments.
The dataset Data with supplementary columns:
PVALUE: The pvalue of the MA test,
Mu: Estimate of the mean marker effect,
Tau: Estimate of the variance of the marker effect, for the Random model only,
the Zscores for each environment if DropZScores = FALSE
.
require(dplyr)
# Import the data
data("metaData")
# Compute the inter-environment correlation matrix
matCorr <- metaGE.cor(metaData, Threshold = 0.8)
# Fixed Effect
FeDF <- metaGE.fit(metaData, matCorr, Method = "Fe")
head(FeDF %>% select(CHR, POS, MARKER, Mu, Tau, PVALUE))
# Random Effect
ReDF <- metaGE.fit(metaData, matCorr, Method = "Re")
head(ReDF %>% select(CHR, POS, MARKER, Mu, Tau, PVALUE))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.