Подгонка и построение гармонической регрессионной модели временного ряда в R - PullRequest
1 голос
/ 30 марта 2019

В R, используя библиотеку HarmonicRegression, я пытаюсь согласовать модель гармонической регрессии с трендом, но я пытаюсь выяснить, как вызвать функцию harmonic.regression и построить построенную модель.

library(TSA)    
library(HarmonicRegression)
data(tempdub)

har_model <- harmonic.regression(tempdub, time(tempdub), Tau = 24, normalize = TRUE, norm.pol = FALSE, norm.pol.degree = 1, trend.eliminate = FALSE, trend.degree = 1)

plot(ts(fitted(har_model), freq=1, start=c(1964,1)), type='l', ylim=range(c(fitted(har_model), tempdub))) 

points(tempdub)

Первая ошибка, с которой я сталкиваюсь:

Error in if (nrow(inputts) != length(inputtime)) stop(paste("Length of time series (inputts):",  : 
  argument is of length zero

Есть идеи, как этого достичь?

...