CommentsWe appreciate the opportunity to reply to Bolstad's (2017) comment on our paper, "Quantifying nonadditive selection caused by indirect ecological effects" (terHorst et al. 2015). We respectfully disagree with Bolstad's argument that our method does not properly quantify nonadditive selection in response to indirect ecological effects, as it certainly does for many of the biological scenarios we envision. However, as with all tests related to ecological indirect effects, the appropriateness of our approach depends on the underlying model that best describes the multispecies interaction. Specifically, the null model we propose applies only to fitness effects of multispecies interactions that are best modeled with additive rather than multiplicative models. The lntransformations Bolstad proposes are appropriate for constructing null models of nonadditive selection when the underlying species interactions are best modeled multiplicatively. The use of additive vs. multiplicative models has a long history in the study of indirect ecological effects (Case and Bender 1981, Billick and Case 1994, Wootton 1994b, which we revisit briefly here and extend to measures of selection. Specifically, we clarify our goals and approach and address Bolstad's two major criticisms of our paper: (1) that we incorrectly infer indirect ecological effects from our measure of the non-additivity of selection, and (2) that log fitness (in contrast to fitness relativized across treatments as we propose) is more appropriate for quantifying selection in response to multiple species over varying timescales. We explain our reasoning and use a simple model to demonstrate when our null model is appropriate and when Bolstad's log fitness model is most appropriate for identifying nonadditive selection in multispecies communities.First, to clarify the goals of our original paper, we provided a method to detect and estimate the strength of nonadditive selection, which occurs when selection on a trait imposed by one selective agent is altered in the presence of a second selective agent. Our goal was not to use evidence for nonadditive selection to infer indirect ecological effects; ecologists have many simpler methods available for detecting indirect ecological effects (Strauss 1991, Wootton 1994a, Menge 1995 that do not require the large sample sizes and additional trait measurements needed to understand natural selection. Most studies estimate the consequences of indirect ecological effects for species abundances, rather than fitness or selection (Miller and terHorst 2012), and rarely measure traits. Indirect ecological effects do not necessarily result in non-additivity of selection, and nonadditive selection can result from processes other than indirect ecological effects. For example, an indirect ecological effect that alters fitness will not necessarily also result in nonadditivity of selection on the traits under consideration. Rather, for any species to alter selection on another it must alter the covariance between relative fitness and a ...