vows (version 0.4)

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)