У меня есть данные временного ряда, которые хранятся в большом фрейме данных. Существует 7 переменных, которые были измерены на нескольких сотнях индивидуумов (таким образом, 7 отдельных временных рядов для каждого индивидуума), и существует около 2000 наблюдений для каждого индивидуума (в 7 раз для каждой переменной). Этот вопрос относится к нескольким аспектам моего анализа, но в настоящее время я делаю EDA, разлагаю каждый временной ряд и пытаюсь определить, подходит ли аддитивная модель или мультипликативная модель для каждой серии. Я хотел бы выполнять эти операции без необходимости поднабора и создания отдельного объекта типа ts для каждой индивидуальной / переменной комбинации.
Это создает фрейм данных, который отформатирован аналогично моему
set.seed(123)
ID = factor(letters[seq(15)])
Time = c(1000,1200,1234,980,1300,1020,1180,1908,1303,
1045,1373,1111,1097,1167,1423)
df <- data.frame(ID = rep(ID, Time), Time = sequence(Time))
df[paste0('Var', c(1:7))] <- rnorm(sum(Time))