gave up on data.table, not sure on fread() arguments - not that many rows anyway.
summarised by mean and standard dev