ANOVA与n维观察？

Question

我开始在R中使用ANOVA进行以下案例：35个国家在120个维度中的得分为[0,1]。通过一定的标准，我先将它们分成8组（1到8组）。这些都存储在data.frame“datos”中，其中包含以下列：维，国家，值，组。 我可以通过以下方式对这120个维度中的每个维度进行ANOVA：

for(i in 1:120){
  L=datos$dimension==rownames(dim1[i,]) #dim1 is a matrix that has the order of the dimensions I use.
  pivot<-datos[L,] #dataset with only the data of dimension i
  aovresult<-aov(value~group, data=pivot) #do ANOVA on dimension i based on groups
  p_value[i]<-summary(aovresult)[[1]][["Pr(>F)"]][[1]] #extract the p_value
}
p_value<-as.data.frame(cbind(rownames(dim1[i,]),p_value))
colnames(p_value)<-c("dimension","p_value")

如您所见，矩阵p_value具有每个维度的ANOVA的p_values。现在，我想知道如何在相同的时间对所有维度进行ANOVA，因为我知道原始分组（1到8）是基于它们的所有维度。

这是一个悬而未决的问题，所以欢迎任何建议。