nice_mod {rempsyc} | R Documentation |
Easy moderations
Description
Easily compute moderation analyses, with effect sizes, and format in publication-ready format.
Usage
nice_mod(
data,
response,
predictor,
moderator,
moderator2 = NULL,
covariates = NULL,
b.label = "b",
standardize = TRUE,
mod.id = TRUE,
ci.alternative = "two.sided",
...
)
Arguments
data |
The data frame |
response |
The dependent variable. |
predictor |
The independent variable. |
moderator |
The moderating variable. |
moderator2 |
The second moderating variable, if applicable. |
covariates |
The desired covariates in the model. |
b.label |
What to rename the default "b" column (e.g.,
to capital B if using standardized data for it to be converted
to the Greek beta symbol in the |
standardize |
Logical, whether to standardize the
data before fitting the model. If |
mod.id |
Logical. Whether to display the model number, when there is more than one model. |
ci.alternative |
Alternative for the confidence interval of the sr2. It can be either "two.sided (the default in this package), "greater", or "less". |
... |
Further arguments to be passed to the |
Details
The effect size, sr2 (semi-partial correlation squared, also
known as delta R2), is computed through effectsize::r2_semipartial.
Please read the documentation for that function, especially regarding
the interpretation of the confidence interval. In rempsyc
, instead
of using the default one-sided alternative ("greater"), we use the
two-sided alternative.
To interpret the sr2, use effectsize::interpret_r2_semipartial()
.
For the easystats equivalent, use report::report()
on the lm()
model object.
Value
A formatted dataframe of the specified lm model, with DV, IV, degrees of freedom, regression coefficient, t-value, p-value, and the effect size, the semi-partial correlation squared, and its confidence interval.
See Also
Checking simple slopes after testing for moderation:
nice_slopes
, nice_lm
,
nice_lm_slopes
. Tutorial:
https://rempsyc.remi-theriault.com/articles/moderation
Examples
# Make the basic table
nice_mod(
data = mtcars,
response = "mpg",
predictor = "gear",
moderator = "wt"
)
# Multiple dependent variables at once
nice_mod(
data = mtcars,
response = c("mpg", "disp", "hp"),
predictor = "gear",
moderator = "wt"
)
# Add covariates
nice_mod(
data = mtcars,
response = "mpg",
predictor = "gear",
moderator = "wt",
covariates = c("am", "vs")
)
# Three-way interaction
x <- nice_mod(
data = mtcars,
response = "mpg",
predictor = "gear",
moderator = "wt",
moderator2 = "am"
)
x
# Get interpretations
cbind(x, Interpretation = effectsize::interpret_r2_semipartial(x$sr2))