¥ Available in PDF (scanned from hand-written pages) for the
truly devoted
H File 3dREMLfit_mathnotes.pdf
¥ Outline of REML estimation methodology
H What is
REML and why do we care?
¥ Matrix algebra for efficient solution of the many linear
systems that must be solved for each voxel
H Sparse
matrix factorizations, multiplications, and solvers
¥ How ARMA(1,1) parameters are estimated in 3dREMLfit
H Optimizing
REML log-likelihood function over a discrete grid of (a,b) values, using
2D binary search
H Must solve
a GLSQ problem for each (a,b) tested, for
each voxel
¥ How statistics are implemented as GLTs
H Testing
null hypothesis Gb=0 for
arbitrary matrix G
¥ Derivation of ARMA(1,1) formulas
H For
completeness, and because we all love equations