Я анализирую данные временных рядов эффективности использования биомассы / дождя растений. У меня есть набор данных с переменной отклика и несколькими факторами (интенсивность выпаса), а также копии для каждого фактора (4 фактора, 5 или 7 повторов).
Мне нужно создать объекты временных рядов из каждого отдельного набора данных, который включает в себя годовые данные с 1981 по 2007 год. У меня есть 24 индивидуальных оценки различных графиков за это время.
Я пытаюсь написать цикл для этого, но с треском проваливаюсь.
Набор данных называется dat.
PMR6yrs - это переменная ответа.
ID - это идентификатор реплики.
каким-то образом он возвращает только один объект TS в конце, поэтому я думаю, что проблема заключается в их хранении.
То, что я до сих пор пробовал, это:
tsPMR<- NULL
for (i in 1:24) {
tsPMR6<- ts(dat$PMR_6yrs[dat$ID == i],
start = c(1981), end = c(2007), frequency = 1)
ts<-rbind(tsPMR, tsPMR6)
}
Может ли кто-нибудь помочь? спасибо см