## Not run:
# mypheno <- read.table("MYPHENO.dat", header = TRUE,
# as.is = TRUE, na.strings = c("NA", "-9"))
# system("plink --bfile MYBED --freq --out MYEXPORT")
# system("plink --bfile MYBED --recodeA --out MYEXPORT")
# mydata <- read.snpdata.plink("MYEXPORT", phenotypes = mypheno)
# ## End(Not run)
Run the code above in your browser using DataLab