these are the 6 scripts needed to do the variance analysis (figure 4)