RESUMO
BACKGROUND: For investigating the individual-environment interplay and individual differences in response to environmental exposures as captured by models of environmental sensitivity including Diathesis-stress, Differential Susceptibility, and Vantage Sensitivity, over the last few years, a series of statistical guidelines have been proposed. However, available solutions suffer of computational problems especially relevant when sample size is not sufficiently large, a common condition in observational and clinical studies. METHOD: In the current contribution, we propose a Bayesian solution for estimating interaction parameters via Monte Carlo Markov Chains (MCMC), adapting Widaman et al. (Psychological Methods, 17, 2012, 615) Nonlinear Least Squares (NLS) approach. RESULTS: Findings from an applied exemplification and a simulation study showed that with relatively big samples both MCMC and NLS estimates converged on the same results. Conversely, MCMC clearly outperformed NLS, resolving estimation problems and providing more accurate estimates, particularly with small samples and greater residual variance. CONCLUSIONS: As the body of research exploring the interplay between individual and environmental variables grows, enabling predictions regarding the form of interaction and the extent of effects, the Bayesian approach could emerge as a feasible and readily applicable solution to numerous computational challenges inherent in existing frequentist methods. This approach holds promise for enhancing the trustworthiness of research outcomes, thereby impacting clinical and applied understanding.