ANOSIM, PERMANOVA, SIMPER, Indicator Value
Pablo E. GutiƩrrez-Fonseca
31 October, 2022
Step 1
Load libraries.
Cargar las librerias que necesitas.
libraries <- c("vegan", "ggplot2", "dplyr", "indicspecies", "labdsv", "MASS",
"labdsv", "cluster")
lapply(libraries, require, character.only = TRUE)## [[1]]
## [1] TRUE
##
## [[2]]
## [1] TRUE
##
## [[3]]
## [1] TRUE
##
## [[4]]
## [1] TRUE
##
## [[5]]
## [1] TRUE
##
## [[6]]
## [1] TRUE
##
## [[7]]
## [1] TRUE
##
## [[8]]
## [1] TRUEStep 2
Data
Cargar los datos.
Step 3
Seleccionar la matrix de especies.
laselva_sp <- dplyr::select(laselva_full, Baetid:Collemb)
str(laselva_sp)
ncol(laselva_sp)
nrow(laselva_sp)ANOSIM, PERMANOVA, SIMPER, Indicator Value
Analysis of Similitude
laselva.dist <- vegdist(laselva_sp, method="bray")
laselva.anosim <- anosim(laselva.dist, laselva_full$month)
summary(laselva.anosim)##
## Call:
## anosim(x = laselva.dist, grouping = laselva_full$month)
## Dissimilarity: bray
##
## ANOSIM statistic R: 0.6113
## Significance: 0.001
##
## Permutation: free
## Number of permutations: 999
##
## Upper quantiles of permutations (null model):
## 90% 95% 97.5% 99%
## 0.0372 0.0498 0.0607 0.0724
##
## Dissimilarity ranks between and within classes:
## 0% 25% 50% 75% 100% N
## Between 6 2364.25 3728.75 4839.75 5886 3974
## four 34 1808.50 3007.75 4229.50 5739 276
## one 5 23.00 50.00 98.00 278 6
## three 1 719.50 1473.00 2479.50 5254 1035
## two 2 362.25 1011.00 1779.00 4392 595plot(laselva.anosim)## Warning in (function (z, notch = FALSE, width = NULL, varwidth = FALSE, : some
## notches went outside hinges ('box'): maybe set notch=FALSE
dis_well <- betadisper(laselva.dist, laselva_full$month)
permutest(dis_well)##
## Permutation test for homogeneity of multivariate dispersions
## Permutation: free
## Number of permutations: 999
##
## Response: Distances
## Df Sum Sq Mean Sq F N.Perm Pr(>F)
## Groups 3 0.39374 0.131248 12.803 999 0.001 ***
## Residuals 105 1.07640 0.010251
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1well.HSD <- TukeyHSD(dis_well)
well.HSD ## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = distances ~ group, data = df)
##
## $group
## diff lwr upr p adj
## one-four -0.26670779 -0.409460808 -0.12395477 0.0000226
## three-four -0.09939867 -0.165957611 -0.03283974 0.0009687
## two-four -0.13582474 -0.205877984 -0.06577149 0.0000105
## three-one 0.16730912 0.029519024 0.30509921 0.0106105
## two-one 0.13088305 -0.008628506 0.27039461 0.0742302
## two-three -0.03642607 -0.095714703 0.02286257 0.3808697par(mfrow=c(1,2))
plot(well.HSD, las=1)
plot(dis_well, las=1)
PERMANOVA
laselva_adonis <- adonis2(laselva_sp ~ month, data = laselva_full, permutations = 999, method="bray")
print(laselva_adonis)## Permutation test for adonis under reduced model
## Terms added sequentially (first to last)
## Permutation: free
## Number of permutations: 999
##
## adonis2(formula = laselva_sp ~ month, data = laselva_full, permutations = 999, method = "bray")
## Df SumOfSqs R2 F Pr(>F)
## month 3 7.2292 0.43757 27.23 0.001 ***
## Residual 105 9.2922 0.56243
## Total 108 16.5213 1.00000
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Paso 5. Comparasion entre grupos.
library(devtools)## Loading required package: usethis##
## Attaching package: 'devtools'## The following object is masked from 'package:permute':
##
## checkinstall_github("pmartinezarbizu/pairwiseAdonis/pairwiseAdonis") ## Skipping install of 'pairwiseAdonis' from a github remote, the SHA1 (ece560d2) has not changed since last install.
## Use `force = TRUE` to force installationlibrary(pairwiseAdonis)
pair.mod <- pairwise.adonis(laselva_sp, factors= laselva_full$month)
pair.mod## pairs Df SumsOfSqs F.Model R2 p.value p.adjusted sig
## 1 one vs two 1 0.4709231 8.484656 0.1865389 0.001 0.006 *
## 2 one vs three 1 1.2988558 17.220251 0.2640323 0.001 0.006 *
## 3 one vs four 1 1.9675994 13.803568 0.3467922 0.001 0.006 *
## 4 two vs three 1 2.1434846 30.313948 0.2773109 0.001 0.006 *
## 5 two vs four 1 5.2546510 52.808573 0.4809148 0.001 0.006 *
## 6 three vs four 1 2.3704514 22.268346 0.2466905 0.001 0.006 *
Paso 6. Vamos a ver quien cuales especies contribuyen a la disimilaridad entre esos grupos.
################ SIMPER ######################
laselva_simper <- with(laselva_full, simper(laselva_sp, laselva_full$month), permutations = 999)
summary(laselva_simper, ordered = TRUE,
digits = max(3,getOption("digits") - 3))##
## Contrast: one_two
##
## average sd ratio ava avb cumsum p
## Tanytar 2.575e-01 7.173e-02 3.590e+00 2.101e+04 5.843e+03 0.585 0.002 **
## Tanypod 5.135e-02 2.953e-02 1.739e+00 8.070e+03 5.833e+03 0.702 1.000
## Chirono 3.847e-02 3.027e-02 1.271e+00 4.768e+03 5.015e+03 0.789 1.000
## Copepod 1.549e-02 1.471e-02 1.053e+00 9.040e+02 3.700e+01 0.824 0.002 **
## Ceratop 1.230e-02 1.411e-02 8.720e-01 1.087e+03 1.424e+03 0.852 1.000
## Oligoch 6.400e-03 4.550e-03 1.406e+00 4.850e+02 1.590e+02 0.867 0.653
## Macrone 6.340e-03 5.830e-03 1.087e+00 3.800e+02 1.500e+02 0.881 0.566
## Ortochl 5.010e-03 7.910e-03 6.340e-01 1.440e+02 3.260e+02 0.892 0.910
## Hexatom 4.970e-03 7.390e-03 6.730e-01 1.050e+02 2.810e+02 0.904 0.854
## Polycha 4.790e-03 4.450e-03 1.077e+00 2.360e+02 1.440e+02 0.915 0.682
## Farrode 3.860e-03 7.140e-03 5.410e-01 1.050e+02 2.690e+02 0.923 0.998
## Nematom 2.510e-03 3.130e-03 8.010e-01 1.310e+02 6.300e+01 0.929 0.510
## Epigomp 2.400e-03 2.100e-03 1.139e+00 1.830e+02 1.840e+02 0.935 1.000
## Xestoch 2.130e-03 1.940e-03 1.100e+00 1.310e+02 1.260e+02 0.939 0.996
## Phylloi 1.970e-03 3.030e-03 6.510e-01 2.600e+01 1.170e+02 0.944 0.940
## Caenis 1.850e-03 1.220e-03 1.521e+00 1.310e+02 4.500e+01 0.948 0.481
## Triplec 1.580e-03 1.540e-03 1.028e+00 7.900e+01 5.200e+01 0.952 0.517
## Hexacyl 1.550e-03 1.440e-03 1.076e+00 7.900e+01 6.300e+01 0.955 0.811
## Nectops 1.540e-03 1.510e-03 1.020e+00 9.200e+01 3.400e+01 0.959 0.235
## Banyall 1.230e-03 2.640e-03 4.660e-01 0.000e+00 8.100e+01 0.962 0.914
## Platyhe 1.150e-03 1.120e-03 1.025e+00 6.600e+01 1.000e+01 0.964 0.267
## Tipulid 1.070e-03 1.770e-03 6.020e-01 0.000e+00 6.100e+01 0.967 1.000
## Ptiloda 9.300e-04 1.110e-03 8.330e-01 5.200e+01 1.300e+01 0.969 0.327
## Acari 9.200e-04 1.500e-03 6.170e-01 1.300e+01 4.900e+01 0.971 0.954
## Oecetis 8.200e-04 1.040e-03 7.840e-01 2.600e+01 3.700e+01 0.973 0.729
## Argia 8.000e-04 1.090e-03 7.370e-01 2.600e+01 3.100e+01 0.974 0.324
## America 7.500e-04 1.720e-03 4.350e-01 0.000e+00 4.800e+01 0.976 0.854
## Chrisop 6.600e-04 9.000e-04 7.280e-01 2.600e+01 1.800e+01 0.978 0.684
## Perisso 6.400e-04 8.100e-04 7.840e-01 2.600e+01 2.400e+01 0.979 0.788
## Calliba 5.800e-04 1.050e-03 5.560e-01 1.300e+01 2.800e+01 0.980 0.768
## Neoelmi 5.600e-04 1.370e-03 4.100e-01 0.000e+00 3.600e+01 0.982 0.882
## Neargyr 5.600e-04 2.090e-03 2.650e-01 0.000e+00 3.000e+01 0.983 0.474
## Coenag 5.400e-04 5.300e-04 1.010e+00 2.600e+01 1.300e+01 0.984 0.249
## Chimarr 4.400e-04 1.900e-03 2.340e-01 0.000e+00 3.400e+01 0.985 0.665
## Libellu 4.400e-04 6.800e-04 6.500e-01 1.300e+01 1.900e+01 0.986 0.584
## Heterag 4.000e-04 8.800e-04 4.540e-01 0.000e+00 2.400e+01 0.987 0.834
## Simuli 3.800e-04 5.900e-04 6.470e-01 1.300e+01 1.500e+01 0.988 0.740
## Collemb 3.600e-04 4.900e-04 7.430e-01 1.300e+01 1.500e+01 0.989 0.775
## Anacro 3.600e-04 8.400e-04 4.280e-01 0.000e+00 2.400e+01 0.990 0.794
## Dixella 3.400e-04 5.400e-04 6.320e-01 1.300e+01 1.000e+01 0.990 0.204
## Coleopt 3.000e-04 7.500e-04 3.990e-01 0.000e+00 1.900e+01 0.991 0.586
## Baetid 2.700e-04 4.300e-04 6.340e-01 1.300e+01 4.000e+00 0.992 0.298
## Pelonom 2.400e-04 4.200e-04 5.700e-01 1.300e+01 0.000e+00 0.992 0.031 *
## Anophel 2.400e-04 4.200e-04 5.700e-01 1.300e+01 0.000e+00 0.993 0.061 .
## Perithe 2.400e-04 4.200e-04 5.700e-01 1.300e+01 0.000e+00 0.993 0.068 .
## Autroti 2.100e-04 6.000e-04 3.520e-01 0.000e+00 1.300e+01 0.994 0.712
## Macrelm 2.100e-04 1.080e-03 1.920e-01 0.000e+00 1.200e+01 0.994 0.428
## Elga 1.900e-04 4.500e-04 4.280e-01 0.000e+00 1.200e+01 0.995 0.325
## Smicrid 1.900e-04 5.000e-04 3.730e-01 0.000e+00 1.200e+01 0.995 0.844
## Clognia 1.800e-04 9.300e-04 1.930e-01 0.000e+00 1.000e+01 0.996 0.624
## Scirtid 1.700e-04 4.400e-04 3.800e-01 0.000e+00 1.200e+01 0.996 0.396
## Lutroch 1.700e-04 4.500e-04 3.780e-01 0.000e+00 1.000e+01 0.996 0.924
## Forcipo 1.600e-04 8.000e-04 2.000e-01 0.000e+00 9.000e+00 0.997 0.835
## Gyrete 1.500e-04 3.300e-04 4.470e-01 0.000e+00 9.000e+00 0.997 0.756
## Phyllog 1.300e-04 4.700e-04 2.880e-01 0.000e+00 7.000e+00 0.997 0.953
## Hyrudi 1.300e-04 4.400e-04 2.920e-01 0.000e+00 7.000e+00 0.998 0.578
## Nehalen 1.000e-04 4.100e-04 2.440e-01 0.000e+00 6.000e+00 0.998 0.616
## Polycen 1.000e-04 4.400e-04 2.220e-01 0.000e+00 7.000e+00 0.998 0.821
## Xenelmi 1.000e-04 3.500e-04 2.840e-01 0.000e+00 6.000e+00 0.998 0.996
## Hemerodr 8.000e-05 3.400e-04 2.260e-01 0.000e+00 6.000e+00 0.998 0.907
## Marilia 7.000e-05 2.200e-04 3.010e-01 0.000e+00 4.000e+00 0.999 0.696
## Palaemn 5.000e-05 2.100e-04 2.440e-01 0.000e+00 3.000e+00 0.999 0.286
## Neotric 5.000e-05 2.100e-04 2.440e-01 0.000e+00 3.000e+00 0.999 0.292
## Limonia 5.000e-05 2.900e-04 1.700e-01 0.000e+00 3.000e+00 0.999 0.555
## Ambrysu 5.000e-05 1.900e-04 2.440e-01 0.000e+00 3.000e+00 0.999 0.911
## Dolichop 4.000e-05 1.800e-04 2.410e-01 0.000e+00 3.000e+00 0.999 0.548
## Heterel 4.000e-05 1.700e-04 2.440e-01 0.000e+00 3.000e+00 0.999 0.215
## Ulmerit 4.000e-05 1.700e-04 2.450e-01 0.000e+00 3.000e+00 0.999 0.192
## Ephyd2 4.000e-05 1.500e-04 2.440e-01 0.000e+00 3.000e+00 0.999 0.167
## Sphaeri 3.000e-05 1.700e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.616
## Pseudol 3.000e-05 1.700e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.822
## Sialis 3.000e-05 1.600e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.363
## Leptohy 3.000e-05 1.500e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.603
## Protopti 3.000e-05 1.500e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.364
## Neocord 3.000e-05 1.500e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.374
## Desmogo 2.000e-05 1.500e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.362
## Cyrnell 2.000e-05 1.300e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.380
## Neoplas 2.000e-05 1.200e-04 1.700e-01 0.000e+00 1.000e+00 1.000 0.886
## Guajir 2.000e-05 1.200e-04 1.710e-01 0.000e+00 1.000e+00 1.000 0.383
## Tricory 2.000e-05 1.000e-04 1.710e-01 0.000e+00 1.000e+00 1.000 0.337
## Hetaeri 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.339
## Macrothe 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.360
## Lestida 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.343
## Agriogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.358
## Rhagove 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.360
## Paracym 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.354
## Leptne 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.346
## Macroste 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.357
## Hydropt 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.358
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Contrast: one_three
##
## average sd ratio ava avb cumsum p
## Tanytar 3.833e-01 6.589e-02 5.816e+00 2.101e+04 2.496e+03 0.582 0.001 ***
## Tanypod 1.156e-01 2.782e-02 4.156e+00 8.070e+03 2.521e+03 0.757 0.391
## Chirono 5.770e-02 3.031e-02 1.904e+00 4.768e+03 2.073e+03 0.845 0.999
## Copepod 1.930e-02 1.820e-02 1.058e+00 9.040e+02 6.800e+00 0.874 0.001 ***
## Ceratop 1.100e-02 7.800e-03 1.416e+00 1.087e+03 9.398e+02 0.891 1.000
## Macrone 1.030e-02 1.772e-02 5.830e-01 3.800e+02 2.301e+02 0.906 0.303
## Oligoch 9.600e-03 5.980e-03 1.605e+00 4.850e+02 1.299e+02 0.921 0.272
## Polycha 5.700e-03 6.050e-03 9.470e-01 2.360e+02 1.253e+02 0.930 0.546
## Hexatom 4.200e-03 6.490e-03 6.500e-01 1.050e+02 1.504e+02 0.936 0.946
## Farrode 3.900e-03 5.220e-03 7.470e-01 1.050e+02 1.834e+02 0.942 0.995
## Ortochl 3.800e-03 3.980e-03 9.520e-01 1.440e+02 1.185e+02 0.948 0.989
## Tipulid 3.000e-03 5.430e-03 5.590e-01 0.000e+00 1.469e+02 0.952 0.992
## Nematom 2.700e-03 2.180e-03 1.230e+00 1.310e+02 3.080e+01 0.956 0.457
## Caenis 2.500e-03 1.480e-03 1.720e+00 1.310e+02 1.820e+01 0.960 0.236
## Epigomp 2.500e-03 1.860e-03 1.360e+00 1.830e+02 1.002e+02 0.964 1.000
## Xestoch 2.200e-03 1.740e-03 1.287e+00 1.310e+02 6.610e+01 0.968 0.991
## Nectops 1.800e-03 1.830e-03 9.680e-01 9.200e+01 1.030e+01 0.970 0.137
## Triplec 1.600e-03 1.220e-03 1.279e+00 7.900e+01 1.820e+01 0.973 0.509
## Hexacyl 1.500e-03 1.770e-03 8.710e-01 7.900e+01 1.940e+01 0.975 0.825
## Platyhe 1.400e-03 1.360e-03 1.034e+00 6.600e+01 1.370e+01 0.977 0.164
## Phylloi 1.400e-03 2.680e-03 5.180e-01 2.600e+01 5.580e+01 0.979 0.996
## Ptiloda 1.100e-03 1.370e-03 8.300e-01 5.200e+01 1.370e+01 0.981 0.203
## Acari 8.000e-04 1.210e-03 6.760e-01 1.300e+01 3.300e+01 0.982 0.986
## Simuli 8.000e-04 2.600e-03 3.130e-01 1.300e+01 2.850e+01 0.983 0.636
## Banyall 7.000e-04 1.370e-03 5.030e-01 0.000e+00 3.530e+01 0.984 0.999
## Chrisop 7.000e-04 1.010e-03 6.740e-01 2.600e+01 9.100e+00 0.986 0.671
## Perisso 7.000e-04 7.400e-04 8.850e-01 2.600e+01 1.590e+01 0.986 0.766
## Argia 6.000e-04 1.030e-03 6.160e-01 2.600e+01 3.400e+00 0.987 0.460
## Coenag 6.000e-04 5.700e-04 9.880e-01 2.600e+01 2.300e+00 0.988 0.228
## Chimarr 6.000e-04 3.060e-03 1.830e-01 0.000e+00 3.190e+01 0.989 0.638
## Oecetis 5.000e-04 7.900e-04 6.790e-01 2.600e+01 5.700e+00 0.990 0.934
## America 5.000e-04 1.670e-03 2.990e-01 0.000e+00 2.510e+01 0.991 0.955
## Neoelmi 5.000e-04 1.530e-03 3.000e-01 0.000e+00 2.390e+01 0.991 0.945
## Calliba 4.000e-04 6.300e-04 6.920e-01 1.300e+01 1.140e+01 0.992 0.923
## Libellu 4.000e-04 5.400e-04 6.710e-01 1.300e+01 6.800e+00 0.993 0.722
## Baetid 3.000e-04 6.100e-04 5.700e-01 1.300e+01 3.400e+00 0.993 0.210
## Collemb 3.000e-04 5.800e-04 5.900e-01 1.300e+01 5.700e+00 0.994 0.804
## Pelonom 3.000e-04 5.200e-04 5.730e-01 1.300e+01 0.000e+00 0.994 0.001 ***
## Anophel 3.000e-04 5.200e-04 5.730e-01 1.300e+01 0.000e+00 0.995 0.064 .
## Perithe 3.000e-04 5.200e-04 5.730e-01 1.300e+01 0.000e+00 0.995 0.064 .
## Xenelmi 3.000e-04 6.000e-04 4.860e-01 0.000e+00 1.370e+01 0.996 0.968
## Dixella 3.000e-04 4.700e-04 5.900e-01 1.300e+01 1.100e+00 0.996 0.240
## Heterag 3.000e-04 9.900e-04 2.750e-01 0.000e+00 1.480e+01 0.996 0.950
## Phyllog 2.000e-04 7.400e-04 3.280e-01 0.000e+00 1.140e+01 0.997 0.862
## Lutroch 2.000e-04 8.100e-04 2.740e-01 0.000e+00 1.140e+01 0.997 0.918
## Hemerodr 2.000e-04 6.700e-04 2.960e-01 0.000e+00 1.030e+01 0.997 0.727
## Anacro 2.000e-04 7.000e-04 2.690e-01 0.000e+00 1.030e+01 0.998 0.962
## Polycen 2.000e-04 5.800e-04 2.760e-01 0.000e+00 8.000e+00 0.998 0.740
## Smicrid 1.000e-04 6.900e-04 2.020e-01 0.000e+00 6.800e+00 0.998 0.945
## Nehalen 1.000e-04 7.600e-04 1.750e-01 0.000e+00 6.800e+00 0.998 0.620
## Autroti 1.000e-04 4.300e-04 2.800e-01 0.000e+00 5.700e+00 0.998 0.921
## Forcipo 1.000e-04 4.000e-04 2.890e-01 0.000e+00 5.700e+00 0.999 0.920
## Neoplas 1.000e-04 4.100e-04 2.430e-01 0.000e+00 4.600e+00 0.999 0.759
## Pseudol 1.000e-04 4.000e-04 2.490e-01 0.000e+00 4.600e+00 0.999 0.687
## Coleopt 1.000e-04 3.000e-04 3.060e-01 0.000e+00 4.600e+00 0.999 0.977
## Ambrysu 1.000e-04 3.000e-04 3.050e-01 0.000e+00 4.600e+00 0.999 0.841
## Gyrete 1.000e-04 2.600e-04 2.610e-01 0.000e+00 3.400e+00 0.999 0.974
## Neargyr 0.000e+00 3.200e-04 1.480e-01 0.000e+00 2.300e+00 0.999 0.892
## Clognia 0.000e+00 2.200e-04 2.120e-01 0.000e+00 2.300e+00 0.999 0.951
## Macrelm 0.000e+00 2.100e-04 2.120e-01 0.000e+00 2.300e+00 1.000 0.813
## Macrothe 0.000e+00 2.800e-04 1.480e-01 0.000e+00 2.300e+00 1.000 0.439
## Dolichop 0.000e+00 2.600e-04 1.480e-01 0.000e+00 2.300e+00 1.000 0.604
## Sphaeri 0.000e+00 1.700e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.714
## Agriogo 0.000e+00 1.700e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.449
## Hydropt 0.000e+00 1.700e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.449
## Hetaeri 0.000e+00 1.600e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.460
## Marilia 0.000e+00 1.600e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.915
## Lestida 0.000e+00 1.600e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.470
## Paracym 0.000e+00 1.600e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.468
## Macroste 0.000e+00 1.600e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.451
## Scirtid 0.000e+00 1.500e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.983
## Limonia 0.000e+00 1.500e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.692
## Rhagove 0.000e+00 1.400e-04 1.480e-01 0.000e+00 1.100e+00 1.000 0.439
## Guajir 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.485
## Ulmerit 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.740
## Tricory 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.472
## Leptohy 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.712
## Elga 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.981
## Neocord 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.482
## Desmogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.427
## Palaemn 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.728
## Heterel 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.748
## Sialis 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.470
## Leptne 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.445
## Neotric 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.717
## Cyrnell 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.453
## Protopti 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.450
## Ephyd2 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.716
## Hyrudi 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.913
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Contrast: one_four
##
## average sd ratio ava avb cumsum p
## Tanytar 4.830e-01 5.269e-02 9.167e+00 2.101e+04 6.266e+02 0.564 0.001 ***
## Tanypod 1.709e-01 2.911e-02 5.872e+00 8.070e+03 9.301e+02 0.764 0.010 **
## Chirono 9.420e-02 2.833e-02 3.325e+00 4.768e+03 7.729e+02 0.874 0.618
## Copepod 2.220e-02 2.098e-02 1.058e+00 9.040e+02 2.200e+00 0.900 0.001 ***
## Ceratop 1.630e-02 9.240e-03 1.760e+00 1.087e+03 4.760e+02 0.919 0.997
## Oligoch 1.050e-02 6.590e-03 1.599e+00 4.850e+02 6.550e+01 0.931 0.234
## Macrone 8.600e-03 7.990e-03 1.072e+00 3.800e+02 2.620e+01 0.941 0.370
## Polycha 5.200e-03 3.850e-03 1.358e+00 2.360e+02 1.310e+01 0.948 0.642
## Tipulid 3.900e-03 7.850e-03 5.020e-01 0.000e+00 1.681e+02 0.952 0.925
## Epigomp 3.300e-03 2.200e-03 1.503e+00 1.830e+02 5.240e+01 0.956 0.975
## Ortochl 3.200e-03 2.470e-03 1.293e+00 1.440e+02 2.400e+01 0.960 0.989
## Nematom 2.900e-03 2.400e-03 1.221e+00 1.310e+02 2.400e+01 0.963 0.411
## Caenis 2.900e-03 1.710e-03 1.702e+00 1.310e+02 1.750e+01 0.967 0.166
## Xestoch 2.700e-03 1.670e-03 1.629e+00 1.310e+02 2.180e+01 0.970 0.919
## Farrode 2.600e-03 1.940e-03 1.358e+00 1.050e+02 7.860e+01 0.973 1.000
## Hexatom 2.600e-03 2.870e-03 8.990e-01 1.050e+02 1.750e+01 0.976 0.966
## Nectops 2.000e-03 2.170e-03 9.240e-01 9.200e+01 2.200e+00 0.978 0.126
## Triplec 1.900e-03 1.470e-03 1.312e+00 7.900e+01 8.700e+00 0.980 0.383
## Hexacyl 1.700e-03 1.960e-03 8.700e-01 7.900e+01 1.960e+01 0.982 0.744
## Platyhe 1.600e-03 1.610e-03 1.002e+00 6.600e+01 6.600e+00 0.984 0.135
## Collemb 1.600e-03 6.180e-03 2.590e-01 1.300e+01 5.460e+01 0.986 0.276
## Ptiloda 1.200e-03 1.500e-03 7.930e-01 5.200e+01 2.200e+00 0.988 0.205
## Phylloi 9.000e-04 1.230e-03 7.290e-01 2.600e+01 1.960e+01 0.989 0.996
## Chrisop 8.000e-04 1.180e-03 6.680e-01 2.600e+01 8.700e+00 0.990 0.570
## Argia 7.000e-04 1.220e-03 5.730e-01 2.600e+01 0.000e+00 0.990 0.442
## Leptne 7.000e-04 3.320e-03 2.060e-01 0.000e+00 3.060e+01 0.991 0.264
## Oecetis 7.000e-04 1.000e-03 6.690e-01 2.600e+01 1.090e+01 0.992 0.799
## Coenag 7.000e-04 6.600e-04 9.890e-01 2.600e+01 0.000e+00 0.993 0.162
## Acari 5.000e-04 1.030e-03 5.070e-01 1.300e+01 8.700e+00 0.993 0.987
## Perisso 5.000e-04 9.100e-04 5.740e-01 2.600e+01 0.000e+00 0.994 0.872
## Calliba 4.000e-04 6.200e-04 6.350e-01 1.300e+01 4.400e+00 0.994 0.886
## Simuli 4.000e-04 6.500e-04 5.910e-01 1.300e+01 4.400e+00 0.995 0.637
## Anophel 4.000e-04 6.200e-04 6.040e-01 1.300e+01 2.200e+00 0.995 0.039 *
## Baetid 4.000e-04 6.100e-04 6.040e-01 1.300e+01 2.200e+00 0.996 0.230
## Perithe 4.000e-04 6.100e-04 6.040e-01 1.300e+01 2.200e+00 0.996 0.042 *
## Banyall 4.000e-04 9.200e-04 3.920e-01 0.000e+00 1.530e+01 0.997 0.994
## Libellu 3.000e-04 6.100e-04 5.730e-01 1.300e+01 0.000e+00 0.997 0.692
## Pelonom 3.000e-04 6.100e-04 5.730e-01 1.300e+01 0.000e+00 0.997 0.001 ***
## Hyrudi 3.000e-04 1.540e-03 2.060e-01 0.000e+00 1.310e+01 0.998 0.369
## Dixella 3.000e-04 5.300e-04 5.740e-01 1.300e+01 0.000e+00 0.998 0.271
## Xenelmi 3.000e-04 6.500e-04 4.130e-01 0.000e+00 1.090e+01 0.998 0.930
## Neoelmi 2.000e-04 1.040e-03 2.060e-01 0.000e+00 8.700e+00 0.999 0.937
## Smicrid 2.000e-04 7.500e-04 2.640e-01 0.000e+00 8.700e+00 0.999 0.770
## Lutroch 1.000e-04 4.000e-04 3.740e-01 0.000e+00 6.600e+00 0.999 0.919
## Forcipo 1.000e-04 3.500e-04 2.980e-01 0.000e+00 4.400e+00 0.999 0.826
## Clognia 1.000e-04 3.500e-04 2.970e-01 0.000e+00 4.400e+00 0.999 0.727
## Anacro 1.000e-04 3.400e-04 2.980e-01 0.000e+00 4.400e+00 0.999 0.963
## Marilia 1.000e-04 4.800e-04 2.060e-01 0.000e+00 4.400e+00 1.000 0.531
## America 1.000e-04 3.300e-04 2.980e-01 0.000e+00 4.400e+00 1.000 0.994
## Phyllog 1.000e-04 2.600e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.974
## Heterag 1.000e-04 2.600e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.991
## Leptohy 1.000e-04 2.400e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.438
## Ambrysu 1.000e-04 2.400e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.870
## Neoplas 1.000e-04 2.400e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.740
## Gyrete 0.000e+00 2.400e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.936
## Hemerodr 0.000e+00 2.400e-04 2.060e-01 0.000e+00 2.200e+00 1.000 0.920
## Guajir 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.262
## Ulmerit 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.461
## Tricory 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.280
## Hetaeri 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.233
## Elga 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.833
## Macrothe 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.236
## Pseudol 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.693
## Neocord 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.227
## Nehalen 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.692
## Lestida 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.253
## Desmogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.254
## Agriogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.251
## Palaemn 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.439
## Rhagove 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.236
## Macrelm 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.705
## Heterel 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.469
## Paracym 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.261
## Scirtid 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.819
## Coleopt 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.972
## Sialis 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.268
## Neargyr 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.782
## Macroste 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.263
## Hydropt 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.251
## Neotric 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.438
## Chimarr 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.907
## Autroti 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.944
## Polycen 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.838
## Cyrnell 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.258
## Protopti 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.255
## Limonia 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.432
## Ephyd2 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.463
## Dolichop 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.582
## Sphaeri 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.455
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Contrast: two_three
##
## average sd ratio ava avb cumsum p
## Tanytar 1.180e-01 8.263e-02 1.428e+00 5.843e+03 2.496e+03 0.248 1.000
## Tanypod 1.075e-01 6.763e-02 1.590e+00 5.833e+03 2.521e+03 0.475 0.733
## Chirono 1.032e-01 7.371e-02 1.399e+00 5.015e+03 2.073e+03 0.692 0.333
## Ceratop 2.814e-02 2.753e-02 1.022e+00 1.424e+03 9.398e+02 0.751 0.997
## Ortochl 1.060e-02 1.603e-02 6.614e-01 3.260e+02 1.185e+02 0.773 0.198
## Hexatom 1.034e-02 1.587e-02 6.517e-01 2.810e+02 1.504e+02 0.795 0.291
## Farrode 9.700e-03 1.311e-02 7.402e-01 2.690e+02 1.834e+02 0.816 0.840
## Macrone 9.550e-03 2.705e-02 3.532e-01 1.500e+02 2.301e+02 0.836 0.714
## Polycha 7.830e-03 1.413e-02 5.541e-01 1.440e+02 1.253e+02 0.852 0.552
## Oligoch 7.050e-03 9.080e-03 7.756e-01 1.590e+02 1.299e+02 0.867 0.902
## Tipulid 5.590e-03 8.140e-03 6.865e-01 6.100e+01 1.469e+02 0.879 1.000
## Epigomp 4.700e-03 4.210e-03 1.115e+00 1.840e+02 1.002e+02 0.888 0.996
## Phylloi 4.470e-03 6.770e-03 6.594e-01 1.170e+02 5.580e+01 0.898 0.476
## Xestoch 4.130e-03 4.530e-03 9.124e-01 1.260e+02 6.610e+01 0.907 0.777
## Banyall 2.900e-03 4.590e-03 6.324e-01 8.100e+01 3.530e+01 0.913 0.375
## Nematom 2.650e-03 6.530e-03 4.062e-01 6.300e+01 3.080e+01 0.918 0.760
## Acari 2.230e-03 3.220e-03 6.911e-01 4.900e+01 3.300e+01 0.923 0.639
## Hexacyl 2.050e-03 2.230e-03 9.156e-01 6.300e+01 1.940e+01 0.927 0.725
## Triplec 1.970e-03 3.450e-03 5.720e-01 5.200e+01 1.820e+01 0.931 0.435
## America 1.960e-03 3.720e-03 5.286e-01 4.800e+01 2.510e+01 0.936 0.309
## Caenis 1.750e-03 2.460e-03 7.093e-01 4.500e+01 1.820e+01 0.939 0.894
## Neoelmi 1.590e-03 3.170e-03 5.003e-01 3.600e+01 2.390e+01 0.943 0.625
## Chimarr 1.500e-03 5.160e-03 2.905e-01 3.400e+01 3.190e+01 0.946 0.523
## Oecetis 1.330e-03 2.140e-03 6.225e-01 3.700e+01 5.700e+00 0.949 0.332
## Simuli 1.260e-03 4.270e-03 2.954e-01 1.500e+01 2.850e+01 0.951 0.612
## Neargyr 1.260e-03 4.570e-03 2.752e-01 3.000e+01 2.300e+00 0.954 0.020 *
## Nectops 1.240e-03 1.930e-03 6.425e-01 3.400e+01 1.030e+01 0.957 0.377
## Copepod 1.160e-03 1.620e-03 7.201e-01 3.700e+01 6.800e+00 0.959 0.936
## Calliba 1.100e-03 2.150e-03 5.111e-01 2.800e+01 1.140e+01 0.961 0.514
## Argia 1.100e-03 2.210e-03 4.967e-01 3.100e+01 3.400e+00 0.964 0.001 ***
## Heterag 1.060e-03 2.030e-03 5.224e-01 2.400e+01 1.480e+01 0.966 0.280
## Perisso 1.000e-03 1.410e-03 7.056e-01 2.400e+01 1.590e+01 0.968 0.471
## Anacro 8.700e-04 1.680e-03 5.208e-01 2.400e+01 1.030e+01 0.970 0.353
## Chrisop 8.500e-04 1.590e-03 5.321e-01 1.800e+01 9.100e+00 0.972 0.889
## Libellu 7.200e-04 1.310e-03 5.512e-01 1.900e+01 6.800e+00 0.973 0.118
## Ptiloda 7.100e-04 1.670e-03 4.259e-01 1.300e+01 1.370e+01 0.975 0.743
## Platyhe 6.700e-04 1.210e-03 5.506e-01 1.000e+01 1.370e+01 0.976 0.966
## Coleopt 6.500e-04 1.420e-03 4.547e-01 1.900e+01 4.600e+00 0.977 0.023 *
## Lutroch 6.200e-04 1.440e-03 4.346e-01 1.000e+01 1.140e+01 0.979 0.789
## Phyllog 6.200e-04 1.440e-03 4.323e-01 7.000e+00 1.140e+01 0.980 0.769
## Collemb 6.200e-04 1.070e-03 5.752e-01 1.500e+01 5.700e+00 0.981 0.802
## Xenelmi 6.000e-04 1.090e-03 5.515e-01 6.000e+00 1.370e+01 0.983 0.995
## Smicrid 5.500e-04 1.360e-03 4.018e-01 1.200e+01 6.800e+00 0.984 0.840
## Autroti 5.400e-04 1.210e-03 4.467e-01 1.300e+01 5.700e+00 0.985 0.258
## Coenag 5.200e-04 1.150e-03 4.515e-01 1.300e+01 2.300e+00 0.986 0.026 *
## Forcipo 4.900e-04 1.720e-03 2.861e-01 9.000e+00 5.700e+00 0.987 0.686
## Macrelm 4.800e-04 2.170e-03 2.194e-01 1.200e+01 2.300e+00 0.988 0.114
## Clognia 4.200e-04 1.890e-03 2.239e-01 1.000e+01 2.300e+00 0.989 0.612
## Hemerodr 4.200e-04 1.170e-03 3.584e-01 6.000e+00 1.030e+01 0.990 0.726
## Polycen 4.000e-04 1.100e-03 3.615e-01 7.000e+00 8.000e+00 0.991 0.521
## Nehalen 3.900e-04 1.390e-03 2.816e-01 6.000e+00 6.800e+00 0.991 0.546
## Elga 3.700e-04 8.400e-04 4.439e-01 1.200e+01 0.000e+00 0.992 0.001 ***
## Dixella 3.700e-04 1.080e-03 3.414e-01 1.000e+01 1.100e+00 0.993 0.052 .
## Gyrete 3.600e-04 7.100e-04 5.030e-01 9.000e+00 3.400e+00 0.994 0.432
## Scirtid 3.200e-04 7.700e-04 4.169e-01 1.200e+01 1.100e+00 0.994 0.011 *
## Hyrudi 2.600e-04 9.100e-04 2.881e-01 7.000e+00 0.000e+00 0.995 0.741
## Baetid 2.300e-04 7.300e-04 3.095e-01 4.000e+00 3.400e+00 0.995 0.736
## Pseudol 2.200e-04 7.300e-04 2.984e-01 1.000e+00 4.600e+00 0.996 0.675
## Ambrysu 2.100e-04 5.600e-04 3.846e-01 3.000e+00 4.600e+00 0.996 0.824
## Neoplas 2.000e-04 7.100e-04 2.748e-01 1.000e+00 4.600e+00 0.997 0.848
## Marilia 1.500e-04 4.500e-04 3.268e-01 4.000e+00 1.100e+00 0.997 0.822
## Dolichop 1.300e-04 4.900e-04 2.660e-01 3.000e+00 2.300e+00 0.997 0.392
## Limonia 1.300e-04 5.900e-04 2.132e-01 3.000e+00 1.100e+00 0.998 0.315
## Sphaeri 1.100e-04 4.900e-04 2.225e-01 1.000e+00 1.100e+00 0.998 0.547
## Neotric 1.000e-04 4.000e-04 2.422e-01 3.000e+00 0.000e+00 0.998 0.003 **
## Palaemn 1.000e-04 4.000e-04 2.445e-01 3.000e+00 0.000e+00 0.998 0.005 **
## Heterel 7.000e-05 3.000e-04 2.436e-01 3.000e+00 0.000e+00 0.998 0.002 **
## Ulmerit 7.000e-05 2.700e-04 2.451e-01 3.000e+00 0.000e+00 0.998 0.002 **
## Macrothe 6.000e-05 4.300e-04 1.471e-01 0.000e+00 2.300e+00 0.999 0.750
## Ephyd2 6.000e-05 2.400e-04 2.446e-01 3.000e+00 0.000e+00 0.999 0.001 ***
## Sialis 6.000e-05 3.400e-04 1.699e-01 1.000e+00 0.000e+00 0.999 0.021 *
## Leptohy 5.000e-05 3.100e-04 1.701e-01 1.000e+00 0.000e+00 0.999 0.735
## Protopti 5.000e-05 3.100e-04 1.702e-01 1.000e+00 0.000e+00 0.999 0.013 *
## Neocord 5.000e-05 2.900e-04 1.703e-01 1.000e+00 0.000e+00 0.999 0.010 **
## Desmogo 5.000e-05 2.800e-04 1.704e-01 1.000e+00 0.000e+00 0.999 0.006 **
## Agriogo 4.000e-05 2.900e-04 1.460e-01 0.000e+00 1.100e+00 0.999 0.746
## Hydropt 4.000e-05 2.900e-04 1.460e-01 0.000e+00 1.100e+00 0.999 0.746
## Hetaeri 4.000e-05 2.700e-04 1.463e-01 0.000e+00 1.100e+00 1.000 0.715
## Cyrnell 4.000e-05 2.300e-04 1.708e-01 1.000e+00 0.000e+00 1.000 0.004 **
## Lestida 4.000e-05 2.500e-04 1.465e-01 0.000e+00 1.100e+00 1.000 0.730
## Paracym 4.000e-05 2.500e-04 1.465e-01 0.000e+00 1.100e+00 1.000 0.763
## Macroste 4.000e-05 2.500e-04 1.466e-01 0.000e+00 1.100e+00 1.000 0.737
## Guajir 3.000e-05 1.900e-04 1.710e-01 1.000e+00 0.000e+00 1.000 0.002 **
## Rhagove 3.000e-05 2.100e-04 1.471e-01 0.000e+00 1.100e+00 1.000 0.750
## Tricory 3.000e-05 1.600e-04 1.711e-01 1.000e+00 0.000e+00 1.000 0.002 **
## Perithe 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.936
## Pelonom 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.730
## Leptne 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.752
## Anophel 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.926
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Contrast: two_four
##
## average sd ratio ava avb cumsum p
## Tanytar 2.174e-01 9.799e-02 2.219e+00 5.843e+03 6.266e+02 0.290 0.001 ***
## Tanypod 2.011e-01 8.402e-02 2.394e+00 5.833e+03 9.301e+02 0.558 0.001 ***
## Chirono 1.688e-01 9.207e-02 1.834e+00 5.015e+03 7.729e+02 0.783 0.001 ***
## Ceratop 3.954e-02 3.556e-02 1.112e+00 1.424e+03 4.760e+02 0.835 0.036 *
## Ortochl 1.166e-02 1.949e-02 5.982e-01 3.260e+02 2.400e+01 0.851 0.135
## Farrode 9.170e-03 1.380e-02 6.641e-01 2.690e+02 7.860e+01 0.863 0.848
## Hexatom 9.090e-03 1.674e-02 5.428e-01 2.810e+02 1.750e+01 0.875 0.607
## Tipulid 8.380e-03 1.296e-02 6.462e-01 6.100e+01 1.681e+02 0.886 0.843
## Oligoch 6.840e-03 7.910e-03 8.653e-01 1.590e+02 6.550e+01 0.895 0.866
## Polycha 6.490e-03 1.423e-02 4.563e-01 1.440e+02 1.310e+01 0.904 0.755
## Epigomp 6.280e-03 5.520e-03 1.138e+00 1.840e+02 5.240e+01 0.912 0.071 .
## Macrone 5.600e-03 9.140e-03 6.130e-01 1.500e+02 2.620e+01 0.920 0.834
## Xestoch 5.030e-03 5.940e-03 8.473e-01 1.260e+02 2.180e+01 0.926 0.115
## Phylloi 4.780e-03 7.980e-03 5.998e-01 1.170e+02 1.960e+01 0.933 0.325
## Banyall 3.070e-03 5.560e-03 5.530e-01 8.100e+01 1.530e+01 0.937 0.304
## Collemb 3.030e-03 1.163e-02 2.609e-01 1.500e+01 5.460e+01 0.941 0.519
## Nematom 3.020e-03 7.770e-03 3.881e-01 6.300e+01 2.400e+01 0.945 0.495
## Hexacyl 2.510e-03 2.690e-03 9.328e-01 6.300e+01 1.960e+01 0.948 0.083 .
## Acari 2.340e-03 4.090e-03 5.720e-01 4.900e+01 8.700e+00 0.952 0.529
## Triplec 2.310e-03 4.500e-03 5.134e-01 5.200e+01 8.700e+00 0.955 0.179
## Caenis 2.200e-03 3.110e-03 7.077e-01 4.500e+01 1.750e+01 0.958 0.292
## Oecetis 1.810e-03 2.790e-03 6.470e-01 3.700e+01 1.090e+01 0.960 0.008 **
## America 1.780e-03 3.820e-03 4.673e-01 4.800e+01 4.400e+00 0.962 0.477
## Neargyr 1.570e-03 6.020e-03 2.607e-01 3.000e+01 0.000e+00 0.964 0.039 *
## Neoelmi 1.560e-03 3.430e-03 4.564e-01 3.600e+01 8.700e+00 0.967 0.573
## Copepod 1.370e-03 2.030e-03 6.758e-01 3.700e+01 2.200e+00 0.968 0.796
## Nectops 1.370e-03 2.490e-03 5.490e-01 3.400e+01 2.200e+00 0.970 0.244
## Argia 1.330e-03 2.890e-03 4.600e-01 3.100e+01 0.000e+00 0.972 0.004 **
## Leptne 1.160e-03 5.720e-03 2.037e-01 0.000e+00 3.060e+01 0.974 0.490
## Calliba 1.150e-03 2.710e-03 4.250e-01 2.800e+01 4.400e+00 0.975 0.387
## Chrisop 1.070e-03 2.040e-03 5.254e-01 1.800e+01 8.700e+00 0.976 0.405
## Heterag 1.010e-03 2.080e-03 4.850e-01 2.400e+01 2.200e+00 0.978 0.452
## Perisso 9.300e-04 1.890e-03 4.904e-01 2.400e+01 0.000e+00 0.979 0.657
## Anacro 9.200e-04 1.820e-03 5.068e-01 2.400e+01 4.400e+00 0.980 0.343
## Hyrudi 8.900e-04 2.950e-03 3.001e-01 7.000e+00 1.310e+01 0.981 0.350
## Chimarr 8.400e-04 3.280e-03 2.561e-01 3.400e+01 0.000e+00 0.983 0.717
## Libellu 7.400e-04 1.610e-03 4.597e-01 1.900e+01 0.000e+00 0.984 0.206
## Smicrid 7.200e-04 1.620e-03 4.466e-01 1.200e+01 8.700e+00 0.985 0.459
## Coleopt 6.700e-04 1.750e-03 3.842e-01 1.900e+01 0.000e+00 0.985 0.120
## Xenelmi 6.500e-04 1.330e-03 4.912e-01 6.000e+00 1.090e+01 0.986 0.937
## Platyhe 6.400e-04 1.370e-03 4.663e-01 1.000e+01 6.600e+00 0.987 0.901
## Clognia 6.200e-04 2.400e-03 2.578e-01 1.000e+01 4.400e+00 0.988 0.269
## Coenag 6.100e-04 1.480e-03 4.122e-01 1.300e+01 0.000e+00 0.989 0.027 *
## Lutroch 6.000e-04 1.220e-03 4.927e-01 1.000e+01 6.600e+00 0.990 0.743
## Simuli 6.000e-04 1.410e-03 4.252e-01 1.500e+01 4.400e+00 0.990 0.825
## Forcipo 5.900e-04 2.160e-03 2.735e-01 9.000e+00 4.400e+00 0.991 0.460
## Macrelm 5.200e-04 2.730e-03 1.917e-01 1.200e+01 0.000e+00 0.992 0.276
## Elga 4.700e-04 1.050e-03 4.484e-01 1.200e+01 0.000e+00 0.993 0.001 ***
## Ptiloda 4.700e-04 1.340e-03 3.511e-01 1.300e+01 2.200e+00 0.993 0.912
## Autroti 4.700e-04 1.310e-03 3.589e-01 1.300e+01 0.000e+00 0.994 0.488
## Phyllog 4.400e-04 1.310e-03 3.359e-01 7.000e+00 2.200e+00 0.994 0.891
## Dixella 4.300e-04 1.390e-03 3.109e-01 1.000e+01 0.000e+00 0.995 0.062 .
## Gyrete 4.100e-04 8.600e-04 4.773e-01 9.000e+00 2.200e+00 0.995 0.245
## Scirtid 3.500e-04 8.900e-04 3.935e-01 1.200e+01 0.000e+00 0.996 0.051 .
## Marilia 3.000e-04 9.100e-04 3.237e-01 4.000e+00 4.400e+00 0.996 0.382
## Baetid 2.600e-04 7.400e-04 3.555e-01 4.000e+00 2.200e+00 0.997 0.564
## Nehalen 2.400e-04 9.800e-04 2.445e-01 6.000e+00 0.000e+00 0.997 0.653
## Hemerodr 2.200e-04 7.000e-04 3.099e-01 6.000e+00 2.200e+00 0.997 0.914
## Polycen 1.900e-04 8.100e-04 2.367e-01 7.000e+00 0.000e+00 0.998 0.898
## Ambrysu 1.800e-04 5.800e-04 3.152e-01 3.000e+00 2.200e+00 0.998 0.837
## Leptohy 1.500e-04 5.700e-04 2.649e-01 1.000e+00 2.200e+00 0.998 0.330
## Neoplas 1.300e-04 4.800e-04 2.599e-01 1.000e+00 2.200e+00 0.998 0.897
## Neotric 1.200e-04 5.000e-04 2.421e-01 3.000e+00 0.000e+00 0.998 0.041 *
## Palaemn 1.200e-04 4.900e-04 2.453e-01 3.000e+00 0.000e+00 0.999 0.043 *
## Limonia 1.200e-04 6.900e-04 1.710e-01 3.000e+00 0.000e+00 0.999 0.367
## Sphaeri 1.000e-04 5.600e-04 1.702e-01 1.000e+00 0.000e+00 0.999 0.551
## Dolichop 9.000e-05 3.900e-04 2.332e-01 3.000e+00 0.000e+00 0.999 0.611
## Perithe 9.000e-05 4.300e-04 2.033e-01 0.000e+00 2.200e+00 0.999 0.559
## Anophel 9.000e-05 4.300e-04 2.033e-01 0.000e+00 2.200e+00 0.999 0.561
## Heterel 8.000e-05 3.500e-04 2.435e-01 3.000e+00 0.000e+00 0.999 0.063 .
## Pseudol 8.000e-05 4.800e-04 1.706e-01 1.000e+00 0.000e+00 0.999 0.905
## Ulmerit 8.000e-05 3.200e-04 2.454e-01 3.000e+00 0.000e+00 0.999 0.080 .
## Sialis 7.000e-05 4.300e-04 1.707e-01 1.000e+00 0.000e+00 1.000 0.236
## Ephyd2 7.000e-05 2.800e-04 2.447e-01 3.000e+00 0.000e+00 1.000 0.137
## Protopti 7.000e-05 3.900e-04 1.708e-01 1.000e+00 0.000e+00 1.000 0.244
## Neocord 6.000e-05 3.600e-04 1.709e-01 1.000e+00 0.000e+00 1.000 0.211
## Desmogo 6.000e-05 3.400e-04 1.710e-01 1.000e+00 0.000e+00 1.000 0.222
## Cyrnell 5.000e-05 2.700e-04 1.711e-01 1.000e+00 0.000e+00 1.000 0.215
## Guajir 4.000e-05 2.200e-04 1.712e-01 1.000e+00 0.000e+00 1.000 0.200
## Tricory 3.000e-05 1.800e-04 1.713e-01 1.000e+00 0.000e+00 1.000 0.218
## Hetaeri 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.529
## Macrothe 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.529
## Lestida 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.545
## Agriogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.559
## Rhagove 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.529
## Paracym 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.532
## Pelonom 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.549
## Macroste 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.569
## Hydropt 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.559
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Contrast: three_four
##
## average sd ratio ava avb cumsum p
## Tanytar 1.437e-01 9.175e-02 1.566e+00 2.496e+03 6.266e+02 0.243 0.424
## Tanypod 1.327e-01 8.957e-02 1.481e+00 2.521e+03 9.301e+02 0.468 0.003 **
## Chirono 1.101e-01 7.628e-02 1.443e+00 2.073e+03 7.729e+02 0.654 0.084 .
## Ceratop 4.607e-02 3.670e-02 1.255e+00 9.398e+02 4.760e+02 0.733 0.001 ***
## Tipulid 1.901e-02 2.563e-02 7.420e-01 1.469e+02 1.681e+02 0.765 0.001 ***
## Farrode 1.417e-02 1.804e-02 7.856e-01 1.834e+02 7.860e+01 0.789 0.017 *
## Macrone 1.251e-02 4.599e-02 2.721e-01 2.301e+02 2.620e+01 0.810 0.250
## Oligoch 1.088e-02 1.807e-02 6.018e-01 1.299e+02 6.550e+01 0.828 0.036 *
## Hexatom 1.015e-02 2.250e-02 4.514e-01 1.504e+02 1.750e+01 0.846 0.410
## Ortochl 9.160e-03 1.630e-02 5.619e-01 1.185e+02 2.400e+01 0.861 0.626
## Polycha 9.120e-03 2.272e-02 4.013e-01 1.253e+02 1.310e+01 0.877 0.256
## Epigomp 7.220e-03 6.760e-03 1.068e+00 1.002e+02 5.240e+01 0.889 0.001 ***
## Xestoch 5.210e-03 6.560e-03 7.934e-01 6.610e+01 2.180e+01 0.898 0.043 *
## Collemb 5.030e-03 2.318e-02 2.172e-01 5.700e+00 5.460e+01 0.906 0.166
## Phylloi 4.740e-03 9.380e-03 5.058e-01 5.580e+01 1.960e+01 0.914 0.326
## Nematom 3.480e-03 5.680e-03 6.127e-01 3.080e+01 2.400e+01 0.920 0.256
## Acari 2.950e-03 5.790e-03 5.101e-01 3.300e+01 8.700e+00 0.925 0.071 .
## Banyall 2.940e-03 4.590e-03 6.413e-01 3.530e+01 1.530e+01 0.930 0.394
## Hexacyl 2.280e-03 2.930e-03 7.783e-01 1.940e+01 1.960e+01 0.934 0.266
## Caenis 2.230e-03 3.870e-03 5.760e-01 1.820e+01 1.750e+01 0.938 0.254
## Simuli 2.220e-03 9.080e-03 2.444e-01 2.850e+01 4.400e+00 0.941 0.145
## Neoelmi 2.040e-03 5.650e-03 3.616e-01 2.390e+01 8.700e+00 0.945 0.226
## Leptne 1.930e-03 9.500e-03 2.032e-01 0.000e+00 3.060e+01 0.948 0.075 .
## America 1.880e-03 5.280e-03 3.565e-01 2.510e+01 4.400e+00 0.951 0.438
## Triplec 1.830e-03 3.820e-03 4.786e-01 1.820e+01 8.700e+00 0.954 0.590
## Xenelmi 1.760e-03 2.920e-03 6.014e-01 1.370e+01 1.090e+01 0.957 0.002 **
## Chimarr 1.460e-03 7.550e-03 1.929e-01 3.190e+01 0.000e+00 0.960 0.446
## Platyhe 1.290e-03 2.670e-03 4.848e-01 1.370e+01 6.600e+00 0.962 0.035 *
## Chrisop 1.280e-03 2.500e-03 5.100e-01 9.100e+00 8.700e+00 0.964 0.094 .
## Perisso 1.260e-03 2.080e-03 6.034e-01 1.590e+01 0.000e+00 0.966 0.035 *
## Calliba 1.210e-03 2.590e-03 4.667e-01 1.140e+01 4.400e+00 0.968 0.317
## Oecetis 1.180e-03 2.430e-03 4.853e-01 5.700e+00 1.090e+01 0.970 0.602
## Phyllog 1.110e-03 3.080e-03 3.618e-01 1.140e+01 2.200e+00 0.972 0.036 *
## Hyrudi 1.070e-03 5.360e-03 1.994e-01 0.000e+00 1.310e+01 0.974 0.203
## Lutroch 1.060e-03 2.670e-03 3.969e-01 1.140e+01 6.600e+00 0.976 0.057 .
## Nectops 1.020e-03 2.430e-03 4.215e-01 1.030e+01 2.200e+00 0.978 0.801
## Smicrid 1.000e-03 3.030e-03 3.291e-01 6.800e+00 8.700e+00 0.979 0.111
## Ptiloda 9.900e-04 2.880e-03 3.422e-01 1.370e+01 2.200e+00 0.981 0.235
## Heterag 9.100e-04 2.630e-03 3.465e-01 1.480e+01 2.200e+00 0.983 0.550
## Anacro 8.200e-04 2.100e-03 3.901e-01 1.030e+01 4.400e+00 0.984 0.528
## Hemerodr 7.900e-04 2.530e-03 3.099e-01 1.030e+01 2.200e+00 0.985 0.053 .
## Forcipo 6.800e-04 1.660e-03 4.082e-01 5.700e+00 4.400e+00 0.986 0.279
## Copepod 6.700e-04 1.620e-03 4.115e-01 6.800e+00 2.200e+00 0.988 0.999
## Neoplas 5.500e-04 1.880e-03 2.910e-01 4.600e+00 2.200e+00 0.989 0.019 *
## Libellu 5.200e-04 1.370e-03 3.765e-01 6.800e+00 0.000e+00 0.989 0.768
## Polycen 5.200e-04 1.850e-03 2.787e-01 8.000e+00 0.000e+00 0.990 0.224
## Autroti 4.700e-04 1.900e-03 2.492e-01 5.700e+00 0.000e+00 0.991 0.511
## Clognia 4.700e-04 1.360e-03 3.459e-01 2.300e+00 4.400e+00 0.992 0.446
## Ambrysu 4.300e-04 1.200e-03 3.589e-01 4.600e+00 2.200e+00 0.993 0.038 *
## Pseudol 4.200e-04 1.640e-03 2.535e-01 4.600e+00 0.000e+00 0.993 0.063 .
## Nehalen 4.100e-04 2.350e-03 1.759e-01 6.800e+00 0.000e+00 0.994 0.367
## Marilia 3.700e-04 1.490e-03 2.450e-01 1.100e+00 4.400e+00 0.995 0.164
## Gyrete 3.600e-04 1.110e-03 3.225e-01 3.400e+00 2.200e+00 0.995 0.443
## Baetid 3.200e-04 1.310e-03 2.472e-01 3.400e+00 2.200e+00 0.996 0.345
## Coleopt 3.000e-04 9.800e-04 3.031e-01 4.600e+00 0.000e+00 0.996 0.899
## Argia 2.000e-04 7.800e-04 2.601e-01 3.400e+00 0.000e+00 0.997 1.000
## Neargyr 1.700e-04 1.150e-03 1.476e-01 2.300e+00 0.000e+00 0.997 0.923
## Leptohy 1.500e-04 7.500e-04 2.021e-01 0.000e+00 2.200e+00 0.997 0.188
## Perithe 1.500e-04 7.500e-04 2.021e-01 0.000e+00 2.200e+00 0.997 0.164
## Anophel 1.500e-04 7.500e-04 2.021e-01 0.000e+00 2.200e+00 0.998 0.162
## Macrelm 1.400e-04 6.400e-04 2.114e-01 2.300e+00 0.000e+00 0.998 0.763
## Coenag 1.400e-04 6.400e-04 2.114e-01 2.300e+00 0.000e+00 0.998 0.990
## Sphaeri 1.200e-04 8.300e-04 1.461e-01 1.100e+00 0.000e+00 0.998 0.412
## Macrothe 1.100e-04 7.700e-04 1.484e-01 2.300e+00 0.000e+00 0.999 0.225
## Agriogo 1.000e-04 7.000e-04 1.469e-01 1.100e+00 0.000e+00 0.999 0.237
## Hydropt 1.000e-04 7.000e-04 1.469e-01 1.100e+00 0.000e+00 0.999 0.237
## Dolichop 1.000e-04 6.400e-04 1.486e-01 2.300e+00 0.000e+00 0.999 0.540
## Hetaeri 9.000e-05 6.000e-04 1.475e-01 1.100e+00 0.000e+00 0.999 0.234
## Lestida 8.000e-05 5.300e-04 1.478e-01 1.100e+00 0.000e+00 0.999 0.234
## Paracym 8.000e-05 5.300e-04 1.478e-01 1.100e+00 0.000e+00 0.999 0.219
## Dixella 8.000e-05 5.300e-04 1.478e-01 1.100e+00 0.000e+00 1.000 0.968
## Macroste 8.000e-05 5.200e-04 1.479e-01 1.100e+00 0.000e+00 1.000 0.233
## Scirtid 7.000e-05 4.900e-04 1.480e-01 1.100e+00 0.000e+00 1.000 0.973
## Limonia 6.000e-05 4.200e-04 1.483e-01 1.100e+00 0.000e+00 1.000 0.726
## Rhagove 6.000e-05 3.900e-04 1.484e-01 1.100e+00 0.000e+00 1.000 0.225
## Guajir 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.634
## Ulmerit 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.877
## Tricory 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.649
## Elga 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.998
## Neocord 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.634
## Desmogo 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.654
## Palaemn 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.868
## Heterel 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.873
## Pelonom 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.625
## Sialis 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.664
## Neotric 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.891
## Cyrnell 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.628
## Protopti 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.639
## Ephyd2 0.000e+00 0.000e+00 NaN 0.000e+00 0.000e+00 1.000 0.870
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## Permutation: free
## Number of permutations: 999
Paso 7. Vamos a determinar las especies que son indecadoras de cada periodo.
max.order = 1 (which is equal to duleg=TRUE) -> singletons #Unicas en un sitio max.order = 2 -> singletons and pairs # Unicas y comparitdas max.order = 3 -> singletons, pairs and triplets
################ Indicator Value ###############
#ind_species <- multipatt(laselva_sp, laselva_full$month, max.order=1,
# func="IndVal",control=how(nperm=999))
#summary(ind_species)