Построение временных рядов с картой () из мур - PullRequest
0 голосов
/ 19 мая 2018

Я пытаюсь сгенерировать сюжет, похожий на этот:

enter image description here

Здесь приведено пошаговое руководство -> https://medium.com/@erickramer/beautiful-data-science-with-functional-programming-and-r-a3f72059500b

Однако код, предоставленный на этом сайте, не создает для меня сюжет, вместо этого я получаю эту ошибку:

> forecasts1 = tsdf %>%
+     map(auto.arima) %>%
+     map(forecast, h=10)
Error in is.constant(x) : 
  (list) object cannot be coerced to type 'double'

Это несмотря на то, что я точно воспроизвел их форматирование данных.Вот наши наборы данных для сравнения:

> str(tsdf)
Classes ‘tbl_df’, ‘tbl’ and 'data.frame':   89 obs. of  1 variable:
 $ time_series:List of 89
  ..$ 1_1  : Time-Series  from 2013 to 2017: 8981338 10707490 11410597 10816217 12263765 ...
  ..$ 1_10 : Time-Series  from 2013 to 2017: 12645212 13510638 13133558 13542970 16074675 ...
  ..$ 1_2  : Time-Series  from 2013 to 2017: 19028892 20626896 19952328 20865263 22547313 ...
  ..$ 1_3  : Time-Series  from 2013 to 2017: 7081624 8317481 8374427 8330653 9643845 ...
  ..$ 1_4  : Time-Series  from 2013 to 2017: 25421637 30934941 30756101 27977317 32417608 ...

И предоставленные примеры данных (на которых код работал, согласно веб-сайту):

> str(time_series)
List of 9
 $ Germany    : Time-Series [1:52] from 1960 to 2011: 684721 716424 749838   ...
 $ Singapore  : Time-Series [1:52] from 1960 to 2011: 7208 7795 8349   ...
 $ Finland    : Time-Series [1:37] from 1975 to 2011: 85842 86137 86344   ...

Я не могу показатьсяпонять это, хотя это может иметь какое-то отношение к тому факту, что у их временных рядов есть одна надежная конечная точка, но у моих временных рядов есть несколько разных месячных конечных точек.

Любая помощь с этим очень приветствуется!

* ОБНОВЛЕНИЕ *

После применения предложения Akruns я сохранил исключительно вектор временных рядов в спискевот так:

tsdf <- akrun %>%
select(time_series)

Я тогда подгоняю модель вот так:

tsdf$time_series %>% map(auto.arima) %>%
map(forecast, h=12)

... а потом сюжет ...

... выглядит ужасно.

enter image description here

Нужно ли преобразовывать шкалу y_axis?Или сделать какое-то различие с данными, прежде чем строить ариму?Очень ценю любые предложения!

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...