Read plink output files

Share:

Description

Read plink output files

Usage

1
readPlinkOut(filename, ...)

Arguments

filename

Filenames of plink output files, see collenv$.plink_out_ext

...

passed to one of readAssoc, readQassoc, readLinear, readLogistic

Value

data.frame

Author(s)

Kaiyin Zhong, Fan Liu

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
## Not run: 
dat1 = readPlinkOut("assoc/mmp13.assoc")
dat2 =    readAssoc("assoc/mmp13.assoc")
all(na.omit(dat1 == dat2))
dat1 = readPlinkOut("assoc/mmp13.assoc", c("CHR", "SNP", "P", "OR"))
dat2 =    readAssoc("assoc/mmp13.assoc", c("CHR", "SNP", "P", "OR"))
all(na.omit(dat1 == dat2))
dat1 = readPlinkOut("assoc/mmp13.qassoc")
dat2 =   readQassoc("assoc/mmp13.qassoc")
all(na.omit(dat1 == dat2))
dat1 = readPlinkOut("assoc/mmp13.qassoc", c("CHR", "SNP", "P", "R2"))
dat2 =   readQassoc("assoc/mmp13.qassoc", c("CHR", "SNP", "P", "R2"))
all(na.omit(dat1 == dat2))

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.