msir (version 1.1)

plot.msir: Plot method for 'msir' objects.

Description

Plots directions and other information from MSIR estimation.

Usage

## S3 method for class 'msir':
plot(x, which, 
     type = c("pairs", "2Dplot", "spinplot", "evalues", "coefficients"), 
     span = NULL, std = TRUE, ylab, xlab, restore.par = TRUE, ...)

Arguments

x

a 'msir' object.

which

a vector of value(s) giving the directions for which the plot should be drawn.

type

the type of plot to be drawn.

span

the span of smoother (only for type = "pairs" | "2Dplot").

std

if TRUE coefficients are standardized (only for type = "coefficients").

ylab

a character string for the y-axis label.

xlab

a character string for the x-axis label.

restore.par

if TRUE the graphical parameters (see par) changed are restored to the previous state. If you want to manipulate the resulting plot you should set restore.par = FALSE.

...

additional arguments.

References

Scrucca, L. (2011) Model-based SIR for dimension reduction. Computational Statistics and Data Analysis, 55(11), 3010-3026.

See Also

msir

Examples

# 2-dimensional response curve
n = 300; p = 5
b1 = c(1, 1, 1, rep(0, p-3))
b2 = c(1,-1,-1, rep(0, p-3))
b = cbind(b1,b2)
x = mvrnorm(n, mu = rep(0,p), Sigma = diag(p))
y = x %*% b1 + (x %*% b1)^3 + 4*(x %*% b2)^2 + rnorm(n)
MSIR = msir(x, y)
summary(MSIR)
plot(MSIR)
plot(MSIR, which = 1:2)
plot(MSIR, type = "2Dplot", which = 1, span = 0.7)
plot(MSIR, type = "2Dplot", which = 2, span = 0.7)
plot(MSIR, type = "spinplot")
plot(MSIR, type = "evalues")
plot(MSIR, type = "coefficients")