草庐IT

在NLME和LME4中拟合相同的型号

数据来自这里library(nlme)dat0在这里,我使用使用的非线性混合效应模型nlme在里面nlme包裹。该模型具有4个固定效应和4个随机效应。我在方差互动矩阵上指定了一个对角线结构,每个矩阵patid形成一个小组。library(lme4)deriv_mod在这里,我想使用lme4包裹。从文档看来formula为了nlmer也必须具有梯度组件,因此我使用了deriv功能首先。但是,我不确定如何指定其余参数?这deriv_mod~list(p1~1,b1~1,p2~1,b2~1)+list(p1~1,b1~1,p2~1,b2~1)是指定4个固定效果(在第一个列表对象中)及其相应的4个随机

关于 r:How to translate glmer() call to lme();并包括 list() 用于随机效果

Howtotranslateglmer()calltolme();andincludinglist()forrandomeffects我之前使用lme4包中的glmer()运行了混合模型分析。我现在想使用nlme包中的lme()来运行相同的分析。这是因为随后使用的函数需要输出或调用lme()混合模型。随后使用的函数尝试使用函数segmented.lme()在数据中查找断点。这个函数的代码可以在这里找到:https://www.researchgate.net/publication/292986444_segmented_mixed_models_in_R_code_and_data之前,我使

关于 r:How to translate glmer() call to lme();并包括 list() 用于随机效果

Howtotranslateglmer()calltolme();andincludinglist()forrandomeffects我之前使用lme4包中的glmer()运行了混合模型分析。我现在想使用nlme包中的lme()来运行相同的分析。这是因为随后使用的函数需要输出或调用lme()混合模型。随后使用的函数尝试使用函数segmented.lme()在数据中查找断点。这个函数的代码可以在这里找到:https://www.researchgate.net/publication/292986444_segmented_mixed_models_in_R_code_and_data之前,我使