This vignette illustrates how to use `cond_effect()`

from
the `stdmod`

package. More about this package can be found in
`vignette("stdmod", package = "stdmod")`

or at https://sfcheung.github.io/stdmod/.

`cond_effect()`

Can DoIt can compute the conditional effect of a predictor (the focal variable) on an outcome variable (dependent variable) for selected levels of the moderator:

```
#> Level conscientiousness emotional_stability Effect S.E. t p Sig
#> High 3.950 0.012 0.117 0.107 0.915
#> Medium 3.343 0.214 0.083 2.560 0.011 *
#> Low 2.736 0.415 0.115 3.601 0.000 ***
```

It can also compute standardized conditional moderation effect of this predictor:

```
#> Level conscientiousness emotional_stability Effect S.E. t p Sig
#> High 1.000 0.007 0.063 0.107 0.915
#> Medium 0.000 0.115 0.045 2.560 0.011 *
#> Low -1.000 0.223 0.062 3.601 0.000 ***
```

Nonparametric bootstrap percentile confidence interval can also be
formed for standardized conditional effect using
`cond_effect_boot()`

.

`cond_effect()`

is not designed to be a versatile tool. It
is designed to be a function “good-enough” for common scenarios.
Nevertheless, it can report some useful information along with the
conditional effects, as demonstrated below.

`x`

), and Moderator
(`w`

)`output`

: The output of`lm()`

,`std_selected()`

, or`std_selected_boot()`

, with at least one interaction term. Bootstrap estimates in`std_selected_boot()`

will be ignored because bootstrapping will be done for each level again.`x`

: The predictor (focal variable).`w`

: The moderator.

These are the only required arguments. Just setting them can generate the graph:

```
library(stdmod)
data(sleep_emo_con)
lm_out <- lm(sleep_duration ~ age + gender +
emotional_stability * conscientiousness,
sleep_emo_con)
cond_out <- cond_effect(output = lm_out,
x = "emotional_stability",
w = "conscientiousness")
cond_out
#> The effects of emotional_stability on sleep_duration, conditional on conscientiousness:
#>
#> Level conscientiousness emotional_stability Effect S.E. t p Sig
#> High 3.950 0.012 0.117 0.107 0.915
#> Medium 3.343 0.214 0.083 2.560 0.011 *
#> Low 2.736 0.415 0.115 3.601 0.000 ***
#>
#>
#> The regression model:
#>
#> sleep_duration ~ age + gender + emotional_stability * conscientiousness
#>
#> Interpreting the levels of conscientiousness:
#>
#> Level conscientiousness % Below From Mean (in SD)
#> High 3.950 83.60 1.00
#> Medium 3.343 49.60 0.00
#> Low 2.736 16.60 -1.00
#>
#> - % Below: The percent of cases equal to or less than a level.
#> - From Mean (in SD): Distance of a level from the mean, in standard
#> deviation (+ve above, -ve below).
```

By default, the print method of `cond_effect()`

output
prints the conditional effects, OLS standard errors, *t*
statistics, *p*-values, and significant test results, along with
other information such as the value of each level of the moderator, its
distance from the mean, the percentage of cases equal to or less than
this level. The regression model is also printed. If only the table of
effects is needed, call `print()`

and set
`table_only`

to `TRUE`

:

```
print(cond_out, table_only = TRUE)
#> Level conscientiousness emotional_stability Effect S.E. t p Sig
#> High 3.950 0.012 0.117 0.107 0.915
#> Medium 3.343 0.214 0.083 2.560 0.011 *
#> Low 2.736 0.415 0.115 3.601 0.000 ***
```

More options in printing the output can be found in the help page of
`print.cond_effect()`

.

If the moderator is a numeric variable, then, by default, the conditional effects for three levels of the moderators will be used: one standard deviation (SD) to the mean (“Low”), mean (“Medium”), and one SD above mean (“High”).

Users can also use percentiles to define “Low”, “Medium”, and “High”
by setting `w_method`

to `"percentile"`

. The
default are 16th percentile, 50th percentile, and 84th percentile, which
corresponds approximately to one SD below mean, mean, and one SD above
mean, respectively, for a normal distribution.

