modalities.fmri.spm.reml¶
Module: modalities.fmri.spm.reml
¶
Functions¶
- nipy.modalities.fmri.spm.reml.orth(X, tol=1e-07)¶
Compute orthonormal basis for the column span of X.
Rank is determined by zeroing all singular values, u, less than or equal to tol*u.max().
- INPUTS:
X – n-by-p matrix
- OUTPUTS:
- B – n-by-rank(X) matrix with orthonormal columns spanning
the column rank of X
- nipy.modalities.fmri.spm.reml.reml(sigma, components, design=None, n=1, niter=128, penalty_cov=1.2664165549094176e-14, penalty_mean=0)¶
Adapted from spm_reml.m
ReML estimation of covariance components from sigma using design matrix.
- INPUTS:
sigma – m-by-m covariance matrix components – q-by-m-by-m array of variance components
mean of sigma is modeled as a some over components[i]
- design – m-by-p design matrix whose effect is to be removed for
ReML. If None, no effect removed (???)
n – degrees of freedom of sigma penalty_cov – quadratic penalty to be applied in Fisher algorithm.
If the value is a float, f, the penalty is f * identity(m). If the value is a 1d array, this is the diagonal of the penalty.
- penalty_mean – mean of quadratic penalty to be applied in Fisher
algorithm. If the value is a float, f, the location is f * np.ones(m).
- OUTPUTS:
C – estimated mean of sigma h – array of length q representing coefficients
of variance components
cov_h – estimated covariance matrix of h