For extracting the genotype counts from the original raw data from PLINK.
1 | tidy_hwdata(file, flip = TRUE)
|
1 2 3 4 5 6 7 8 9 10 11 12 | ## Not run:
# this is how tidy versions were created (run in package dir)
files <- sapply(RAWFILES, eve102_data, include_raw=TRUE)
for (i in seq_along(files)) {
outfile <- file.path("inst", "extdata", paste0(names(files)[i], ".txt.gz"))
outcon <- gzfile(outfile, "w")
write.table(tidy_hwdata(file.path("inst", "extdata", basename(files[i]))),
file=outcon, sep="\t", quote=FALSE, row.names=FALSE)
close(outcon)
}
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.