now, uses prepare_beds_and_contrast.R