我正在尝试在 R 中运行方差分析模型。我有一个包含 3 行和 12 列的数据文件。每行都是解释变量特定水平的数据。单元格 [i,j] 是级别 i 的第 j 个响应。该文件的扩展名为“.dat”。我正在运行以下 R 代码来尝试获取 36 x 2 数据框来运行 anova 模型,而不是 3 x 12 原始数据框:
data <- read.table("usedcar.dat", row.names = 1)
young <- data[1,]
med <- data[2,]
old <- data[3,]
Price <- c(young, med, old)
Age <- as.factor(c(rep(1,12), rep(2,12), rep(3,12)))
data <- cbind(Age, Price)
data <- as.data.frame(data)
但是当我尝试从中获取方差分析模型时,我收到无效列表类型错误:
m1 <- aov(Price ~ Age, data = data)
Error in model.frame.default(formula = Price ~ Age, data = data, drop.unused.levels = TRUE) : invalid type (list) for variable 'Price'
我在这里做错了什么?
这是一个随机矩阵,如果有帮助的话:
replicate(12, rnorm(3))
这是 str(data) 结果:
str(data)
'data.frame': 36 obs. of 2 variables:
$ Age :List of 36
..$ 1 : int 1
..$ 2 : int 1
..$ 3 : int 1
...
..$ 36: int 3
$ Price:List of 36
..$ 1 : int 2300
...
..$ 36: int 2075