Use fread() instead of read.delim() for speed