Learn R Programming

OmicKriging (version 1.0)

computeGRM: Compute the genetic relatedness matrix

Description

This function calls the software GCTA to calculate the genetic relationship matrix from all autosomal SNP genotype data in PLINK binary PED file format.

Usage

computeGRM(bedfullheader, grmfullheader = "genomic", 
gctaname = "gcta", idfile = NULL, snpfile = NULL)

Arguments

bedfullheader
path to PLINK binary .bed/.bim/.fam files
grmfullheader
output path and file name
gctaname
path to gcta executable
idfile
path to a file listing the subset of individuals to include in GRM calculation (2 columns labeled: FID, IID)
snpfile
path to a file listing a subset of SNPs, one per line and no column header, to included in the GRM calculation, default is to include all autosomal SNPs in the .bim file

Value

    References

    Wheeler HE, et al. (2013) Poly-Omic Prediction of Complex Traits: OmicKriging. arXiv:1303.1788 http://arxiv.org/abs/1303.1788

    Yang J, Lee S, Goddard M, Visscher P (2011) GCTA: a tool for genome-wide complex trait analysis. American Journal of Human Genetics 88: 76-82. http://www.complextraitgenomics.com/software/gcta/

    Purcell, S. et al. (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics 81:559-75. http://pngu.mgh.harvard.edu/~purcell/plink/