```
data(sleep_emo_con)
lm_out <- lm(sleep_duration ~ age + gender +
emotional_stability * conscientiousness,
sleep_emo_con)
cond_out <- cond_effect(output = lm_out,
x = "emotional_stability",
w = "conscientiousness",
w_method = "percentile")
print(cond_out, title = FALSE, model = FALSE)
#> Level conscientiousness emotional_stability Effect S.E. t p Sig
#> High 4.000 -0.004 0.122 -0.034 0.973
#> Medium 3.400 0.195 0.084 2.322 0.021 *
#> Low 2.700 0.427 0.119 3.600 0.000 ***
#>
#> Interpreting the levels of conscientiousness:
#>
#> Level conscientiousness % Below From Mean (in SD)
#> High 4.000 87.20 1.08
#> Medium 3.400 57.00 0.09
#> Low 2.700 16.60 -1.06
#>
#> - % Below: The percent of cases equal to or less than a level.
#> - From Mean (in SD): Distance of a level from the mean, in standard
#> deviation (+ve above, -ve below).
```

Note that the empirical percentage of cases equal to or less than a level may not be exactly equal to that for the requested percentile if the number of cases is small and/or the number of unique values of the moderator is small.

If the moderator is a categorical variable (a string variable or a factor), then the conditional effect of the moderator for each value of this categorical moderator will be printed:

```
set.seed(61452)
sleep_emo_con$city <- sample(c("Alpha", "Beta", "Gamma"),
nrow(sleep_emo_con), replace = TRUE)
lm_cat <- lm(sleep_duration ~ age + gender + emotional_stability*city,
sleep_emo_con)
cond_out <- cond_effect(lm_cat,
x = "emotional_stability",
w = "city")
print(cond_out, title = FALSE, model = FALSE)
#> Level city emotional_stability Effect S.E. t p Sig
#> Alpha Alpha 0.408 0.135 3.027 0.003 **
#> Beta Beta 0.351 0.147 2.388 0.017 *
#> Gamma Gamma 0.020 0.149 0.131 0.896
```

If one or more variables are standardized, the OLS confidence
intervals are not appropriate (Cheung, Cheung, Lau, Hui,
& Vong, 2022; Yuan & Chan,
2011). Users can call `cond_effect_boot()`

to use
nonparametric bootstrapping to form the percentile confidence interval
for each conditional effect.

`conf`

: The level of confidence, expressed as a proportion. Default is .95, requesting a 95% confidence interval.`nboot`

The number of bootstrap samples to drawn. Should be at least 2000 but 5000 is preferable.

To make the results reproducible, call `set.seed()`

before
calling `cond_effect_boot()`

, as illustrated below.

```
lm_out <- lm(sleep_duration ~ age + gender +
emotional_stability * conscientiousness,
sleep_emo_con)
# Standardize all variables and do the moderated regression again
# Use to_standardize as a shortcut to to_center and to_scale
lm_std <- std_selected(lm_out,
to_standardize = ~ .)
set.seed(897043)
cond_std <- cond_effect_boot(output = lm_std,
x = "emotional_stability",
w = "conscientiousness",
nboot = 2000)
print(cond_std, model = FALSE, title = FALSE, level_info = FALSE)
#> Level conscientiousness emotional_stability Effect CI Lower CI Upper S.E.
#> High 1.000 0.007 -0.107 0.117 0.063
#> Medium 0.000 0.115 0.029 0.202 0.045
#> Low -1.000 0.223 0.071 0.364 0.062
#> t p Sig
#> 0.107 0.915
#> 2.560 0.011 *
#> 3.601 0.000 ***
#>
#> [CI Lower, CI Upper] shows the 95% nonparametric bootstrap confidence
#> interval(s) (based on 2000 bootstrap samples).
#>
#> Note:
#>
#> - The variable(s) sleep_duration, emotional_stability,
#> conscientiousness is/are standardized.
#> - The conditional effects are the standardized effects of
#> emotional_stability on sleep_duration.
```

Please refer to the help page of `cond_effect()`

and
`cond_effect_boot()`

for other options available, such as
defining the number of SDs from mean to define “Low” and “High”, the
percentiles to be used, or using parallel processing to speed up
bootstrapping.

Cheung, S. F., Cheung, S.-H., Lau, E. Y. Y., Hui, C. H., & Vong,
W. N. (2022) Improving an old way to measure moderation effect in
standardized units. *Health Psychology*, *41*(7), 502-505.
https://doi.org/10.1037/hea0001188.

Yuan, K.-H., & Chan, W. (2011). Biases and standard errors of
standardized regression coefficients. *Psychometrika, 76*(4),
670-690. https://doi.org/10.1007/s11336-011-9224-6