comparing branches
corrected old-bw: had not been using varset.interacted for calculating degrees of freedom when using interacted data