In methexp/HMMTreeR: An R Interface to HMMTree

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
library("HMMTreeR")

# Example 1: HMMTree example files for PD model
res <- lc(model_file = "pd.eqn", data_file = "pd.dat")

# extract fit statistics
fit_statistics(res)

# extract parameter estimates
parameter_estimates(res)

# compute weighted means for parameter estimates
weighted_means(res)

# compute weighted means for parameter estimates from winning model
weighted_means(res[[length(res)]])