screen.vox: Screen voxels for a voxelwise smoothing object

In vows: Voxelwise Semiparametrics

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.

Author(s)

Lei Huang [e-mail] and Philip Reiss [e-mail]

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)))

vows documentation built on May 2, 2019, 9:26 a.m.