>Anhead(An)>bartlett.test(Migrate_rates~Ecosystem,data = An)>aov(formula = Migrate_rates~Ecosystem,data = An)>m1summary(m1)>TukeyHSD(m1)Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Migrate_rates ~ Ecosystem, data = An)
$Ecosystem
diff lwr upr p adj
DA-WL-DA-DG -0.24876286 -0.275197906 -0.22232782 0.0000000
DA-WS-DA-DG -0.07545086 -0.101885906 -0.04901582 0.0000000
DG-DA-DA-DG 0.28863293 0.262197888 0.31506797 0.0000000
DG-WL-DA-DG 0.13369423 0.107259185 0.16012927 0.0000000
DG-WS-DA-DG 0.02940491 0.002969868 0.05583995 0.0146976
WL-DA-DA-DG 0.06237151 0.035936467 0.08880655 0.0000000
WL-DG-DA-DG -0.06469621 -0.091131255 -0.03826117 0.0000000
WL-WS-DA-DG 1.08520482 1.058769783 1.11163987 0.0000000
WS-DA-DA-DG 0.59251876 0.566083718 0.61895380 0.0000000
WS-DG-DA-DG 2.44088868 2.414453641 2.46732372 0.0000000
WS-WL-DA-DG 0.23277454 0.206339499 0.25920958 0.0000000
DA-WS-DA-WL 0.17331200 0.146876959 0.19974704 0.0000000
DG-DA-DA-WL 0.53739579 0.510960752 0.56383084 0.0000000
DG-WL-DA-WL 0.38245709 0.356022050 0.40889213 0.0000000
DG-WS-DA-WL 0.27816777 0.251732732 0.30460282 0.0000000
WL-DA-DA-WL 0.31113437 0.284699332 0.33756941 0.0000000
WL-DG-DA-WL 0.18406665 0.157631609 0.21050169 0.0000000
WL-WS-DA-WL 1.33396769 1.307532648 1.36040273 0.0000000
WS-DA-DA-WL 0.84128162 0.814846582 0.86771667 0.0000000
WS-DG-DA-WL 2.68965155 2.663216506 2.71608659 0.0000000
WS-WL-DA-WL 0.48153741 0.455102364 0.50797245 0.0000000
DG-DA-DA-WS 0.36408379 0.337648752 0.39051883 0.0000000
DG-WL-DA-WS 0.20914509 0.182710049 0.23558013 0.0000000
DG-WS-DA-WS 0.10485577 0.078420732 0.13129081 0.0000000
WL-DA-DA-WS 0.13782237 0.111387331 0.16425741 0.0000000
WL-DG-DA-WS 0.01075465 -0.015680391 0.03718969 0.9754251
WL-WS-DA-WS 1.16065569 1.134220647 1.18709073 0.0000000
WS-DA-DA-WS 0.66796962 0.641534582 0.69440466 0.0000000
WS-DG-DA-WS 2.51633955 2.489904505 2.54277459 0.0000000
WS-WL-DA-WS 0.30822540 0.281790363 0.33466045 0.0000000
DG-WL-DG-DA -0.15493870 -0.181373744 -0.12850366 0.0000000
DG-WS-DG-DA -0.25922802 -0.285663062 -0.23279298 0.0000000
WL-DA-DG-DA -0.22626142 -0.252696462 -0.19982638 0.0000000
WL-DG-DG-DA -0.35332914 -0.379764185 -0.32689410 0.0000000
WL-WS-DG-DA 0.79657190 0.770136854 0.82300694 0.0000000
WS-DA-DG-DA 0.30388583 0.277450788 0.33032087 0.0000000
WS-DG-DG-DA 2.15225575 2.125820712 2.17869079 0.0000000
WS-WL-DG-DA -0.05585839 -0.082293430 -0.02942335 0.0000000
DG-WS-DG-WL -0.10428932 -0.130724359 -0.07785428 0.0000000
WL-DA-DG-WL -0.07132272 -0.097757760 -0.04488768 0.0000000
WL-DG-DG-WL -0.19839044 -0.224825482 -0.17195540 0.0000000
WL-WS-DG-WL 0.95151060 0.925075556 0.97794564 0.0000000
WS-DA-DG-WL 0.45882453 0.432389491 0.48525957 0.0000000
WS-DG-DG-WL 2.30719446 2.280759414 2.33362950 0.0000000
WS-WL-DG-WL 0.09908031 0.072645273 0.12551536 0.0000000
WL-DA-DG-WS 0.03296660 0.006531558 0.05940164 0.0026917
WL-DG-DG-WS -0.09410112 -0.120536165 -0.06766608 0.0000000
WL-WS-DG-WS 1.05579992 1.029364874 1.08223496 0.0000000
WS-DA-DG-WS 0.56311385 0.536678808 0.58954889 0.0000000
WS-DG-DG-WS 2.41148377 2.385048732 2.43791881 0.0000000
WS-WL-DG-WS 0.20336963 0.176934590 0.22980467 0.0000000
WL-DG-WL-DA -0.12706772 -0.153502764 -0.10063268 0.0000000
WL-WS-WL-DA 1.02283332 0.996398275 1.04926836 0.0000000
WS-DA-WL-DA 0.53014725 0.503712209 0.55658229 0.0000000
WS-DG-WL-DA 2.37851717 2.352082133 2.40495222 0.0000000
WS-WL-WL-DA 0.17040303 0.143967991 0.19683807 0.0000000
WL-WS-WL-DG 1.14990104 1.123465997 1.17633608 0.0000000
WS-DA-WL-DG 0.65721497 0.630779931 0.68365001 0.0000000
WS-DG-WL-DG 2.50558490 2.479149855 2.53201994 0.0000000
WS-WL-WL-DG 0.29747075 0.271035713 0.32390580 0.0000000
WS-DA-WL-WS -0.49268607 -0.519121107 -0.46625102 0.0000000
WS-DG-WL-WS 1.35568386 1.329248817 1.38211890 0.0000000
WS-WL-WL-WS -0.85243028 -0.878865325 -0.82599524 0.0000000
WS-DG-WS-DA 1.84836992 1.821934882 1.87480497 0.0000000
WS-WL-WS-DA -0.35974422 -0.386179260 -0.33330918 0.0000000
WS-WL-WS-DG -2.20811414 -2.234549183 -2.18167910 0.0000000
其中除了WL-DG和DA-WS之间差异不显著外,其余都显著。
R语言绘制密度图(Density)>library(ggplot2)>library(ggpubr)>Denhead(Den)>spsp
显著性检验
>bartlett.test(Migrate.rate~Ecosystem,data = Den)>m1summary(m1)>TukeyHSD(m1)Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Migrate.rate ~ Ecosystem, data = Den)
$Ecosystem
diff lwr upr p adj
D-W-D-D -0.1845951 -0.1989791 -0.1702111 0
W-D-D-D 0.6134542 0.5990702 0.6278382 0
W-W-D-D 0.5146732 0.5002892 0.5290572 0
W-D-D-W 0.7980493 0.7836653 0.8124333 0
W-W-D-W 0.6992683 0.6848843 0.7136523 0
W-W-W-D -0.0987810 -0.1131650 -0.0843970 0