我正在使用自定义对比度进行单因素方差分析(每个物种)。
[,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
我将强度0.5与5、5与12.5进行比较,依此类推。这些是我正在处理的数据
具有以下结果
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”的全球平均值。在强度1中,我正在将强度0.5的平均值与5进行比较。
如果我理解正确的话,强度1的系数2.2应该是强度水平0.5和5的均值之差的一半。但是,我的手工计算与摘要的计算不匹配。任何人都可以弄清楚我在做什么错吗?
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
能否提供您用来估算的R中的lm()函数。您如何使用对比功能?
—
菲利普(Philippe)2012年
btw
—
飞行
geom_points(position=position_dodge(width=0.75))
将解决您的绘图中的点不与框对齐的方式。
自从我提出问题以来,@ flies一直在介绍
—
RomanLuštrik18年
geom_jitter
,这是所有抖动的geom_point()参数的快捷方式。
我没有注意到那里的抖动。做
—
飞行
geom_jitter(position_dodge)
工作?我一直在使用geom_points(position_jitterdodge)
闪避将点添加到箱形图中。
@flies在
—
RomanLuštrik18年
geom_jitter
此处查看文档。根据以上答案,根据我的经验,我发现没有必要使用箱形图。曾经 如果我有很多点,我会使用小提琴图来显示点的密度,其细节要比箱形图精细得多。绘制许多点或密度不方便时,又发明了箱形图。也许是时候开始考虑删除此(残障人士)可视化了。