require(fda)
# Load Montana data
montana_data <- loadMontanaData()
# Extract variables
Temp <- montana_data$Temp
NDVI <- montana_data$NDVI
# Create a list for Montana data
Montana_Data <- list(Temp / sd(Temp), NDVI)
# Define basis functions
bs1 <- create.bspline.basis(c(0, 23), 11)
bs2 <- create.bspline.basis(c(1, 33), 13)
bs2d <- list(bs2, bs2)
bsmv <- list(bs1, bs2d)
# Convert to funts object
Y <- funts(X = Montana_Data, basisobj = bsmv,
start = as.Date("2008-01-01"),
end = as.Date("2013-09-30"),
vnames = c("Normalized Temperature (\u00B0C)" , "NDVI"),
dnames = list("Time", c("Latitude", "Longitude")),
tname = "Date"
)
Run the code above in your browser using DataLab