/usr/lib/R/site-library/vegan/ONEWS is in r-cran-vegan 2.4-6-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453 454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497 498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516 517 518 519 520 521 522 523 524 525 526 527 528 529 530 531 532 533 534 535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565 566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607 608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625 626 627 628 629 630 631 632 633 634 635 636 637 638 639 640 641 642 643 644 645 646 647 648 649 650 651 652 653 654 655 656 657 658 659 660 661 662 663 664 665 666 667 668 669 670 671 672 673 674 675 676 677 678 679 680 681 682 683 684 685 686 687 688 689 690 691 692 693 694 695 696 697 698 699 700 701 702 703 704 705 706 707 708 709 710 711 712 713 714 715 716 717 718 719 720 721 722 723 724 725 726 727 728 729 730 731 732 733 734 735 736 737 738 739 740 741 742 743 744 745 746 747 748 749 750 751 752 753 754 755 756 757 758 759 760 761 762 763 764 765 766 767 768 769 770 771 772 773 774 775 776 777 778 779 780 781 782 783 784 785 786 787 788 789 790 791 792 793 794 795 796 797 798 799 800 801 802 803 804 805 806 807 808 809 810 811 812 813 814 815 816 817 818 819 820 821 822 823 824 825 826 827 828 829 830 831 832 833 834 835 836 837 838 839 840 841 842 843 844 845 846 847 848 849 850 851 852 853 854 855 856 857 858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893 894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941 942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991 992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041 1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091 1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140 1141 1142 1143 1144 1145 1146 1147 1148 1149 1150 1151 1152 1153 1154 1155 1156 1157 1158 1159 1160 1161 1162 1163 1164 1165 1166 1167 1168 1169 1170 1171 1172 1173 1174 1175 1176 1177 1178 1179 1180 1181 1182 1183 1184 1185 1186 1187 1188 1189 1190 1191 1192 1193 1194 1195 1196 1197 1198 1199 1200 1201 1202 1203 1204 1205 1206 1207 1208 1209 1210 1211 1212 1213 1214 1215 1216 1217 1218 1219 1220 1221 1222 1223 1224 1225 1226 1227 1228 1229 1230 1231 1232 1233 1234 1235 1236 1237 1238 1239 1240 1241 1242 1243 1244 1245 1246 1247 1248 1249 1250 1251 1252 1253 1254 1255 1256 1257 1258 1259 1260 1261 1262 1263 1264 1265 1266 1267 1268 1269 1270 1271 1272 1273 1274 1275 1276 1277 1278 1279 1280 1281 1282 1283 1284 1285 1286 1287 1288 1289 1290 1291 1292 1293 1294 1295 1296 1297 1298 1299 1300 1301 1302 1303 1304 1305 1306 1307 1308 1309 1310 1311 1312 1313 1314 1315 1316 1317 1318 1319 1320 1321 1322 1323 1324 1325 1326 1327 1328 1329 1330 1331 1332 1333 1334 1335 1336 1337 1338 1339 1340 1341 1342 1343 1344 1345 1346 1347 1348 1349 1350 1351 1352 1353 1354 1355 1356 1357 1358 1359 1360 1361 1362 1363 1364 1365 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 1393 1394 1395 1396 1397 1398 1399 1400 1401 1402 1403 1404 1405 1406 1407 1408 1409 1410 1411 1412 1413 1414 1415 1416 1417 1418 1419 1420 1421 1422 1423 1424 1425 1426 1427 1428 1429 1430 1431 1432 1433 1434 1435 1436 1437 1438 1439 1440 1441 1442 1443 1444 1445 1446 1447 1448 1449 1450 1451 1452 1453 1454 1455 1456 1457 1458 1459 1460 1461 1462 1463 1464 1465 1466 1467 1468 1469 1470 1471 1472 1473 1474 1475 1476 1477 1478 1479 1480 1481 1482 1483 1484 1485 1486 1487 1488 1489 1490 1491 1492 1493 1494 1495 1496 1497 1498 1499 1500 1501 1502 1503 1504 1505 1506 1507 1508 1509 1510 1511 1512 1513 1514 1515 1516 1517 1518 1519 1520 1521 1522 1523 1524 1525 1526 1527 1528 1529 1530 1531 1532 1533 1534 1535 1536 1537 1538 1539 1540 1541 1542 1543 1544 1545 1546 1547 1548 1549 1550 1551 1552 1553 1554 1555 1556 1557 1558 1559 1560 1561 1562 1563 1564 1565 1566 1567 1568 1569 1570 1571 1572 1573 1574 1575 1576 1577 1578 1579 1580 1581 1582 1583 1584 1585 1586 1587 1588 1589 1590 1591 1592 1593 1594 1595 1596 1597 1598 1599 1600 1601 1602 1603 1604 1605 1606 1607 1608 1609 1610 1611 1612 1613 1614 1615 1616 1617 1618 1619 1620 1621 1622 1623 1624 1625 1626 1627 1628 1629 1630 1631 1632 1633 1634 1635 1636 1637 1638 1639 1640 1641 1642 1643 1644 1645 1646 1647 1648 1649 1650 1651 1652 1653 1654 1655 1656 1657 1658 1659 1660 1661 1662 1663 1664 1665 1666 1667 1668 1669 1670 1671 1672 1673 1674 1675 1676 1677 1678 1679 1680 1681 1682 1683 1684 1685 1686 1687 1688 1689 1690 1691 1692 1693 1694 1695 1696 1697 1698 1699 1700 1701 1702 1703 1704 1705 1706 1707 1708 1709 1710 1711 1712 1713 1714 1715 1716 1717 1718 1719 1720 1721 1722 1723 1724 1725 1726 1727 1728 1729 1730 1731 1732 1733 1734 1735 1736 1737 1738 1739 1740 1741 1742 1743 1744 1745 1746 1747 1748 1749 1750 1751 1752 1753 1754 1755 1756 1757 1758 1759 1760 1761 1762 1763 1764 1765 1766 1767 1768 1769 1770 1771 1772 1773 1774 1775 1776 1777 1778 1779 1780 1781 1782 1783 1784 1785 1786 1787 1788 1789 1790 1791 1792 1793 1794 1795 1796 1797 1798 1799 1800 1801 1802 1803 1804 1805 1806 1807 1808 1809 1810 1811 1812 1813 1814 1815 1816 1817 1818 1819 1820 1821 1822 1823 1824 1825 1826 1827 1828 1829 1830 1831 1832 1833 1834 1835 1836 1837 1838 1839 1840 1841 1842 1843 1844 1845 1846 1847 1848 1849 1850 1851 1852 1853 1854 1855 1856 1857 1858 1859 1860 1861 1862 1863 1864 1865 1866 1867 1868 1869 1870 1871 1872 1873 1874 1875 1876 1877 1878 1879 1880 1881 1882 1883 1884 1885 1886 1887 1888 1889 1890 1891 1892 1893 1894 1895 1896 1897 1898 1899 1900 1901 1902 1903 1904 1905 1906 1907 1908 1909 1910 1911 1912 1913 1914 1915 1916 1917 1918 1919 1920 1921 1922 1923 1924 1925 1926 1927 1928 1929 1930 1931 1932 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 1944 1945 1946 1947 1948 1949 1950 1951 1952 1953 1954 1955 1956 1957 1958 1959 1960 1961 1962 1963 1964 1965 1966 1967 1968 1969 1970 1971 1972 1973 1974 1975 1976 1977 1978 1979 1980 1981 1982 1983 1984 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 2226 2227 2228 2229 2230 2231 2232 2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245 2246 2247 2248 2249 2250 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265 2266 2267 2268 2269 2270 2271 2272 2273 2274 2275 2276 2277 2278 2279 2280 2281 2282 2283 2284 2285 2286 2287 2288 2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323 2324 2325 2326 2327 2328 2329 2330 2331 2332 2333 2334 2335 2336 2337 2338 | -*-Text-*-
VEGAN RELEASE VERSIONS
======================
CHANGES IN VEGAN 1.17-12
- This is a maintenance release which improves robustness of
several functions. A major release is expected soon.
- tolerance.cca: a new function to find the species response
widths (a.k.a. tolerances) and sample heteregeneities from a
cca() result.
- adonis: much faster.
- betadiver: argument 'index' was renamed to 'method' so that the
function can be used similarly as dist(), vegdist() and other
distance functions. This allows using betadiver() as a distance
function in metaMDS().
- cca/rda/capscale: handling of aliased and other zero-rank
components changed in vegan 1.17-11, but not all support
functions were adapted to these changes in that release. Now the
following functions cope with the changes and are more robust:
capscale, anova.cca, bstick.cca, goodness.cca, predict.cca,
screeplot.cca, calibrate.cca, deviance.cca/rda, ordiplot3d,
ordiresids, ordirgl, ordixytplot.
- isomap: text or points are plotted more cleanly, and text uses
ordilabel().
- make.cepnames: a bit more flexible and robust, with a new
argument that allows selecting first and second name component
instead of first and last.
- metaMDSrotate: rotation could be slightly off with more than two
dimensions when there are NA scores.
- ordiellipse, ordihull, ordispider, orditorp: accept NA scores.
- ordilabel: gained argument 'select'.
- ordiplot: better handling of graphical arguments, arg 'cex' can
be set by the user.
- swan: gained new argument 'maxit' that allows setting maximum
number of iterations. Defaults 'maxit = Inf' or the current
behaviour of continuing as long as there are zeros that can be
filled.
CHANGES IN VEGAN 1.17-11
- This is a maintenance release which improves the robustness of
several functions but introduces no new important features.
- cca/rda/capscale: handling of aliased and other zero-rank
components changed. They are no longer made NULL and left
undisplayed, but now they are shown as zero rank. The general
output also changed so that proportions of inertia are only
shown if there are conditions or constraints, but not for
unconstrained analysis. The capscale() function no longer shows
species scores if these are unavailable. Several support
functions more robust.
- nobs: R 2.13.0 introduced generic function nobs() to find the
number of observations, and this release provides nobs() for
several vegan results. (This does not make vegan dependent on R
2.13.0, but vegan works with older R as well.)
- prc: allows anova(..., by = "axis") and other 'by' cases for
prc() results. Some rda() support functions still fail, but now
they stop informatively.
- specaccum: gained argument 'groups' which can be used to find
the number of species in subsets of the data.
CHANGES IN VEGAN 1.17-10
- This is minor revision that mainly fixes vegan with respect to
changes in the currently released R 2.13.0. Most importantly,
cmdscale() output changed in R 2.13.0 and because of this
capscale() could fail in some rare situations with argument 'add
= TRUE'. This vegan bug made BiodiversityR package fail its
tests in R 2.13.0.
- metaMDSrotate: gained argument na.rm = TRUE.
CHANGES IN VEGAN 1.17-9
- anova of cca/rda/capscale results gave wrong results in partial
models. The bug was introduced in vegan 1.17-7.
- diversity and related functions rarefy, rrarefy and specnumber
now accept vector input. Earlier a single site had to be
analysed either as a single-row matrix or using the non-default
setting MARGIN = 2.
- drarefy: new function that returns a matrix of probabilities
that a species occurs in a rarefied sample of a given size.
- metaMDS: it is possible to supply a starting configuration with
argument 'previous.best'. A previous metaMDS or isoMDS result can
also be given as a starting configuration. If the starting
configuration has a higher number of dimensions than requested,
the extra ones are dropped, and if the starting configuration
has fewer dimensions, random scores for extra dimensions will be
added. This may help in running metaMDS over a range of
dimensionalities.
- metaMDSrotate: can now rotate metaMDS solutions with any number
of dimensions so that the first axis is parallel to a fitted
environmental vector. Previously, only two dimensional solutions
worked.
- ordilabel: gained argument 'xpd' that allows labels outside the
plot area. This allows labels above axes, for instance.
- ordisurf: gained several new arguments to control the mgcv::gam
fitting. Also gained an argument to suppress plotting, and a new
plot method. The fitted model can be specified with a formula.
- prestonfit and friends: default is now 'tiesplit = TRUE' (which
was a new feature introduced in vegan 1.17-8).
CHANGES IN VEGAN 1.17-8
- prestonfit: fixed a bug in as.preston(): the largest octave
could be missing with 'tiesplit = TRUE'.
- decorana: spurious axes scores and eigenvalues could be reported
when the eigenvalues actually were zero. This was rarely a
problem with real data, but occurred only in arbitrary examples.
- procrustes: checks input.
CHANGES IN VEGAN 1.17-7
- anova.cca: more robust when models were fitted without 'data'
argument, or when 'na.action' or 'subset' was used.
- capscale: implemented 'subset' in model definition. Additive
constant with 'add = TRUE' is taken into account in predict()
and fitted(). Implemented simulate() which returns a
dissimilarity matrix with random error about predicted values.
- eigenvals: can now extract eigenvalues of some ade4 and labdsv
result objects.
- nestednodf: did not use binary data when weighted = FALSE was
used together with order = FALSE. Reported by Daniel Spitale.
- prestonfit: gained option to split tied frequencies (1, 2, 4, 8,
etc.) between adjacent octaves.
- specaccum: implemented choice of using either "individuals" or
"sites" as x-axis in plot(). Corrected a typo in the result
object: now returns indeed "individuals".
CHANGES IN VEGAN 1.17-6
- capscale: vegan 1.17-5 defined total inertia as the sum of
absolute values of eigenvalues which changed the total inertia
from previous versions. This is changed so that the total
inertia is the sum of all eigenvalues, i.e., negative
eigenvalues are subtracted from the total which was the
definition in 1.17-4 and earlier. The proportions of inertia
components are now expressed for non-negative eigenvalues only,
and a new item "Real Total" (sum of positive eigenvalues) is
added to the output if there are negative eigenvalues
("Imaginary" component). The function no longer returns the
negative eigenvalues.
- CCorA: scaling of object scores changed: they are no longer
divided with sqrt(n-1). Minor fixes for improved robustness.
- metaMDS: automatically changes some arguments to non-default
values if input data contains negative entries. Help page gives
advice to do so for non-community data.
- wascores: checks that input weights are non-negative.
NEW FEATURES AND BUG FIXES IN VEGAN 1.17-5
- anova.cca: Empty models with no constrained component have
correct degrees of freedom. Tests 'by = "axis"' are really
marginal for all axes (first axis used to be non-marginal to
latter axes).
- as.mlm.cca/rda, intersetcor, vif.cca: avoid bug in qr.X in R
2.12.0 which use wrong variable names in aliased models (fixed
also in R 2.12.1).
- betadisper: plot did not allow selecting axes. Reported by Sarah
Goslee.
- capscale: failed with NA action (reported by Nevil Amos). Total
inertia is defined as the sum of absolute values of all
eigenvalues consistently with cmdscale (from R 2.12.1) and other
vegan functions.
- cca/rda/capscale: print proportions of inertia components.
- decorana: Fortran code prints warnings of failed convergence.
Default plot() uses points if items have no name labels to be
plotted.
- deviance.cca/rda: return 0 (instead of NA) for over-paremetrized
models with no unconstrained component.
- eigenvals: uses sum of absolute values of eigenvalues in models
that may have negative eigenvalues (cmdscale, wcmdscale,
capscale). Can extract eigenvalues of pcnm() and cmdscale() of R
2.12.1. Uses zapsmall() to print near-zero
eigenvalues.
- mantel.correlog: P-values could slightly off. The function uses
internally mantel() to evaluate the statistic and its
significance. The mantel() function adds one both to the
denominator and numerator, but mantel.correlog() did not notice
this and made the addition for the second time. When the
mantel.correlog() reported a positive value of the statistic, it
reversed the direction of the one-sided test of mantel(), but
did not handle tied values correctly in this reversal.
- nestednodf: Gustavo Carvalho's new version that also implements
a new quantitative method of Almeida-Neto & Ulrich (Env Mod
Software, 26, 173-178; 2011).
- oecosimu: takes care that the statistic is evaluated with binary
data when null models are binary (regression introduced in vegan
1.17-0). Can handle NA values in permutations (as na.rm = TRUE).
- ordilabel: gained new argument 'col' to set the text colour of
the labels separately from 'border'. Function ordiellipse() uses
this in filled polygons to set the colour of labels to that of
borders or foreground instead of the colour of the fill.
- ordiR2step: new function to model selection in rda() and
capscale() based on adjusted R2 following the recommendations of
Blanchet, Legendre & Borcard (Ecology 89, 2623-2632; 2008).
- ordistep: correct a name clash when the fitted model had term or
item called 'mod'. Reported by Richard Telford (Bergen, Norway).
- pcnm: returns the truncated distanced matrix with argument
'dist.ret'.
- prc: rewritten by Cajo ter Braak with a fix to scaling of
coefficients for full compatibility with Canoco.
- procrustes: Fixed translation following the report and fix by
Christian Dudel (Bochum, Germany). New predict() method to
rotate 'newdata' points to the target. Graphics with plot()
select the direction of arrows with new argument 'to.target',
and use ordilabel() to display the labels with 'type = "t"'.
- rankindex: can take a list of dissimilarity functions as the
argument to allow the use of other dissimilarity functions than
vegdist() of vegan.
- rda/capscale: the scores() can take two scaling constants
('const'). One is used for species, second for WA and LC scores
of sites. This allows using scalings of prcomp(), princomp() or
Canoco (like documented in the design decisions vignette).
- swan: disconnected data caused infinite loop. Now zeros are left
to the result with disconnected data.
- treedist: can handle empty or one-node trees. Gained new
argument 'relative' (defaults TRUE): if FALSE, finds raw
dissimilarities of tree heights. Help page tells that relative
tree distances are in range 0..2 instead of 0..1, since
combining two trees may add a new common root.
- treedive: handles trivial cases of zero (diversity NA) or one
species (diversity 0).
- vegdist: help page gives now binary variants of the indices
using designdist() notation.
- vif.cca: can handle models with aliased terms.
- wcmdscale: returns goodness of fit statistic (item GOF) and
handles negative eigenvalues consistently with cmdscale() and
capscale().
NEW FEATURES AND BUG FIXES IN VEGAN 1.17-4
- MOStest: a new set of functions to implement generalized
Mitchell-Olds & Shaw test for the location of the quadratic
extreme in a given interval. The test can also be used for the
location of the optimum of the Gaussian response model. In
addition to the basic test, there are Fieller and deviance
profile methods for the confidence interval of the location of
the quadratic extreme or Gaussian optimum.
- mrpp & meandist: mrpp() does not evaluate the Classification
Strength (CS) of Van Sickle & Hughes (J. N. Am. Benthol. Soc.,
19: 370-384; 2000) any longer. The old version misinterpreted
the weighting scheme of CS, and with correct scheme there is no
exact relationship between CS and the corresponding MRPP
statistic and therefore permutation tests are not available. CS
is still evaluated in meandist(), but with corrected weighting
scheme, and no tests are performed. Reported by Dr John Van
Sickle (Corvallis OR).
- msoplot: legend text was in a wrong order. Reported by Daniel
Borcard.
- vegdist: Marti Anderson variant of the Gower distance (method =
"altGower") should be without range standardization. Reported by
Sergio Garcia.
- ordiellipse: labels disappeared in their background with 'draw =
"polygon"'.
- predict.cca and predict.rda: gained argument 'newdata' for 'type
= "response"' and 'type = "working"'. These return estimates
of the response data, or dissimilarities in capscale().
NEW FEATURES AND BUG FIXES IN VEGAN 1.17-3
- adonis: handles ties correctly.
- anova.cca and permutest.cca: permutations of cca result is
faster and now nearly equally fast (or slow) in rda and cca.
- betadisper: 'type = "median"' (the default) was not computing
the spatial median on the real and imaginary axes separately.
Reported by Marek Omelka.
- cca, rda and capscale failed when Condition() was a factor, but
the constraints had only continuous variables.
- envfit: defaults to use 999 permutations instead of skipping
permutations.
- mantel and mantel.partial: faster permutation.
- mantel.correlog: upgraded and faster.
- ordiarrows, ordisegments and ordispider gained argument 'label'
to label the groups corresponding to drawn arrows or lines.
- ordiresids: de-weights residuals and fitted values in CCA so
that these are identical to the values shuffled in simulate.cca.
- permutest.cca: adds observed statistic among permutations when
printing the result (does not influence anova.cca or
calculations).
- RsquareAdj: The 'rda' method used wrong number of degrees of
freedom in rank deficit models (number of dependent variables
was lower than the rank of constraints and conditions). Default
method handles vector input.
- scores: default method works with one row of scores. The scores
of "cca" and "rda" methods always have names, even if there are
no names in the input data.
- Permutation tests: The available permutation tests are described
in a new help files accessed via ?permutations.
NEW FEATURES IN VEGAN 1.17-2
- permutest.betadisper: printed P-values in a wrong order in
displaying the the matrix of pairwise tests (the values were
correct, but formatting failed). Reported by Dan O'Shea.
- nesteddisc: failed if the most species rich sites were tied. The
function is now much faster (though still slow), but the price
is that it does not try as hard to find the optimal ordering.
- screeplot.cca & friends: new argument 'legend' for all methods
(except 'decorana') to draw a legend if the observed and broken
stick distribution are both plotted.
- ordistep: adds an 'anova' item to the final model similarly as
the standard step(). You can suppress the tracing ('trace =
FALSE'), and find the model build history save in 'anova' in a
compact form.
- densityplot.oecosimu: gave warnings when there was only one
statistic and hence one lattice panel.
- predict functions for cca and rda objects match 'newdata' by
dimension names instead of index.
- simulate.rda and simulate.cca have new argument 'rank' which
allows using lower rank presentation of fitted values
(including rank = 0).
- treedist: new function to find dissimilarities of species
property trees of communities. The property trees can be, e.g.,
functional diversity trees, taxonomies or phylogenies. Similar
in spirit to UniFrac distance (C. Lozupone & R. Knight, Appl
Environ Microbiol 71:8225-8235; 2005), but completely different
in design and works only with binary data.
NEW FEATURES AND FIXES IN VEGAN 1.17-1
- multipart: new functions for multiplicative partionining of
gamma diversity into alpha and beta diversity components.
- CCorA: Fixed choice of scores in biplots -- Used PC scores
instead of observed scores in the right-hand-side biplot
panel. The biplot function has several new and enhanched
alternatives of plots.
- envfit: ignored weights in cca() results for factors or a single
continuous variable. The bug was introduced with NA handling
upgrade in vegan 1.17-0. The problem with fitted vectors was
reported by Richard Telford (Bergen, Norway).
CHANGES AND NEW FEATURES IN VEGAN 1.17-0
- Guillaume Blanchet joined the vegan team.
- New function to partition data-set diversity (gamma) into
within-plot (alpha) and between-plot (beta) diversity
components. Function adipart performs additive
partitioning (gamma = alpha + beta). Function hiersimu performs
hierarchical null model testing similar to adipart but by using
custom function to calculate statistics for levels of a hierarchy.
- Subsets and missing value handling added to constrained
ordination methods cca(), rda() and capscale(). The missing
values (NA) can be handled with setting na.action (defaults
na.fail). With na.action = na.omit, observations with missing
values are removed, and with na.action = na.exclude they are
kept, but scores may be NA. However, the WA scores for sites are
available in non-partial models with na.action = na.exclude.
The 'subset' can be defined using any variable in the
constraining data set or species in the dependent data.
- Functions for fitting environmental variables onto ordination
(envfit, ordisurf) are aware of missing values in constraints or
NA values in scores. The ordination plot functions also can
handle NA scores.
- New and upgraded quantitative null models. In particular,
quantitative swap models (function abuswap) allows generetating
null matrices where marginal totals and fill are fixed, or row
and column fills are fixed, or row and column fills and either
row or column totals are fixed.
- oecosimu: rewritten to handle quantitative null models. Gained
keyword 'alternative' for "two-sided", or "less" and "greater"
one-sided tests. More robust with degenerate solutions. The
'method' can now be a user-supplied function. New support
functions as.ts() and as.mcmc which transform sequential models
into form that can be analysed using tools for time-series
(as.ts()) or MCMC sequences of the 'coda' package (as.mcmc).
- calibrate: calibrate is now a generic function with a new method
calibrate.ordisurf() in addition to the old calibrate.cca. These
find the estimates of environmental variables from ordination.
- ordistep: stepwise selection of terms in constrained ordination
(cca, rda and capscale) using permutation tests instead of
pseudo-AIC that is used by the standard step() function.
- pcnm: new function to find weighted principal coordinates of
neighbour matrix (PCNM) from distances between points. These are
typically used for spatial filtering in constrained
ordination. The function uses weighted analysis and can
therefore produce PCNM for correspondence analysis in addition
to PCA and RDA.
- betadisper: can use spatial medians which are now the default
method. Preliminary tests indicate that spatial medians correct
the anti-conservative Type I errors reported by Stewart
Schultz.
- decostand & vegdist: new transformation 'log' which implements
Marti Anderson's scaling log(x) + 1 for x>0 (which is not at all
the same as log(x+1)), and vegdist has alternative Gower
function that skips double zeros. Together these implement the
"modified Gower" distance of Anderson et al. (Ecology Letters 9,
683-693; 2006). Feature request #473 by Etienne Laliberte in
vegan.r-forge.r-project.org.
- model.matrix.cca & model.frame.cca: new functions to reconstruct
the model frame and model matrix (model matrices in partial
models) of constraints used in ordination methods cca(), rda()
and capscale().
- simulate.cca & simulate.cca: simulate the community (response)
data from the result of cca() or rda() under alternative
hypothesis. Error is added to the fitted values from
ordination. The function uses either Gaussian error or permutes
residuals and adds these to fitted values.
- spandepth: a new function to find the depth of each node in the
minimum spanning tree produced by spantree(). Feature request by
W. E. Sharp.
- alias.cca: gained argument names.only.
- metaMDSrotate: a new function to rotate metaMDS so that first
axis is parallel to an environmental variable.
- msoplot: uses standard legend.
- nesteddisc: new method that tries to find the smallest possible
value of the statistic in tied data. Slow, but fixes the
problems of the published method. The problem with tied values
in nesteddisc was found with Carsten Dormann.
- contribdiv: gained a plot method.
- rarefy: failed with one site and many sample sizes like
rarefy(rpois(10, 2), sample=2:3).
- meandist: plot can draw histograms as an alternative to
dendrograms.
- plot functions for 'decorana' and 'cca' and friends and
'ordiplot' use 'linestack' if only one axis was requested.
CHANGES IN VEGAN 1.15-4
- Changed package dependence: vegan does not depend on 'ellipse'.
- anosim: user interface identical to 'mrpp'. Accepts now data
matrix and finds dissimilarities internally.
- betadisper: fix removal of zero eigenvalues with non-Euclidean
distances. This may change the results slightly, but in most
cases the effects are minor or none.
- capscale: has now 'fitted' and 'residual' methods, and 'predict'
works with 'type = "response"'. These return dissimilarities
that produce the same ordination as the original data.
- indpower: new function to find the indicator power of species to
predict presence of other species (Halme et al., Conservation
Biology 23, 1008-1016; 2009). Closely related to the the power
to predict probabilities in beals().
- mantel.correlog: new functions to produce multivariate Mantel
correlograms (Legendre & Legendre, Numerical Ecology, section
13.1.5; 1998).
- metaMDS: accepts user-supplied dissimilarities. Species scores,
data transformation, step-across and half-change scaling are
unavailable with user-supplied dissimilarities, but random
starts, PC rotation and scaling to original range of input
dissimilarities work.
- nestedtemp: row and column labels can be turned on/off
independently in plots.
- ordihull & ordiellipse: new argument 'label' to plot group names
for hulls or ellipses. Return (invisibly) the data for plotted
convex hulls or ellipses. In 'ordihull' this is a list of hull
vertices, and in 'ordiellipse' a list of (scaled) covariance and
centre data for ellipses. New 'summary' methods find the centres
and areas of plotted hulls or ellipses. Argument 'draw' has new
option "none" that suppresses all drawing so that only the data
for summary can be extracted without plotting.
- orditorp: works with reversed axes, like with 'xlim = c(3, -3)'.
- ordixyplot: has a panel function for arrows: a lattice variant
of ordiarrows.
- poolaccum, estaccumR: new functions to find the specpool() or
estimateR() estimates of extrapolated species richness in random
accumulations of sites. These have 'plot' and 'summary' methods.
- scores: functions biplot, points, text, ordilabel, ordiplot3d,
ordixyplot and spantreee did not pass all arguments to scores()
function. In particular, this concerned rda where 'scaling' and
'const' arguments could not be used within these functions.
- radfit: new function 'radlattice' for a lattice plots of fitted
models for a single site.
- rrarefy: a new function to generate random rarefied communities.
- RsquareAdj: now a generic function to find adjusted R squared
with special cases to 'rda', 'cca', 'lm' and 'glm'.
- summary.cca: cleaner output.
- wcmdscale: returns scaled scores for axes with negative
eigenvalues.
CHANGES IN VEGAN 1.15-3
- anosim, mantel, mrpp, protest, envfit: did not include the
observed statistic among permutations. Functions adonis and
permutest.cca did this correctly, but did not print results
neatly.
- anova.cca: name clash if data were indexed with 'i'.
- capscale: fixed handling of negative eigenvalues with
non-Euclidean distances. The total inertia is the sum of all
eigenvalues so that negative eigenvalues are subtracted from the
total. The total inertia of negative components and their rank
(number) is given as 'Imaginary' component, and the negative
eigenvalues are listed after unconstrained positive eigenvalues.
The procedure is based on Gower, Linear Algebra and its
Applications 67, 81-97 (1985). New argument 'sqrt.dist' takes
square root of the internally calculated dissimilarities and
avoids negative eigenvalues with some indices, such as vegan
Jaccard and Bray-Curtis. The adjustment is corrected for indices
with upper limit of one, and the reported eigenvalues and
inertia components are reduced by a factor of sqrt(n-1) and are
similar to those reported by 'cmdscale' or 'wcmdscale'.
- eigenvals: a new function to extract eigenvalues from rda, cca,
capscale, wcmdscale, prcomp, princomp, svd or eigen. If the
result object contains squareroots of eigenvalues, they are
squared. The summary method also finds proportions and
cumulative proportions explained. Function summary.cca now uses
this to display eigenvalues.
- kendall.global: could get wrong counts of ties in large data
sets.
- meandist: new sister function to 'mrpp' that finds mean
within-group and between-group dissimilarities. The summary
function finds overall averages of these, and returns all three
MRPP variants plus classification strength. The plot method
draws a dendrogram based on the mean dissimilarity matrix, with
leaves hanging to within-groups dissimilarity. The functions
follow Sickle, J. Agric. Biol. Envir. Stat. 2, 370-388 (1997).
- ordisurf: fits now linear or quadratic trend surfaces if
'knots' argument is set to 0, 1, or 2.
- orditkplot: copes with NA and NaN scores.
- ordixyplot: mixed x and y axes for biplot arrows and class
centroids. Function failed in constrained ordination.
- tsallis: gained argument 'hill' to find results analogous to
Hill numbers in renyi(), or the species number equivalents of
indices.
- wcmdscale: removes now zero eigenvalues instead of the last
eigenvalue. The bug was copied from cmdscale(), which still has
the bug in R 2.9.0 (plus another that was not copied to
wcmdscale).
CHANGES IN VEGAN 1.15-2
- adonis: adds one to numerator and denominator of permutation
tests, the default number of iterations was raised from 5 to
999, and the result object got a 'terms' component. Uses much
less memory allowing analysis of larger problems.
- anova.cca: anova(..., by = "axis") gained new keyword 'cutoff'
to stop permutation tests after exceeding the given cut off
level of significance. The second term of anova(..., by =
"margin") used different random numbers than other terms.
- beals: Completely rewritten by Miquel de Caceres. Knows now
also the cross validated version of Beals smoothing and other
choices described by De Caceres & Legendre, Oecologia 156,
657-669; 2008.
- betadisper: handles missing values both in dissimilarities and
grouping.
- cca/rda: cleaner output of summary() in unconstrained models.
- commsimulator: simulated result retains original row and column
names.
- contribdiv: new function for contribution diversity (Lu et al.,
Basic and Applied Ecology 8, 1-12; 2007).
- decostand: gained dots in argument list, so that
stressplot(metaMDS(x, dist = "gower", trymax = 40)) works.
- dispindmorisita: new function for the Morisita index of
dispersion. See Krebs, Ecological Methodology; 1999.
- kendall.global and kendall.post: new functions for Kendall's
coefficient of concordance. In ecology these can be used to
identify significant species associations (Legendre, J Agric
Biol Environm Stat 10, 226-245; 2005).
- metaMDS: more robust with distances like Euclidean and Manhattan
which have no upper limit. The stepacross works correctly for
these, but gives a warning that its use may not be sensible.
There is a better heuristics to avoid half-change scaling with
these indices. The 'halfchange' argument is now honoured when
given in metaMDS() call.
- mrpp: returns within-class dissimilarities. Evaluates
classification strength (Van Sickle, Biological and
Environmental Statistics, 2, 370-388; 1997) if weight.type = 3
was used.
- nestednodf: new nestedness function of overlap and decreasing
fill (Almeida-Neto et al., Oikos 117, 1227-1239; 2008). Coding
by Gustava Carvalho.
- ordirgl: no more superfluous warnings with type = "t".
- ordisurf: gained an argument 'bubble' to draw bubble plots, or
vary point sizes according the value of the observed variable.
The (invisible) return object now has an item 'grid' for the
fitted value and grid values.
- orditkplot: can produce TIFF graphics, if installed R has this
capability.
- procrustes: new text() function.
- radfit: broken-stick model (function rad.null) failed with
Gaussian and Gamma error families. The plot.radfit command
gains argument log = "y", which allows using arithmetic scales
or log-log scales where Zipf model is a straight line.
- spantree: retains names and uses labels in plot.
- datasets: Oribatid mites (data 'mite') got taxon names.
- help files updated, and do not raise errors or warnings with the
Rd parser version 2 in R 2.9.0.
CHANGES IN VEGAN 1.15-1
- betadisper and related functions now work when the 'group'
argument represents a single group/level.
- decorana: 'text' function failed.
- cca/rda/capscale: cleaner handling of results with missing
(NULL) row or column names, also in plots.
- cca/rda/capscale: new functions 'head' and 'tail' for 'summary'
of cca results to show only some rows of the scores.
- nestedness: method nesteddisc was found to be strongly dependent
on the ordering of tied column (species) frequencies. A warning
is now issued, and development version of vegan in
http://vegan.r-forge.r-project.org/ has an experimental function
that tries to find the optimal ordering.
- permatfull: 'row'/'column' arguments were mixed.
- radfit: all methods work consistently for communities of 0, 1 or
2 species.
- wcmdscale: typo in the result object items corrected.
CHANGES IN VEGAN 1.15-0
GENERAL
- Peter Solymos joined the vegan development team.
NEW FUNCTIONS
- add1.cca and drop1.cca: functions to implement permutation tests
with argument test = "permutation". Function drop1.cca uses
anova.cca(..., by = "margin") and add1.cca implements a new type
of analysis for single term additions.
- ordilabel: new alternative for cluttered ordination plots. Text
is written on opaque labels. The texts still cover each other,
but at least the uppermost are readable, and the ordering can be
controlled by 'priority' argument. Similar to 's.label' in ade4
package.
- ordipointlabel: new alternative for cluttered ordination plots.
Points are in fixed positions, but their text labels are
located to avoid overplotting. The optimization uses simulated
annealing of 'optim' function. Returns an "orditkplot" object so
that the results can be edited with 'orditkplot'. Similar to
'pointLabel' function in maptools package.
- permatfull, permatswap: functions to permute quantitative count
data.
- treedive: estimation of functional diversity defined as the
height of the dendrogram of species properties (Petchey & Gaston,
Ecology Letters 9, 741-758; 2006). With a helper function
'treeheight' to find the height of a 'hclust' dendrogram.
- tsallis: Tsallis entropy family as an alternative to 'renyi'.
- wcmdscale: weighted metric scaling a.k.a. weighted principal
coordinates analysis; only uses row weights.
NEW FEATURES AND FIXES
- anova.cca: default permutation model changed from "direct" to
"reduced" after Pierre Legendre demonstrated that model
"reduced" was superior in term of Type I error.
- anova.cca: handles smoothly models where constrained or
unconstrained models are NULL and the tests are impossible (used
to stop with error).
- anova.cca(..., by = "margin") was handling wrongly 'x' in
'~Condition(x) + x + z', or model formulae where same variables
where used both as conditions and (aliased) constraints.
- cca/rda/capscale: improved robustness in the formula
interface. Partial handling of "cca" objects produced by
ade4:::cca.
- commsimulator: swap, trial swap and quasiswap written in C and
*much* faster (100x in some tests).
- oecosimu: accepts now a vector of statistics, and the user can
give the name of the statistic in the function call.
- ordiplot: did not use partial matching for "sites" and
"species".
- orditkplot: Improved user interface. Improved zooming into
graphs. Imitates R plotting characters ('pch'). Label font
family, size and font type can be vectors.
- permuted.index2 and associated functions allow for restricted
permutations of strata (i.e., restricted shuffling of blocks).
- specaccum: removes missing (all zero) species which gave sd = NA
with method = "exact".
CHANGES IN VEGAN 1.13-2
- anova.cca (and permutest.cca) now calculate the residual df as
(number of rows) - (rank of constraints and conditions) - 1
instead of using the rank of the residual ordination. With
single response variable in 'rda' the degrees of freedom and
F-values are now identical to those in 'lm' (linear
model). However, the change does not influence significances
(P-values), because these have been always found by permutation,
and this change does not influence the order statistic used in
permutations.
- oecosimu: always estimates the original statistic with binary
data even when the user supplies quantitative data.
- orditkplot: no superfluous pointer lines when moving labels.
- procrustes plot failed if two configurations were exactly
identical.
- rda (and capscale) use internal scaling constant so that the
returned site and species scores with scalings 1, 2 or 3
together provide a biplot approximation of the original
data. This scaling constant is calculated internally, but now
its numerical value is returned as an attribute of 'scores.rda',
and 'summary' displays its value. Vignette on "Design decisions
and implementation" in vegan explains the calculation of the
internal scaling constant.
- summary of capscale identical to the summary of rda.
- scores of rda and capscale with scaling = 0 really return the
unmodified scores without the scaling constant.
CHANGES IN VEGAN 1.13-1
GENERAL
- Helene Wagner joined the vegan team.
FIXES FOR HANDLING RANKS IN CONSTRAINED ORDINATION
- cca, rda and capscale had only known one kind of rank: the rank
of the ordination result (= number of axes). In fact there are
two other types of ranks: the rank of constraints after removing
conditions, and the rank of constraints + conditions, where
conditions refer to terms "partialled out". This hit those that
tried to use rda with single response variable instead of
community matrix, and also some partial models were handled
wrongly. The changes mainly concern cases where the rank of
constraints is higher than the rank of the ordination (more
constraints than ordination axes). The changes have visible
effects in following support functions:
- alias: no superfluous aliasing of terms, partial models aliased
correctly.
- anova: always uses rank of the constraints which also fixes
anova(..., by = "terms") when single response variable was used.
- calibrate.cca: identifies cases when the rank of constraints is
higher than the rank of ordination and refuses to analyse these.
- extractAIC: uses rank of constraints for degrees of freedom.
- predict: predict(..., type = "lc", newdata = somedata) works in
partial analysis.
OTHER FIXES
- adonis: faster and improved documentation.
- intersetcor and inertcomp now check that input really is from
constrained ordination instead of giving obscure error messages.
- lines.spantree knows again graphical arguments such as 'col',
'lty'.
- mite.xy: new data set on the spatial coordinates of sample sites
of Oribatid mites.
- mso: plot.mso replaced with msoplot, and plot function now
displays the ordination scatter plot. Function msoplot collapses
distances larger than half of the maximum distance into a single
distance class. Printed result shows the variogram data.
- ordicluster, ordiellipse, ordispider, orglspider, ordisurf,
factorfit and vectorfit could fail with non-vegan ordination
objects in R 2.7.0 whose weights.default function gave error if
object had no weights (these worked in R 2.6.2 and earlier).
- taxa2dist: has a "method" name for distances.
CHANGES IN VEGAN 1.13-0
GENERAL
- Based on development version 1.12-15 (revision 354 at
http://vegan.r-forge.r-project.org)
NEW FUNCTIONS
- betadiver: beta diversity functions as reviewed by Patricia
Koleff et al. (J. Anim. Ecol., 72, 367-382; 2003), with a plot
function to produce triangular plots.
- mso: Helene Wagner's multiscale ordination or spatial
partitioning of cca and rda. This is taken from the Ecological
Archives with minimal edition with the permission of Helene
Wagner.
- nestedtemp: matrix temperature method for oecosimu following
Rodriguez-Girones & Santamaria (J. Biogeogr. 33, 924-935; 2006),
but still without iterative optimization of row and column
ordering.
- TukeyHSD.betadisper: pairwise comparisons for betadisper.
NEW FEATURES AND FIXES
- adonis: returns both species and site scores.
- betadisper: was not calculating distance to centroid correctly
for observations where the imaginary distance to centroid was
greater than the real distance. Now takes the absolute value of
the combined distance before taking the square root. This is
in-line with Marti Anderson's PERMDISP2.
- BCI: example has spatial coordinates of plots.
- biplot.rda: argument for arrow col.
- capscale: accepts other distance functions than vegdist, and can
use metaMDSdist for extended dissimilarites & "metaPCoA".
- designdist: knows 2x2 contingency table notation with a, b, c,
and d.
- metaMDS: transforms data like with distances for species WA.
- orditkplot: allows editing labels and zooming into plot.
- permDisper: renamed to permutest.betadisper.
- permutest: now a generic function, currently with methods
'cca' and 'permDisper'.
- rarefy: accepts vectors of sample sizes.
- rgl.isomap: dynamic rgl plots for isomap.
- screeplot.cca etc: return invisibly xycoords.
- specaccum: returns numbers of individuals with method = "rarefaction".
- summary.cca: returns more statistics on "variance explained by axes".
- zzz: vegan got startup message.
CHANGES IN VEGAN VERSION 1.11-4
- A critical bug fix in adonis: there was a critical bug in adonis
code, and inconsistent statistics were used in permutations so
that P-values were grossly wrong in multi-variables models
(single variable models were OK). In addition, df were wrong in
deficit rank models, and unused factor levels were not dropped.
All adonis users should upgrade and rerun their analyses.
CHANGES IN VEGAN VERSION 1.11-3
GENERAL
- Bug fixes from Rev. 305 on http://r-forge.r-project.org/.
FIXES
- anova.cca: number of permutations could exceed perm.max.
- permuted.index2: updated to the version in devel branch. This
means bugfixes in permCheck, numPerms and permuted.index2. It
also adds allPerms to get all possible permutations when
complete enumeration is feasible. In addition, there is now a
function permuplot to graphically show the current permutation
scheme.
- plot.cca, plot.envfit and associates: automatic scaling of
biplot arrows and fitted vectors was wrong when axes were
reversed (like 'xlim = c(1,-1)') or the origin was shifted in
plot.envfit (like 'at = c(1,1)'). Added internal function
'ordiArrowMul'.
- plot.procrustes: failed if two configurations were identical.
- varpart4: sum of squares was wrong if called directly instead of
being called via varpart(). Reported by Guillaume Blanchet.
CHANGES IN VEGAN VERSION 1.11-2
GENERAL
- minor bug fixes and documentation updates from the devel trunk.
- version 1.11-1 was made but never released to CRAN, and this is
the first minor release of the 1.11 series.
FIXES AND UPDATES
- bstick.princomp: works now.
- numPerms: was returning incorrect number of permutations
when type = "strata" selected.
- permuted.index2: was permuting samples within levels of strata
as well as permuting the levels themselves if type = "strata"
selected.
- Documentation: diversity-vegan gains discussion on taxonomic
diversity and using designdist for analysing beta diversity.
Proof-reading and updates in diversity-vegan and FAQ-vegan.
CHANGES IN VEGAN VERSION 1.11-0
GENERAL
- Based on devel version 1.10-13 (rev. 205 at R-Forge).
- Gavin Simpson joined the vegan team.
- Suggests now 'tcltk' (for orditkplot).
NEW FUNCTIONS
- anova.cca gained a new support function to analyse marginal
effects of individual terms (which are similar to Type III
effects). Defined with argument 'by = "margin"'.
- betadisper: new functions for Marti Anderson's analysis of
homogeneity of multivariate dispersions.
- biplot.rda: biplot function for PCA run with rda. Arrows are
used instead of points.
- CCorA: Canonical correlation analysis with a robust algorithm,
with permutation test and plot function.
- oecosimu: functions to analyse nestedness of communities (such
as on islands or patches). Function oecosimu is a general
wrapper, and commsimulator generates null-communities of various
types (r00, r0, r1, r2, c0, swap, trial swap, backtracking,
quasiswap). The nestedness can be analysed with functions like
nestedchecker (number of checkerboard units), nestedn0 (measure
N0), nesteddisc (discrepancy), but users can supply their own
functions or even use some standard R functions such as
chisq.test.
- ordiresids: similar diagnostic plots as in plot.lm for
constrained ordination: Residuals ~ Fitted, sqrt(abs(Residuals))
~ Fitted, and qqmath(~ Residuals) using Lattice graphics.
- orditkplot: interactive and editable plotting function. Function
displays one set of points (species, sites) using both points
and labels (text). The points are fixed, but labels can be
dragged to better positions with mouse. The edited plots can be
saved as EPS, exported to various graphical formats (EPS, PDF,
PNG, JPEG, BMP, XFIG depending on the system) or dumped back to
the R session for plotting and further processing.
- ordixyplot: a set of functions for Lattice graphics of
ordination results. Includes ordixyplot for 2D graphics,
ordisplom for pairs plots, and ordicloud for 3D graphics. All
can be subsetted and formatted in the usual Lattice way.
- permuted.index2: New version of permuted.index() that now allows
restricted permutations. Can produce permutations for
time-series or line transects and for spatial grids. These can
also be nested within 'strata'. permuted.series() and
permuted.grid() are the relevant workhorse
functions. Permutation options are set by new function
permControl(). Currently used only in betadisper, but we plan to
migrate vegan functions to permuted.index2() in the devel
version, and will eventually replace the current
permuted.index(). With support function permCheck for checking
permutation schemes.
NEW DATA SETS
- sipoo: birds in the Sipoo archipelago (Finland, too close to
Helsinki).
NEW FEATURES AND FIXES
- adonis: accepts any 'dist' object as input.
- as.mlm.cca, as.mlm.rda use now correct names for variables in
aliased models. The data were pivoted correctly in R, but the
labels were not.
- anova.cca assesses now P value as (hits+1)/(tries+1).
- anova.cca: anova(..., by = "axis") failed when fitted model had
terms like poly(x,2) or log(x).
- bgdispersal uses now a more powerful statistic for the McNemar
test (in terms of Type 1 error rate).
- calibrate.cca does correct pivoting in aliased models.
- capscale: negative scaling in plot works similarly as in rda.
- decorana does not crash R when called with NULL row data, such
as decorana(dune[FALSE,]). Method predict(..., type="sites")
works correctly with downweighted analysis.
- fitted.cca, fitted.rda gained argument type = "working" to get
the fitted values and residuals used internally in calculation
(in cca() these are weigthed and Chi-square standardized
values).
- isomap checks that input data are dissimilarities or can be
changed into dissimilarities without warnings.
- metaMDS gains argument wascores (defaults TRUE) to suppress
calculation of species scores.
- orditorp now handles "..." more cleanly.
- scores.cca, scores.rda accept display = c("species", "sites").
- summary.prc honours argument 'axis'.
- taxa2dist issues now a warning if called with 'check = FALSE'
and some distances == 0, typically meaning that basal taxa
(species) were not coded.
- varpart failed if there were unused levels in factors.
- wascores returns now NA for missing (all zero) species instead
of failing.
DOCUMENTS
- new documents: FAQ, a simple introduction to ordination in
vegan, a detailed explanation of diversity methods. New
formatting.
- Added these NEWS.
OLDER VEGAN VERSIONS
Version 1.8-8 (Oct 2, 2007)
* Minor bugfix release for upcoming R-2.6.0. Based on the
http://r-forge.r-project.org/projects/vegan/ revision 17 (= 1.8-7)
with ported bug fix revisions (see below for revision numbers).
* anova.cca: by = "term" failed in partial model. This was broken
in 1.8-6 by introducing a test against deficit rank models (r47).
* cascadeKM: Calinski index works now when the input data is a
data.frame (r57:58).
* flush.console: metaMDSiter and bioenv use now flush.console()
so that Windows people also see the trace (r56).
* ordispantree: made defunct, was deprecated in 1.8-1 (r38).
* scores: handles now numeric data frames (r25).
* summary.cca: failed if only one type of scores was requested
(r50, 52).
* taxondive: Fixed dim checking and matching species names in
community data and taxonomic distance data (r21).
* tweaks to pass --pedantic R CMD check, mainly in formatting
source files, unused variables in source files and superfluous
braces in help files (r39, r46, r62).
* Updated FAQ-vegan.pdf to the current version at R-Forge.
Version 1.8-7 (August 24, 2007)
* Based on devel version 1.9-34.
* DESCRIPTION: M. Henry H. Stevens (Miami University, Oxford,
Ohio) joined the vegan team.
* adonis: new function for nonparametric MANOVA that is
appropriate for even extremely wide matrices sometimes associated
with gene data and with diverse ecological communities. Author
Hank Stevens.
* taxondive: a new function for indices of taxonomic diversity and
distinctness after Clarke & Warwick (Mar Ecol Prog Ser 216,
265--278, 2001 and other papers). With a helper function
'taxa2dist' to turn taxonomies into distances with an option for
variable step length (Clarke & Warwick, Mar Ecol Prog Ser 184,
21--29, 1999), and a toy data set on the taxonomy of dune meadow
species ('dune.taxon'). With a help, testing and pressure from
Mike Cappo, James Cook Uni, Qld.
* bgdispersal: previous version was partly garbled (by me), and
P. Legendre provided a corrected one.
* designdist: keeps 'dist' attributes even when the 'method'
function drops them. Swapped the order to (terms, methods) in the
default name.
* metaMDS: issues a warning if data are disconnected. Passes extra
arguments to other 'distfun' than 'vegdist' so that you can set
'terms' in 'designdist' etc (metaMDSdist). Can now do trymax=0 or
skip random starts and give you enhanced 'isoMDS' result
(metaMDSiter).
* ordiplot: failed if number of species was equal to number of
sites (and so did plot.metaMDS and plot.isomap using this).
* plot.profile.fisherfit: corrected a harmless error detected by
checkUsagePackage().
* predict.rda: removed some dead (but heavy) code from type =
"response".
Version 1.8-6 (May 9, 2007)
* Based on devel version 1.9-23.
* as.mlm.cca, as.mlm.rda: new functions to refit constrained
ordination result (cca, rda, capscale) as a multiple response
linear model. You can find influence statistics (Cook's distance,
hat values) from the refitted model. You also can find t-values
etc., but these have the same bias as in other software and should
not be used.
* bgdispersal: a new function for dispersal direction in
biogeography (Legendre & Legendre 1998, section 13.3.4). Author
Pierre Legendre.
* designdist: a new function for defining your own dissimilarity
index or for estimating beta diversity (Koleff et al.,
J. Ecol. 72, 367-382; 2003).
* isomap: a new function for isometric feature mapping of
Tenenbaum et al. (Science 290, 2319-2323; 2000).
* screeplot, bstick: new functions to draw screeplots of vegan
ordination results with brokenstick lines, and alternative
screeplot functions for prcomp and princomp with brokenstick.
Author Gavin L. Simpson.
* swan: a new function for the degree of absence (Swan 1970,
Ecology 51, 89-102).
* anova.cca: now refuses to do 'by = "terms"' if the rank of
constraints is higher than the rank of the community matrix.
* bioenv: gains argument 'partial' to perform partial bioenv.
* cca, rda, capscale: can now handle longer expression within
'Condition()' (ordiParseFormula). Used to drop observations with
missing values in unused variables (ordiGetData).
* goodness.cca, goodness.rda. 'statistic = "distance"' was wrongly
implemented. Now refuses to find "distance" in constrained
analysis: distances of constrained and unconstrained components do
not add up to to distances in unconstrained ordination.
* metaMDS (metaMDSdist): gains argument 'distfun' to use other
dissimilarity functions than vegdist.
* renyiaccum: used a variable that was not defined as an argument.
Added support functions persp.renyiaccum and rgl.renyiaccum
(Roeland Kindt).
* stressplot: R2's renamed to 'non-metric fit' and 'linear fit'.
* Doc: Corrected reference to Hurlbert in diversity (thanks to
Ralph Grundel). Updated references (varpart, renyiaccum). Removed
discussion on t-values in cca from vignettes, because as.mlm.cca
now implements those. General cleanup and better utf-8 encoding.
Version 1.8-5 (January 11, 2007)
* Based on devel version 1.9-12.
* no.shared (manifest in metaMDS): prints thousands of lines of
debugging info that I forgot to deactive in release. Not fatal,
but extremely annoying.
* capscale: inertia name as "unknown" if the dissimilarity object
does not have a "method" name. Suggested by Roeland Kindt.
* DESCRIPTION: license is now explicitly GPL v2 (but not later).
Version 1.8-4 (January 8, 2007)
* Based on devel version 1.9-10.
* cascadeKM: a new function to wrap kmeans and optimality criteria
for classification (Sebastien Durand, Pierre Legendre & Marie
Helene Ouellette).
* renyiaccum: a new function for Renyi (and Hill) accumulation
curves (Roeland Kindt).
* bioenv: bioenv.formula uses now "na.action = NULL" in
'model.frame', and bioenv.default passes arguments to 'cor' which
means that you can set NA treatment in 'cor' using argument "use".
* cca, rda: added "..." to formula versions to satisfy tests in
R-DEVEL.
* cca, rda, capscale: used to fail if called within other
functions. Now data always evaluated in the environment of
formula using new internal function ordiGetData.
* anova.cca: checks that the model has both residual and
constrained components or stops with understandable error message
(used to stop with incomprehensible error message).
* print.summary.cca, print.summary.decorana: have now arguments
'head' and 'tail' to print only a part of species and site
scores. Suggested by Gavin Simpson.
* metaMDS: checks now that the input data ('comm') is not a 'dist'
object (like many users have had).
* ordisurf: Does not depend on package 'akima' any longer, but
directly finds fitted values in a regular grid using
'predict.gam'. Added pnpoly.c to find which of these values are
within the convex hull defined by data. Results also look neater
with sparse data now. Added argument 'labcex' passed to 'contour'
for changing size of contour labels. Setting 'labcex = 0' will
suppress drawing labels (by setting drawlabels = FALSE in
'contour()').
* orditorp: handles now vector arguments of 'col', 'pcol', 'cex',
and 'pch'.
* rad.zipfbrot: less likely to overflow to NA coefficients during
iteration.
* renyi: added a plot function, and documented together with
renyiaccum() instead of diversity().
* scores.default: Knows now about ade4 objects. Primarily looks
for scores scaled by eigenvalues both for sites and species.
* specaccum: Added new conditioned method of Colwell et
al. with estimated sd based on extrapolated richness (Roeland
Kindt).
* vegdist.c: More informative warnings with 'method' name (useful
with rankindex).
* DESCRIPTION: listed 'require()d' packages in the "Suggests:"
field to satisfy more anal tests in R-DEVEL. Mention diversity
analysis in the "Description:".
Version 1.8-3 (Sept 29, 2006)
* Based on devel version 1.9-2.
* varespec.rda, varechem.rda: saved in binary form, because old
ascii form gave warning in R-2.4.0-rc.
* vegdist: added Chao index (of Jaccard type) that should take
into account missing pairs of species. Checks that Binomial index
is non-negative. Identical sites could have dissimilarity of
magnitude 1e-17 after some standardizations, but now values <1e-15
are zapped to zero.
* estimateR: uses now standard unbiased formulation of Chao.
* renyi: should work now (really!).
* metaMDS: with zero = "add", zeros now replaced with
min(dis[dis>0])/2 (used to be 1E-4) (metaMDSdist). Sets number of
tries also when this was not set previously (metaMDSiter)
Version 1.8-2 (June 13, 2006)
* version 1.8-1 failed test in Windows because of a wrong encoding
name. Explicit \enc added for non-ascii words. Kurt Hornik and
Uwe Ligges diagnosed this and led me to see the light.
* similar to devel version 1.7-97.
Version 1.8-1 (June 12, 2006)
* Based on devel version 1.7-96.
* Pierre Legendre joined the vegan team.
* beals: a new function for Beals smoothing.
* bioenv: added 'trace' argument.
* cca/rda/capscale: accept several 'Condition' elements in the
formula.
* capscale: capscale(y ~ ., data=...) or expansion of "." on the
rhs works now. Documentation recognizes now db-RDA as the real
mother method.
* scores.cca, summary.cca etc: rewrite so that is cleaner and
easier to maintain. User visible changes are scaling by species
standard deviation (negative scaling) for 'rda', scaling=0
(no scaling) for all methods and slightly changed output and
improved user control in summary. These scalings actually were
documented in 1.6-10, although I dropped them just before the
release. predict.cca, predict.rda: work now with 'newdata' even
when not called with formula.
* anova.cca: new argument 'by' for tests of single terms or axes:
with by = "terms" performs individual test for constraints, and
with by = "axis" a separate sequential test for each axis. New
argument 'first' to analyse only the first axis instead of all
constrained variation.
* intersetcor: new function for the interset correlation or the
(weighted) correlation between individual constraints (contrasts)
and invidual axes in cca/rda/capscale. (Not recommended.)
* decostand: does not automatically convert matrix to a
data.frame. NA handling more consistent now (thanks to Tyler Smith
for diagnosis). Adds attribute 'decostand' giving the "method".
* linestack: now really uses median as the midpoint with the odd
number of cases, and does not give superfluous warnings with three
or less items. New argument 'labels' to replace the default text
in plot. The old argument 'label' renamed to 'side'. However, the
function still works with the old syntax, but gives a warning if
the old argument 'label' is used for 'side'. Returns invisibly the
shifted positions of labels.
* metaMDS: 'postMDS' sets now attributes similarly when called
independently or within 'metaMDS'. 'metaMDS' forwards arguments
to 'postMDS' except 'halfchange'. Change of phrasing in
'print'. Added handling of zero dissimilarities into
'metaMDSdist': either "fail" or "add" 1E-4 into zeros.
* mrpp: new function for the multiresponse permutation procedure
(MRPP). Code by Henry Stevens (Miami Univ, Oxford, Ohio).
* ordiarrows: draws arrowhead only in the last segment. New
argument 'startmark' for marking the starting point of the arrow.
* ordisurf: new arguments 'main' for the title and 'nlevels' and
'levels' for the number of contour or their values.
* orditorp: arguments for text colour and text character expansion
changed to standard 'col' and 'cex' from previous 'tcol' and
'tcex'.
* procrustes: 'summary' prints now rotation matrix, translation
and scale, and honours 'digits'.
* prc: new functions for Principal Response Curves (PRC) of van
den Brink and ter Braak (Envir. Toxicol. Chem, 18, 138-148; 1999).
This is a special rda() model with dedicated summary and plot
functions. New data set 'pyrifos' to demonstrate 'prc'.
* radfit: added brokenstick as a null model (rad.null), removed
rad.veil (as it was a bad idea originally), corrected minor bugs
in rad.preempt (which did not fail gracefully). Line colours
start from the point colour in plot.radfit.frame. 'print' uses
"g" format for coefficients and adds 'digits' argument. Added
'summary.radfit.frame' that simply prints each model.
* rankindex: uses now cluster:::daisy when 'grad' includes
factors.
* spantree: now a method function with 'plot', 'lines' and
'cophenetic' methods. 'lines' replaces 'ordispantree'. The plot
has a weird, new way of finding configuration for a spanning tree
from cophenetic distances (unpublished). Documented separately.
* specaccum: 'plot' honours now 'ylim'.
* specpool: Chao richness was wrongly defined, but now uses the
biased formula (error introduced in 1.6-5, correct earlier).
Failed with zero species or if there were no species that occurred
exactly one in the 'pool' (thanks to Emmanuel Castella, Geneve CH,
for the bug report).
* varpart: new functions for unbiased partitioning of variation by
two to four explanatory tables in RDA or linear regression. The
author of these functions is Pierre Legendre & co (Univ
Montreal). New data sets 'mite', 'mite.env' and 'mite.pcnm' to
demonstrate the functions.
* vegandocs: new function to display *all* pdf documentation and
ChangeLog. This really should be in base R -- this is a kluge to
fill the hole.
* vegemite: added argument to 'select' a subset of sites. Drops
missing species from the table. Prints number of species and
number of sites and the used cover scale at the end of the table.
Passes arguments (i.e., 'maxabund') to coverscale().
* coverscale: added argument 'maxabund' to 'scale = "log"'.
Returns the name of the cover scale as an attribute for vegemite()
to print.
* vegdist: now first checks input and then transforms (if
needed). Thanks to Tyler Smith,
* Internal changes: permuted.index acceptes NULL strata as an
alternative to missing strata. ordispantreee deprecated.
spider.cca removed.
* Documentation: general cleanup in help files. New chapter on
t-values in cca/rda/capscale in vegan-FAQ. New pdf document on
partioning with varpart by Pierre Legendre & co. Non-latin
characters now use UTF-8 in documentation. R manual says that you
should not use non-latin characters in help files, but that was
written by Englishmen. However, this seems to cause distress to
some users of a US West Coast OS (Windows), but OK with mainstream
OS's (Linux, MacOS). .
Version 1.6-10 (September 26, 2005)
* Based on devel version 1.7-77. Checked with R 2.1.1 (stable)
and R 2.2.0 (alpha).
* rda: negative scalings explicitly ignored and treated as
corresponding positive values. Function summary.rda used to fail
with NA centroids.
* permutest.cca & anova.cca: permutation of 'cca' result now
re-weights environmental data properly with permuted community
weights (this breaks compliance with popular proprietary
software). New default method 'direct' that always permutes the
data instead of residuals. Now clearly faster basic routines, but
re-weighting in 'cca' is costly, and permutations may even be
slower than earlier in small data sets. The permutest.cca returns
more data: constrained and residual total inertia, degrees of
freedom in the input model, and .Random.seed used in iterations.
Thanks to Pierre Legendre for pushing me to make this faster.
* ordination plot functions: obey now xlim and ylim.
* ordination text and points functions: a new argument 'select'
that can be used to select a subset of items, and a new argument
'labels' for text used instead of the default row names.
* points.cca, text.cca: biplot arrows will be scaled automatically
to fit the current graph if 'arrow.mul' is not given. The new
behaviour is similar to the default in plot.cca and in
plot.envfit.
* orditorp: a new "ordination text or points" function to add text
or points to an existing plot: adds text if this can be done
without overwriting other text labels, and points otherwise.
* linestack: a new function to draw labelled one-dimensional
diagrams without overwriting the labels (a primitive one, and
could easily be improved: submissions are welcome).
* ordirgl, orgltext: adapted to changes in rgl package version
0.65. Workaround for older rgl packages, too.
* decostand: Added Hellinger transformation. Empty columns and
rows become now 0 instead of NaN in most methods, except when the
input data contains negative values. Warns on input with negative
entries or on output containing NaN. Has now argument
'range.global' for method 'range' based on the code supplied by
Tyler Smith.
* plot.radfit: puts now legend "topright" in R >= 2.1.0 (which has
this keyword).
* read.cep: issues a warning if vegan was compiled with gfortran,
which has a bug that may corrupt the result. The bug concerns
multiline input with T format modifier, and it was corrected in
http://gcc.gnu.org/ml/gcc-patches/2005-09/msg00126.html, but still
bugs most released versions of gcc.
* vegdist: new indices 'raup' and 'binomial'. Method 'raup'
implements probabilistic Raup-Crick index and is based on the code
submitted by Michael Bedward. Method 'binomial' implements
Millar's index, officially published as "Binomial deviance as a
dissimilarity measure" (the C code has been in vegan for about two
years, but now I tell about it and add the public interface).
Method 'mountford' will give NA for any comparison involving an
empty site. Warns about empty sites or negative entries with all
methods except 'euclidean' and 'manhattan'.
* Documentation: added documentation of the cca/rda/capscale
result object. Dontruns used more neatly.
Version 1.6-9 (April 22, 2005)
* Maintenance release: 1.6-8 failed in R 2.1.0 patched and R 2.2.0
devel (works in R 2.1.0 release) due to problems in
as.preston. Based on devel version 1.7-62.
* as.fisher, as.preston: used table() in a way that failed in R
2.1.0 patched
* calibrate.cca: new function to predict or calibrate or
bioindicate the values of environmental constraints from
community composition (ordination).
* decostand: new argument na.rm (defaults FALSE) for ignoring
missing values in row, column or matrix standardizations.
* vegdist: new argument na.rm (defaults FALSE) for pairwise
deletion of missing vaues in dissimilarity calculation.
Version 1.6-8 (April 18, 2005)
* Based on devel version 1.7-59. Adapted to R 2.1.0 beta.
* DESCRIPTION: gives due credit to Roeland Kindt and Bob O'Hara as
co-authors.
* documentation: updates in capscale, vegdist. vegan-FAQ adapted
to changes in Sweave in R 2.1.0.
* several methods assumed that input is count data, but silently
accepted floating point numbers and gave wrong results. Now they
stop with error with non-integer input: fisherfit, prestonfit,
prestondistr, rarefy, fisher.alpha, estimateR.
* bioenv: uses now 'cor' instead of 'cor.test', and does not give
so many superfluous warnings and is marginally faster. Changed
printed output, so that gives 'call' instead of names of community
and environmental data, since bioenv.formula could not handle
these cleanly.
* capscale: has now argument 'add' to use an additive constant to
all dissimilarities so that all computed eigenvalues are
non-negative. This is an argument of underlying 'cmdscale'
function, which implements the "Correction method 2" of Legendre &
Legendre (1998), p. 434.
* decorana: checks now that there are no negative data entries.
* dune: cleaner site names.
* envfit: 'plot' now automatically scales arrows similarly as
'plot.cca' if 'arrow.mul' is not specified and arrows are added to
an old plot. Has now a 'scores' function.
* goodness.cca, predict.cca: documented separately.
* goodness.metaMDS: new function to assess pointwise goodness of
fit in metaMDS or isoMDS.
* humpfit: user can now give starting values of parameters.
'summary.humpfit' returns 'cov.unscaled' so that users can apply
'ellipse' function of 'ellipse' package to display approximate
confidence ellipses for pairs of parameters (with Normal
assumptions).
* make.cepnames: No longer duplicates one component names, but
'abbreviate's them to eight characters.
* metaMDS: Split to independent metaunits 'metaMDSdist' and
'metaMDSiter'. New function 'metaMDSredist' tries to reconstruct
dissimilarity matrices and transformations. 'postMDS' skips
halfchange scaling with too few points (with a warning) instead of
stopping with an error. Prints now info about 'postMDS'
operations. New 'scores' function. Improved 'plot' function with
'display' argument and labelling of axes. Argument 'shrink' to
undo expansion of species scores in plot or scores functions.
Workaround for a future bug in 'isoMDS' which drops names of
points in R 2.1.0. Updates number of 'tries' with
'previous.best'.
* ordiplot: has now 'display' argument, so that only species or
sites can be plotted.
* ordiplot3d: a new function for 3D static ordination diagrams.
* ordirgl: new function for dynamic 3D graphics of ordination
results. With support functions orglpoints, orgltext,
orglsegments and orglspider to add graphical items to dynamic
plots. Needs 'rgl' package.
* predict.cca, predict.rda: type = "wa" with 'newdata' works now
even when some species were removed from the original fit because
they were all-zero (missing). Default now to 'model = "CA"' in
unconstrained analysis without "CCA" component.
* predict.decorana: a new function for 'decorana' results. Similar
to 'predict.cca' (but refuses to give impossible results).
* procrustes: obeys now 'choices' even when input is a matrix.
Plots now projections of rotated axes. New argument 'scores' so
that can handle other than "site" scores.
* protest: passes now arguments to 'scores' so that now can handle
other than "site" scores on more than two dimensions.
* rda: handles now species with constant values (typically
missing or all zeros).
* spenvcor: new function to find the "species -- environment
correlation" in contrained ordination (cca, rda, capscale).
* stressplot: a new function to plot Shepard diagram for 'metaMDS'
or 'isoMDS'.
* summary.cca/summary.rda: print now call.
* vegemite: can now 'use' 'dendrogram' objects.
Version 1.6-7 (Jan 24, 2005)
* Based on devel version 1.7-42.
* plot.envfit did not plot vectors. Thanks to Roeland Kindt and
Ron E. VanNimwegen for bug reports and fixes.
* ordisurf obeys now keyword 'display', and '...' will transfer
arguments to 'scores'.
* ordisegments: obeys now keyword 'display'.
* ordigrid: works now.
* ordihull, ordiarrows, ordisegments, ordispider, ordiellipse:
have a new keyword 'show.groups' to show only specified groups.
With the help of this argument it is possible to use different
colours and linetypes for each group. Further, it makes possible
to plot results of logical comparisons (such as A1 > 4) by setting
show.groups = T. Two user requests.
* new diagnostic and helper functions for 'cca', 'rda' and
'capscale': 'goodness' to estimate the proportion of inertia
accounted for or residuals for sites or species; 'inertcomp' to
decompose species and site inertia for conditioned, contrained and
residual componets; 'vif.cca' to estimate the variance inflation
factors for constraints and conditions; 'fitted' and 'residuals'
to approximate data by ordination scores; 'predict' to approximate
data or find site or species scores, possibly with 'newdata';
'coef' to find the regression coefficients. 'alias.cca' is now
similar to 'alias.lm' (simplified version of the latter).
Version 1.6-6 (Jan 7, 2005)
* Based on version 1.7-34, but without still experimental
functions predict, fitted, vif & goodness for cca and rda
objects.
* mantel: implemented partial mantel test (function
'mantel.partial'). Both 'mantel' and 'mantel.partial' should be
marginally faster and not give so many superfluous warnings.
* plot.envfit: should be more reliable now. A new keyword 'at' to
position the bunch of arrows in other places than the origin.
Thanks to Roeland Kindt for fixes.
* rankindex: Default correlation now "spearman" since the older
default ("kendall") was far too slow in large data sets. Should
not give so many superfluous warnings.
* summary.rda: Site scores were wrong with 'scaling = 3'.
Influences 'plot' and 'scores' commands for 'rda' and 'capscale'
results with 'scaling = 3'.
* vegdist: Issues a warning if 'method = "morisita"' is used with
non-integer data.
Version 1.6-5 (Oct 12, 2004)
* Based on version 1.7-27. Checked with R-1.9.1 (Linux, MacOS X)
and R-2.0.0 (Linux, patched version in Windows XP). Passed other
tests, but some examples in 'humpfit' failed in Windows XP, and
are not run on that platfrom. General cleanup in documentation.
Does not 'require(mva)' any longer.
* anosim: corrected the equation in docs (function was
correct). Thanks to Yong Cao for notifying this.
* bioenv.formula: finds now variables from attach'ed data.frames
and 'data' need not be specified.
* capscale & deviance.rda: capscale was modified so that its
result could be handled with step (added a terms component).
Function deviance.capscale was deleted so that deviance.rda will
be used. A literature reference to AIC in CCA/RDA was added to the
documentation.
* cca/rda/capscale: Small eigenvalues are made to zero and rank
reduced accordingly. Most often influences 'capscale'.
* cca/rda/capscale: Now alias constraints that are collinear with
other constraints or conditions in partial analysis. New function
'alias.cca' to print aliased contrasts. NB this does not influence
the proper ordination results, but only to the selection of biplot
scores and centroids displayed. 'summary' no longer gives a
redundant warning.
* decorana: removes empty species with a warning instead of
stopping with error. Still error with empty sites. The behaviour
is similar to 'cca'.
* envfit (vectorfit, factorfit): Has now a formula interface.
Bottlenecks in permutation changed now into C functions. Much
faster, in particular 'factorfit' which was 25 times faster in my
tests. plot.envfit could drop names (dimensions). Now honours
'scaling' argument in cca/rda/capscale (used to 'scaling = 2'
always).
* envfit, ordisurf: have now weights parameter 'w'. Weights are
used for the results of 'cca' and 'decorana' as default, or they
can be supplied by the user. In this way, envfit gives equal
results to biplot scores and centroids in 'cca', and ordisurf is
consistent with envfit. Both can now access 'lc' scores of 'cca'
etc.
* estimateR: Abundance-based estimators of extrapolated richness,
written by Bob O'Hara <bob.ohara@helsinki.fi>.
* humpfit: Asesses now the dispersion parameter in non-Poisson
models (summary.humpfit). Added a profile method that inherits
from profile.glm, so that you can use plot.profile.glm,
pairs.profile.glm and confint.profile.glm of MASS for displaying
profile and finding the confidence intervals from the profile
likelihood.
* fisherfit: prints now standard error of the estimate. Added
'profile' and 'confint' methods.
* metaMDS: A new function to collect all Peter Minchin's
recommendations to run nonmetric multidimensional scaling: (1)
adequate data transformation with Wisconsin double and sqrt if
needed, (2) use of adequate dissimilarity measure, (3) possible
stepacross if needed (this Minchin didn't have), (4) run NMDS with
several random starts and stop after finding two similar
solutions, (5) scale results with postMDS, and (6) add species
scores with wascores. Function has print and plot methods.
Function scores.default changed to understand metaMDS.
* specpool: Chao is now based on unbiased equation. Standard
errors now estimated for Chao, 1st-order jackknife and bootstrap
richness, but SE of 2nd-order jackknife still missing.
* vegdist: an option for binary indices, since some users believed
these are not in vegan, although you can get them with
'decostand'.
Version 1.6-4 (June 10, 2004)
* Based on 1.7-12.
* anova.cca: Changed defaults to beta=0.01 and perm.max=10000.
* bioenv: A new function implementing the BIO-ENV procedure of
Clarke & Ainsworth (Mar. Ecol. Prog. Ser. 92, 205-219; 1993). The
function finds the subset of scaled environmental variables with
maximum (rank) correlation with community data.
* cca: added ter Braak scaling, a.k.a. Hill scaling, with
negative values of `scaling'. This scaling should approximate
Hill's non-linear rescaling of `decorana', but it the scaling is
linear and does not rescale. The scaling is very simple: the
results from corresponding positive scaling are multipiled with
sqrt(1/(1-lambda)).
* plot.cca: Better heuristics used for the length of biplot
arrows, and they are now longer in general. Axes are drawn for
biplot arrows also with `text.cca' and `points.cca' functions.
* deviance.cca and extractAIC.cca: auxilliary functions which
allow automatic model building with `step' or `stepAIC' (MASS)
functions in cca and rda. Unfortunately, cca and rda do not have
deviance or AIC, and these functions are certainly wrong and
dangerous. However, with continuous candidate variables they
select the variables in the same order as Canoco.
* humpfit: A new family of functions to fit the no-interaction
model (Oksanen, J., J. Ecol. 84, 293-295; 1996) to species
richness -- biomass data.
* vegdist: corrected a bug in Gower index -- range standardization
was made to sites, but it should be done to species. Thanks to
Brett Melbourne <bamelbourne@ucdavis.edu> for the bug report.
Version 1.6-3 (Mar 22, 2004)
* Still based on 1.7-3, but ported some changes from version
1.7-7 (envfit, procrustes). This version was produced to prepare
for incompatible changes in coming R-1.9.0.
* vegan-FAQ: Does not show the vegan and R versions the
documentations were built, because of stupid incompatible change
in R-1.9.0 of the future. package.description() function was
unnecessarily replaced with packageDescription, and to accomodate
recent, present and future versions of R seemed to be too much
hassle.
* plot.envfit: Added option 'p.max' to display only environmental
variables assessed to have P values less or equal to the given
limit.
* plot.procrustes: added option kind = 0 to draw only axes and
functions points.procrustes and lines.procrustes to add points and
line segments or arrows to the plot.
Version 1.6-2: Based on devel version 1.7-3.
* renyi: should work now (again?).
* prestondistr: The truncation level was badly chosen. Now default
level is set to -1, or log2(-1) = 0.5 which might be the left
limit of the first octave. This is now a parameter, leaving the
choice as the responsibility of the user.
* cca.default: Removes missing species and issues a warning
instead of stopping with error.
Version 1.6-1: Based on devel version 1.7-1.
* Minor maintanenance release to satisfy tests in R-devel. There
was a buglet in ordisegments that caused a warning. I had
introduced a trick to handle printCoefmat before it was invented
to replace print.coefmat, but this failed in R-devel.
* Vignettes included in this release. These were only in devel
versions of my web pages previously.
Version 1.6-0 (Fisher)
* Based on devel version 1.5-58. Passes checks with R-1.8.0 and
R-1.7.0 and compiles into a Windows binary. Many functions were
contributed by Roeland Kindt (Nairobi, Kenya).
* BCI: A new data set on tree counts on 1-ha plots in the Barro
Colorado Island. Thanks to Roeland Kindt and Richard S. Condit.
* capscale: A new function for [partial] constrained analysis of
principal coordinates (CAP).
* cca, rda: Can now handle cases where the number of constraints
is higher than the number of species.
* cca, decorana: Find now if data have empty rows or columns.
* decostand: Added Chi-square transformation.
* distconnected: A new function to find if vegetation data are
disconnected below given threshold dissimilarity level. Utility
function no.shared gives a logical dissimilarity matrix with
values TRUE for cases with no shared species.
* fisher.alpha: Should work with large sample sizes now (bug
report thanks to Ariel Bergamini).
* fisherfit, prestonfit, prestondistr: New functions to fit
Fisher's logarithmic series and Preston's log-normal model.
* make.cepnames: New utility function to change Latin names into
4+4 character CEP names.
* ordicluster, ordispantree: New functions in the ordihull family
to overlay a cluster dendrogram or a spanning tree onto an
ordination. Weights are now really used with cca and decorana
plots.
* radfit: New function to fit Ranked Abundance Dominance or
Dominance/ Diversity models with maximum likelihood: Pre-emption,
log-normal, veiled log-normal, Zipf and Zipf--Mandelbrot models.
* rankindex: Can now use stepacross and new dissimilarity indices
of vegdist.
* rarefy: Can now optionally find SE of rarefied richness.
* renyi: A new function to find Rényi diversities or Hill numbers
with any scale (thanks to Roeland Kindt).
* scores.ordiplot: should be more robust now.
* spantree: A new function to find a minimum spanning tree using
only dissimilarities below a threshold or disregarding NA
dissimilarities.
* specaccum: A new function for species accumulation
curves. Thanks to Roeland Kindt.
* specnumber: A very simple utility function to find species
richness.
* specpool: New function for extrapolated species richness in a
species pool, or for estimating the number of unseen species.
* stepacross: New function for estimation of dissimilarities
between sites that have no shared species. Implements both
flexible shortest paths and their approximation known as extended
dissimilarities.
* vegdist: New dissimilarity indices Jaccard (finally), Morisita,
Horn--Morisita and Mountford.
Version 1.4-4
* Based on devel version 1.5-35.
* decorana: Finds now eigenvalues from the final solution as ratio
of biased weighted variances of site and species scores. The
values returned by the legacy Fortran function are called
``decorana values''.
* downweight: passes the downweighting fraction as an attribute,
and decorana catches and prints the fraction.
* wascores: Uses now biased variances for expading WAs and returns
the shrinkage factors as an attribute "shrinkage". Shrinkage
factors are equal to eigenvalues in CCA when only this one
variable is used as constraint. Function `eigengrad' returns only
these eigenvalues.
* rda summary/plot bugfix: Failed with non-default scaling (1 or
3) when there were no factor constraints. Corrected a statement
giving a harmless warning in plot.cca.
* ordiplot knows now option type = "text". New simpleminded
functions points.ordiplot and text.ordiplot for adding items in an
ordiplot graph.
* ordiarrows, ordisegments could drop dimensions (when will I
learn!) and fail if only one arrow or one segment should be
drawn. One argument was missing in a `gl' command in `ordigrid'.
Version 1.4-3 (Luova)
* Based on devel version 1.5-30
* `cca' and `rda' use canonical expansion of the formula and
return the `terms' component. Functions like `terms', `formula'
and `update' magically started to work with cca and rda (and eight
lines of code would allow `step' and `stepAIC' to work magically
in model building... but it's so much magic that I don't trust
it). Some of the allowed things are now `mod <- cca(dune ~ . -
Use, dune.env)'
(all variables but `Use' in dune.env), and `update(mod, . ~ . -
Manure)' (remove `Manure' from the previous).
* `cca' and `rda' find centroids of factor levels with factor
constraints (only with formula interface). These can be accessed
with `display="cn"' in `scores', `summary' and `plot'. If
`display="cn"' is requested in `plot', biplot arrows are still
drawn unless there is a class centroid with the same name:
continuous variables are still shown as arrows and ordered factors
both as arrows as centroids.
* Enhanced user control of many ordination plots. Some accept more
graphical parameters (`ordiplot', `plot.procrustes'). New
functions `points.cca' and `text.cca' to add these items into
CCA/RDA plots (documented separately).
* `identify.ordiplot' knows more alternatives, among them
`plot.procrustes'. This was helped with new, very generic
`scores.ordiplot'.
* New functions to add graphical items in ordination diagrams:
`ordihull' draws convex hulls for groups, `ordiarrows' draws
arrows, `ordisegments' segments, `ordigrid' combines points in to
a grid, `ordispider' combines points to their (weighted) centroids
or WA scores to the corresponding LC score in cca/rda,
`ordiellipse' draws dispersion or confidence ellipses for
points. All these take a groups argument for selecting subsets of
points. `ordispider' obsoletes `spider.cca' (introduced in the
previous release).
* Implemented simple ``deshrinking'' of weighted averages in
`wascores'.
* A new diversity function `fisher.alpha' that estimates the
diversity as the alpha parameter of Fisher's log-series. Thanks to
Bob O'Hara <bob.ohara@helsinki.fi> for the main code.
* `summary.decorana' now in canonical form so that printing is
done by `print.summary.decorana' and user can intercept and catch
the result. Several users had been surprised of the earlier
non-canonical behaviour. `summary' knows again the prior weights
used in downweighting of rare species (these evaporated in 1.4-0).
* rda summary/scores/plot bugfix: Scaling of site scores was wrong
with scaling alternatives 1 and 3.
* envfit (vectorfit, factorfit) bug fixes: Can now handle single
vectors (labelling, permutation). Recognize now `choices' so that
the user can select ordination axes.
* Checked with released R-1.7.0.
Version 1.4-2
Based on devel version 1.5-18
* New functions and data sets:
* dune: A new data set -- the classic Dune Meadow data with the
environmental data.
* ordiplot, identify.ordiplot: New functions to plot "any"
ordination result and identify plotted points. Intended to provide
similar functionality as `plot', `plotid' and `specid' functions
in `labdsv' without name clashes. Functions `plot.cca' and
`plot.decorana' return (invisibly) `ordiplot' objects as well, and
these can be used by `identify' to label species and sites.
* rda: Redundancy Analysis, or optionally, yet another PCA. This
is a spin-off from cca(), and it is documented together with
cca. The only new functions are rda, rda.default, rda.formula and
summary.rda. Otherwise rda uses cca methods which were changed to
be aware of rda.
* read.cep: A function to read Cornell Ecology Package or CEP
formatted files into R. This has been in devel versions since
1.1-1, but it never made to the releases, since this worked only
in Linux, but crashed R in Windows. It seemed to work in R-1.6.1
with gcc 3.2
(MinGW) in Windows, so it is included. Treat with caution.
* spider.cca: A tiny function that joins each WA score to the
corresponding LC score in cca/rda plots.
Minor changes:
* Checked with R-1.7.0 (unstable devel version, snapshot
2003-02-05) and corrected as needed.
* Added symmetric scaling (=3) to alternitves in cca and rda.
* Minor improvements and bug fixes in vegemite.
Version 1.4-1 (Logan)
Based on devel version 1.5-10.
New features:
* Permutation tests added to envfit (vectorfit, factorfit) and
procrustes (called protest).
* All permutation tests take now argument strata: If given,
permutations made only within strats. Concerns anosim, anova.cca
(permutest.cca), mantel, envfit (factorfit, envfit) and protest
(in procrustes).
* fitted.procrustes(): A new function.
* rarefy(): A new function for rarefaction species richness.
Other changes:
* cca.default(): Handles now NULL matrices X and Z: skips them.
* cca.formula(): Knows now cca(X ~ 1) and permforms unconstrained
CA, and cca(X ~ ., data) and perfors CCA using all variables in
`data' as constraints. Has now na.action=na.fail so that cca
crashes more informatively (used to crash mysteriously). A more
graceful na.action may come. Assignment "=" corrected to the true
blue "<-" (used to crash in old R).
* envfit (factorfit, vectorfit): Use `scores' now. Various
bugfixes and should work now with single factors or vectors.
* plot.envfit(): Options `choices' works now like documented.
* ordisurf(): A new name to surf() to avoid name clash with
labdsv.
* procrustes(): added option `symmetric' for symmetric rotation
and goodness of fit statistic.
* plot.procrustes(): Keyword `axes' changed to `choices' to be
consistent with other functions.
* vegemite(): A new name to vegetab. The name was chosen because
the output is so compact (and to avoid confusion with function
vegtab in labdsv). Function checks now that the cover codes fit in
one column and splits the output if it does not fit into used page
width.
Version 1.4-0
Based on devel version 1.3-19.
New functions:
* `anosim': Analysis of Similarities.
* `cca': [Partial] [Constrained] Correspondence Analysis. This is
an alternative to Dray's `CAIV' (package `CoCoAn'). Differences
include formula interface, WA scores in addition to LC scores,
partial CCA, residual CA after constraints, algorithm based on svd
(and much faster), more standard support functions (`plot',
`summary', `print').
* `anova.cca' and `permutest.cca': Permutation tests for `cca'.
* `envfit': A new wrapper function which calls either `vectorfit'
or new `factorfit' function, and `plot.envfit' for easier display
of results in graphs.
* `mantel': Mantel test for two
dissimilarity matrices.
* `scores': A generic function to extract scores from `cca',
`decorana' or even from some common ordination methods in R.
* `surf': Surface fitting for ordination, intended as an
alternative to `vectorfit'.
* `vegetab': Prints a compact, ordered vegetation table, together
with `coverscale' to transform data into compact format.
Other changes:
* `downweight': directly callable instead of being embedded in
`decorana'. Now it can be used with other downweighting thresholds
or with other CA's than `decorana', e.g., `cca', `ca', `CAIV'.
* `plot.decorana': new keyword `type', removed keyword `tol'.
* `vectorfit': centres now ordination before fitting, so it will
work with other ordinations than MDS.
* Mainly user invisible changes in `plot.decorana',
`plot.procrustes', `postMDS', `print.summary.procrustes',
`vectorfit', `procrustes'. One of the main changes was that
`eqscplot' of `MASS' library was replaced with `plot(..., asp=1)'
of standard R (but this is not S-plus compliant).
* Edited help files. More examples run from help files.
* Tried to improve LabDSV compliance, see
http://labdsv.nr.usu.edu/.
* Checked with R-1.5.0. However, `log' in diversity not changed
into `logb' (yet), although base is specified, because `logb'
won't work in R pre-1.5.0.
Version 1.2-1
* Checked with R-1.4-0 (frozen snapshot) tools: Documentation and
method consistency corrected (summary.decorana.R, decorana.Rd),
and file permissions changed.
* decorana.Rd: removed dependence on package `multiv' in the
example, since it is not a `recommended package' -- replaced with
`decorana(..., ira=1)'.
Version 1.2-0
Based on experimental version 1.1-6. Doesn't include read.cep() of the
experimental version, because the function is certainly unportable to
all platforms (works on Red Hat Linux 7.1, but perhaps nowhere else).
Major new feature:
* decorana(): R port of Mark Hill's DECORANA for Detrended
correspondence analysis, with methods print(), summary() and
plot().
Minor changes:
* postMDS(): finds isoMDS scores and doesn't destroy list
structure.
* vectorfit: finds isoMDS scores and decorana row scores (which
are uncentred, though).
* Proof reading in documentation.
Version 1.0-3
*`plot.procrustes': selection of axes plotted, improved scaling in
graphics, control of axis labels.
* `postMDS': added plot to demonstrate half-change scaling.
* `procrustes': target can now have lower dimensionality than
rotated matrix.
Version 1.0-2
* Fixed bugs in documentation.
* `plot.procrustes': Fixed passing `...' to graphics functions.
Version 1.0-1
* The first public version.
|