function_group_warp_bayes: Bayesian Group Warping

In fdasrvf: Elastic Functional Data Analysis

Bayesian Group Warping

Description

This function aligns a set of functions using Bayesian SRSF framework

Usage

function_group_warp_bayes(
  f,
  time,
  iter = 50000,
  powera = 1,
  times = 5,
  tau = ceiling(times * 0.04),
  gp = seq(dim(f)[2]),
  showplot = TRUE
)

Arguments

f	matrix (N x M) of M functions with N samples
time	sample points of functions
iter	number of iterations (default = 150000)
powera	Dirichlet prior parameter (default 1)
times	factor of length of subsample points to look at (default = 5)
tau	standard deviation of Normal prior for increment (default ceil(times*.4))
gp	number of colors in plots (defaults seq(dim(f)[2]))
showplot	shows plots of functions (default = T)

Value

Returns a list containing

f0	original functions
f_q	f aligned quotient space
gam_q	warping functions quotient space
f_a	f aligned ambient space
gam_a	warping ambient space
qmn	mean srsf

References

Cheng, W., Dryden, I. L., and Huang, X. (2016). Bayesian registration of functions and curves. Bayesian Analysis, 11(2), 447-475.

Examples

## Not run: 
  out <- function_group_warp_bayes(simu_data$f, simu_data$time)

## End(Not run)

fdasrvf documentation built on Nov. 19, 2023, 1:09 a.m.