0

我是 R 的新手，而且（我承认！）对统计数据不太好，所以如果我的问题太琐碎，我很抱歉，但我真的很感激有关此事的一些提示。

对于我们拥有的 2 个不同的种植系统（农林业和农业）中的每一个，我有 9 个点（图）土壤湿度测量值超过 2 个月（每周测量值）。我们还有最近的树（大于 5cm DBH）和每个地块中的精确测量点之间的距离（以米为单位）（在农林业中为 4.2 至 12m，在农业中为 50m）。因此，对于 18 个地块（农林业中的 9 个和农业中的 9 个）中的每一个，我都有一个湿度 (y) 随时间 (x) 的曲线（表现相似但由于天气波动而变化）。我需要知道的是：

  1. 这些测量点之间的湿度随时间的变化是否取决于（或受其影响）树木的距离？意思是，如果树木更靠近测量点（在种植园中间？

  2. 这些曲线（湿度 x 时间）是否彼此显着不同？

我首先考虑将农林业系统的每 3 个树木测量点（距树木的较小距离、中等距离和较高距离）和来自农林业的所有 9 个点进行分组，并将它们用作复制，因为它们的行为更相似。然而这让我有点困惑。

所以......我一直在考虑使用ez包装中的重复测量方差分析。所以在这种情况下，我有：

str(SanPedro)

data.frame':   450 obs. of  6 variables:
 Parcel  : Factor w/ 2 levels "Forest","Agriculture": 1 1 1 1 1 1 1 1 1 1 ...
 Distance: Factor w/ 4 levels "A","B","C","D": 1 1 1 1 1 1 1 1 1 1 ...
 Plot    : num  1 1 1 1 1 1 1 1 1 1 ...
 Date    : Date, format: "0011-07-20" "0011-07-24" ...
 Humidity: num  0.217 0.205 0.199 0.2 0.192 0.181 0.184 0.18 0.179 0.178 ...
 Number  : num  1 2 3 4 5 6 7 8 9 10 ..

当我尝试运行ezANOVAas

ezANOVA(data=SanPedro, dv=Humidity, wid=Number, within=Parcel, between=Plot, type=1, return_aov=TRUE)

我懂了：

 Warning: Converting "Number" to factor for ANOVA.
 Warning: "Plot" will be treated as numeric.
 Error in ezANOVA_main(data = data, dv = dv, wid = wid, within = within,  : 
  One or more cells is missing data. Try using ezDesign() to check your data.

如果我检查ezDesign(SanPedro)，我得到：

ezDesign(SanPedro)

Error in as.list(c(x, y, row, col)) : 
 argument "x" is missing, with no default

最后，我并不真正理解数据的问题，我什至不确定这ezANOVA是否真的适合我的案例......我非常感谢任何关于此事的提示和想法！感谢一个loooot！=)

4

0 回答 0