Hi All, I have some genetic data and I want to resample it changing the number of individuals and/or the number of markers to show the effects of smaller sample sizes on the pairwise-IBD calculations. The “standard” I need to follow is PLINK, specifically the --genome option. I could create multiple subsets of the data, export it an have PLINK do the needed calculations, and then reimport the results in R, but this is obviously pretty silly, especially because it would require creating, exporting and importing a few hundreds files. I could write the needed functions to do what PLINK does, but I’d rather not reinvent the wheel so it might be worth asking: is there any package/function that would do (in R) what PLINK does with --genome?
Bets wishes F -- Federico Calboli Ecological Genetics Research Unit Department of Biosciences PO Box 65 (Biocenter 3, Viikinkaari 1) FIN-00014 University of Helsinki Finland federico.calb...@helsinki.fi ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
