plot.semipar.mp: Plot massively parallel semiparametric models

In vows: Voxelwise Semiparametrics

Description

Given a massively parallel smoothing object produced by semipar.mp, the function plots the fitted smooth(s) for a given point (e.g., at a given voxel).

Usage

## S3 method for class 'semipar.mp'
plot(x, Y, arr.ind = NULL, which.vox = NULL,
  which.smooth = NULL, coverage = 0.95, length.new = 100, ylim = NULL,
  ylab = NULL, ...)

Arguments

x	an object of class "semipar.mp".
Y	an n \times V outcome matrix.
arr.ind	a 3-element vector specifying the element of the 3-dimensional array of locations (e.g., voxels) for which plotting is desired. If NULL, which.vox must be specified.
which.vox	the index of the voxel to be plotted. If NULL, arr.ind must be specified.
which.smooth	the index of the smooth term of which the confidence interval plot is to be displayed. The default value is NULL which refers to displaying the plots for all the smooth terms in the model.
coverage	the confidence level of the pointwise confidence intervals in the plot.
length.new	length of the vector of ordered variables with which to predict.
ylim, ylab,	arguments to be passed to plot.
...	arguments to be passed to plot.

Author(s)

Yin-Hsiu Chen [e-mail], Philip Reiss [e-mail], and Lan Huo

Examples

n<-32
Ys <- matrix(0, n, 5)
for(i in 1:n) Ys[i,]<--2:2+rnorm(5, i^2, i^0.5)+sin(i)
x1 <- rnorm(n,0,5)
x2 <- 1:n+runif(n, 1, 20)
semipar.obj <- semipar.mp(~x1+sf(x2,k=10),Y=Ys,lsp=seq(5,50,,30))
plot(semipar.obj, Y=Ys, which.vox=2)

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