-
Notifications
You must be signed in to change notification settings - Fork 0
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
1 parent
7db9225
commit 8df0ca7
Showing
1 changed file
with
25 additions
and
0 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,25 @@ | ||
Fit a spectral neural mass model, specified using the DCM conventional structure (DCM with DCM.M.pE, DCM.M.pC,DCM.M.f, DCM.M.IS and DCM.xY specified). | ||
|
||
Parameter estimation using log-likelihood (and/or Free Energy) optimisation with Levenberg-Marquardt. | ||
|
||
The functions estimate model parameters by maximising the log-likelihood (or free energy) through iterative optimisation. The algorithms employ the Levenberg-Marquardt | ||
method with dynamic updates of hyperparameters (observation variance and damping factor). | ||
|
||
Usage: | ||
|
||
Given a fully specified DCM, do: | ||
|
||
''' | ||
M = aFitDCM(DCM); | ||
M.aloglik(num_iter) ; <-- log likelihood estimation or | ||
M.aloglikFE(num_iter) ; <-- free enery estimation | ||
|
||
% to re-run / add more iterations: | ||
M.update_parameters(M.Ep) | ||
M.aloglikFE(num_iter) | ||
|
||
% and to access posteriors: | ||
M.Ep | ||
M.CP | ||
M.F | ||
''' |