创建样本数据:
norm_power <- rnorm(2298)
phoneme <- c("t", "v", "s", "p")
intensity <- c("loud", "soft")
unique_ID <- as.factor(rep(c("elec1", "elec2", "elec3", "elec4", "elec5", "elec6", "elec7",
"elec8", "elec9", "elec10"), length.out = 2298))
Subject <- as.factor(rep(c("DBS1", "DBS2", "DBS3", "DBS4"), length.out = 2298))
sample_data <- data.frame(norm_power, unique_ID, Subject, stringsAsFactors = FALSE)
sample_data$intensity[sample(1:nrow(sample_data), nrow(sample_data), FALSE)] <- rep(intensity, c(1149,1149))
sample_data$phoneme[sample(1:nrow(sample_data), nrow(sample_data), FALSE)] <- rep(phoneme, c(574, 574, 575, 575))
我正在尝试执行方差分析并提取样本数据的多个子集的系数-参见上文。当我将数据集过滤到单独的电极(即“ elec1”)中时,我可以获得结果(请参见输出):
sample_data_elec1 <- filter(sample_data, unique_ID == "elec1")
output <-ezANOVA(dv = norm_power, within=.(phoneme, intensity), wid=Subject, data=sample_data_elec1)
但是,当我尝试一次对sample_data的多个“子集”完成类似的操作时,我只会为所有操作生成NA。请参阅下面的“ anova_data”:
get_coefs <-function(d){
output <-ezANOVA(dv = norm_power, within=.(phoneme, intensity), wid=Subject, data=d)
tibble('phoneme_F' = output[["ANOVA"]][[4]][1], "phoneme_pval" = output[["ANOVA"]][[5]][1],
"intensity_F" = output[["ANOVA"]][[4]][2], "intensity_pval" = output[["ANOVA"]][[5]][2],
"interact_F" = output[["ANOVA"]][[4]][3], "interact_pval"= output[["ANOVA"]][[5]][3])
}
anova_data <- sample_data %>%
group_by(unique_ID) %>%
do({get_coefs(.)}) %>%
ungroup()
感谢您的帮助-如果可以进一步说明,请告诉我!
您收到错误,因为您正在运行重复测量方差分析的测试。使用this website中的示例,它具有非常严格的设计:
library(ez) bilingual=read.delim("http://coltekin.net/cagri/R/data/bilingual.txt") # this works ezANOVA(data=bilingual,dv=mlu,wid=.(subj), within=.(language, age), between=gender)
并且您可以看到,对于每个主题,本质上每个元素之前/之后或它们的组合都只有一个测量值。
table(bilingual$age,bilingual$language,bilingual$subj) , , = 1 home.only school firstgrade 1 1 preschool 1 1 secondgrade 1 1 , , = 2 home.only school firstgrade 1 1 preschool 1 1 secondgrade 1 1
返回到示例,您需要使用
within_full
,因为您有多个度量,某些示例数据集可以使用,但是您会看到一些抱怨,这是正确的,因为设计并不意味着对于此测试:
output <-ezANOVA(dv = norm_power,between=.(phoneme, intensity), within_full=.(phoneme, intensity), wid=.(Subject), data=droplevels(sample_data_elec1))
而且您也错过了之间的影响,根据反复测量,我已经在上面为您补充了这些影响。现在,如果您的研究或测试不适合repeated measure anova,则只需将要测试的所有变量放在以下之间:
output <-ezANOVA(dv = norm_power,between=.(phoneme, intensity), wid=.(Subject), data=droplevels(sample_data_elec1))