Я хотел установить ковариационную структуру Тёплица, используя GLS. Week
- моя переменная времени, а group
- индикатор двоичной группы. Я понимаю, что должен указать дополнительный аргумент в функции CorARMA
(p
и q
), но я не понимаю, что это такое? Может кто-нибудь объяснить?
gls_par_top <- gls(data=df,weight ~ week + group2 + group2*week,corr=corARMA(form=~1 | id,p=?,q=?),weights = varIdent(form=~1|week))