Create a random effects dataset object
This function is intended for advanced use cases in which users require detailed control of sampling algorithms and data structures.
Minimal input validation and error checks are performed – users are responsible for providing the correct inputs.
For tutorials on the "proper" usage of the stochtree's advanced workflow, we provide several vignettes at stochtree.ai
Usage
createRandomEffectsDataset(group_labels, basis, variance_weights = NULL)
Arguments
- group_labels
Vector of group labels
- basis
Matrix of bases used to define the random effects regression (for an intercept-only model, pass an array of ones)
- variance_weights
(Optional) Vector of observation-specific variance weights
Value
RandomEffectsDataset object
Examples
rfx_group_ids <- sample(1:2, size = 100, replace = TRUE)
rfx_basis <- matrix(rnorm(3*100), ncol = 3)
weight_vector <- rnorm(100)
rfx_dataset <- createRandomEffectsDataset(rfx_group_ids, rfx_basis)
rfx_dataset <- createRandomEffectsDataset(rfx_group_ids, rfx_basis, weight_vector)