function_group_warp_bayes {fdasrvf}

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)