vows (version 0.5)

semipar4d: Massively parallel semiparametric regression for 4-dimensional data

Description

This is a wrapper function for semipar.mp to handle 3D image responses.

Usage

semipar4d(arr4d, formula, lsp, data, range.basis = NULL, knots = "quantile", rm.constr = FALSE, random = NULL, store.reml = FALSE, store.fitted = FALSE)

Arguments

arr4d

a 4-dimensional response array, where the first 3 dimensions refer to spatial coordinates and the last dimension corresponds to different images.

formula, lsp, data, range.basis, knots, rm.constr, random, store.reml, store.fitted

see semipar.mp.

Value

An object of class "semipar.mp", with two changes. (1) If store.fitted = TRUE, the fitted values are given as a 4-dimensional array. (2) A call component is included.

See Also

semipar.mp

Examples

data(test)
d4 = test$d4
x = test$x
semi.obj = semipar4d(d4, ~sf(x), lsp=-5:5, data=data.frame(x = x))
plot(semi.obj, which.vox = 4)