Расчет времени для присвоения атрибута tsp - PullRequest
0 голосов
/ 30 октября 2019

Я хочу присвоить атрибут tsp для вектора

#dummy data
a <- c(600L, 200L, 0L, 100L, 200L, 100L, 500L, 1000L, 200L, 800L, 
       350L, 1050L, 450L, 750L, 250L, 100L, 100L, 300L, 600L, 650L, 
       450L, 350L, 250L, 450L, 200L, 150L, 250L, 12250L, 650L, 50L, 
       50L, 300L, 350L, 400L, 550L, 1300L, 250L, 100L, 150L, 150L, 150L, 
       150L)

Временные ряды начинаются в апреле 2016 года и заканчиваются в сентябре 2019 года, первого числа каждого месяца. Частота ежемесячная, поэтому 12. Я хочу, чтобы правильные десятичные значения использовались для каждого year_month. Как я могу вычислить это, не округляя результаты, поскольку это вызывает ошибку?

library(lubridate)
start <- decimal_date(as.Date("2016-04-01"))  #2016.249
end   <- decimal_date(as.Date("2019-09-01"))  #2019.666

tsp(a) <- c(start, end, 12) - not working
tsp(a) <- c(2016.25, 2019.66666666667, 12) #using correct tsp info    
...