# NOT RUN {
# define the harmonic acceleration operator used in the
# analysis of the daily temperature data
Lcoef <- c(0,(2*pi/365)^2,0)
harmaccelLfd <- vec2Lfd(Lcoef, c(0,365))
hmat <- vec2Lfd(matrix(Lcoef, 1), c(0, 365))
# }
# NOT RUN {
all.equal(harmaccelLfd, hmat)
# }
# NOT RUN {
# }
Run the code above in your browser using DataCamp Workspace