not needed; sumstats.R is sufficient for present purposes
updated and now works - takes iterations, n, s as input and makes a nested list that is iterations long with the appropriate number of means and variances
moved unused files