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)
``` |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.