如何在 R 中为 manova 准备数据答案

【问题标题】：How to prepare data in R for manova如何在 R 中为 manova 准备数据
【发布时间】：2012-01-20 20:20:21
【问题描述】：

文本中的数据如下所示：

                       Initial Speed            
Pedal Rotation  10 mph  25 mph  40 mph  55 mph
  5 degrees         0.35    0.19    0.14    0.10
  8 degrees         0.37    0.28    0.19    0.19
  10 degrees        0.42    0.30    0.29    0.23

表格中的数字代表加速度。

我想知道如何准备好这些信息以在 R 中运行 manova 测试。

我已经从列和行创建了变量：

mph10 <- c(0.35, 0.37, 0.42)
mph25 <- c(0.19, 0.28, 0.30)
mph40 <- c(0.14, 0.19, 0.29)
mph55 <- c(0.10, 0.19, 0.23)

和

d5 <- c(0.35, 0.19, 0.14, 0.10)
d8 <- c(0.37, 0.28, 0.19, 0.23)
d10 <- c(0.42, 0.30, 0.29, 0.23)

我在 Youtube 上找到了一些教程，这些教程展示了如何进行 manova，它们都使用 .csv 文件并直接使用 manova，而没有解释数据在文件中的排列方式或如何制作自己的文件。

所以，我不知道下一步该做什么......

【问题讨论】：

标签： r analysis variance manova

【解决方案1】：

通常，您希望数据设置在各自的列中包含多变量响应，在各自的列中包含预测变量。在我看来（不确定你测量的是什么）你的文本文件已经设置好了......但我会删除第一列中重复的“度数”这个词。

我检查了 youtube（r 中的 manova），如果您使用头骨示例（我发现的第二件事）并意识到 manova 命令的 cbind 部分中的所有内容都只是单独的列将 data.frame 组合在一起作为多元响应变量，你应该没问题。

【讨论】：

我真的不明白。我不知道头骨数据与上表有何关系。我上面的表格中的“年份”是什么？我应该 cbind 踏板旋转列，然后我会使用 MPH 行之一作为“年”？
恭喜你在后续问题中得到它 :)

【解决方案2】：

您应该在分析之前将数据转换为数据框，您上面的代码只会创建未链接的列。建议的代码（可能是更有效的方法，这是可行的）：

degrees <- data.frame(Degrees = c("d5","d8","d10"), MPH10=c(0.35, 0.37, 0.32), 
MPH25=c(0.19, 0.28, 0.30), MPH40=c(0.14, 0.19, 0.29), MPH55=c(0.10, 0.19, 0.23))

数据测试：

degrees

  Degrees MPH10 MPH25 MPH40 MPH55
1      d5  0.35  0.19  0.14  0.10
2      d8  0.37  0.28  0.19  0.19
3     d10  0.32  0.30  0.29  0.23

您现在应该能够分析 manova 的度数了。

当您为 MPH 变量执行 c(...) 时，您将它们设置为观察值，而不是变量。如果要输入列而不是行，则使用rbind() 命令。因此命令mph10 <- c(0.35, 0.37, 0.42) 将给mph10 提供以下输出：

[1] 0.35 0.37 0.42

而mph10 <- rbind(0.35, 0.37, 0.42) 给了

     [,1]
[1,] 0.35
[2,] 0.37
[3,] 0.42

HTH。

【讨论】：