View source: R/long-grmtree-latentpar.R
| latentpar_longitudinal_grmtree | R Documentation |
Extracts the latent trait distribution parameters from each terminal node of
a longitudinal GRM tree: the T2 latent mean shift (\mu_{T2}), the T2
latent variance (\sigma^2_{T2}), and the T1-T2 correlation.
latentpar_longitudinal_grmtree(object, node = NULL, ...)
object |
A |
node |
Optional vector of node IDs. If NULL, all terminal nodes. |
... |
Additional arguments (currently unused). |
A data.frame with columns:
Terminal node ID
Sample size in the node
Latent mean at T2 (positive = improvement)
Latent variance at T2
Correlation between T1 and T2 latent traits
longitudinal_grmtree for Phase 1 (tree fitting),
fscores_longitudinal_grmtree, for computing factor scores for
longitudinal GRMTree
library(grmtree)
# Load the synthetic longitudinal data
data("grmtree_long_data", package = "grmtree")
# Prepare the wide-format response matrix
items_t1 <- c("MOS_Listen", "MOS_Info", "MOS_Advice_Crisis", "MOS_Confide",
"MOS_Advice_Want", "MOS_Fears", "MOS_Personal", "MOS_Understand")
ld <- prepare_longitudinal_data(
data = grmtree_long_data,
items_t1 = items_t1,
items_t2 = paste0(items_t1, "_year1"),
covariates = c("sex", "age", "residency", "job",
"education", "comorbidity_count", "ever_smoker")
)
# Phase 1: fit the longitudinal GRM tree
ltree <- longitudinal_grmtree(
resp_wide ~ sex + age + residency + job +
education + comorbidity_count + ever_smoker,
data = ld, n_items = 8,
control = grmtree.control(minbucket = 200)
)
# Print latent trait summary
latent <- latentpar_longitudinal_grmtree(ltree)
print(latent)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.