时间戳

首页 > 解决方案 > anova_test 在主题 ANOVA 中没有以三种方式返回 Mauchly

anova - anova_test 在主题 ANOVA 中没有以三种方式返回 Mauchly

问题描述

我正在使用一个名为 sleep 的数据集（在此处找到：https ://drive.google.com/file/d/15ZnsWtzbPpUBQN9qr-KZCnyX-0CYJHL5/view ）在受试者 ANOVA 中运行三种方式，比较基于刺激、剥夺、和时间。在使用来自 rstatix 的 anova_test 之前，我已经成功地做到了这一点。我想查看球形输出，但它没有出现在输出中。我已经在主题数据集中想出了其他三种方式，所以我不确定为什么会这样。这是我的代码：

anova_test(data = sleep, dv = Performance, wid = Subject, within = c(Stimulation, Deprivation, Time))

我还尝试将它保存到一个对象并使用 get_anova_table，但这看起来并没有什么不同。

sleep_aov <- anova_test(data = sleep, dv = Performance, wid = Subject, within = c(Stimulation, Deprivation, Time))
get_anova_table(sleep_aov, correction = "GG")

这是我从互联网上提取的理想数据集，因此我开始认为数据的 W 为 1（完美的球形度），因此 rstatix 正在跳过此输出。这是 anova_test 做的吗？

这也是我使用返回 Mauchly 的数据集的代码：

weight_loss_long <- pivot_longer(data = weightloss, cols = c(t1, t2, t3), names_to = "time", values_to = "loss")
weight_loss_long$time <- factor(weight_loss_long$time)
anova_test(data = weight_loss_long, dv = loss, wid = id, within = c(diet, exercises, time))

标签： anovarstatix

解决方案

根本不是专家，但可能是因为您的因素只有两个级别。

从 anova_summary() 帮助：

“值返回类 anova_test 的对象，包含用于独立测量 ANOVA 的 ANOVA 表的数据框。但是，对于重复/混合测量 ANOVA，它是一个包含以下组件的列表返回：

ANOVA：包含 ANOVA 结果的数据框

Mauchly 的球形检验：如果存在任何具有 2 个以上水平的 Ss 内变量，则包含 Mauchly 球形检验结果的数据框。仅报告具有超过 2 个级别的效果，因为球形度必然适用于仅具有 2 个级别的效果。

球形校正：如果存在任何 Ss 内变量，则数据框包含 Greenhouse-Geisser 和 Huynh-Feldt epsilon 值，以及相应的校正 p 值。"

推荐阅读