I have very large matrices of genomic data (counts of nucleotides across dozens of individuals for millions of sites) that we process with an R script, and then write to fil