read.bed: Reading data from binary PLINK files

In zhangh12/SCAT: Summary Based Conditional Association Test

Description

Loads genotype data from PLINK format files .bed, .bim, and .fam.

Usage

read.bed(bed, bim, fam, sel.snps = NULL, sel.subs = NULL, encode012 = TRUE)

Arguments

bed	the name of the bed file.
bim	the name of the bim file.
fam	the name of the fam file.
sel.snps	a character vector of SNPs to be extracted from the plink files. The default is NULL, i.e., all SNPs are extracted.
sel.subs	an optional character vector specifying a subset of subject IDs to be extracted from the plink files. These IDs should be matched with the second column of fam files. The default is NULL, i.e., all subjects are extracted.
encode012	logical. Encoding the genotypes using 0/1/2 if TRUE, or using symbols of the reference and effect alleles if FALSE. The default is TRUE.

Value

A data frame of genotypes of specified subjects in the plink files.

Examples

library(SCAT)

# Load the sample data

bed <- system.file("extdata", package = 'SCAT', 'chr5.bed')
bim <- system.file("extdata", package = 'SCAT', 'chr5.bim')
fam <- system.file("extdata", package = 'SCAT', 'chr5.fam')

## first five SNPs
b <- read.table(bim, header = FALSE, as.is = TRUE, nrows = 5)
## first 50 subjects
f <- read.table(fam, header = FALSE, as.is = TRUE, nrows = 50)
geno <- read.bed(bed, bim, fam, sel.snps = b[, 2], sel.subs = f[, 2])

dim(geno) # 50 x 5

zhangh12/SCAT documentation built on May 14, 2019, 5:02 a.m.