Skip to contents

Perform Single Marker Testing (SMT) on a particular variant

Usage

TestMarker(x, g, add.noise = FALSE)

Arguments

x

a null model object as returned by FitNull

g

a genotype vector

add.noise

if 'raw' return raw residuals, if FALSE simply rank normalize residuals, if TRUE add structured noise along the subspace spanned by the background covs and genotype so that the residuals are 'full rank': truly from N(0,I), not a curved Gaussian.

Value

a list with three elements

  • `p-value`: single marker test p-value

  • `y`: the residual phenotype vector after testing the genotype

  • `Q`: an orthogonal matrix with columns spanning the column space of the tested genotype and the background covariates in the null model, useful for further testing