The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
DynForest is a R package aiming to predict an outcome
using multivariate longitudinal predictors. The method is based on
random forest principle where the longitudinal predictors are modeled
through the random forest. DynForest currently supports
continuous, categorical and survival outcome. The methodology is fully
described for a survival outcome in the paper:
Devaux A., Helmer C., Genuer R. and Proust-Lima C. (2023). Random survival forests with multivariate longitudinal endogenous covariates. Statistical Methods in Medical Research. <doi:10.1177/09622802231206477>
DynForest user guide is also available in the paper:
Devaux A., Proust-Lima C. and Genuer R. (2023). Random Forests for time-fixed and time-dependent predictors: The DynForest R package. arXiv. <doi:10.48550/arXiv.2302.02670>
DynForest package version 1.1.3 could be install from
the CRAN
with:
install.packages("DynForest")Development version of DynForest is also available from
GitHub
with:
# install.packages("devtools")
devtools::install_github("anthonydevaux/DynForest")library(DynForest)
#> Registered S3 method overwritten by 'cmprsk':
#> method from
#> plot.cuminc lcmm
data(pbc2)
# Get Gaussian distribution for longitudinal predictors
pbc2$serBilir <- log(pbc2$serBilir)
pbc2$SGOT <- log(pbc2$SGOT)
pbc2$albumin <- log(pbc2$albumin)
pbc2$alkaline <- log(pbc2$alkaline)# Build longitudinal data
timeData <- pbc2[,c("id","time",
"serBilir","SGOT",
"albumin","alkaline")]
# Create object with longitudinal association for each predictor
timeVarModel <- list(serBilir = list(fixed = serBilir ~ time,
random = ~ time),
SGOT = list(fixed = SGOT ~ time + I(time^2),
random = ~ time + I(time^2)),
albumin = list(fixed = albumin ~ time,
random = ~ time),
alkaline = list(fixed = alkaline ~ time,
random = ~ time))
# Build fixed data
fixedData <- unique(pbc2[,c("id","age","drug","sex")])
# Build outcome data
Y <- list(type = "surv",
Y = unique(pbc2[,c("id","years","event")]))dynforest()
function# Run DynForest function
res_dyn <- dynforest(timeData = timeData, fixedData = fixedData,
timeVar = "time", idVar = "id",
timeVarModel = timeVarModel, Y = Y,
ntree = 50, nodesize = 5, minsplit = 5,
cause = 2, ncores = 15, seed = 1234)summary(res_dyn)
#> dynforest executed for survival (competing risk) outcome
#> Splitting rule: Fine & Gray statistic test
#> Out-of-bag error type: Integrated Brier Score
#> Leaf statistic: Cumulative incidence function
#> ----------------
#> Input
#> Number of subjects: 312
#> Longitudinal: 4 predictor(s)
#> Numeric: 1 predictor(s)
#> Factor: 2 predictor(s)
#> ----------------
#> Tuning parameters
#> mtry: 3
#> nodesize: 5
#> minsplit: 5
#> ntree: 50
#> ----------------
#> ----------------
#> dynforest summary
#> Average depth per tree: 5.94
#> Average number of leaves per tree: 20.44
#> Average number of subjects per leaf: 9.67
#> Average number of events of interest per leaf: 4.34
#> ----------------
#> Computation time
#> Number of cores used: 15
#> Time difference of 37.15092 secs
#> ----------------These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.