# NOT RUN {
data(Balt) # Baltimore mean sea level record
ts1 <- ts(Balt[2], start = Balt[1, 1]) # create time series input object
g <- gap.fillview(ts1, station_name = "Baltimore", fillgaps = 1) # SSA gap fill
# }
# NOT RUN {
t <- custom.trend(g, station_name = "Baltimore (USA)", iter = 500, trend = c(1,2),
vlm = 0.6)
# }
# NOT RUN {
data(t)
str(t) # check structure of object
# }
Run the code above in your browser using DataLab