Changed SNP dosage input to read.csv(file, row.names = 1) because we expect a matrix as input here with rownames. load_file function doesn't work properly here yet