0

我正在使用一个名为 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))
4

1 回答 1

0

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

从 anova_summary() 帮助:

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

ANOVA:包含 ANOVA 结果的数据框

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

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

于 2021-07-16T17:32:36.430 回答