Я пытаюсь выполнить многомерный анализ прерывистых временных рядов в R, используя пакет "its.analysis".
Я получаю следующую ошибку:
"Ошибка в заказе (c (1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), c (1L,1L, 1L,: невыполненный тип 'list' в 'listgreater'
Дополнительно: Предупреждающие сообщения:
1: В cbind (x, time): количество строк результата не являетсякратное длине вектора (arg 1)
2: в мин. (temp $ depvar): нет не пропущенных аргументов в min; возвращается Inf
3: в max (temp $ depvar):нет не пропущенных аргументов для max; возвращает -Inf "
Код:
r_1 = c(-23,-27,-5,-14,-21,4,-2,-6,-6,0)
cov1 = c(0.6,-0.3,0.1,-0.2,-0.8,-0.4,-0.5,-0.2,-0.1,-1.5)
cov2 = c(0,0,0.5,-0.4,-0.1,-1.2,0.8,-0.9,0,0)
time = 1:10
intn = ts(c(rep(0,4),rep(1,6)),start = 1, end = 10)
#intn = ts(as.factor(c(rep(0,4),rep(1,6))),start = 1, end =10 )
#use ts.union to return ts.object
ts_r = ts.union(r_1,cov1,cov2,time,intn)
#use ts.object as data in itsa.model
model.ym <- itsa.model(data = ts_r,time = "time",depvar = "r_1",interrupt_var = "intn",
covariates = c("cov1","cov2"),bootstrap = FALSE)
#use ts.union to return dataframe
df_r = ts.union(r_1,cov1,cov2,time,intn,dframe = TRUE)
#use dataframe as data in itsa.model
model.ym <- itsa.model(data = df_r,time = "time",depvar = "r_1",interrupt_var = "intn",
covariates = c("cov1","cov2"),bootstrap = FALSE)