2014-09-10 53 views
0

我想考慮一下我的巢式研究設計在eHOF包的HOF功能。 我們採用嵌套在樣本塊和嵌套在村莊(隨機)中的點進行採樣。嵌套設計和F-,對與R 2值在eHOF包

到目前爲止我用GLMM描述上物種豐富度生境類型(生境,說明,固定的)的效果(SPR,響應變量,固定)。因此,我使用的LME(NLME-包)其中,I可以結合嵌套研究設計在功能:

model1<-lme(SpR~Habitat,random=~1|Village/Block,data=birds) # full model 

現在我想描述沿一梯度上SPR的效果。我的HOF功能如下所示:

hof.blocks<-HOF(SpR,Indexc,M=max(SpR),bootstrap=100, 
       family=poisson) 

有沒有辦法將嵌套研究設計合併到HOF中?

此外,我無法獲得的P值,也不F-或R2-值。我忽略了什麼?

套餐:eHOF 版本:1.4 [R版本3.0.2

回答

0

不,不是的時刻。你可以嘗試使用nlme()適合用手非線性模型,也許看在eHOF包代碼提示如何指定公式等,但包裝直接使用非線性優化器,而不是nlme()所以你需要做翻譯工作才能得到nlme可以使用的模型公式,並制定出其他細節。

+0

感謝您對這個想法!但不幸的是,這種方法對於我所有不同的物種進行測試都非常耗時。我正在尋找其他可能性來整合嵌套設計。 eHOF軟件包的作者之後運行了GLM,其中包括嵌套研究設計和不包括嵌套研究設計,並比較了模型的AIC。儘管如此,如何通過eHOF包獲得p值和F值呢? – Lula 2014-09-22 11:14:10