有人可以为R中的gamm
函数推荐一种包含自相关的方法(如gamm(...,correlation=corAR1())
中的mgcv
函数),但也建议用于处理二进制响应数据? gamm()
帮助文件有关于将其用于二进制数据的明确警告。
" gamm对二进制数据表现不佳,因为它使用PQL
。最好将gam
与s(...,bs="re")
条款或gamm4
一起使用。"
最好,我可以说,gamm4
目前还没有实现自相关。
希望我只是遗漏了一些明显的东西。
谢谢!