[Vegan-commits] r2115 - in pkg/vegan: R man tests/Examples

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Fri Mar 2 11:48:35 CET 2012


Author: jarioksa
Date: 2012-03-02 11:48:35 +0100 (Fri, 02 Mar 2012)
New Revision: 2115

Modified:
   pkg/vegan/R/simper.R
   pkg/vegan/man/simper.Rd
   pkg/vegan/tests/Examples/vegan-Ex.Rout.save
Log:
Merge simper updates from pull request #6 at github.com

Modified: pkg/vegan/R/simper.R
===================================================================
--- pkg/vegan/R/simper.R	2012-03-01 14:34:17 UTC (rev 2114)
+++ pkg/vegan/R/simper.R	2012-03-02 10:48:35 UTC (rev 2115)
@@ -34,7 +34,7 @@
                 contr[(j-1)*n.a+k, ] <- md / sum(me)	
             }
         }
-        average <- colMeans(contr) * 100
+        average <- colMeans(contr)
         
         if(nperm > 0){
             if (trace)
@@ -51,7 +51,7 @@
                         contrp[(j-1)*n.a+k, ] <- mdp / sum(mep)  
                     }
                 }
-                perm.contr[ ,p] <- colMeans(contrp) * 100
+                perm.contr[ ,p] <- colMeans(contrp)
             }
         p <- (apply(apply(perm.contr, 2, function(x) x >= average), 1, sum) + 1) / (nperm + 1)
         } 
@@ -62,13 +62,13 @@
         overall <- sum(average)
         sdi <- apply(contr, 2, sd)
         ratio <- average / sdi
-        av.a <- colMeans(group.a)
-        av.b <- colMeans(group.b) 
+        ava <- colMeans(group.a)
+        avb <- colMeans(group.b) 
         ord <- order(average, decreasing = TRUE)
-        cusum <- cumsum(average[ord] / overall * 100)
+        cusum <- cumsum(average[ord] / overall)
         out <- list(species = colnames(comm), average = average,
-                    overall = overall, sd = sdi, ratio = ratio, ava = av.a,
-                    avb = av.b, ord = ord, cusum = cusum, p = p)
+                    overall = overall, sd = sdi, ratio = ratio, ava = ava,
+                    avb = avb, ord = ord, cusum = cusum, p = p)
         outlist[[paste(comp[i,1], "_", comp[i,2], sep = "")]] <- out
     }
     attr(outlist, "permutations") <- nperm
@@ -86,7 +86,7 @@
         names(cusum[[i]]) <- spec[[i]]
     }
     ## this probably fails with empty or identical groups that have 0/0 = NaN
-    out <- lapply(cusum, function(z) z[seq_len(min(which(z >= 70)))])
+    out <- lapply(cusum, function(z) z[seq_len(min(which(z >= 0.7)))])
     print(out)
     invisible(x)
 }
@@ -105,7 +105,7 @@
         } 
     } 
     else {
-        out <- lapply(object, function(z) data.frame(cbind(contr = z$average, sd = z$sd, 'contr/sd' = z$ratio, av.a = z$ava, av.b = z$avb, p = z$p)))
+        out <- lapply(object, function(z) data.frame(cbind(contr = z$average, sd = z$sd, 'contr/sd' = z$ratio, ava = z$ava, avb = z$avb, p = z$p)))
     }
     attr(out, "digits") <- digits
     attr(out, "permutations") <- attr(object, "permutations")
@@ -138,4 +138,4 @@
     if ((np <- attr(x, "permutations")) > 0)
         cat("P-values based on", np, "permutations\n")
     invisible(x)
-}
+}
\ No newline at end of file

Modified: pkg/vegan/man/simper.Rd
===================================================================
--- pkg/vegan/man/simper.Rd	2012-03-01 14:34:17 UTC (rev 2114)
+++ pkg/vegan/man/simper.Rd	2012-03-02 10:48:35 UTC (rev 2115)
@@ -48,7 +48,7 @@
   \eqn{d_{jk}=\sum_{i=1}^S d_{ijk}}{d[jk] = sum(i=1..S) d[ijk]}. 
   
   The \code{simper} functions performs pairwise comparisons of groups
-  of sampling units and finds the average percentage contributions
+  of sampling units and finds the average contributions
   of each species to the average overall Bray-Curtis dissimilarity.
 
   The function displays most important species for each pair of
