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?
(PS this was first posted on R-help but I was recommended I post here as well).
