This function computes a statistical map and residuals which are the objects necessary to perform the parametric bootstrap joint (PBJ) inference procedure.

bootStats(images, coefficients = NULL, X, Xred,
  statistic = function(stat, ...) stat, ...)

Arguments

images

An n by V matrix of voxels inside the brain.

coefficients

An n by V matrix of parameter estimates to subtract from the bootstrapped values. If computing p-values this should be the observed parameter estimates. If computing confidence intervals then this should be NULL (the default).

X

Design matrix for the full model.

Xred

Design matrix for the reduced model.

statistic

A function of the form statistic(stat, ...) that returns a statistic of interest.

...

Arguments passed to statistic function.

Value

Returns a vector of the resulting estimate