ebayes_EM {ScoreEB}

Empirical Bayes for multi-locus selection

Description

Empirical Bayes using expectation–maximization algorithm.

Usage

ebayes_EM(x,z,y,EMB.tau,EMB.omega)

Arguments

x

fixed effect vector or matrix.

z

genotype data.

y

phenotype data.

EMB.tau

one of hyperparameters in inverse chi-square distribution.

EMB.omega

one of hyperparameters in inverse chi-square distribution.

Value

u

The effect values of markers, and their absolute values are used as the basis for further screening.

Examples

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)