readPlinkToMatrixByIndex {seqminer} | R Documentation |
Read a gene from BGEN file and return a genotype matrix
readPlinkToMatrixByIndex(plinkFileObject, sampleIndex, markerIndex)
plinkFileObject |
a PlinkFileObject obtained by openPlink() |
sampleIndex |
integer, 1-basd, index of samples to be extracted |
markerIndex |
integer, 1-basd, index of markers to be extracted |
genotype matrix, marker by sample
http://zhanxw.com/seqminer/ for online manual and examples
## these indice are nonsynonymous markers for 1:196621007-196716634", ## refer to the readVCFToMatrixByRange() fileName = system.file("plink/all.anno.filtered.extract.bed", package = "seqminer") fileName = sub(fileName, pattern = ".bed", replacement = "") sampleIndex = seq(3) markerIndex =c(14, 36) cfh <- readPlinkToMatrixByIndex(fileName, sampleIndex, markerIndex)