test.peak: Likelihood ratio tests for QTL peaks

In steibelj/gwaR: Functions for performing GWA from GBLUP

Description

function for performing likelihood ratio test h_peak^2=0 of QTL peak using a gblup model

Usage

test.peak(gb, x, peak_pos = NULL)

Arguments

gb	a gblup model
x	A matrix of standardized genotypes with animals in columns and markers in rows
peak_pos	an index vector used to select the rows (markers) of x that capture the QTL peak. if NULL: the whole x matrix is used.

Value

a list with elements

• pvalue the pvalue of the LRT

• llik Log-likelihood of full model, reduced model and their difference

• VARS a dataframe with variance component estimates under reduced and full model

See Also

gblup lrt

steibelj/gwaR documentation built on May 30, 2019, 10:45 a.m.