plot.SIR_threshold: Graphical output of SIR_threshold

In SIRthresholded: Sliced Inverse Regression with Thresholding

Graphical output of SIR_threshold

Description

Display the 10 first eigen values and the estimated index versus Y of the thresholded SIR model.

Usage

## S3 method for class 'SIR_threshold'
plot(x, choice = "", ...)

Arguments

x	A SIR_threshold object
choice	the graph to plot: "eigvals" Plot the eigen values of the matrix of interest. "estim_ind" Plot the estimated index by the SIR model versus Y. "" Plot every graphs (default).
...	arguments to be passed to methods, such as graphical parameters (not used here).

Value

No return value

Examples

# Generate Data
set.seed(10)
n <- 500
beta <- c(1,1,rep(0,8))
X <- mvtnorm::rmvnorm(n,sigma=diag(1,10))
eps <- rnorm(n)
Y <- (X%*%beta)**3+eps

# Apply SIR with hard thresholding
res = SIR_threshold(Y, X, H = 10, lambda = 0.2, thresholding = "hard")

# Eigen values
plot(res,choice="eigvals")

# Estimated index versus Y
plot(res,choice="estim_ind")

SIRthresholded documentation built on July 10, 2023, 2:03 a.m.