У меня ежедневные данные за 6 лет. Первоначальные исследования данных показали сезонность дня недели (два дня) и зимних месяцев. Я хотел бы создать временной ряд, включающий эти две сезонности.
Ниже приведен код, который я придумал, но я не уверен, как включить сезонность:
inds<-seq(as.Date("2011-01-01"), as.Date("2016-12-31"),by="day")
myts<-ts(mydataframe,start=c(2011,as.numeric(format(inds[1],"%j"))),frequency=365.25)
Я тоже пробовал msts, но не уверен, прав ли я:
msts <- msts(mydataframe, seasonal.periods = c(7,90),
start = decimal_date(as.Date("2011-01-01")))
Я ценю ваш совет.