對R中的所有變量重復測量方差分析

Question

我執行重復測量方差分析，我有任何變量（4個或更多），當我這樣做時，我必須分別對每個變量進行方差分析。 這對我來說rep.anova方便，例如，這里有兩個變量rep.anova的代碼-tember_body和SAD（我有8個度量）

library("car")
tabl20=read.csv("path to dataset",sep=";",dec=",")
ageLevels <- c(1, 2, 3,4,5,6,7,8)
ageFactor <- as.factor(ageLevels)
ageFrame <- data.frame(ageFactor)
ageBind <- cbind(tabl20$temper_body_v0,
                 tabl20$temper_body_v1,
                 tabl20$temper_body_v2,
                 tabl20$temper_body_v3,
                 tabl20$temper_body_v4,
                 tabl20$temper_body_v5,
                 tabl20$temper_body_v6,
                 tabl20$temper_body_v7)
ageModel <- lm(ageBind ~ 1)
analysis <- Anova(ageModel, idata = ageFrame, idesign = ~ageFactor)
summary(analysis)
######
tabl20=read.csv("path to data",sep=";",dec=",")
ageLevels <- c(1, 2, 3,4,5,6,7,8)
ageFactor <- as.factor(ageLevels)
ageFrame <- data.frame(ageFactor)
ageBind <- cbind(tabl20$SAD_v0,
                 tabl20$SAD_v1,
                 tabl20$SAD_v2,
                 tabl20$SAD_v3,
                 tabl20$SAD_v4,
                 tabl20$SAD_v5,
                 tabl20$SAD_v6,
                 tabl20$SAD_v7)

ageModel <- lm(ageBind ~ 1)
analysis <- Anova(ageModel, idata = ageFrame, idesign = ~ageFactor)
summary(analysis)

如何一次對所有變量執行重復測量方差分析，而不是分別對每個變量進行分析。

Answer 1

我找到答案了，也許對某人有用

library(car) 
tabl20=read.csv("c:/path/,sep=";",dec=",") 
ageLevels <- c(1, 2, 3,4,5,6,7,8) 
ageFactor <- as.factor(ageLevels) 
ageFrame <- data.frame(ageFactor) 

tabl20 <- tabl20[, order(names(tabl20))] 

measures <- function(data = tabl20, n = 4) { 
  list <- list(n) 
  for(i in 0:3) {list[[i+1]] <- as.matrix(cbind(data[, 8*i + 1:8])) 
  } 
  list 
} 

measures_list <- measures() 

models <- lapply( 
  measures_list, function(x) { 
    ageModel <- lm(x ~ 1) 
    Anova(ageModel, idata = ageFrame, idesign = ~ageFactor) 
  } )

但我不知道如何對所有變量進行繪圖

如果我這樣做情節（模型）

我得到這個錯誤

Error in xy.coords(x, y, xlabel, ylabel, log) : 
  'x' is a list, but does not have components 'x' and 'y'

也許有人告訴我，如何為所有變量創建圖