Plot optimizer trace from a ferx fit

Description

Produces a two-panel diagnostic plot from the per-iteration trace written when optimizer_trace = TRUE was passed to [ferx_fit](ferx_fit.qmd). The top panel shows OFV over iterations; the bottom panel shows a method-specific convergence metric (gradient norm for gradient-based methods, MH accept rate for SAEM, LM lambda for Gauss-Newton). Phase boundaries are drawn as vertical dashed lines.

Usage

ferx_plot_trace(fit, log_ofv = FALSE)

Arguments

  • fit: A ferx_fit object or path to a trace CSV file (see [ferx_trace](ferx_trace.qmd)).
  • log_ofv: Logical; plot OFV on a log scale relative to the final value (\(OFV - OFV_{final}\)). Default FALSE.

Seealso

Other diagnostics: [check_diagnostics](check_diagnostics.qmd)(), [ferx_cor_matrix](ferx_cor_matrix.qmd)(), [ferx_cov_screen](ferx_cov_screen.qmd)(), [ferx_estimates](ferx_estimates.qmd)(), [ferx_eta_cov](ferx_eta_cov.qmd)(), [ferx_runlog_iters](ferx_runlog_iters.qmd)(), [ferx_trace](ferx_trace.qmd)(), [ferx_warnings](ferx_warnings.qmd)(), [summary.ferx_fit](summary.ferx_fit.qmd)()

Concept

diagnostics

Value

Invisibly returns the trace data frame (from [ferx_trace](ferx_trace.qmd)).

Examples

ex  <- ferx_example("warfarin")
fit <- ferx_fit(ex$model, ex$data, method = "gn", covariance = FALSE,
                optimizer_trace = TRUE)
tr  <- ferx_trace(fit)
head(tr)