- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
我有一个混合模型,数据如下所示:
> head(pce.ddply)
subject Condition errorType errors
1 j202 G O 0.00000000
2 j202 G P 0.00000000
3 j203 G O 0.08333333
4 j203 G P 0.00000000
5 j205 G O 0.16666667
6 j205 G P 0.00000000
每个受试者都为错误类型(O 或 P)提供两个数据点,并且每个受试者处于条件 G (N=30) 或 N (N=33) 中。 errorType 是一个重复变量,Condition 是一个中间变量。我对主要影响和交互作用都感兴趣。所以,首先是方差分析:
> summary(aov(errors ~ Condition * errorType + Error(subject/(errorType)),
data = pce.ddply))
Error: subject
Df Sum Sq Mean Sq F value Pr(>F)
Condition 1 0.00507 0.005065 2.465 0.122
Residuals 61 0.12534 0.002055
Error: subject:errorType
Df Sum Sq Mean Sq F value Pr(>F)
errorType 1 0.03199 0.03199 10.52 0.001919 **
Condition:errorType 1 0.04010 0.04010 13.19 0.000579 ***
Residuals 61 0.18552 0.00304
条件不重要,但错误类型和交互很重要。
但是,当我使用 lmer 时,我得到了一组完全不同的结果:
> lmer(errors ~ Condition * errorType + (1 | subject),
data = pce.ddply)
Linear mixed model fit by REML
Formula: errors ~ Condition * errorType + (1 | subject)
Data: pce.ddply
AIC BIC logLik deviance REMLdev
-356.6 -339.6 184.3 -399 -368.6
Random effects:
Groups Name Variance Std.Dev.
subject (Intercept) 0.000000 0.000000
Residual 0.002548 0.050477
Number of obs: 126, groups: subject, 63
Fixed effects:
Estimate Std. Error t value
(Intercept) 0.028030 0.009216 3.042
ConditionN 0.048416 0.012734 3.802
errorTypeP 0.005556 0.013033 0.426
ConditionN:errorTypeP -0.071442 0.018008 -3.967
Correlation of Fixed Effects:
(Intr) CndtnN errrTP
ConditionN -0.724
errorTypeP -0.707 0.512
CndtnN:rrTP 0.512 -0.707 -0.724
所以对于 lmer,Condition 和交互很重要,但 errorType 不重要。
此外,lmer 的结果与 glm 的结果完全相同,这让我相信出了什么问题。
谁能帮我理解为什么它们如此不同?我怀疑我使用的 lmer 不正确(尽管我已经尝试了许多其他版本,例如 (errorType | subject) 并得到类似的结果。
最佳答案
我认为答案是两种方法都采用不同的方法来处理方差。 ANOVA 划分方差,组(重复测量)变量只是这样做的另一个考虑因素。 ANOVA 假设组方差或同方差的同质性,如果明显违反此假设,则 ANOVA 可能不是正确的方法。
另一方面,Lmer 本质上是一个用于多级建模的函数。在多级框架内,您可以显式建模方差,引入固定效应和随机效应(基本上是方差)之间的区别。异方差在这里不是问题。
另一种看待它的方式是,ANOVA 采用无池化方法(每条线都是分开的),Lmer 采用部分池化方法(线共享一些信息)。
此外,ANOVA 使用 OLS 估计,而 lmer 使用 ML 版本(在您的情况下为 REML)。
这是我能解释的最好的,这至少应该足以让您在自己的问题研究中朝着正确的方向前进。但是,要获得更详尽的答案,您可能确实想在 CrossValidated 上提问。
关于r - 为什么 aov 和 lmer 之间存在巨大差异?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/20959054/
我有一个包含 142 个数据条目的数据集:两次测量 121 个人(两年,治疗前和治疗后,年份 = 0 或 1),第二年有 46 个人在处理区,其余在控制区绘图(处理 = 0 或 1)。以下是一些示例数
我正在尝试编写一个函数来收集我在脚本中经常使用的一些调用 我在我的例子中使用了 lme4 包的 sleepstudy 数据 这是我开始使用的功能(的简化版本): trimModel1 <- funct
我经常遇到这个问题:我想用约束拟合多级回归。我不知道该怎么做。我通常最终会使用 lavaan,因为它允许对回归系数设置约束。但它仍然不能有随机斜率模型(只有随机截距,事实是我也不知道如何在 lavaa
我有一个具有以下结构的数据框: > t str(t) 'data.frame': 699 obs. of 7 variables: $ Awns : int 0
我正在尝试使用 lmer 函数运行混合效果模型。我的实验包括使用一些相同的个体在不同温度下的代谢率(一些缺失数据)。文本文件的结构如下所示: > str(data.by.animal) 'data.f
lme4 包中的函数lmer 默认使用minqa 包中的bobyqa 作为优化算法。 根据以下帖子https://stat.ethz.ch/pipermail/r-sig-mixed-models/2
我正在使用 lme4 运行混合模型在 R: full_mod3=lmer(logcptplus1 ~ logdepth*logcobb + (1|fyear) + (1 |flocation), da
我有一些纵向数据,我想从中获得指定时间的预测均值。该模型包括 2 个项,它们的交互作用和时间变量的样条项。当我尝试获取预测均值时,我得到“mm %*% fixef(m4) 中的错误:不一致的参数” 我
我目前正在尝试帮助一位同事,但我根本找不到解决方案。所以我希望其他人可以帮助我们。 我有一个数据集,其中包含使用不同研究设计评估的权重数据,针对不同研究中的不同物种(一项研究包括多种设计和多种物种)。
我正在为具有 4 个级别的预测器 root.type 上的单个响应变量运行线性混合模型;当我运行模型时,我只想要有关整个因素的信息,但它一直将其拆分为多个级别。有什么想法吗? Ca.auto |t|)
当您有一个包含大量因素和相互作用的多级模型时,固定效应矩阵的相关性大小可能会变得非常大且不清楚。 我可以使用 symbolic.cor=T打印方法中的参数以更清晰地打印摘要,如下所示: ratbrai
我想使用回归模型而不是“方差分析”(AOV)函数在 R 中运行重复测量方差分析。 这是我的 3 个主题内因素的 AOV 代码示例: m.aov<-aov(measure~(task*region*ac
我对混合模型使用了以下语法,然后 step 但它不起作用。 它通常是这样工作的还是我实际上不能使用 lmer 使用反向消除?谢谢! fullmodel<-lmer(Eeff~NDF+ADF+CP+NE
我有四个位置和四个基板的移植实验(取自每个位置)。我已经确定了每个位置和基质组合中每个种群的存活率。本实验重复 3 次。 我创建了一个 lmm 如下: Survival.model <- lmer(S
我需要提取 standard error来自 lmer 输出的方差分量. library(lme4) model <- lmer(Reaction ~ Days + (1|Subject), slee
我在 R 装了一个模型与 lmer() -功能来自 lme4包裹。我缩放了因变量: mod fixef(mod) (Intercept) X1 X2
我已经拟合了一个 lmer 模型,现在我正在尝试根据实际系数而不是缩放系数来解释系数。 我的顶级模特是: lmer(logcptplus1~scale.t6+scale.logdepth+(1|lo
我刚刚将 lme4 更新到 1.0-4 版,当我运行 lmer() 时,我之前收敛的混合效果模型现在会打印此警告: Warning message: In (function (fn, par, lo
我曾经使用下面的代码来计算 lmer 模型的标准化系数。但是,随着 lme 的新版本,返回对象的结构发生了变化。 如何调整函数 stdCoef.lmer 以使其与新的 lme4 版本一起使用? # I
我正在尝试在 drake 计划中安装一些 lme4::lmer 模型,但出现错误 'data' not found, and some variables missing from formula e
我是一名优秀的程序员,十分优秀!