https://github.com/cran/fda
Tip revision: 877347f65135b34c238ac4f27de992aa38d15347 authored by J. O. Ramsay on 03 November 2009, 00:00:00 UTC
version 2.4.0
version 2.4.0
Tip revision: 877347f
deriv.fd.R
deriv <- function(expr, ...)UseMethod('deriv')
deriv.fd <- function(expr, Lfdobj=int2Lfd(1), returnMatrix=FALSE, ...)
{
# Applies linear differential operator LFD to functional data object FD
# and returns the result as functional data object DERIVFD.
# Last modified 8 May 2012 by Jim Ramsay
fdobj <- expr
if (!inherits(fdobj, "fd")) stop(
"Argument FD not a functional data object.")
Lfdobj <- int2Lfd(Lfdobj)
basisobj <- fdobj$basis
nbasis <- basisobj$nbasis
rangeval <- basisobj$rangeval
evalarg <- seq(rangeval[1], rangeval[2], len=10*nbasis+1)
Lfdmat <- eval.fd(evalarg, fdobj, Lfdobj, returnMatrix)
Lfdcoef <- project.basis(Lfdmat, evalarg, basisobj, returnMatrix)
Dfdnames <- fdobj$fdnames
Dfdnames[[3]] <- paste("D",Dfdnames[[3]])
Dfdobj <- fd(Lfdcoef, basisobj, Dfdnames)
return(Dfdobj)
}