【问题标题】:Nested design and F-, p- & R²-values in eHOF packageeHOF 封装中的嵌套设计和 F、p 和 R² 值
【发布时间】:2014-11-04 14:08:19
【问题描述】:

我想在 eHOF 包的 HOF 函数中考虑我的嵌套研究设计。 我们采样了嵌套在样本块中的点和嵌套在村庄中的点(随机)。

到目前为止,我使用 GLMM 来描述栖息地类型(栖息地、解释性、固定)对物种丰富度(SpR、响应变量、固定)的影响。因此我使用了 lme (of nlme-package) 我可以将嵌套学习设计合并到函数中:

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

【问题讨论】:

    标签: r nested hierarchical logistic-regression random-effects


    【解决方案1】:

    不,目前没有。您可以尝试使用nlme() 手动拟合非线性模型,也许在 eHOF 包代码中查找有关如何指定公式等的提示,但该包直接使用非线性优化器,而不是nlme(),因此您需要进行翻译工作以获得nlme 可以使用并计算出其他细节的模型公式。

    【讨论】:

    • 感谢您提出这个想法!但不幸的是,这种方法对我所有不同的物种进行测试会非常耗时。我一直在寻找其他可能性来整合嵌套设计。 eHOF 软件包的作者随后运行了包含和不包含嵌套研究设计的 GLM,并比较了模型的 AIC。然而,如何通过 eHOF 包获得 p 值和 F 值?
    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 2018-12-19
    • 2019-12-09
    • 2019-04-13
    • 2019-07-13
    • 1970-01-01
    • 2011-03-14
    • 2015-10-17
    相关资源
    最近更新 更多