## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----eval=FALSE---------------------------------------------------------------
# # install.packages("remotes")
# remotes::install_github("traitecoevo/hmde")

## -----------------------------------------------------------------------------
library(dplyr)
library(ggplot2)
library(hmde)

## -----------------------------------------------------------------------------
prior_pars(hmde_model("canham_multi_ind"))

## ----eval=FALSE---------------------------------------------------------------
# # Build fit and extract estimates
# set.seed(2026)
# canham_multi_ind_fit <- hmde_data_template("canham_multi_ind",
#                                            obs_data = Tree_Size_Data) |>
#   hmde_run(chains = 1, cores = 1, iter = 1000)
# 
# Tree_Size_Ests <- hmde_estimates(fit = canham_multi_ind_fit,
#                                  obs_data = Tree_Size_Data)

## -----------------------------------------------------------------------------
summary(Tree_Size_Ests)

# Plot fitted growth function pieces
hmde_plot_de_pieces(Tree_Size_Ests)

## -----------------------------------------------------------------------------
#Plots of size over time for a sample of 5 individuals
hmde_plot_obs_est_inds(Tree_Size_Ests, n_ind_to_plot = 5)