I am interested in understanding if two predictors variables are associated with treatment outcome (DV). I included two covariates that 1) differed across my levels of IV (age) and 2) another that is strongly associated with the DV (baseline functioning). I ran a simultaneous multiple regression with the covariates and both IVs entered into the same time to estimate the unique relationships between the IVs with the DV.
A reviewer implied that hierarchical regression was needed (enter covariates into step 1 and IVs into step 2. My understanding is this would be asking a different question and that it is totally appropriate to include covariates in a single step in simultaneous regression.
Another reviewer suggested that only ANCOVA was appropriate to address this question. I'm pretty confident that this is a consequence of the reviewer being trained in ANOVA but not regression models, and is a misunderstanding.
Welcoming feedback and any academic citations to support the statistical approach.