telesca_warp: Telesca Hierarchical Warping

In wzhorton/warptk: Warping Toolkit

Description

Implements a method of bayesian hierarchical curve registration described in a paper by Telesca (date).

Usage

telesca_warp(
  y_list,
  niter = 1000,
  nburn = 1000,
  int_q = 5,
  int_p = 20,
  asig = 0.1,
  bsig = 0.1,
  at = 0.1,
  bt = 0.1,
  al = 0.1,
  bl = 0.1,
  progress = TRUE,
  debug = FALSE
)

Arguments

y_list	list of curves with equal length
niter, nburn	MCMC iterations and burns. Total run time is the sum.
int_q	number of internal knots specified for the warping function
int_p	number of internal knots specified for the data curves
asig, bsig, at, bt, al, bl	hyperparameters for inverse gamma priors on model variances
progress	logical; indicates whether or not to print percent progress.
debug	logical; stops function right before return and enters debug mode.

Value

a list of warped curves, the estimated mean curve, and a vector of MCMC acceptance rates. Curve list elements are vectors with length matching the elements of y_list.

wzhorton/warptk documentation built on Nov. 14, 2020, 10:01 p.m.