Minimum Detectable Effect Size for multiple linear regression I am trying to calculate the minimum detectable effect size (MDE) of an explanatory variable in a multiple linear regression.
The regression looks like the following:
$$
y_i = \beta_0 + \beta_1X_1 + \beta_2*X_2 + \beta_3*X_3 + \beta_4*X_4 + \epsilon_i
$$
whereby I am only interested in the MDE of $X_1$ in this multiple linear regression.
I have found code for linear regression in R: pwr.f2.test().
However, I am not sure if this can eliminate the minimum detectable effect size for one variable. Does anyone know how to get the minimum detectable effect size of one explanatory variable in a multiple linear regression? How can I implement this in R?

Here is a useful formula to remember:
$$\pm \beta_{j}^{a}=\frac{\left(z_{1-\alpha / 2}+z_{\gamma}\right) \sigma_{y \mid \mathbf{x}}}{\sigma_{x_{j}} \sqrt{n\left(1-\rho_{j}^{2}\right)}}$$
Here, $\beta_j^\alpha$ is the minimal detectable effect for a covariate assuming:

*

*Our FPR is $\alpha$ and our desired power is $\gamma$.  If you use the popular $\alpha=0.05$ and $\gamma=0.8$ then $z_{a-\alpha/2} \approx 1.96$ and $z_\gamma \approx 0.84$.


*The residual variation is $\sigma_{y\vert x}$.  This is akin to the standard deviation of the residuals assuming the linear model is the true model.


*The standard deviation of the covariate is $\sigma_{x_j}$.  If the covariate is continuous, this can always be rescaled to 1 through standardization.  If the covariate is binary then you can use $\sigma_{x_j} = \sqrt{f(1-f)}$ where $f$ is the fraction of your sample with the covariate set to 1.


*$n$ is the sample size


*$1/(1-\rho_j^2)$ is the variance inflation factor for the covariate (which can be hard to estimate a priori).
