register.newfd {fda} | R Documentation |
This function registers a new functional data object to pre-computed warping functions.
register.newfd(yfd, Wfd,type=c('direct','monotone','periodic'))
yfd |
a multivariate functional data object defining the functions to be
registered with Wfd .
|
Wfd |
a functional data object defining the registration functions to be used
to register yfd . This can be the result of either landmarkreg
or register.fd .
|
type |
indicates the type of registration function.
|
Only shift registration is considered for the periodic case.
a functional data object defining the registered curves.
Ramsay, James O., and Silverman, Bernard W. (2006), Functional Data Analysis, 2nd ed., Springer, New York.
Ramsay, James O., and Silverman, Bernard W. (2002), Applied Functional Data Analysis, Springer, New York, ch. 6 & 7.
# Register the lip data with landmark registration, then register the first # derivatives with the pre-computed warping functions. # Lip data: lipfd <- smooth.basisPar(liptime, lip, 6, Lfdobj=int2Lfd(4), lambda=1e-12)$fd names(lipfd$fdnames) <- c("time(seconds)", "replications", "mm") # Landmark Registration: lipmeanmarks <- apply(lipmarks,2,mean) wnbasis <- 6 wnorder <- 4 wbreaks <- c(0,lipmeanmarks,0.35) warpbasis <- create.bspline.basis(nbasis=wnbasis, norder=wnorder, breaks=wbreaks); WfdPar <- fdPar(fd(basisobj=warpbasis), 2, 1e-4) lipreglist <- landmarkreg(lipfd, as.matrix(lipmarks), lipmeanmarks, WfdPar) # And warp: Dlipfd = deriv.fd(lipfd,Lfdobj=1) Dlipregfd = register.newfd(Dlipfd,lipreglist$warpfd,type='direct')