MBNMAtime (version 0.2.4)

fitplot: Plot fitted values from MBNMA model

Description

Plot fitted values from MBNMA model

Usage

fitplot(
  mbnma,
  treat.labs = NULL,
  disp.obs = TRUE,
  n.iter = round(mbnma$BUGSoutput$n.iter/4),
  n.thin = mbnma$BUGSoutput$n.thin,
  ...
)

Value

Generates a plot of fitted values from the MBNMA model and returns a list containing the plot (as an object of class c("gg", "ggplot")), and a data.frame of posterior mean fitted values for each observation.

Arguments

mbnma

An S3 object of class "mbnma" generated by running a time-course MBNMA model

treat.labs

A character vector of treatment labels with which to name graph panels. Can use mb.network()[["treatments"]] with original dataset if in doubt.

disp.obs

A boolean object to indicate whether raw data responses should be plotted as points on the graph

n.iter

number of total iterations per chain (including burn in; default: 2000)

n.thin

thinning rate. Must be a positive integer. Set n.thin > 1 to save memory and computation time if n.iter is large. Default is max(1, floor(n.chains * (n.iter-n.burnin) / 1000)) which will only thin if there are at least 2000 simulations.

...

Arguments to be sent to ggplot2::ggplot()

Details

Fitted values should only be plotted for models that have converged successfully. If fitted values (theta) have not been monitored in mbnma$parameters.to.save then additional iterations will have to be run to get results for these.

Examples

# \donttest{
# Make network
painnet <- mb.network(osteopain)

# Run MBNMA
mbnma <- mb.run(painnet,
  fun=temax(pool.emax="rel", method.emax="common",
    pool.et50="abs", method.et50="random"))

# Plot fitted values from the model
# Monitor fitted values for 500 additional iterations
fitplot(mbnma, n.iter=500)
# }