我比较了在不同栖息地观察到的不同行为的比例。我通过绘制 R 中使用的数据来可视化这些差异。我还运行了方差分析和 Tukey 的 HSD 测试。如果我正确地解释了 Tukey 输出,它会一直告诉我在图中看起来非常不同的组比较实际上根本没有区别!?!?如果有人可以解释我的误解或可以建议一种不同的方法来测试栖息地之间行为的显着差异,我将不胜感激。
一些 Tukey 输出:
| comparison | diff | lwr | upr | p adj|
|---------------------|---------|------------|------|--------|
|Urban-Edge | 31.900000 | 9.670532 |54.12947 |0.0011390|
|Edge-Orchard | 2.190909 |-19.527468| 23.90929 |0.9996715|
在分析中绘制和使用的一些数据:
|habitat |behavior | frequency|
|--------|---------|----------|
|Urban |Flying| 0.0119146119|
|Urban |Flying| 0.0143968228|
|Urban |Flying| 0.0143968228|
|Urban |Flying| 0.0054608638|
|Urban |Flying| 0.0049644216|
|Urban |Flying| 0.0047989409|
|Urban |Flying| 0.0038060566|
|Urban |Flying| 0.0049644216|
|Urban |Flying| 0.0105907662|
|Urban |Flying| 0.0157206685|
|Edge |Flying| 0.0013238458|
|Edge |Flying| 0.0039715373|
|Edge |Flying| 0.0021512494|
|Edge |Flying| 0.0082740361|
|Edge |Flying| 0.0031441337|
|Edge |Flying| 0.0033096144|
|Edge |Flying| 0.0028131723|
|Edge |Flying| 0.0044679795|
|Edge |Flying| 0.0049644216|
|Edge |Flying| 0.0038060566|
|Orchard |Flying| 0.0026476915|
|Orchard| Flying| 0.0021512494|
|Orchard |Flying| 0.0026476915|
|Orchard| Flying| 0.0009928843|
|Orchard |Flying| 0.0033096144|
|Orchard |Flying| 0.0077775939|
|Orchard| Flying| 0.0043024988|
|Orchard |Flying| 0.0034750952|
|Orchard |Flying| 0.0041370180|
|Orchard |Flying| 0.0033096144|
|Orchard |Flying| 0.0033096144|
Simple ANOVA and Tukey code:
oneFly <- aov(Freq ~ Habitat, data = Fly)
summary(oneFly)
tukey.oneFly <-TukeyHSD(oneFly)
tukey.oneFly