SKAT_NULL_emmaX: Get parameters and residuals from the null model with kinship...

In ttimbers/SKAT: SNP-Set (Sequence) Kernel Association Test

Description

Compute model parameters and residuals for SKAT with incorporating kinship structure.

Usage

SKAT_NULL_emmaX (formula, data=NULL, K=NULL, 
Kin.File=NULL, ngrids=100, llim=-10, ulim=10, esp=1e-10) 
  
 

Arguments

formula	an object of class “formula”: a symbolic description of the NULL model to be fitted.
data	an optional data frame containing the variables in the model (default=NULL). If it is NULL, the variables are taken from 'environment(formula)'
K	a kinship matrix. If K=NULL, the function reads Kin.File to get a kinship matrix.
Kin.File	an emmax-kin output file name. If K=NULL, the function reads this file.
ngrids	Number of grids to search optimal variance component
llim	Lower bound of log ratio of two variance components
ulim	Upper bound of log ratio of two variance components
esp	Tolerance of numerical precision error

Details

The Emma package code was used to implement this function.

Resampling is not implemented.

Value

This function returns an object that has model parameters and residuals of the NULL model of no association between genetic variables and outcome phenotypes. After obtaining it, use SKAT function to carry out the association test.

Author(s)

Seunggeun Lee

Examples

data(SKAT.fam.example)
attach(SKAT.fam.example)


obj<-SKAT_NULL_emmaX(y ~ X, K=K)
SKAT(Z, obj)$p.value
	

ttimbers/SKAT documentation built on June 1, 2019, 2:48 a.m.