allDatafullRuns.R control file, row 30 where marginals is defined, should probably be marginals[,-1] instead of marginals[1, -1].