nlme

    0热度

    1回答

    R对于1级和2级变量的LME协变量的排序？

    我有R中1级和2个变量我dataframe (df)纵向数据： ID Year Gender Race MathScore DepressionScore MemoryScore 1 1999 M C 80 15 80 1 2000 M C 81 25 60 1 2001 M C 70 50 75 2 1999 F C 65 15
    r panel-data nlme longitudinal 2017-03-22

    0热度

    1回答

    如何报告nlme混合效应模型的总体结果

    我想报告nlme包中的一个因子lme的结果。我想知道A对y的整体影响。要做到这一点我将与空模型的模型对比： m1 <- lme(y~A,random=~1|B/C,data=data,weights=varIdent(form = ~1|A),method="ML") m0 <- lme(y~1,random=~1|B/C,data=data,weights=varIdent(form = ~
    r lme4 mixed-models nlme 2016-12-14

    0热度

    1回答

    错误使用nlme;畸形因子

    我正在使用nlme软件包来学习多层次模型，并遵循教科书“Discovering Statistics Using R”的示例。 Mixed Models Code 的数据集是蜜月Period.dat，在他们的同伴网站也可下载。 Data Set - Multilevel Models require(nlme) require(reshape2) satisfactionData = read
    r nlme code-statistics 2016-12-20

    0热度

    1回答

    Wald Chi-Squared在一个混合效应模型中的两个变量之间的测试r

    我试图完成纵向数据分析中的作业。 现在的问题是比较模型中年龄的横截面和纵向效应（基线横截面：基底，纵向年龄：年龄变化）的差异。 模型予代码等： 适合< -lme（logfev1〜baseage + agechange +高度，随机=〜1 | ID，相关性= corAR1（形式=〜访问| ID），logfev1） 在Stata的，我们只需要如下代码：测试baseage = agechange，那么答
    r mixed-models nlme longitudinal 2017-02-15

    0热度

    1回答

    如何将glmer（）调用转换为lme（）;并包括用于随机效应的列表（）

    我以前在软件包lme4中使用glmer（）运行混合模型分析。我现在想在包nlme中使用lme（）运行相同的分析。这是因为随后使用的函数需要输出或调用lme（）混合模型。 随后使用的函数尝试使用函数segmented.lme（）在数据中查找断点。这个函数的代码可以在这里找到：https://www.researchgate.net/publication/292986444_segmented_mi
    r lme4 mixed-models nlme 2017-02-20

    0热度

    1回答

    nlme：使用CSH协方差模型拟合混合模型

    我正尝试使用nlme软件包在R中使用重复测量（MMRM）模型拟合混合模型。 数据的结构如下： 每个患者属于三个组（grp）之一，并被分配到一个治疗组（trt）。 患者结果（y）在6次访问（访问）期间测量。 我想在不同访问中使用具有异构差异的复合对称模型（如SAS的PROC MIXED的CSH类型，https://support.sas.com/documentation/cdl/en/statug
    r nlme 2017-04-13

    3热度

    1回答

    lme模型的用户定义函数拟合：错误

    我开始编写一个函数，它可以与nlme建立线性混合模型。我遇到一个错误：Error in eval(expr, envir, enclos) : object 'value' not found，我认为这是由于R不知道在哪里找到数据帧变量（例如，value）。如果这实际上是错误发生的原因，那么我如何知道value和timepoint属于Dat中变量（在下面的可重现代码中）的函数？ require(n
    r nlme 2016-09-30

    2热度

    1回答

    错误，同时运行与NLME

    我试图运行GLS功能，但我不断收到此错误信息GLS： “错误model.frame.default（公式=〜+变种FS，数据=列表（MINBIO15 = C（37L，： 对象不是矩阵” 我的数据看起来像这样： MINBIO15 MAXBIO15 FS Achyranthes_aspera 37 117 0 Achyropsis_avicularis 28 86 0 A
    r nlme 2016-10-03

    1热度

    1回答

    在R中指定随机效果之间的非线性关系

    我正在使用多级模型尝试描述纵向变化中的不同模式。当随机效应完全相关时，Dingemanse et al (2010)描述了“扇出”模式。然而，我发现当随机效应之间的关系是非线性的但在观察到的时间间隔内单调递增时会出现类似的模式。在这种情况下，随机效应并不完全相关，而是由函数描述。 请参阅下面的示例以获取此示例。这个例子仍然具有很高的截距 - 斜率相关性（> .9），但是可以得到低于.7的相关性，同
    r mixed-models multi-level nlme r2jags 2017-03-02

    4热度

    1回答

    随机效应变量是否会自动作为R中lmer（或lme）中的因子？

    据我所知，具有混合效果模型作为随机效果连续或数字变量没有多大意义（例如，请参阅here）。 但我不知道是如果lme4::lmer或R中nlme::lme有意阻止你这样做...... 具体来说，就是我要问的是：如果我提供lmer（或lme）任何非因子（非分类）变量作为随机效应，函数是否会自动将其作为一个因子来处理？ 插入factor()直接进入11聚物（如使用时lm通常的方法）产生以下错误：虽然上述
    r lme4 categorical-data mixed-models nlme 2016-10-10
最新问题