Description Usage Arguments Details Value Author(s) References See Also Examples
resample.sprinter
is a wrapper, which subsamples the data and performs the function sprinter on each subsample. It evaluates the interaction terms by an improvement in prediction performance. The importance of an interaction is judged by variable inclusion frequencies and mean coefficients for each interaction term.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | resample.sprinter(x, time,
status,
fold = 5,
oob.rel = 0.632,
mandatory,
repetitions = 25,
n.inter.candidates = 1000,
screen.main,
screen.inter = fit.rf,
fit.final,
args.screen.main = list(),
args.screen.inter = list(),
args.fit.final = args.screen.main,
orthogonalize = TRUE,
parallel = FALSE,
mc.cores = detectCores(), ...)
|
x |
n * p matrix of covariates. |
time |
vector of length |
status |
censoring indicator, i.e., vector of length |
fold |
number of subsamples. Default is 5. |
oob.rel |
frequency of subsamples. |
mandatory |
vector with variable names of mandatory covariates, where parameter estimation should be performed unpenalized. |
repetitions |
number of repetitions of the interaction screening approach. Repetitions are performed by creating subsamples and applying the interaction screening approach on each subsample dataset separately. |
n.inter.candidates |
minimal number of potential interaction candidates, which are considered in the final model building step. |
screen.main |
function for screening potential main effects. See more in details. |
screen.inter |
function for detecting potential interaction candidates. See more in details. |
fit.final |
function for building final Cox proportional hazards model. Default is the function set in screen.main. See more in details. |
args.screen.main |
list of arguments which should be used in the main effects detection step. |
args.screen.inter |
list of arguments which should be used in the interaction screening step. |
args.fit.final |
list of arguments which should be used in the final model building step. |
orthogonalize |
logical value. If true all variables are made orthogonal to those that are assessed as main effects by |
parallel |
logical value indicating whether the interaction screening step should be performed parallel. |
mc.cores |
the numbers of cores to use. |
... |
additional arguments. |
resample.sprinter
provides the possibility of evaluating the relevance of an interaction by using resampling techniques and the resultant variable inclusion frequency.
This function subsamples the data and applies the function sprinter
on each subsample.
As a result the user can print the interactions with the largest inclusion frequencies and their mean coefficients.
Returns a list of length fold
of objects produced by sprinter
.
Written by Isabell Hoffmann isabell.hoffmann@uni-mainz.de.
Sariyar, M., Hoffmann, I. Binder, H. (2014). Combining techniques for screening and evaluating interaction terms on high-dimensional time-to-event data. BMC Bioinformatics 15:58.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 | ############
# Example for resampling and screening main effects by CoxBoost,
# screening interactions by random forest and
# building a final Cox proportional hazards model by CoxBoost:
# system.time with one core:
# user system elapsed
# 399.78 2.40 403.02
# For a faster run set repetitions down!
############
## Not run:
# Create survival data with interactions:
data <- simul.int(518,n = 200, p = 500,
beta.int = 1.0,
beta.main = 0.9,
censparam = 1/20,
lambda = 1/20)[[1]]
set.seed(123)
resamcb <- resample.sprinter( x=data[,1:500],
time = data$obs.time,
status= data$obs.status,
fold = 10,
repetitions = 5,
mandatory = c("ID1","ID2"),
n.inter.candidates = 1000,
screen.main = fit.CoxBoost,
fit.final = fit.CoxBoost,
args.screen.main = list(stepno = 5),
parallel = F)
summary(resamcb)
summary(resamcb, optional = T)
## End(Not run)
# true coefficients:
# Clin.cov1 Clin.cov2 ID5:ID6 ID7:ID8
# 0.9 -0.9 1 -1
## Example for resampling and screening main effects by CoxBoost,
## screening interactions by random forest and
## building a final Cox proportional hazards model by CoxBoost:
##
##system.time with one core:
## user system elapsed
## 479.55 2.39 482.40
## For a faster run set repetitions down!
## Not run:
# Create survival data with interactions:
data <- simul.int(518,n = 200, p = 500,
beta.int = 1.0,
beta.main = 0.9,
censparam = 1/20,
lambda = 1/20)[[1]]
resamunicox <- resample.sprinter( x=data[,1:500],
time = data$obs.time,
status= data$obs.status,
fold = 10,
repetitions = 10,
mandatory = c("ID1","ID2"),
n.inter.candidates = 1000,
screen.main = fit.uniCox,
fit.final = fit.uniCox,
parallel = FALSE)
summary(resamunicox)
summary(resamunicox, optional = TRUE)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.