1
GLS
有誰知道如何與這個錯誤解決問題:錯誤R更多有關corAR1,不是唯一的
library(nlme)
gls(Number.of.Fish~Julian.Date+Temperature+Size.Class,na.action=na.omit,data=mydata,correlation=corAR1(form=~Julian.Date))
Error in Initialize.corAR1(X[[1L]], ...) :
covariate must have unique values within groups for "corAR1" objects
感謝您的幫助。不幸的是,它是計數數據 – Francesca 2014-09-10 19:43:52
** mgcv **包中有這樣做的方法。一個是'bam()',它可以適應具有AR(1)誤差模型的GAM。不知道這將如何處理重複的觀察。否則,您可以使用'gamm()',它通過MASS包中的'glmmPQL()'適用,並使用** nlme **中的關聯結構。最後,'魔術'頁面(IIRC)包括一個使用'gam()'和一些手工操作來擬合AR(1)的例子。 – 2014-09-10 19:59:42