srcpcc {sensitivity} | R Documentation |
srcpcc
computes the standardized regression coefficients (SRC)
and the partial correlation coefficients (PCC). Analysis can be done
on the ranks; then the indices are the standardized rank regression
coefficients (SRRC) and the partial rank correlation coefficients
(PRCC).
srcpcc(model = NULL, x, pcc = TRUE, rank = FALSE, nboot = 0, conf = 0.95, ...)
model |
the model |
x |
the input sample |
pcc |
logical. If TRUE , the P(R)CCs are computed |
rank |
logical. If TRUE , the analysis is done on the
ranks |
nboot |
the number of bootstrap replicates |
conf |
the confidence level for bootstrap confidence intervals |
... |
any other arguments for model which are passed
unchanged each time it is called |
model
is a function or a predictor (a class with a
predict
method) computing the response y
based on the
sample given by x
. If no model is specified, the indices will be
computed when one gives the response.
srcpcc
returns an object of class "srcpcc"
.
An object of class "srcpcc"
is a list containing the following
components:
y |
the response |
src |
the estimations of the SRC indices (or SRRC if rank analysis is requested) |
pcc |
if requested, the estimations of the PCC indices (or PRCC if rank analysis is requested) |
The number of model evaluations is n where n is the
size of the sample x
.
Saltelli, A., Chan, K. and Scott, E. M., 2000, Sensitivity analysis, Wiley.
# linear model : Y = X1 + X2 + X3 model1 <- function(x) x[, 1] + x[, 2] + x[, 3] # a 100-sample with X1 ~ U(0.5, 1.5) # X2 ~ U(1.5, 4.5) # X3 ~ U(4.5, 13.5) n <- 100 x <- data.frame(X1 = runif(n, 0.5, 1.5), X2 = runif(n, 1.5, 4.5), X3 = runif(n, 4.5, 13.5)) # sensitivity analysis sa <- srcpcc(model = model1, x = x, nboot = 100) print(sa) par(mfrow = c(1,2)) plot(sa, ask = FALSE)