spm_LAP_iS.m 780 B

1234567891011121314151617181920
  1. function [iS] = spm_LAP_iS(p,R)
  2. % default precision function for LAP models (hidden states)
  3. % FORMAT [iS] = spm_LAP_iS(p,R)
  4. %
  5. % p{1} - vector of precisions for causal states (v)
  6. % p{2} - vector of precisions for hidden states (v)
  7. % R - generalised precision matrix
  8. %
  9. % iS - precision matrix for generalised states (causal and then hidden)
  10. %__________________________________________________________________________
  11. % Copyright (C) 2008 Wellcome Trust Centre for Neuroimaging
  12. % Karl Friston
  13. % $Id: spm_LAP_iS.m 3694 2010-01-22 14:16:51Z karl $
  14. % precisions for generalised errors on causal (V) and hidden states (W)
  15. %==========================================================================
  16. V = kron(R,diag(exp(p.h)));
  17. W = kron(R,diag(exp(p.g)));
  18. iS = blkdiag(V,W);