Runs the variational Laplace (Fisher scoring) loop on the structure built by
dcm_peb_prepare, and returns the group-level PEB together with
the first-level DCMs updated under the empirical priors. Mirrors the
estimation half of SPM25's spm_dcm_peb.
dcm_peb_run(prep, verbose = TRUE)The list returned by dcm_peb_prepare.
Logical. Report free energy per iteration via
message(). Silence with suppressMessages() or
verbose = FALSE.
List with two elements: PEB, the group-level result
(Ep group parameter expectations, Cp their covariance,
Eh/Ch log-precision estimates, F free energy, and the
Pnames/Xnames/Snames labels), and P, the input
DCMs with Ep, Cp, M$pE, M$pC and F
updated under the empirical priors.
Progress is reported as VL Iter k: F=... dF=... [t=...], where
t is the log step size of the Fisher scoring update. The loop stops
when the step size collapses or the free-energy increase falls below 1e-4.