brmViolin {pcvr} | R Documentation |
Function to visualize hypotheses tested on brms models similar to those made using growthSS outputs.
Description
Function to visualize hypotheses tested on brms models similar to those made using growthSS outputs.
Usage
brmViolin(
fit,
params = NULL,
hyp = "num/denom>1.05",
compareX = "a",
againstY = "b",
group_sep = "[.]",
groups_into = c("group2"),
x = NULL,
facet = NULL,
cores = getOption("mc.cores", 1),
returnData = FALSE
)
Arguments
fit |
A brmsfit object or a list of brmsfit objects |
params |
A list of parameters to use from the fit. Defaults to NULL in which case all growth model parameters are used. |
hyp |
A character string defining the hypothesis to be tested. Defaults to "num/denom > 1.05". The "num" and "denom" names should be kept, but the direction and magnitude can be changed. |
compareX |
Which groups in the model should be compared as numerator options? Defaults to NULL in which case a plot will not be made but the data will be returned. |
againstY |
Which group in the model should be used as the denominator (typically a control group to compare against)? Defaults to NULL in which case a plot will not be made but the data will be returned. |
group_sep |
A regex pattern to match the separation of grouping
terms in the models group term (see the formula argument of |
groups_into |
A vector of column names to make after groups are split by group_sep. If this or groups_sep are NULL then no groups are assumed. |
x |
The variable to be plotted on the x axis (should be from groups_into). |
facet |
The variable to be used to facet the ggplot (should be another option from groups_into). If left NULL then the plot will only be faceted by params. Note that with the nature of againstY this faceting is often redundant but it does add labels which are helpful for keeping results organized.. |
cores |
Optional number of cores to run hypotheses in parallel. Defaults to 1 unless the "mc.cores" option is set. |
returnData |
Logical, should data be returned? This is treated as TRUE if a plot will not be generated but otherwise defaults to FALSE. |
Value
Returns a ggplot showing a brms model's posterior distributions as violins and filled by posterior probability of some hypothesis.
Examples
set.seed(123)
simdf <- growthSim(
"logistic", n = 20, t = 25,
params = list("A" = c(200, 160), "B" = c(13, 11), "C" = c(3, 3.5))
)
ss <- growthSS(
model = "logistic", form = y ~ time | id / group, sigma = "spline",
list("A" = 130, "B" = 10, "C" = 3),
df = simdf, type = "brms"
)
fit <- fitGrowth(ss, backend = "cmdstanr", iter = 500, chains = 1, cores = 1)
brmViolin(fit, hyp = "num/denom>1.05",
compareX = "a",
againstY = "b", returnData = TRUE)