data("park_visitation")
data("flickr_userdays")
park <- "YELL"
pud_ts <- ts(park_visitation[park_visitation$park == park,]$pud, start = 2005, freq = 12)
pud_ts <- log(pud_ts)
nps_ts <- ts(park_visitation[park_visitation$park == park,]$nps, start = 2005, freq = 12)
nps_ts <- log(nps_ts)
nps_decomp <- auto_decompose(nps_ts)
trend_proxy <- log(flickr_userdays)
vm <- visitation_model(pud_ts,trend_proxy,ref_series = nps_ts)
plot(vm)
Run the code above in your browser using DataLab