mixed-models

    1热度

    1回答

    我试图从在其他软件中指定的模型,HLM,在R上的输出,这里匹配:http://justpaste.it/q10n 到目前为止我试图将模型(带有相关截距的随机坡度)时,不匹配的开固定效应: m1 <- lmer(formula = mathach ~ 1 + freered + (1 + hrs | school), data = dat) 这里是将R summary(m1) Linear mix

    0热度

    1回答

    我在R中有一个线性混合效应模型(lme4),我想使用MixedModels包在Julia中运行它。我的问题是,Julia的MixedModels不允许在lme4中使用偏移量变量。有什么建议么? 例子:Y ~ X1 + X2 + X3 + X2|Group + X3|Group + offset(X4)

    0热度

    1回答

    我注意到,当使用包含因子类型预测变量的lme4包中的lmer函数指定模型时,后缀指示预测器是因子水平的字符串,如这里用于治疗的情况下: library(afex) data(obk.long) m1 <- lmer(value ~ treatment + (1|id), obk.long) summary(m1) Fixed effects: Estimate Std.

    0热度

    1回答

    我正在组织使用混合模型的数据。 关于nr个观测数据是不平衡的。 它看起来是这样的: Id <- c("A","A","A","A", "B", "B", "B", "C", "C", "D","D","D","D","D","D") v1 <- c(5,5,5,5,8,8,8,15,15,4,4,4,4,4,4) v2 <- c(3,3,3,3,5,5,5,10,10,8,8,8,8,8,8)

    -1热度

    1回答

    我有具有以下参数的混合模型: 组1 为组2 随机不同的斜率和截距项A的斜率和截距术语这是由组索引效果/组内的对象 有一种方法用PROC MIXED建模此?我似乎无法弄清楚如何为两组获得不同的斜坡/截距。

    0热度

    2回答

    我有一个变量yi,代表了一系列不同研究(Site)随时间推移的治疗效果nyears。还有两个分组因子,分别为两个级别:N(Nhigh/Nlow)和Myc(AM/ECM)。我需要知道yi随着时间的推移nyears显示了显着的正向还是负向趋势,并且如果趋势在子组N x Myc之间变化。 混合效应模型显示显著三重互动nyears * N * Myc library(lme4) library(car)

    2热度

    2回答

    因此,我对R相当陌生。我正在研究一个项目,其中有一个数据集用于评估鸟类年龄。我们有来自95个人的> 400,000个观测值。我的任务是这样的: “如果图上有一些数据点,那么这个图会更有说服力,你可以从lme4的getcall中获取ranef的数据点,因此,运行模型时不需要年龄,然后,抽出随机效应项(每个个体的值),然后绘制它们(y =随机效应,x =年龄)。 这是有问题的图(它是在Excel中使用

    2热度

    1回答

    我试图运行一个混合效果模型,该模型包含三个带有交互作用和随机截距和斜率的固定效果。我想在glmmadmb指定的模型是: > fit_zipoiss_ambig<-glmmadmb(AmbigCount~Posn.c*mood.c*Valence.c + offset(InputAmbig) + (1+Valence.c|mood.c/Chain), data = Data, zeroInflati

    0热度

    1回答

    我装修的R.混合模型我的模型是: model <- lmer(provision_rate~breeding_type+nestling_age+time+sex:nestling_age)+(1|nest)+(1|individual), data = provision) 性爱是父母的性别。我想绘制交互术语(sex:nestling_age)和provision_rate之间的关系。我不知

    1热度

    1回答

    我在R中的stats::model.matrix函数上贴了标签。在描述中,它会创建一个设计矩阵。它给了我足够多的行数,这并不符合我的数据中的观察数量,也不符合我的模型中的参数数量。 什么是设计矩阵/模型矩阵? 这是我如何使用它: M03b <- glmer(APMs ~ PrePost + Gf + eyeFRF + (1|content) + (eyeFRF|ID), data=mlmData,