profile_plot: Likelihood Profile Plots

In EMC2: Bayesian Hierarchical Analysis of Cognitive Models of Choice

Likelihood Profile Plots

Description

Creates likelihood profile plots from a design and the experimental data by varying one model parameter while holding all others constant.

Usage

profile_plot(
  data,
  design,
  p_vector,
  range = 0.5,
  layout = NA,
  p_min = NULL,
  p_max = NULL,
  use_par = NULL,
  n_point = 100,
  n_cores = 1,
  round = 3,
  true_args = list(),
  ...
)

Arguments

data	A dataframe. Experimental data used, needed for the design mapping
design	A design list. Created using design.
p_vector	Named vector of parameter values (typically created with sampled_pars(design))
range	Numeric. The max and min will be p_vector + range/2 and p_vector - range/2, unless specified in p_min or p_max.
layout	A vector indicating which layout to use as in par(mfrow = layout). If NA, will automatically generate an appropriate layout.
p_min	Named vector. If specified will instead use these values for minimum range of the selected parameters.
p_max	Named vector. If specified will instead use these values for maximum range of the selected parameters.
use_par	Character vector. If specified will only plot the profiles for the specified parameters.
n_point	Integer. Number of evenly spaced points at which to calculate likelihood
n_cores	Number of likelihood points evenly spaced between the minimum and maximum likelihood range.
round	Integer. To how many digits will the output be rounded.
true_args	A list. Optional additional arguments that can be passed to plot.default for the plotting of the true vertical line.
...	Optional additional arguments that can be passed to plot.default.

Value

Vector with highest likelihood point, input and mismatch between true and highest point

Examples

# First create a design
design_DDMaE <- design(data = forstmann,model=DDM,
                      formula =list(v~0+S,a~E, t0~1, s~1, Z~1, sv~1, SZ~1),
                      constants=c(s=log(1)))
# Then create a p_vector:
p_vector=c(v_Sleft=-2,v_Sright=2,a=log(.95),a_Eneutral=log(1.5),a_Eaccuracy=log(2),
          t0=log(.25),Z=qnorm(.5),sv=log(.5),SZ=qnorm(.5))
# Make a profile plot for some parameters. Specifying a custom range for t0.
profile_plot(p_vector = p_vector, p_min = c(t0 = -1.35),
             p_max = c(t0 = -1.45), use_par = c("a", "t0", "SZ"),
             data = forstmann, design = design_DDMaE, n_point = 10)

EMC2 documentation built on April 11, 2025, 5:50 p.m.