fdasrvf (version 1.9.4)

curve_principal_directions: Curve PCA

Description

Calculate principal directions of a set of curves

Usage

curve_principal_directions(betamean, mu, K, mode = "O", no = 3, N = 5)

Arguments

betamean

array (n,T) of mean curve

mu

array (n,T) of mean srvf

K

array (2*T,2*T) covariance matrix

mode

Open ("O") or Closed ("C") curves

no

number of components

N

number of samples on each side of mean

Value

pd list describing principal directions

References

Srivastava, A., Klassen, E., Joshi, S., Jermyn, I., (2011). Shape analysis of elastic curves in euclidean spaces. Pattern Analysis and Machine Intelligence, IEEE Transactions on 33 (7), 1415-1428.

Examples

# NOT RUN {
data("mpeg7")
out = curve_srvf_align(beta[,,1,1:2],maxit=2) # note: use more shapes, small for speed
K = curve_karcher_cov(out$betamean, beta[,,1,1:2])
pd = curve_principal_directions(out$betamean, out$q_mu, K)
# }