Özel kontrastlı tek yönlü bir ANOVA (tür başına) yapıyorum.
[,1] [,2] [,3] [,4]
0.5 -1 0 0 0
5 1 -1 0 0
12.5 0 1 -1 0
25 0 0 1 -1
50 0 0 0 1
burada yoğunluk 0.5'i 5'e karşı 5, 12'ye karşı 5 vb. Bunlar üzerinde çalıştığım veriler
aşağıdaki sonuçlarla
Generalized least squares fit by REML
Model: dark ~ intensity
Data: skofijski.diurnal[skofijski.diurnal$species == "niphargus", ]
AIC BIC logLik
63.41333 67.66163 -25.70667
Coefficients:
Value Std.Error t-value p-value
(Intercept) 16.95 0.2140872 79.17334 0.0000
intensity1 2.20 0.4281744 5.13809 0.0001
intensity2 1.40 0.5244044 2.66970 0.0175
intensity3 2.10 0.5244044 4.00454 0.0011
intensity4 1.80 0.4281744 4.20389 0.0008
Correlation:
(Intr) intns1 intns2 intns3
intensity1 0.000
intensity2 0.000 0.612
intensity3 0.000 0.408 0.667
intensity4 0.000 0.250 0.408 0.612
Standardized residuals:
Min Q1 Med Q3 Max
-2.3500484 -0.7833495 0.2611165 0.7833495 1.3055824
Residual standard error: 0.9574271
Degrees of freedom: 20 total; 15 residual
16.95 "niphargus" için küresel ortalamadır. Yoğunluk1'de, yoğunluk 0,5 ile 5 arasındaki araçları karşılaştırıyorum.
Bu hakkı anladıysam, 2.2 yoğunluğu1 katsayısı, 0.5 ve 5 yoğunluk seviyeleri arasındaki farkın yarısı olmalıdır. Ancak, el hesaplamalarım özetle eşleşmiyor. Yanlış yaptığım birileri çipleyebilir mi?
ce1 <- skofijski.diurnal$intensity
levels(ce1) <- c("0.5", "5", "0", "0", "0")
ce1 <- as.factor(as.character(ce1))
tapply(skofijski.diurnal$dark, ce1, mean)
0 0.5 5
14.500 11.875 13.000
diff(tapply(skofijski.diurnal$dark, ce1, mean))/2
0.5 5
-1.3125 0.5625
geom_points(position=position_dodge(width=0.75))
, planınızdaki noktaların kutularla hizalanma biçimini düzeltir.
geom_jitter
, bu da tüm geom_point () parametreleri için bir kısayol.
geom_jitter(position_dodge)
işi? Ben kullanıyorum geom_points(position_jitterdodge)
kaçmak ile kutudiyagramlar için noktaları ekleyin.
geom_jitter
burada . Yukarıdaki cevabımdan bu yana yaşadığım deneyime göre, kutu grafikleri kullanmanın gereksiz olduğunu düşünüyorum. Hiç. Eğer çok puanım varsa, nokta yoğunluğunu kutu grafiklerden çok daha ince ayrıntılarla gösteren keman grafikleri kullanırım. Boxplotlar birçok noktayı çizerken icat edildi veya yoğunlukları uygun değildi. Belki de bu (özürlü) görselleştirmeyi bırakmayı düşünmenin zamanı gelmiştir.