我一直使用SPSS进行重复测量分析,但是现在我想切换到R。我发现'ez'软件包非常容易理解。但是,我不确定某些问题。
这是一个示例数据集:
install.packages('ez')
library(ez)
bla <- c(1,2,5,8,88,99,102,150)
cat <- rep(c(1,2),4) #group var
ID <- 1:8
RT1 <- c(1500,1600,1700,1333,1444,1555,1666,1888)
RT2 <- c(500,600,800,433,544,655,666,888)
dat <- data_frame(ID, cat,bla, RT1, RT2)
#change to long format
dat2<- gather(dat, var_name,RT,-ID, -cat)
dat2 <- arrange(dat2 , ID, cat, var_name)
rt_anova = ezANOVA(
data = dat2
, dv = .(RT)
, wid = .(ID)
, within = .(var_name)
, between = .(cat)
)
这是我的输出:
> rt_anova
$`ANOVA`
Effect DFn DFd F p p<.05 ges
2 cat 1 6 1.849095e-02 8.962836e-01 2.988056e-03
3 var_name 1 6 2.166000e+03 6.594391e-09 * 9.085505e-01
4 cat:var_name 1 6 3.710802e-31 1.000000e+00 1.702068e-33
是否可以使用e ^ x舍入输出值而不是精确值?这样我得到F值(F <0.01)的两位小数,p值(p <.001)的两位小数,eta²(eta²<.01)的两位小数。
如何添加事后分析?我的首选方法是对所有重大影响进行sidak调整比较。
最好有一种简单的方法已包含在ezANOVA中。还是有一个简短而简单的解决方法?