## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
  echo = TRUE, message = FALSE, warning = FALSE,
  collapse = TRUE, comment = "#>",
  eval = FALSE
)

## ----load---------------------------------------------------------------------
# library(gdpar)

## ----data---------------------------------------------------------------------
# set.seed(20260526)
# n <- 200L
# x <- rnorm(n)
# theta_ref_true <- 1.5
# a_true <- 0.8
# sigma_true <- 0.3
# y <- theta_ref_true + a_true * x + rnorm(n, sd = sigma_true)
# dat <- data.frame(y = y, x = x)
# head(dat)

## ----fit----------------------------------------------------------------------
# fit <- gdpar(
#   formula = y ~ x,
#   data    = dat,
#   family  = gdpar_family("gaussian"),
#   path    = "bayes"
# )

## ----print--------------------------------------------------------------------
# print(fit)

## ----coef---------------------------------------------------------------------
# coefs <- coef(fit)
# coefs

## ----predict------------------------------------------------------------------
# new_x <- data.frame(x = seq(-2, 2, length.out = 11))
# preds <- predict(fit, newdata = new_x, level = 0.9)
# preds

## ----diagnostics--------------------------------------------------------------
# diag <- diagnostics(fit)
# print(diag)