likelihood {ScoreEB}

Carry out likelihood ratio test

Description

Snps selected via EM-Bayes to further identified by likelihood ratio test.

Usage

likelihood(xxn,xxx,yn,bbo)

Arguments

xxn

fixed effect vector or matrix.

xxx

snp matrix which are selected by EM-Bayes.

yn

phenotype data.

bbo

effect value of snp estimated by EM-Bayes.

Value

lod

Odds of logarithm vector of markers.

Examples

data(geno)
data(pheno)
z <- t(geno[,-c(1:4)])
y <- as.matrix(pheno)
n.sample <- dim(z)[1]
m.marker <- dim(z)[2]
x <- as.matrix(rep(1,n.sample))
beta <- rnorm(m.marker)
likelihood(x,z,y,beta)