重复测量方差分析或时间序列分析?

时间:2013-11-14 17:59:52

标签: r time-series anova

我是R的新人,(我承认!)统计数据不太好,所以如果我的问题太琐碎,我很抱歉,但我真的很感激有关这个问题的一些提示。

对于我们拥有的2个不同种植系统(农林业和农业)的2个月(每周测量),我有9点(图)土壤湿度测量值。我们还有最近的树(大于5厘米胸径)和每个地块中的精确测量点之间的距离(在农林中为4.2到12米之间,在农业中为50米)。因此,对于18个地块中的每一个(农业中9个和农业中的9个),我具有随时间变化的湿度(y)(x)(其表现类似但由于天气波动而变化)。我需要知道的是:

  1. 测量点之间的湿度随时间的变化是否取决于树木的距离(或受其影响)?这意味着,如果树木更接近测量点(位于种植园中间,那么树木是否能容纳更多的水或从土壤中取出更多的水?

  2. 这些曲线(湿度x时间)是否显着不同?

  3. 我首先考虑的是对农林业系统的每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 ..
    

    当我尝试将ezANOVA作为

    运行时
    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! =)

0 个答案:

没有答案