Интерполировать высокочастотный временной ряд

Question

У меня есть физические временные ряды в диапазоне данных выборки за 2 года с частотой 30 минут, но есть многократные и широкие интервалы потерянных данных, как вы можете видеть там:

Я пытался с помощью функции na.interp из forecast пакета с плохим результатом (как показано выше):

sapply(dataframeTS[2:10], na.interp)

Я ищу более полезный метод.

ОБНОВЛЕНИЕ:

Вот дополнительная информация о шаблоне, который я хочу записать, конкретно строка data.Этот подвыбор относится к маю.

Answer 1

Возможно, вы захотите попробовать пакет **imputeTS**.Это пакет R, предназначенный для вменения пропущенных значений временных рядов.

Методы na.seadec(), na.seasplit(), na.kalman() могут быть интересны здесь

Есть еще много вариантов алгоритма - вы можетенайдите в этом документе о пакете.

В этом конкретном случае я бы попробовал:

na.seasplit(yourData)

или

na.kalman(yourData)

или

na.seadec(yourData)

Имейте в виду, что, возможно, вам необходимо правильно указывать информацию о сезонности во временном ряду.(вам нужно создать временной ряд (объект ts) и установить параметр frequency)

Тем не менее, возможно, что он вообще не сработает, вам придется попробовать.(если вы можете предоставить данные, я также попробую)