1

文本中的数据如下所示：

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

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

4

2 回答 2

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 的度数。

当你c(...)为你的 MPH 变量做你的时，你将它们设置为观察，而不是变量。如果要输入列而不是行，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。

于 2011-12-15T05:58:15.587 回答
1

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

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

于 2011-12-15T04:55:04.960 回答