I am trying to simulate the residuals of a model that includes orthogonal polynomial terms with DHARMa:
cpr<-glmmTMB(cp ~ poly(zrc,2)*poly(zoc,2) + site +