2014-09-10 46 views
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 

回答

2

您已經採取同一Julian.Date兩個或兩個以上的樣品,你不能有這樣的在corAR1()corStruct()對象中。

你如何解決這個問題將取決於你的採樣設計。這些價值來自兩個不同的網站嗎?如果是這樣site也許應該在corAR1()

corAR1(form ~ Julian.Date | site) 

除此之外,我疲於應付發生在同一天兩個採樣想象等距的觀察時間。

如果您的回答是一個計數,您可能需要再次考慮使用gls(),而不是使用離散計數數據的模型。

+0

感謝您的幫助。不幸的是,它是計數數據 – Francesca 2014-09-10 19:43:52

+0

** mgcv **包中有這樣做的方法。一個是'bam()',它可以適應具有AR(1)誤差模型的GAM。不知道這將如何處理重複的觀察。否則,您可以使用'gamm()',它通過MASS包中的'glmmPQL()'適用,並使用** nlme **中的關聯結構。最後,'魔術'頁面(IIRC)包括一個使用'gam()'和一些手工操作來擬合AR(1)的例子。 – 2014-09-10 19:59:42

相關問題