重複mesuare Anova所有變量在R

我執行重複mesuare Anova，我有任何變量（4或更多），當我這樣做時，我必須分別對每個變量進行anova。它的不方便對我來說，例如這裏的代碼rep.anova 2個變量 - tember_body和SAD（我有8項措施）重複mesuare Anova所有變量在R

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)

如何在執行一次rep.measure ANOVA所有變量，所有的，而不是做分別進行每個變量分析。

來源

2016-01-22 psysky

看起來，這將是構建一個動態公式以在任何想要的變量組合上運行* ANOVA *的問題。我建議你發佈一些示例數據或使用'mtcars'數據重現你的例子，這樣可以更容易地貢獻出潛在的解決方案。 – Konrad

謝謝，我如何附加我的數據集，它很小？ – psysky

看看[這個討論]（http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example），它充滿了關於如何產生一個可重複問題的例子。 – Konrad

我找到答案，也許是有人

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'

也許任何人promt我，如何創建plo t對於所有變量

來源

2016-01-23 15:12:07 psysky

重複mesuare Anova所有變量在R

回答

相關問題