powered by
Empirical Bayes using expectation<U+2013>maximization algorithm.
ebayes_EM(x,z,y,EMB.tau,EMB.omega)
fixed effect vector or matrix.
genotype data.
phenotype data.
one of hyperparameters in inverse chi-square distribution.
The effect values of markers, and their absolute values are used as the basis for further screening.
# NOT RUN { data(geno) data(pheno) EMB.tau <- 0 EMB.omega <- 0 z <- t(geno[,-c(1:4)]) y <- as.matrix(pheno) nsample <- dim(z)[1] x <- as.matrix(rep(1,nsample)) ebayes_EM(x,z,y,EMB.tau,EMB.omega) # }
Run the code above in your browser using DataLab