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.2.0 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") Quick example with survival outcome
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")]))
# 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.05251 secs #> ----------------
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4