DynForest {DynForest} | R Documentation |

Build a random forest using multivariate longitudinal endogenous covariates

```
DynForest(
timeData = NULL,
fixedData = NULL,
idVar = NULL,
timeVar = NULL,
timeVarModel = NULL,
Y = NULL,
ntree = 200,
mtry = NULL,
nodesize = 1,
minsplit = 2,
cause = 1,
nsplit_option = "quantile",
ncores = NULL,
seed = round(runif(1, 0, 10000)),
verbose = TRUE
)
```

`timeData` |
A data.frame containing the id and time measurements variables and the time-dependent predictors. |

`fixedData` |
A data.frame containing the id variable and the time-fixed predictors. Categorical variables should be characterized as factor. |

`idVar` |
A character indicating the name of variable to identify the subjects |

`timeVar` |
A character indicating the name of time variable |

`timeVarModel` |
A list for each time-dependent predictors containing a list of formula for fixed and random part from the mixed model |

`Y` |
A list of output which should contain: |

`ntree` |
Number of trees to grow. Default value set to 200. |

`mtry` |
Number of candidate variables randomly drawn at each node of the trees. This parameter should be tuned by minimizing the OOB error. Default is defined as the square root of the number of predictors. |

`nodesize` |
Minimal number of subjects required in both child nodes to split. Cannot be smaller than 1. |

`minsplit` |
(Only with survival outcome) Minimal number of events required to split the node. Cannot be smaller than 2. |

`cause` |
(Only with competing events) Number indicates the event of interest. |

`nsplit_option` |
A character indicates how the values are chosen to build the two groups for the splitting rule (only for continuous predictors). Values are chosen using deciles ( |

`ncores` |
Number of cores used to grow trees in parallel. Default value is the number of cores of the computer-1. |

`seed` |
Seed to replicate results |

`verbose` |
A logical controlling the function progress. Default is |

The function currently supports survival (competing or single event), continuous or categorical outcome.

FUTUR IMPLEMENTATIONS:

Continuous longitudinal outcome

DynForest function return a list with the following elements:

`data` | A list containing the data used to grow the trees |

`rf` | A table with each tree in column. Provide multiple characteristics about the tree building |

`type` | Outcome type |

`times` | A numeric vector containing the time-to-event for all subjects |

`cause` | Indicating the cause of interest |

`causes` | A numeric vector containing the causes indicator |

`Inputs` | A list of 3 elements: `Longitudinal` , `Numeric` and `Factor` . Each element contains the names of the predictors |

`Longitudinal.model` | A list of longitudinal markers containing the formula used for modeling in the random forest |

`param` | A list containing the hyperparameters |

`comput.time` | Computation time |

Anthony Devaux (anthony.devaux@u-bordeaux.fr)

Devaux A., Helmer C., Dufouil C., Genuer R., Proust-Lima C. (2022). Random survival forests for competing risks with multivariate longitudinal endogenous covariates. arXiv <doi: 10.48550/arXiv.2208.05801>

`summary.DynForest compute_OOBerror compute_VIMP compute_gVIMP predict.DynForest plot.DynForest`

```
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)
# Sample 100 subjects
set.seed(1234)
id <- unique(pbc2$id)
id_sample <- sample(id, 100)
id_row <- which(pbc2$id%in%id_sample)
pbc2_train <- pbc2[id_row,]
timeData_train <- pbc2_train[,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_train <- unique(pbc2_train[,c("id","age","drug","sex")])
# Build outcome data
Y <- list(type = "surv",
Y = unique(pbc2_train[,c("id","years","event")]))
# Run DynForest function
res_dyn <- DynForest(timeData = timeData_train, fixedData = fixedData_train,
timeVar = "time", idVar = "id",
timeVarModel = timeVarModel, Y = Y,
ntree = 50, nodesize = 5, minsplit = 5,
cause = 2, ncores = 2, seed = 1234)
```

[Package *DynForest* version 1.1.0 Index]