@drjulie_b So one thing is to use `group=interaction(dph,otherthing)`, so that'll make combinations of the dph and otherthing. So that's splt into number of groups of #dph*#otherthing

Maybe inspecting `otherthing` to see what values are in there? Is a factor with levels that are now filtered out?