Compute PLS-PCE Sensitivity Indexes

Share:

Description

Compute the optimal number of components, the PLS-PCE sensitivity indexes and related results.

Usage

1
calcPLSPCE(pce, nc = 2)

Arguments

pce

an object of class PCEpoly. Design to analyze.

nc

integer. Required number of components.

Value

An object of class PLSPCE.

Examples

1
2
3
4
5
6
7
8
9
### Load the dataset
load(system.file("extdata",  "ishigami200.Rda", package="plspolychaos"))
X <- ishi200[, -ncol(ishi200)] # inputs
Y <- ishi200[,  ncol(ishi200)] # output
### Creation of the full polynomials
degree <- 6 # polynomial degree
pce <- polyLeg(X, Y, degree) 
### Compute the PLS-PCE sensitivity indexes for 25 components
ret <-  calcPLSPCE(pce, nc=25)