如何在R中为manova准备数据

Question

文本中的数据如下所示：

                       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，而没有解释数据是如何在文件中排列的，也没有解释如何制作自己的文件。

所以，我不知道下一步该怎么做……

Answer 1

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

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。

Answer 2

通常，您希望在数据设置中使用多变量响应，每个响应位于各自的列中，而预测变量则位于各自的列中。在我看来(不确定你测量了什么)，你的文本文件已经设置好了……但我会去掉在第一列中重复出现的单词“度”。

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