@@ -65,15 +65,15 @@
 \value{
   A list of class \code{"simper"} with following items:
   \item{species}{The species names.}
-  \item{average}{average contribution to overall dissimilarity.}
-  \item{overall}{The overall between-group dissimilarity (times 100).} 
-  \item{sd}{standard deviation of contribution.} 
-  \item{ratio}{mean to sd ratio.}
-  \item{ava, avb}{average abundances per group.}
+  \item{average}{Average contribution to overall dissimilarity.}
+  \item{overall}{The overall between-group dissimilarity.} 
+  \item{sd}{Standard deviation of contribution.} 
+  \item{ratio}{Average to sd ratio.}
+  \item{ava, avb}{Average abundances per group.}
   \item{ord}{An index vector to order vectors by their contribution or
     order \code{cusum} back to the original data order.}
-  \item{cusum}{Ordered cumulative per cent contribution.}
-  \item{p}{permutation \eqn{p}-value. Probability of getting a larger
+  \item{cusum}{Ordered cumulative contribution.}
+  \item{p}{Permutation \eqn{p}-value. Probability of getting a larger
     or equal average contribution in random permutation of the group
     factor.}
 }

Modified: pkg/vegan/tests/Examples/vegan-Ex.Rout.save
===================================================================
--- pkg/vegan/tests/Examples/vegan-Ex.Rout.save	2012-03-01 14:34:17 UTC (rev 2114)
+++ pkg/vegan/tests/Examples/vegan-Ex.Rout.save	2012-03-02 10:48:35 UTC (rev 2115)
@@ -5998,246 +5998,240 @@
 cumulative contributions of most influential species:
 
 $BF_SF
-   Agrsto    Alogen    Lolper    Trirep    Poatri    Leoaut    Brohor    Achmil 
- 9.824271 18.254830 25.956958 33.367870 40.734444 47.729205 53.120026 57.946526 
-   Plalan    Elyrep    Brarut 
-62.522255 67.016196 71.098133 
+    Agrsto     Alogen     Lolper     Trirep     Poatri     Leoaut     Brohor     Achmil     Plalan     Elyrep     Brarut 
+0.09824271 0.18254830 0.25956958 0.33367870 0.40734444 0.47729205 0.53120026 0.57946526 0.62522255 0.67016196 0.71098133 
 
 $BF_HF
-   Rumace    Poatri    Plalan    Brohor    Lolper    Elyrep    Trirep    Antodo 
- 8.163219 15.193797 21.918333 27.967181 33.969561 39.843338 45.298204 50.276849 
-   Achmil    Brarut    Alogen    Tripra    Junart 
-55.222648 60.021994 64.584333 69.126471 73.366621 
+    Rumace     Poatri     Plalan     Brohor     Lolper     Elyrep     Trirep     Antodo     Achmil     Brarut     Alogen     Tripra 
+0.08163219 0.15193797 0.21918333 0.27967181 0.33969561 0.39843338 0.45298204 0.50276849 0.55222648 0.60021994 0.64584333 0.69126471 
+    Junart 
+0.73366621 
 
 $BF_NM
-  Lolper   Poatri   Poapra   Trirep   Brohor   Brarut   Elepal   Agrsto 
-12.42718 19.92126 27.11756 34.14609 39.58520 44.48077 49.10724 53.69083 
-  Achmil   Leoaut   Antodo   Plalan 
-58.23926 62.53645 66.38182 70.12577 
+   Lolper    Poatri    Poapra    Trirep    Brohor    Brarut    Elepal    Agrsto    Achmil    Leoaut    Antodo    Plalan 
+0.1242718 0.1992126 0.2711756 0.3414609 0.3958520 0.4448077 0.4910724 0.5369083 0.5823926 0.6253645 0.6638182 0.7012577 
 
 $SF_HF
-   Agrsto    Alogen    Lolper    Plalan    Rumace    Elyrep    Poatri    Brarut 
- 8.350879 16.534834 23.934930 30.843624 37.716139 43.334492 48.351753 52.804045 
-   Elepal    Poapra    Antodo    Sagpro    Tripra 
-57.205850 61.423981 65.549838 69.628951 73.696831 
+    Agrsto     Alogen     Lolper     Plalan     Rumace     Elyrep     Poatri     Brarut     Elepal     Poapra     Antodo     Sagpro 
+0.08350879 0.16534834 0.23934930 0.30843624 0.37716139 0.43334492 0.48351753 0.52804045 0.57205850 0.61423981 0.65549838 0.69628951 
+    Tripra 
+0.73696831 
 
 $SF_NM
-  Poatri   Alogen   Agrsto   Lolper   Elepal   Poapra   Brarut   Elyrep 
-10.13601 19.35731 26.67383 33.77578 39.99419 45.26707 50.44725 55.05643 
-  Leoaut   Trirep   Sagpro   Salrep 
-59.26117 63.20111 67.12478 70.91528 
+   Poatri    Alogen    Agrsto    Lolper    Elepal    Poapra    Brarut    Elyrep    Leoaut    Trirep    Sagpro    Salrep 
+0.1013601 0.1935731 0.2667383 0.3377578 0.3999419 0.4526707 0.5044725 0.5505643 0.5926117 0.6320111 0.6712478 0.7091528 
 
 $HF_NM
-   Poatri    Lolper    Rumace    Poapra    Plalan    Brarut    Elepal    Agrsto 
- 9.913221 17.468460 23.917190 29.701331 35.469313 40.365488 44.804851 49.226546 
-   Trirep    Elyrep    Antodo    Junart    Tripra    Salrep 
-53.434466 57.564661 61.543243 65.341300 68.921695 72.432408 
+    Poatri     Lolper     Rumace     Poapra     Plalan     Brarut     Elepal     Agrsto     Trirep     Elyrep     Antodo     Junart 
+0.09913221 0.17468460 0.23917190 0.29701331 0.35469313 0.40365488 0.44804851 0.49226546 0.53434466 0.57564661 0.61543243 0.65341300 
+    Tripra     Salrep 
+0.68921695 0.72432408 
 
 > summary(sim)
 
 Contrast: BF_SF 
 
-        contr       sd  ratio   av.a   av.b  cumsum
-Agrsto 6.1374 0.034193 179.49 0.0000 4.6667   9.824
-Alogen 5.2667 0.036476 144.39 0.6667 4.3333  18.255
-Lolper 4.8116 0.039445 121.98 6.0000 3.0000  25.957
-Trirep 4.6297 0.025525 181.38 4.6667 1.3333  33.368
-Poatri 4.6020 0.033801 136.15 3.6667 4.6667  40.734
-Leoaut 4.3697 0.024922 175.34 4.3333 1.3333  47.729
-Brohor 3.3677 0.025860 130.23 2.6667 0.5000  53.120
-Achmil 3.0152 0.020821 144.82 2.3333 0.1667  57.947
-Plalan 2.8585 0.021549 132.65 2.0000 0.0000  62.522
-Elyrep 2.8074 0.029778  94.28 1.3333 2.0000  67.016
-Brarut 2.5501 0.023902 106.69 2.0000 2.0000  71.098
-Poapra 2.5129 0.023967 104.85 4.0000 2.5000  75.121
-Sagpro 2.4326 0.022149 109.83 0.6667 1.8333  79.014
-Belper 1.9859 0.017088 116.22 1.6667 0.6667  82.193
-Elepal 1.8611 0.042958  43.33 0.0000 1.3333  85.172
-Antodo 1.7543 0.025804  67.98 1.3333 0.0000  87.981
-Junbuf 1.6031 0.023708  67.62 0.0000 1.1667  90.547
-Viclat 1.4671 0.013306 110.26 1.0000 0.0000  92.895
-Hyprad 1.0286 0.015198  67.68 0.6667 0.0000  94.542
-Ranfla 0.9306 0.013595  68.45 0.0000 0.6667  96.031
-Junart 0.6979 0.016109  43.33 0.0000 0.5000  97.148
-Calcus 0.6979 0.016109  43.33 0.0000 0.5000  98.266
-Rumace 0.4526 0.010444  43.33 0.0000 0.3333  98.990
-Cirarv 0.3983 0.009185  43.36 0.0000 0.3333  99.628
-Chealb 0.2326 0.005370  43.33 0.0000 0.1667 100.000
-Empnig 0.0000 0.000000    NaN 0.0000 0.0000 100.000
-Airpra 0.0000 0.000000    NaN 0.0000 0.0000 100.000
-Potpal 0.0000 0.000000    NaN 0.0000 0.0000 100.000
-Tripra 0.0000 0.000000    NaN 0.0000 0.0000 100.000
-Salrep 0.0000 0.000000    NaN 0.0000 0.0000 100.000
+          contr       sd  ratio   av.a   av.b  cumsum
+Agrsto 0.061374 0.034193 1.7949 0.0000 4.6667 0.09824
+Alogen 0.052667 0.036476 1.4439 0.6667 4.3333 0.18255
+Lolper 0.048116 0.039445 1.2198 6.0000 3.0000 0.25957
+Trirep 0.046297 0.025525 1.8138 4.6667 1.3333 0.33368
+Poatri 0.046020 0.033801 1.3615 3.6667 4.6667 0.40734
+Leoaut 0.043697 0.024922 1.7534 4.3333 1.3333 0.47729
+Brohor 0.033677 0.025860 1.3023 2.6667 0.5000 0.53120
+Achmil 0.030152 0.020821 1.4482 2.3333 0.1667 0.57947
+Plalan 0.028585 0.021549 1.3265 2.0000 0.0000 0.62522
+Elyrep 0.028074 0.029778 0.9428 1.3333 2.0000 0.67016
+Brarut 0.025501 0.023902 1.0669 2.0000 2.0000 0.71098
+Poapra 0.025129 0.023967 1.0485 4.0000 2.5000 0.75121
+Sagpro 0.024326 0.022149 1.0983 0.6667 1.8333 0.79014
+Belper 0.019859 0.017088 1.1622 1.6667 0.6667 0.82193
+Elepal 0.018611 0.042958 0.4333 0.0000 1.3333 0.85172
+Antodo 0.017543 0.025804 0.6798 1.3333 0.0000 0.87981
+Junbuf 0.016031 0.023708 0.6762 0.0000 1.1667 0.90547
+Viclat 0.014671 0.013306 1.1026 1.0000 0.0000 0.92895
+Hyprad 0.010286 0.015198 0.6768 0.6667 0.0000 0.94542
+Ranfla 0.009306 0.013595 0.6845 0.0000 0.6667 0.96031
+Junart 0.006979 0.016109 0.4333 0.0000 0.5000 0.97148
+Calcus 0.006979 0.016109 0.4333 0.0000 0.5000 0.98266
+Rumace 0.004526 0.010444 0.4333 0.0000 0.3333 0.98990
+Cirarv 0.003983 0.009185 0.4336 0.0000 0.3333 0.99628
+Chealb 0.002326 0.005370 0.4333 0.0000 0.1667 1.00000
+Empnig 0.000000 0.000000    NaN 0.0000 0.0000 1.00000
+Airpra 0.000000 0.000000    NaN 0.0000 0.0000 1.00000
+Potpal 0.000000 0.000000    NaN 0.0000 0.0000 1.00000
+Tripra 0.000000 0.000000    NaN 0.0000 0.0000 1.00000
+Salrep 0.000000 0.000000    NaN 0.0000 0.0000 1.00000
 
 Contrast: BF_HF 
 
-        contr      sd  ratio   av.a av.b  cumsum
-Rumace 3.8666 0.02606 148.38 0.0000  3.2   8.163
-Poatri 3.3301 0.02579 129.11 3.6667  4.8  15.194
-Plalan 3.1852 0.01830 174.01 2.0000  3.0  21.918
-Brohor 2.8651 0.01799 159.26 2.6667  0.8  27.967
-Lolper 2.8431 0.02215 128.34 6.0000  4.0  33.970
-Elyrep 2.7822 0.02959  94.04 1.3333  2.0  39.843
-Trirep 2.5838 0.01656 156.03 4.6667  2.8  45.298
-Antodo 2.3582 0.02042 115.47 1.3333  1.8  50.277
-Achmil 2.3426 0.01474 158.93 2.3333  1.2  55.223
-Brarut 2.2733 0.01802 126.17 2.0000  2.8  60.022
-Alogen 2.1610 0.02308  93.63 0.6667  1.6  64.584
-Tripra 2.1514 0.02207  97.47 0.0000  1.8  69.126
-Junart 2.0084 0.02555  78.60 0.0000  1.6  73.367
-Leoaut 1.9318 0.01356 142.41 4.3333  2.8  77.445
-Belper 1.8290 0.01486 123.05 1.6667  0.4  81.306
-Agrsto 1.7605 0.02284  77.08 0.0000  1.4  85.023
-Junbuf 1.5000 0.02066  72.60 0.0000  1.2  88.190
-Viclat 1.2848 0.01140 112.74 1.0000  0.0  90.903
-Sagpro 1.1685 0.01297  90.08 0.6667  0.8  93.369
-Elepal 1.0169 0.02111  48.17 0.0000  0.8  95.516
-Hyprad 0.8950 0.01312  68.24 0.6667  0.0  97.406
-Poapra 0.7203 0.01010  71.33 4.0000  3.4  98.927
-Ranfla 0.5084 0.01055  48.17 0.0000  0.4 100.000
-Empnig 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Airpra 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Cirarv 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Potpal 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Calcus 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Salrep 0.0000 0.00000    NaN 0.0000  0.0 100.000
-Chealb 0.0000 0.00000    NaN 0.0000  0.0 100.000
+          contr      sd  ratio   av.a av.b  cumsum
+Rumace 0.038666 0.02606 1.4838 0.0000  3.2 0.08163
+Poatri 0.033301 0.02579 1.2911 3.6667  4.8 0.15194
+Plalan 0.031852 0.01830 1.7401 2.0000  3.0 0.21918
+Brohor 0.028651 0.01799 1.5926 2.6667  0.8 0.27967
+Lolper 0.028431 0.02215 1.2834 6.0000  4.0 0.33970
+Elyrep 0.027822 0.02959 0.9404 1.3333  2.0 0.39843
+Trirep 0.025838 0.01656 1.5603 4.6667  2.8 0.45298
+Antodo 0.023582 0.02042 1.1547 1.3333  1.8 0.50277
+Achmil 0.023426 0.01474 1.5893 2.3333  1.2 0.55223
+Brarut 0.022733 0.01802 1.2617 2.0000  2.8 0.60022
+Alogen 0.021610 0.02308 0.9363 0.6667  1.6 0.64584
+Tripra 0.021514 0.02207 0.9747 0.0000  1.8 0.69126
+Junart 0.020084 0.02555 0.7860 0.0000  1.6 0.73367
+Leoaut 0.019318 0.01356 1.4241 4.3333  2.8 0.77445
+Belper 0.018290 0.01486 1.2305 1.6667  0.4 0.81306
+Agrsto 0.017605 0.02284 0.7708 0.0000  1.4 0.85023
+Junbuf 0.015000 0.02066 0.7260 0.0000  1.2 0.88190
+Viclat 0.012848 0.01140 1.1274 1.0000  0.0 0.90903
+Sagpro 0.011685 0.01297 0.9008 0.6667  0.8 0.93369
+Elepal 0.010169 0.02111 0.4817 0.0000  0.8 0.95516
+Hyprad 0.008950 0.01312 0.6824 0.6667  0.0 0.97406
+Poapra 0.007203 0.01010 0.7133 4.0000  3.4 0.98927
+Ranfla 0.005084 0.01055 0.4817 0.0000  0.4 1.00000
+Empnig 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Airpra 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Cirarv 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Potpal 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Calcus 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Salrep 0.000000 0.00000    NaN 0.0000  0.0 1.00000
+Chealb 0.000000 0.00000    NaN 0.0000  0.0 1.00000
 
 Contrast: BF_NM 
 
-        contr      sd  ratio   av.a   av.b cumsum
-Lolper 9.0681 0.02644 342.92 6.0000 0.3333  12.43
-Poatri 5.4684 0.04465 122.48 3.6667 0.0000  19.92
-Poapra 5.2511 0.01813 289.66 4.0000 0.6667  27.12
-Trirep 5.1287 0.02756 186.11 4.6667 1.8333  34.15
-Brohor 3.9689 0.02920 135.90 2.6667 0.0000  39.59
-Brarut 3.5723 0.02869 124.52 2.0000 2.8333  44.48
-Elepal 3.3759 0.03573  94.49 0.0000 2.1667  49.11
-Agrsto 3.3446 0.03473  96.30 0.0000 2.1667  53.69
-Achmil 3.3190 0.02338 141.98 2.3333 0.3333  58.24
-Leoaut 3.1356 0.02026 154.80 4.3333 3.1667  62.54
-Antodo 2.8060 0.03295  85.17 1.3333 1.3333  66.38
-Plalan 2.7319 0.02193 124.58 2.0000 0.8333  70.13
-Salrep 2.6770 0.02927  91.45 0.0000 1.8333  73.79
-Belper 2.3529 0.01909 123.22 1.6667 0.3333  77.02
-Hyprad 2.1721 0.02450  88.64 0.6667 1.1667  80.00
-Ranfla 2.0314 0.02275  89.28 0.0000 1.3333  82.78
-Elyrep 1.9993 0.02926  68.33 1.3333 0.0000  85.52
-Calcus 1.7833 0.02681  66.53 0.0000 1.1667  87.96
-Junart 1.7694 0.02600  68.06 0.0000 1.1667  90.39
-Viclat 1.5773 0.01447 109.02 1.0000 0.1667  92.55
-Sagpro 1.5432 0.01857  83.10 0.6667 0.5000  94.66
-Airpra 1.3410 0.01969  68.09 0.0000 0.8333  96.50
-Potpal 1.0739 0.01571  68.35 0.0000 0.6667  97.97
-Alogen 0.9997 0.01463  68.33 0.6667 0.0000  99.34
-Empnig 0.4787 0.01105  43.32 0.0000 0.3333 100.00
-Junbuf 0.0000 0.00000    NaN 0.0000 0.0000 100.00
-Rumace 0.0000 0.00000    NaN 0.0000 0.0000 100.00
-Cirarv 0.0000 0.00000    NaN 0.0000 0.0000 100.00
-Tripra 0.0000 0.00000    NaN 0.0000 0.0000 100.00
-Chealb 0.0000 0.00000    NaN 0.0000 0.0000 100.00
+          contr      sd  ratio   av.a   av.b cumsum
+Lolper 0.090681 0.02644 3.4292 6.0000 0.3333 0.1243
+Poatri 0.054684 0.04465 1.2248 3.6667 0.0000 0.1992
+Poapra 0.052511 0.01813 2.8966 4.0000 0.6667 0.2712
+Trirep 0.051287 0.02756 1.8611 4.6667 1.8333 0.3415
+Brohor 0.039689 0.02920 1.3590 2.6667 0.0000 0.3959
+Brarut 0.035723 0.02869 1.2452 2.0000 2.8333 0.4448
+Elepal 0.033759 0.03573 0.9449 0.0000 2.1667 0.4911
+Agrsto 0.033446 0.03473 0.9630 0.0000 2.1667 0.5369
+Achmil 0.033190 0.02338 1.4198 2.3333 0.3333 0.5824
+Leoaut 0.031356 0.02026 1.5480 4.3333 3.1667 0.6254
+Antodo 0.028060 0.03295 0.8517 1.3333 1.3333 0.6638
+Plalan 0.027319 0.02193 1.2458 2.0000 0.8333 0.7013
+Salrep 0.026770 0.02927 0.9145 0.0000 1.8333 0.7379
+Belper 0.023529 0.01909 1.2322 1.6667 0.3333 0.7702
+Hyprad 0.021721 0.02450 0.8864 0.6667 1.1667 0.8000
+Ranfla 0.020314 0.02275 0.8928 0.0000 1.3333 0.8278
+Elyrep 0.019993 0.02926 0.6833 1.3333 0.0000 0.8552
+Calcus 0.017833 0.02681 0.6653 0.0000 1.1667 0.8796
+Junart 0.017694 0.02600 0.6806 0.0000 1.1667 0.9039
+Viclat 0.015773 0.01447 1.0902 1.0000 0.1667 0.9255
+Sagpro 0.015432 0.01857 0.8310 0.6667 0.5000 0.9466
+Airpra 0.013410 0.01969 0.6809 0.0000 0.8333 0.9650
+Potpal 0.010739 0.01571 0.6835 0.0000 0.6667 0.9797
+Alogen 0.009997 0.01463 0.6833 0.6667 0.0000 0.9934
+Empnig 0.004787 0.01105 0.4332 0.0000 0.3333 1.0000
+Junbuf 0.000000 0.00000    NaN 0.0000 0.0000 1.0000
+Rumace 0.000000 0.00000    NaN 0.0000 0.0000 1.0000
+Cirarv 0.000000 0.00000    NaN 0.0000 0.0000 1.0000
+Tripra 0.000000 0.00000    NaN 0.0000 0.0000 1.0000
+Chealb 0.000000 0.00000    NaN 0.0000 0.0000 1.0000
 
 Contrast: SF_HF 
 
-        contr       sd  ratio   av.a av.b  cumsum
-Agrsto 4.7380 0.031273 151.51 4.6667  1.4   8.351
-Alogen 4.6433 0.032897 141.15 4.3333  1.6  16.535
-Lolper 4.1986 0.027007 155.46 3.0000  4.0  23.935
-Plalan 3.9198 0.033208 118.04 0.0000  3.0  30.844
-Rumace 3.8992 0.027369 142.47 0.3333  3.2  37.716
-Elyrep 3.1877 0.029550 107.87 2.0000  2.0  43.334
-Poatri 2.8466 0.021522 132.27 4.6667  4.8  48.352
-Brarut 2.5261 0.021044 120.04 2.0000  2.8  52.804
-Elepal 2.4974 0.038877  64.24 1.3333  0.8  57.206
-Poapra 2.3932 0.019180 124.78 2.5000  3.4  61.424
-Antodo 2.3409 0.021430 109.23 0.0000  1.8  65.550
-Sagpro 2.3144 0.020479 113.01 1.8333  0.8  69.629
-Tripra 2.3080 0.023432  98.50 0.0000  1.8  73.697
-Junart 2.2850 0.025677  88.99 0.5000  1.6  77.724
-Trirep 2.2383 0.019487 114.86 1.3333  2.8  81.669
-Junbuf 2.1643 0.022237  97.33 1.1667  1.2  85.484
-Leoaut 2.0509 0.016422 124.89 1.3333  2.8  89.099
-Achmil 1.5183 0.011393 133.26 0.1667  1.2  91.775
-Brohor 1.3375 0.014504  92.22 0.5000  0.8  94.132
-Ranfla 1.0661 0.013387  79.64 0.6667  0.4  96.011
-Belper 0.9991 0.012571  79.48 0.6667  0.4  97.772
-Calcus 0.6623 0.015076  43.93 0.5000  0.0  98.939
-Cirarv 0.3809 0.008669  43.94 0.3333  0.0  99.611
-Chealb 0.2208 0.005025  43.93 0.1667  0.0 100.000
-Empnig 0.0000 0.000000    NaN 0.0000  0.0 100.000
-Airpra 0.0000 0.000000    NaN 0.0000  0.0 100.000
-Viclat 0.0000 0.000000    NaN 0.0000  0.0 100.000
-Hyprad 0.0000 0.000000    NaN 0.0000  0.0 100.000
-Potpal 0.0000 0.000000    NaN 0.0000  0.0 100.000
-Salrep 0.0000 0.000000    NaN 0.0000  0.0 100.000
+          contr       sd  ratio   av.a av.b  cumsum
+Agrsto 0.047380 0.031273 1.5151 4.6667  1.4 0.08351
+Alogen 0.046433 0.032897 1.4115 4.3333  1.6 0.16535
+Lolper 0.041986 0.027007 1.5546 3.0000  4.0 0.23935
+Plalan 0.039198 0.033208 1.1804 0.0000  3.0 0.30844
+Rumace 0.038992 0.027369 1.4247 0.3333  3.2 0.37716
+Elyrep 0.031877 0.029550 1.0787 2.0000  2.0 0.43334
+Poatri 0.028466 0.021522 1.3227 4.6667  4.8 0.48352
+Brarut 0.025261 0.021044 1.2004 2.0000  2.8 0.52804
+Elepal 0.024974 0.038877 0.6424 1.3333  0.8 0.57206
+Poapra 0.023932 0.019180 1.2478 2.5000  3.4 0.61424
+Antodo 0.023409 0.021430 1.0923 0.0000  1.8 0.65550
+Sagpro 0.023144 0.020479 1.1301 1.8333  0.8 0.69629
+Tripra 0.023080 0.023432 0.9850 0.0000  1.8 0.73697
+Junart 0.022850 0.025677 0.8899 0.5000  1.6 0.77724
+Trirep 0.022383 0.019487 1.1486 1.3333  2.8 0.81669
+Junbuf 0.021643 0.022237 0.9733 1.1667  1.2 0.85484
+Leoaut 0.020509 0.016422 1.2489 1.3333  2.8 0.89099
+Achmil 0.015183 0.011393 1.3326 0.1667  1.2 0.91775
+Brohor 0.013375 0.014504 0.9222 0.5000  0.8 0.94132
+Ranfla 0.010661 0.013387 0.7964 0.6667  0.4 0.96011
+Belper 0.009991 0.012571 0.7948 0.6667  0.4 0.97772
+Calcus 0.006623 0.015076 0.4393 0.5000  0.0 0.98939
+Cirarv 0.003809 0.008669 0.4394 0.3333  0.0 0.99611
+Chealb 0.002208 0.005025 0.4393 0.1667  0.0 1.00000
+Empnig 0.000000 0.000000    NaN 0.0000  0.0 1.00000
+Airpra 0.000000 0.000000    NaN 0.0000  0.0 1.00000
+Viclat 0.000000 0.000000    NaN 0.0000  0.0 1.00000
+Hyprad 0.000000 0.000000    NaN 0.0000  0.0 1.00000
+Potpal 0.000000 0.000000    NaN 0.0000  0.0 1.00000
+Salrep 0.000000 0.000000    NaN 0.0000  0.0 1.00000
 
 Contrast: SF_NM 
 
-        contr       sd  ratio   av.a   av.b cumsum
-Poatri 7.8284 0.040947 191.18 4.6667 0.0000  10.14
-Alogen 7.1219 0.046958 151.67 4.3333 0.0000  19.36
-Agrsto 5.6508 0.044176 127.92 4.6667 2.1667  26.67
-Lolper 5.4851 0.059914  91.55 3.0000 0.3333  33.78
-Elepal 4.8027 0.047168 101.82 1.3333 2.1667  39.99
-Poapra 4.0724 0.031790 128.10 2.5000 0.6667  45.27
-Brarut 4.0008 0.034398 116.31 2.0000 2.8333  50.45
-Elyrep 3.5598 0.038515  92.43 2.0000 0.0000  55.06
-Leoaut 3.2475 0.034813  93.28 1.3333 3.1667  59.26
-Trirep 3.0430 0.031634  96.19 1.3333 1.8333  63.20
-Sagpro 3.0304 0.030477  99.43 1.8333 0.5000  67.12
-Salrep 2.9275 0.032014  91.44 0.0000 1.8333  70.92
-Antodo 2.4541 0.036694  66.88 0.0000 1.3333  74.09
-Calcus 2.2763 0.029443  77.31 0.5000 1.1667  77.04
-Ranfla 2.2566 0.022819  98.89 0.6667 1.3333  79.96
-Junart 2.2543 0.028598  78.83 0.5000 1.1667  82.88
-Hyprad 2.0108 0.031291  64.26 0.0000 1.1667  85.48
-Junbuf 1.9860 0.029034  68.40 1.1667 0.0000  88.06
-Plalan 1.5420 0.022769  67.72 0.0000 0.8333  90.05
-Airpra 1.4883 0.021881  68.02 0.0000 0.8333  91.98
-Belper 1.2317 0.015921  77.37 0.6667 0.3333  93.57
-Potpal 1.1883 0.017407  68.26 0.0000 0.6667  95.11
-Achmil 0.9294 0.014931  62.24 0.1667 0.3333  96.32
-Brohor 0.7172 0.016333  43.91 0.5000 0.0000  97.24
-Rumace 0.5590 0.012751  43.84 0.3333 0.0000  97.97
-Empnig 0.5225 0.012001  43.54 0.0000 0.3333  98.64
-Cirarv 0.4782 0.010889  43.91 0.3333 0.0000  99.26
-Chealb 0.2893 0.006602  43.82 0.1667 0.0000  99.64
-Viclat 0.2792 0.006425  43.45 0.0000 0.1667 100.00
-Tripra 0.0000 0.000000    NaN 0.0000 0.0000 100.00
+          contr       sd  ratio   av.a   av.b cumsum
+Poatri 0.078284 0.040947 1.9118 4.6667 0.0000 0.1014
+Alogen 0.071219 0.046958 1.5167 4.3333 0.0000 0.1936
+Agrsto 0.056508 0.044176 1.2792 4.6667 2.1667 0.2667
+Lolper 0.054851 0.059914 0.9155 3.0000 0.3333 0.3378
+Elepal 0.048027 0.047168 1.0182 1.3333 2.1667 0.3999
+Poapra 0.040724 0.031790 1.2810 2.5000 0.6667 0.4527
+Brarut 0.040008 0.034398 1.1631 2.0000 2.8333 0.5045
+Elyrep 0.035598 0.038515 0.9243 2.0000 0.0000 0.5506
+Leoaut 0.032475 0.034813 0.9328 1.3333 3.1667 0.5926
+Trirep 0.030430 0.031634 0.9619 1.3333 1.8333 0.6320
+Sagpro 0.030304 0.030477 0.9943 1.8333 0.5000 0.6712
+Salrep 0.029275 0.032014 0.9144 0.0000 1.8333 0.7092
+Antodo 0.024541 0.036694 0.6688 0.0000 1.3333 0.7409
+Calcus 0.022763 0.029443 0.7731 0.5000 1.1667 0.7704
+Ranfla 0.022566 0.022819 0.9889 0.6667 1.3333 0.7996
+Junart 0.022543 0.028598 0.7883 0.5000 1.1667 0.8288
+Hyprad 0.020108 0.031291 0.6426 0.0000 1.1667 0.8548
+Junbuf 0.019860 0.029034 0.6840 1.1667 0.0000 0.8806
+Plalan 0.015420 0.022769 0.6772 0.0000 0.8333 0.9005
+Airpra 0.014883 0.021881 0.6802 0.0000 0.8333 0.9198
+Belper 0.012317 0.015921 0.7737 0.6667 0.3333 0.9357
+Potpal 0.011883 0.017407 0.6826 0.0000 0.6667 0.9511
+Achmil 0.009294 0.014931 0.6224 0.1667 0.3333 0.9632
+Brohor 0.007172 0.016333 0.4391 0.5000 0.0000 0.9724
+Rumace 0.005590 0.012751 0.4384 0.3333 0.0000 0.9797
+Empnig 0.005225 0.012001 0.4354 0.0000 0.3333 0.9864
+Cirarv 0.004782 0.010889 0.4391 0.3333 0.0000 0.9926
+Chealb 0.002893 0.006602 0.4382 0.1667 0.0000 0.9964
+Viclat 0.002792 0.006425 0.4345 0.0000 0.1667 1.0000
+Tripra 0.000000 0.000000    NaN 0.0000 0.0000 1.0000
 
 Contrast: HF_NM 
 
-        contr       sd  ratio av.a   av.b  cumsum
-Poatri 7.1553 0.013681 523.02  4.8 0.0000   9.913
-Lolper 5.4533 0.029625 184.08  4.0 0.3333  17.468
-Rumace 4.6546 0.030806 151.09  3.2 0.0000  23.917
-Poapra 4.1750 0.018852 221.46  3.4 0.6667  29.701
-Plalan 4.1633 0.029560 140.84  3.0 0.8333  35.469
-Brarut 3.5340 0.020104 175.79  2.8 2.8333  40.365
-Elepal 3.2043 0.032315  99.16  0.8 2.1667  44.805
-Agrsto 3.1915 0.028887 110.48  1.4 2.1667  49.227
-Trirep 3.0372 0.022871 132.80  2.8 1.8333  53.434
-Elyrep 2.9811 0.038676  77.08  2.0 0.0000  57.565
-Antodo 2.8717 0.024799 115.80  1.8 1.3333  61.543
-Junart 2.7414 0.028537  96.07  1.6 1.1667  65.341
-Tripra 2.5843 0.025972  99.51  1.8 0.0000  68.922
-Salrep 2.5340 0.027291  92.85  0.0 1.8333  72.432
-Alogen 2.4459 0.032399  75.49  1.6 0.0000  75.821
-Leoaut 2.0703 0.014125 146.58  2.8 3.1667  78.689
-Ranfla 1.9285 0.019939  96.72  0.4 1.3333  81.361
-Junbuf 1.8181 0.024648  73.76  1.2 0.0000  83.880
-Hyprad 1.7141 0.026548  64.57  0.0 1.1667  86.255
-Calcus 1.6833 0.024901  67.60  0.0 1.1667  88.587
-Achmil 1.6555 0.014900 111.11  1.2 0.3333  90.881
-Sagpro 1.5282 0.016535  92.43  0.8 0.5000  92.998
-Airpra 1.2605 0.018243  69.10  0.0 0.8333  94.744
-Brohor 1.2094 0.015169  79.73  0.8 0.0000  96.420
-Potpal 1.0105 0.014556  69.42  0.0 0.6667  97.820
-Belper 0.8801 0.013732  64.09  0.4 0.3333  99.039
-Empnig 0.4536 0.010325  43.93  0.0 0.3333  99.668
-Viclat 0.2399 0.005461  43.93  0.0 0.1667 100.000
-Cirarv 0.0000 0.000000    NaN  0.0 0.0000 100.000
-Chealb 0.0000 0.000000    NaN  0.0 0.0000 100.000
+          contr       sd  ratio av.a   av.b  cumsum
+Poatri 0.071553 0.013681 5.2302  4.8 0.0000 0.09913
+Lolper 0.054533 0.029625 1.8408  4.0 0.3333 0.17468
+Rumace 0.046546 0.030806 1.5109  3.2 0.0000 0.23917
+Poapra 0.041750 0.018852 2.2146  3.4 0.6667 0.29701
+Plalan 0.041633 0.029560 1.4084  3.0 0.8333 0.35469
+Brarut 0.035340 0.020104 1.7579  2.8 2.8333 0.40365
+Elepal 0.032043 0.032315 0.9916  0.8 2.1667 0.44805
+Agrsto 0.031915 0.028887 1.1048  1.4 2.1667 0.49227
+Trirep 0.030372 0.022871 1.3280  2.8 1.8333 0.53434
+Elyrep 0.029811 0.038676 0.7708  2.0 0.0000 0.57565
+Antodo 0.028717 0.024799 1.1580  1.8 1.3333 0.61543
+Junart 0.027414 0.028537 0.9607  1.6 1.1667 0.65341
+Tripra 0.025843 0.025972 0.9951  1.8 0.0000 0.68922
+Salrep 0.025340 0.027291 0.9285  0.0 1.8333 0.72432
+Alogen 0.024459 0.032399 0.7549  1.6 0.0000 0.75821
+Leoaut 0.020703 0.014125 1.4658  2.8 3.1667 0.78689
+Ranfla 0.019285 0.019939 0.9672  0.4 1.3333 0.81361
+Junbuf 0.018181 0.024648 0.7376  1.2 0.0000 0.83880
+Hyprad 0.017141 0.026548 0.6457  0.0 1.1667 0.86255
+Calcus 0.016833 0.024901 0.6760  0.0 1.1667 0.88587
+Achmil 0.016555 0.014900 1.1111  1.2 0.3333 0.90881
+Sagpro 0.015282 0.016535 0.9243  0.8 0.5000 0.92998
+Airpra 0.012605 0.018243 0.6910  0.0 0.8333 0.94744
+Brohor 0.012094 0.015169 0.7973  0.8 0.0000 0.96420
+Potpal 0.010105 0.014556 0.6942  0.0 0.6667 0.97820
+Belper 0.008801 0.013732 0.6409  0.4 0.3333 0.99039
+Empnig 0.004536 0.010325 0.4393  0.0 0.3333 0.99668
+Viclat 0.002399 0.005461 0.4393  0.0 0.1667 1.00000
+Cirarv 0.000000 0.000000    NaN  0.0 0.0000 1.00000
+Chealb 0.000000 0.000000    NaN  0.0 0.0000 1.00000
 > 
 > 
 > 



More information about the Vegan-commits mailing list