有谁知道为什么SPSS和R包，rstatix和DescTools之间的部分eta平方值有差异？

问题描述

我正在受试者之间进行2（性别）x 2（建议）ANOVA，R和SPSS均报告了相同的ANOVA统计数据：建议：F = 372.012，效果df = 1，错误df = 661；对于“性别x建议”：F = 45.449，效果df = 1，错误df = 2。

在计算部分eta平方时，R（在多个R包中：rstatix和DescTools）报告的Advice为0.221，性别x Advice为0.031。但是，SPSS和使用Lakens的效果大小电子表格（https://osf.io/ixgcd/）都为Advice带来了0.360，为Gender x Advice带来了0.064。

有人知道R是否以不同于标准的方式计算部分eta平方值吗？

这是一个示例数据集，其中仅包含测试问题所需的变量： https://docs.google.com/spreadsheets/d/15AIyIfTi9YgMWM5FTl163uddPx1E19xfaU-vMDRlJuI/edit?usp=sharing

这是我在RStudio中使用的代码：

# load packages
library(haven)
library(rstatix)
library(DescTools)

# read in data
sample_data <- read_sav([insert file location])

# gather Perception1 and Perception2 into 2 groups
sample_data <- sample_data %>% 
gather(key = "Advice",value = "MaleDom",Perception1,Perception2) %>% 
convert_as_factor(ResponseId,Advice)

# rstatix
# compute anova
anova <- aov(MaleDom ~ Gender*Advice,data = sample_data)

# partial eta squared
partial_eta_squared(anova)

# DescTools
# partial eta squared
EtaSq(anova,type = 2,anova = FALSE)

这是我在SPSS中使用的语法：

GLM Perception1 Perception2 BY Gender
/WSFACTOR=advice 2 Polynomial 
/METHOD=SSTYPE(3)
/POSTHOC=Gender(BTUKEY) 
/PLOT=PROFILE(Gender*advice)
/PRINT=DESCRIPTIVE ETASQ HOMOGENEITY 
/CRITERIA=ALPHA(.05)
/WSDESIGN=advice 
/DESIGN=Gender.

注意：我正在使用SPSS 26版和R 3.6.3版。我有一个带有64位操作系统的Windows 10。

解决方法

使用它来定义数据：

 d.dat <- structure(list(ResponseId = c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30),Gender = c("Woman","Woman","Man","Man"),Perception1 = c(3.33,3.67,1.33,5.33,6.67,3.33,4.33,1,6.33,4.67,3.33),Perception2 = c(6,5.67,6.33)),class = "data.frame",row.names = c(NA,-30L)) 

,

感谢您提供示例数据。出于三种可能的原因，您将获得不同的偏eta平方值：

1. 在SPSS中，您指定Advice是重复测量因子，而在R中，您将其视为对象间因子。
2. 您在SPSS中将多项式对比用于对象内因子Advice，而在R中，您正在使用处理对比。
3. 您在SPSS中使用III型平方和，而在R中使用I型平方和。

您可以像这样在R中指定重复测量方差分析（split-plot design）：

anova <- aov(MaleDom ~ Gender*Advice + Error(ResponseId/Advice),data=sample_data)

然后

EtaSq(anova,type=1,anova=FALSE)    # note type=1,not 3!

对于Advice（样本数据）的部分eta平方的结果为0.2659604。这等于使用以下语法的SPSS输出：

GLM Perception1 Perception2 BY Gender
  /WSFACTOR=Advice 2 Simple
  /MEASURE=MaleDom
  /CONTRAST(Gender)=Deviation(1)
  /METHOD=SSTYPE(1)
  /PRINT=ETASQ
  /CRITERIA=ALPHA(.05)
  /WSDESIGN=Advice
  /DESIGN=Gender.

请注意，来自软件包partial_eta_squared()的{​​{1}}不会处理重复测量rstatix对象，但是来自软件包aov()的{​​{1}}与{{ 1}}来自eta_squared(anova)$Eta_Sq_partial。