# \donttest{
pedfile <- system.file("extdata",
"snp3kvars-CHR8-25947258-25951166-plink.ped",
package = "geneHapR")
mapfile <- system.file("extdata",
"snp3kvars-CHR8-25947258-25951166-plink.map",
package = "geneHapR")
p.link <- import_plink.pedmap(pedfile = pedfile, mapfile = mapfile,
sep_map = "\t", sep_ped = "\t")
p.link <- filter_plink.pedmap(p.link, mode = "POS",
Chr = "chr08", start = 25948004,
end = 25949944)
hapResult <- plink.pedmap2hap(p.link, hapPrefix = "H",
hetero_remove = TRUE,
na_drop = TRUE)
# }
Run the code above in your browser using DataLab