vows (version 0.5)

screen.vox: Screen voxels for a voxelwise smoothing object

Description

Inputs a voxelwise smoothing object as produced by semipar4d, and outputs an object containing the results for a subset of the voxels.

Usage

screen.vox(semi.obj, arr4d, include)

Arguments

semi.obj

an object of class semipar.mp.

arr4d

the 4-dimensional array used to generate the object.

include

a logical matrix indicating which points (or voxels) should be included.

Value

a modified version of semipar.obj, with pointwise coefficients (coef component), pointwise degrees of freedom (pwdf), pointwise log smoothing parameter (pwlsp), and pointwise variance estimate (sigma2) for the points specified by include only.

See Also

semipar.mp

Examples

data(test)
d4 = test$d4
x = test$x
vw.obj = semipar4d(d4, formula = ~sf(x), data = data.frame(x = x), lsp=-5:5)

# Include only the first 600 voxels
sv = screen.vox(vw.obj, d4, rep(1:0, c(600,400)))