I propose adding a function, frequentist_to_bayesian(), to bayestestR to enable seamless conversion of frequentist models (e.g., lm, glm, lme4, gam) into their Bayesian equivalents using packages like brms. This function would complement the existing bayesian_as_frequentist() by providing a two-way workflow between frequentist and Bayesian methods. It should extract the formula, data, and family from the input model and fit a Bayesian counterpart with default or user-specified priors. This would simplify Bayesian exploration for users and allow direct comparison of frequentist and Bayesian results, enhancing the versatility of the package.
Please take into account that I can help in development if needed.