# NOT RUN {
# Here is the code that sets up the single forcing term for the
# refinery data examnple:
TimeData <- RefineryData[,1]
TrayData <- RefineryData[,2]
ValvData <- RefineryData[,3]
conbasis <- create.constant.basis(c(0,193))
alphafdPar <- fdPar(conbasis)
Valvbreaks <- c(0,67,193)
Valvnbasis <- 2
Valvnorder <- 1
Valvbasis <- create.bspline.basis(c(0,193), Valvnbasis, Valvnorder, Valvbreaks)
Valvfd <- smooth.basis(TimeData, ValvData, Valvbasis)$fd
# Fterm Fields: funobj parvec estimate Ufd factor
FTerm <- make.Fterm(alphafdPar, 1.0, TRUE, Valvfd, 1)
FList <- vector("list", 1)
FList[[1]] <- FTerm
#
# For other examples of the use of this function,
# see these examples in the description of the function \code{make.Variable}.
#
# The single forcing term of the head impact data
#
# The cosine and constant forcing terms for the average temperature in
# Montreal.
#
# The 20-step forcing function and its coefficient for the "fda" script data.
#
# The set point forcing term for control variable of the cruise control data.
#
# }
Run the code above in your browser using DataLab