]> git.donarmstrong.com Git - ape.git/blob - ChangeLog
new files for trait simulations
[ape.git] / ChangeLog
1                 CHANGES IN APE VERSION 2.4-2
2
3
4 NEW FEATURES
5
6     o The new functions rTraitCont and rTraitDisc simulate continuous and
7       discrete traits under a wide range of evolutionary models.
8
9     o add.scale.bar() has a new option 'ask' to draw interactively.
10
11
12 BUG FIXES
13
14     o add.scale.bar() drew the bar outside the plotting region with the
15       default options with unrooted or radial trees.
16
17
18
19                 CHANGES IN APE VERSION 2.4-1
20
21
22 NEW FEATURES
23
24     o rtree() and rcoal() now accept a numeric vector for the 'br'
25       argument.
26
27     o vcv() is a new generic function with methods for the classes "phylo"
28       and "corPhyl" so that it is possible to calculate the var-cov matrix
29       for "transformation models". vcv.phylo() can still be used for trees
30       of class "phylo"; its argument 'cor' has been renamed 'corr'.
31
32
33 BUG FIXES
34
35     o bind.tree() failed when 'y' had no root edge.
36
37     o read.nexus() shuffled tip labels when the trees have no branch
38       lengths and there is a TRANSLATE block.
39
40     o read.nexus() does not try to translate node labels if there is a
41       translation table in the NEXUS file. See ?read.nexus for a
42       clarification on this behaviour.
43
44     o plot.multiPhylo() crashed R when plotting a list of trees with
45       compressed tip labels.
46
47     o write.nexus() did not translate the taxa names when asked for.
48
49     o plot.phylo(type = "fan") did not rotate the tip labels correctly
50       when the tree has branch lengths.
51
52     o ace(type = "continuous", method = "ML") now avoids sigma² being
53       negative (which resulted in an error).
54
55     o nj() crashed with NA/NaN in the distance matrix: an error in now
56       returned.
57
58
59
60                 CHANGES IN APE VERSION 2.4
61
62
63 NEW FEATURES
64
65     o base.freq() has a new option 'freq' to return the counts; the
66       default is still to return the proportions.
67
68
69 BUG FIXES
70
71     o seg.sites() did not handle ambiguous nucleotides correctly: they
72       are now ignored.
73
74     o plot(phy, root.edge = TRUE) failed if there was no $root.edge in
75       the tree: the argument is now ignored.
76
77     o add.scale.bar() failed when 'x' and 'y' were given (thanks to Janet
78       Young for the fix).
79
80
81 OTHER CHANGES
82
83     o Trying to plot a tree with a single tip now returns NULL with a
84       warning (it returned an error previously).
85
86     o The way lines representing nodes are coloured in phylograms has
87       been modified (as well as their widths and types) following some
88       users' request; this is only for dichotomous nodes.
89
90     o The argument 'adj' in [node][tip][edge]labels() now works when
91       using 'pie' or 'thermo'.
92
93     o A more informative message error is now returned by dist.dna() when
94       'model' is badly specified (partial matching of this argument is
95       done now).
96
97     o Deprecated functions are now listed in a help page: see
98       help("ape-defunct") with the quotes.
99
100
101 DEPRECATED & DEFUNCT
102
103     o The functions heterozygosity, nuc.div, theta.h, theta.k and
104       theta.s have been moved from ape to pegas.
105
106     o The functions mlphylo, DNAmodel and sh.test have been removed.
107
108
109
110                 CHANGES IN APE VERSION 2.3-3
111
112
113 BUG FIXES
114
115     o add.scale.bar() always drew a horizontal bar.
116
117     o zoom() shuffled tips with unrooted trees.
118
119     o write.nexus() failed to write correctly trees with a "TipLabel"
120       attribute.
121
122     o rcoal() failed to compute branch lengths with very large n.
123
124     o A small bug was fixed in compar.cheverud() (thanks to Michael
125       Phelan for the fix).
126
127     o seg.sites() failed when passing a vector.
128
129     o drop.tip() sometimes shuffled tip labels.
130
131     o root() shuffled node labels with 'resolve.root = TRUE'.
132
133
134
135                 CHANGES IN APE VERSION 2.3-2
136
137
138 BUG FIXES
139
140     o all.equal.phylo() did not compare unrooted trees correctly.
141
142     o dist.topo(... method = "PH85") did not treat unrooted trees
143       correctly (thanks to Tim Wallstrom for the fix).
144
145     o root() sometimes failed to test for the monophyly of the
146       outgroup correctly.
147
148     o extract.clade() sometimes included too many edges.
149
150     o vcv.phylo() did not work correctly when the tree is in
151       "pruningwise" order.
152
153     o nj() did not handle correctly distance matrices with many 0's.
154       The code has also been significantly improved: 7, 70, 160 times
155       faster with n = 100, 500, 1000, respectively.
156
157
158
159                 CHANGES IN APE VERSION 2.3-1
160
161
162 NEW FEATURES
163
164     o The new function is.monophyletic tests the monophyly of a group.
165
166     o There is now a c() method for lists of class "DNAbin".
167
168     o yule.cov() now fits the null model, and its help page has been
169       corrected with respect to this change.
170
171     o drop.tip() has a new option 'rooted' to force (or not) a tree
172       to be treated as (un)rooted.
173
174
175 BUG FIXES
176
177     o dist.gene() failed on most occasions with the default
178       pairwise.deletion = FALSE.
179
180     o read.tree() failed to read correctly the tree name(s).
181
182     o boot.phylo() now treats correctly data frames.
183
184     o del.gaps() did not copy the rownames of a matrix.
185
186     o A small bug was fixed in CDAM.global().
187
188     o ace() failed with large data sets. Thanks to Rich FitzJohn for
189       the fix. With other improvements, this function is now about 6
190       times faster.
191
192     o write.tree() failed with objects of class "multiPhylo".
193
194     o drop.tip(, subtree = TRUE) sometimes shuffled tip labels.
195
196
197 OTHER CHANGES
198
199     o [.multiPhylo and [.DNAbin now respect the original class.
200
201     o Instances of the form class(phy) == "phylo" have been replaced
202       by inherits(phy, "phylo").
203
204     o rcoal() is now faster.
205
206
207 DEPRECATED & DEFUNCT
208
209     o klastorin() has been removed.
210
211
212
213                 CHANGES IN APE VERSION 2.3
214
215
216 NEW FEATURES
217
218     o The new functions CADM.global and CADM.post, contributed by
219       Pierre Legendre, test the congruence among several distance
220       matrices.
221
222     o The new function yule.time fits a user-defined time-dependent
223       Yule model by maximum likelihood.
224
225     o The new function makeNodeLabel creates and/or modifies node
226       labels in a flexible way.
227
228     o read.tree() and write.tree() have been modified so that they can
229       handle individual tree names.
230
231     o plot.phylo() has a new argument 'edge.lty' that specifies the
232       types of lines used for the edges (plain, dotted, dashed, ...)
233
234     o phymltest() has been updated to work with PhyML 3.0.1.
235
236
237 BUG FIXES
238
239     o drop.tip() shuffled tip labels in some cases.
240
241     o drop.tip() did not handle node.label correctly.
242
243     o is.ultrametric() now checks the ordering of the edge matrix.
244
245     o ace() sometimes returned negative values of likelihoods of
246       ancestral states (thanks to Dan Rabosky for solving this long
247       lasting bug).
248
249
250 OTHER CHANGES
251
252     o The data set xenarthra has been removed.
253
254
255
256                 CHANGES IN APE VERSION 2.2-4
257
258 BUG FIXES
259
260     o The bug fix in read.nexus() in version 2.2-3 was wrong: this is
261       now fixed. (Thanks to Peter Wragg for the fix!)
262
263     o A warning message occurred for no reason with ace(method="GLS").
264
265
266 OTHER CHANGES
267
268     o There is now a general help page displayed with '?ape'.
269
270
271
272                 CHANGES IN APE VERSION 2.2-3
273
274
275 NEW FEATURES
276
277     o The new function extract.clade extracts a clade from a tree by
278       specifying a node number or label.
279
280     o fastme.bal() has two new options 'spr' and 'tbr' to perform tree
281       operations of the same names.
282
283     o dist.dna() can now return the number of site differences by
284       specifying model="N".
285
286
287 BUG FIXES
288
289     o chronopl() did not work with CV = TRUE.
290
291     o read.nexus() did not work correctly in some situations (trees on
292       multiple lines with different numbers of lines and/or with
293       comments inserted within the trees).
294
295     o ltt.plot(), ltt.lines(), and mltt.plot() did not count correctly
296       the number of lineages with non-binary trees.
297
298
299 OTHER CHANGES
300
301     o ape has now a namespace.
302
303     o drop.tip() has been improved: it should be much faster and work
304       better in some cases (e.g., see the example in ?zoom).
305
306
307
308                 CHANGES IN APE VERSION 2.2-2
309
310
311 NEW FEATURES
312
313     o dist.gene() has been substantially improved and gains an option
314       'pairwise.deletion'.
315
316     o cbind.DNAbin() has a new option 'fill.with.gaps' and is now
317       more flexible.
318
319
320 BUG FIXES
321
322     o prop.part() failed with a single tree with the default option
323      'check.labels = TRUE'.
324
325    o summary.DNAbin() failed to display correctly the summary of
326      sequence lengths with lists of sequences of 10,000 bases or more
327      (because summary.default uses 4 significant digits by default).
328
329    o read.nexus() failed to read a file with a single tree with line
330      breaks in the Newick string.
331
332    o del.gaps() returned a list of empty sequences when there were no
333      gaps.
334
335
336 OTHER CHANGES
337
338     o phymltest() has been updated for PhyML 3.0 and gains an option
339       'append', whereas the option 'path2exec' has been removed.
340
341     o rbind.DNAbin() and cbind.DNAbin() now accept a single matrix
342       which is returned unchanged (instead of an error).
343
344     o The data sets bird.orders and bird.families are now stored as
345       Newick strings; i.e., the command data(bird.orders) calls
346       read.tree().
347
348
349
350                 CHANGES IN APE VERSION 2.2-1
351
352
353 NEW FEATURES
354
355     o The new function makeLabel() helps to modify labels of trees,
356       lists of trees, or DNA sequences, with several utilities to
357       truncate and/or make them unique, substituting some
358       characters, and so on.
359
360     o The new function del.gaps() removes insertion gaps ("-") in a
361       set of DNA sequences.
362
363     o read.dna() can now read Clustal files (*.aln).
364
365
366 BUG FIXES
367
368     o root() failed with 'resolve.root = TRUE' when the root was
369       already the specified root.
370
371     o Several bugs were fixed in mlphylo().
372
373     o collapsed.singles() did not propagate the 'Nnode' and
374       'node.labels' elements (thanks to Elizabeth Purdom for the fix).
375
376     o read.nexus() failed to remove correctly the comments within
377       trees.
378
379     o read.nexus() failed to read a file with a single tree and no
380       translation of tip labels.
381
382     o read.nexus() failed to place correctly tip labels when reading
383       a single tree with no edge lengths.
384
385     o A bug was fixed in sh.test().
386
387
388 OTHER CHANGES
389
390     o unique.multiPhylo() is faster thanks to a suggestion by Vladimir
391       Minin.
392
393     o The option 'check.labels' of consensus() and prop.part() is now
394       TRUE by default.
395
396     o write.dna() now does not truncate names to 10 characters with
397       the Phylip formats.
398
399
400
401                 CHANGES IN APE VERSION 2.2
402
403
404 NEW FEATURES
405
406     o Four new functions have been written by Damien de Vienne for the
407       graphical exploration of large trees (cophyloplot, subtrees,
408       subtreeplot), and to return the graphical coordinates of tree
409       (without plotting).
410
411     o The new functions corPagel and corBlomberg implement the Pagel's
412       "lambda" and Blomberg et al.'s "ACDC" correlation structures.
413
414     o chronopl() has been improved and gains several options: see its
415       help page for details.
416
417     o boot.phylo() has now an option 'trees' to possibly return the
418       bootstraped trees (the default is FALSE).
419
420     o prop.part() has been improved and should now be faster in all
421       situations.
422
423
424 BUG FIXES
425
426     o read.dna() failed if "?" occurred in the first 10 sites of the
427       first sequence.
428
429     o The x/y aspect of the plot is now respected when plotting a
430       circular tree (type = "r" or "f").
431
432     o Drawing the tip labels sometimes failed when plotting circular
433       trees.
434
435     o zoom() failed when tip labels were used instead of their numbers
436       (thanks to Yan Wong for the fix).
437
438     o drop.tip() failed with some trees (fixed by Yan Wong).
439
440     o seg.sites() failed with a list.
441
442     o consensus() failed in some cases. The function has been improved
443       as well and is faster.
444
445
446
447                 CHANGES IN APE VERSION 2.1-3
448
449
450 BUG FIXES
451
452     o A bug in read.nexus() made the Windows R-GUI crash.
453
454     o An error was fixed in the computation of ancestral character
455       states by generalized least squares in ace().
456
457     o di2multi() did not modify node labels correctly.
458
459     o multi2di() failed if the tree had its attribute "order" set to
460       "cladewise".
461
462
463
464                 CHANGES IN APE VERSION 2.1-2
465
466
467 NEW FEATURES
468
469     o There three new methods for the "multiPhylo" class: str, $,
470       and [[.
471
472     o root() gains the options 'node' and 'resolve.root'
473       (FALSE by default) as well as its code being improved.
474
475     o mltt.plot() has now an option 'log' used in the same way
476       than in plot.default().
477
478
479 BUG FIXES
480
481     o mltt.plot() failed to display the legend with an unnamed
482       list of trees.
483
484     o nodelabels() with pies now correcly uses the argument
485       'cex' to draw symbols of different sizes (which has
486       worked already for thermometers).
487
488     o read.nexus() generally failed to read very big files.
489
490
491 OTHER CHANGES
492
493     o The argument 'family' of compar.gee() can now be a function
494       as well as a character string.
495
496     o read.tree() and read.nexus() now return an unnamed list if
497       'tree.names = NULL'.
498
499     o read.nexus() now returns a modified object of class "multiPhylo"
500       when there is a TRANSLATE block in the NEXUS file: the individual
501       trees have no 'tip.label' vector, but the list has a 'TipLabel'
502       attribute. The new methods '$' and '[[' set these elements
503       correctly when extracting trees.
504
505
506
507                 CHANGES IN APE VERSION 2.1-1
508
509
510 NEW FEATURES
511
512     o The new function rmtree generates lists of random trees.
513
514     o rcoal() now generates a genuine coalescent tree by default
515       (thanks to Vladimir Minin for the code).
516
517
518 BUG FIXES
519
520     o nuc.div() returned an incorrect value with the default
521       pairwise.deletion = FALSE.
522
523
524 OTHER CHANGES
525
526     o The internal codes of bionj(), fastme.bal(), and fastme.ols()
527       have been improved so that they are stabler and faster.
528
529     o R packages used by ape are now loaded silently; lattice and gee
530       are loaded only when needed.
531
532
533
534                 CHANGES IN APE VERSION 2.1
535
536
537 NEW FEATURES
538
539     o The new function identify.phylo identifies clades on a plotted
540       tree using the mouse.
541
542     o It is now possible to subset a list of trees (object of class
543       "multiPhylo") with "[" while keeping its class correct.
544
545     o The new function as.DNAbin.alignment converts DNA sequences
546       stored in the "alignment" format of the package seqinr into
547       an object of class "DNAbin".
548
549     o The new function weight.taxo2 helps to build similarity matrices
550       given two taxonomic levels (usually called by other functions).
551
552     o write.tree() can now take a list of trees (class "multiPhylo")
553       as its main argument.
554
555     o plot.correlogram() and plot.correlogramList() have been
556       improved, and gain several options (see the help page for
557       details). A legend is now plotted by default.
558
559
560 BUG FIXES
561
562     o dist.dna() returned some incorrect values with `model = "JC69"'
563       and `pairwise.deletion = TRUE'. This affected only the
564       distances involving sequences with missing values. (Thanks
565       to Bruno Toupance for digging this bug out.)
566
567     o write.tree() failed with some trees: this is fixed by removing
568       the `multi.line' option (trees are now always printed on a
569       single line).
570
571     o read.nexus() did not correctly detect trees with multiple root
572       edges (see OTHER CHANGES).
573
574
575 OTHER CHANGES
576
577     o The code of mlphylo() has been almost entirely rewritten, and
578       should be much stabler. The options have been also greatly
579       simplified (see ?mlphylo and ?DNAmodel for details).
580
581     o The internal function nTips has been renamed klastorin_nTips.
582
583     o The code of is.ultrametric() contained redundancies and has
584       been cleaned-up.
585
586     o The code of Moran.I() and of correlogram.formula() have been
587       improved.
588
589     o read.tree() and read.nexus() now return an error when trying to
590       read a tree with multiple root edges (see BUG FIXES). The
591       correction applied in previous version did not work in all
592       situations.
593
594     o The class c("multi.tree", "phylo") has been renamed
595       "multiPhylo".
596
597
598 DOCUMENTATION
599
600     o There is now a vignette in ape: see vignette("MoranI", "ape").
601
602
603 DEPRECATED & DEFUNCT
604
605     o as.matching() and as.phylo.matching() do not support branch
606       lengths.
607
608     o correlogram.phylo() and discrete.dist() have been removed.
609
610
611
612                 CHANGES IN APE VERSION 2.0-2
613
614
615 NEW FEATURES
616
617     o The new function matexpo computes the exponential of a square
618       matrix.
619
620     o The new function unique.multi.tree removes duplicate trees from
621       a list.
622
623     o yule() has a new option `use.root.edge = FALSE' that specifies
624       to ignore, by default, the root edge of the tree if it exists.
625
626
627 BUG FIXES
628
629     o which.edge() failed when the index of a single terminal edge was
630       looked for.
631
632     o In diversi.time(), the values returned for model C were
633       incorrect.
634
635     o A bug was fixed in yule() that affected the calculation of the
636       likelihood in the presence of ties in the branching times.
637
638     o There was a bug in the C function mat_expo4x4 affecting the
639       calculations of the transition probabilities for models HKY and
640       GTR in mlphylo().
641
642     o A small bug was fixed in as.matrix.DNAbin (thanks to James
643       Bullard).
644
645     o rtree() did not `shuffle' the tip labels by default, so only a
646       limited number of labelled topologies could be generated.
647
648
649
650                 CHANGES IN APE VERSION 2.0-1
651
652
653 NEW FEATURES
654
655     o The three new functions bionj, fastme.ols, and fastme.bal
656       perform phylogeny estimation by the BIONJ and fastME methods in
657       OLS and balanced versions. This is a port to R of previous
658       previous programs done by Vincent Lefort.
659
660     o The new function chronoMPL performs molecular dating with the
661       mean path lengths method of Britton et al. (2002, Mol. Phyl.
662       Evol. 24: 58).
663
664     o The new function rotate, contributed by Christoph Heibl, swaps
665       two clades connected to the same node. It works also with
666       multichotomous nodes.
667
668     o The new `method' as.matrix.DNAbin() may be used to convert
669       easily DNA sequences stored in a list into a matrix while
670       keeping the names and the class.
671
672
673 BUG FIXES
674
675     o chronopl() failed when some branch lengths were equal to zero:
676       an error message is now returned.
677
678     o di2multi() failed when there was a series of consecutive edges
679       to remove.
680
681
682
683                 CHANGES IN APE VERSION 1.10-2
684
685
686 NEW FEATURES
687
688     o plot.phylo() can now plot circular trees: the option is type =
689       "fan" or type = "f" (to avoid the ambiguity with type = "c").
690
691     o prop.part() has a new option `check.labels = FALSE' which allows
692       to considerably speed-up the calculations of bipartitions. As a
693       consequence, calculations of bootstrap values with boot.phylo()
694       should be much faster.
695
696
697 BUG FIXES
698
699     o read.GenBank() did not return correctly the list of species as
700       from ape 1.10: this is fixed in this version
701
702     o Applying as.phylo() on a tree of class "phylo" failed: the
703       object is now returned unchanged.
704
705
706
707                 CHANGES IN APE VERSION 1.10-1
708
709
710 NEW FEATURES
711
712     o The three new functions Ntip, Nnode, and Nedge return, for a
713       given tree, the number of tips, nodes, or edges, respectively.
714
715
716 BUG FIXES
717
718     o read.nexus() did not set correctly the class of the returned
719       object when reading multiple trees.
720
721     o mllt.plot() failed with objects of class c("multi.tree",
722       "phylo").
723
724     o unroot() did not work correctly in most cases.
725
726     o reorder.phylo() made R freeze in some occasions.
727
728     o Plotting a tree in pruningwise order failed.
729
730     o When plotting an unrooted tree, the tip labels where not all
731       correctly positioned if the option `cex' was used.
732
733
734
735                 CHANGES IN APE VERSION 1.10
736
737
738 NEW FEATURES
739
740     o Five new `method' functions have been introduced to manipulate
741       DNA sequences in binary format (see below).
742
743     o Three new functions have been introduced to convert between the
744       new binary and the character formats.
745
746     o The new function as.alignment converts DNA sequences stored as
747       single characters into the class "alignment" used by the package
748       seqinr.
749
750     o read.dna() and read.GenBank() have a new argument `as.character'
751       controlling whether the sequences are returned in binary format
752       or as character.
753
754
755 BUG FIXES
756
757     o root() failed when the tree had node labels: this is fixed.
758
759     o plot.phylo() did not correctly set the limits on the y-axis with
760       the default setting: this is fixed.
761
762     o dist.dna() returned a wrong result for the LogDet, paralinear,
763       and BH87 models with `pairwise.deletion = TRUE'.
764
765
766 OTHER CHANGES
767
768     o DNA sequences are now internally stored in a binary format. See
769       the document "A Bit-Level Coding Scheme for Nucleotides" for the
770       details. Most functions analyzing DNA functions have been
771       modified accordingly and are now much faster (dist.dna is now
772       ca. 60 times faster).
773
774
775
776                 CHANGES IN APE VERSION 1.9-4
777
778
779 BUG FIXES
780
781     o A bug was fixed in edgelabels().
782
783     o as.phylo.hclust() did not work correctly when the object of
784       class "hclust" has its labels set to NULL: the returned tree has
785       now its tip labels set to "1", "2", ...
786
787     o consensus could fail if some tip labels are a subset of others
788       (e.g., "a" and "a_1"): this is now fixed.
789
790     o mlphylo() failed in most cases if some branch lengths of the
791       initial tree were greater than one: an error message is now
792       issued.
793
794     o mlphylo() failed in most cases when estimating the proportion of
795       invariants: this is fixed.
796
797
798
799                 CHANGES IN APE VERSION 1.9-3
800
801
802 NEW FEATURES
803
804     o The new function edgelabels adds labels on the edge of the tree
805       in the same way than nodelabels or tiplabels.
806
807
808 BUG FIXES
809
810     o multi2di() did not handle correctly branch lengths with the
811       default option `random = TRUE': this is now fixed.
812
813     o A bug was fixed in nuc.div() when using pairwise deletions.
814
815     o A bug occurred in the analysis of bipartitions with large
816       numbers of large trees, with consequences on prop.part,
817       prop.clades, and boot.phylo.
818
819     o The calculation of the Billera-Holmes-Vogtmann distance in
820       dist.topo was wrong: this has been fixed.
821
822
823
824                 CHANGES IN APE VERSION 1.9-2
825
826
827 NEW FEATURES
828
829     o The new function ladderize reorganizes the internal structure of
830       a tree to plot them left- or right-ladderized.
831
832     o The new function dist.nodes computes the patristic distances
833       between all nodes, internal and terminal, of a tree. It replaces
834       the option `full = TRUE' of cophenetic.phylo (see below).
835
836
837 BUG FIXES
838
839     o A bug was fixed in old2new.phylo().
840
841     o Some bugs were fixed in chronopl().
842
843     o The edge colours were not correctly displayed by plot.phylo
844       (thank you to Li-San Wang for the fix).
845
846     o cophenetic.phylo() failed with multichotomous trees: this is
847       fixed.
848
849
850 OTHER CHANGES
851
852     o read.dna() now returns the sequences in a matrix if they are
853       aligned (interleaved or sequential format). Sequences in FASTA
854       format are still returned in a list.
855
856     o The option `full' of cophenetic.phylo() has been removed because
857       it could not be used from the generic.
858
859
860 DEPRECATED & DEFUNCT
861
862     o rotate() has been removed; this function did not work correctly
863       since ape 1.9.
864
865
866
867                 CHANGES IN APE VERSION 1.9-1
868
869
870 BUG FIXES
871
872     o Trees with a single tip were not read correctly in R as the
873       element `Nnode' was not set: this is fixed.
874
875     o unroot() did not set correctly the number of nodes of the
876       unrooted tree in most cases.
877
878     o read.GenBank() failed when fetching very long sequences,
879       particularly of the BX-series.
880
881     o A bug was introduced in read.tree() with ape 1.9: it has been
882       fixed
883
884
885
886                 CHANGES IN APE VERSION 1.9
887
888
889 NEW FEATURES
890
891     o There are two new print `methods' for trees of class "phylo" and
892       lists of trees of class "multi.tree", so that they are now
893       displayed in a compact and informative way.
894
895     o There are two new functions, old2new.phylo and new2old.phylo,
896       for converting between the old and new coding of the class
897       "phylo".
898
899     o dist.dna() has three new models: Barry and Hartigan ("BH87"),
900       LogDet ("logdet"), and paralinear ("paralin").
901
902     o compute.brlen() has been extended: several methods are now
903       available to compute branch lengths.
904
905     o write.dna() can now handle matrices as well as lists.
906
907
908 BUG FIXES
909
910     o cophenetic.phylo() sometimes returned a wrong result with
911       multichotomous trees: this is fixed.
912
913     o rotate() failed when a single tip was specified: the tree is now
914       returned unchanged.
915
916     o ace() did not return the correct index matrix with custom
917       models: this is fixed.
918
919     o multi2di() did not work correctly when resolving multichotomies
920       randomly: the topology was always the same, only the arrangement
921       of clades was randomized: this is fixed. This function now
922       accepts trees with no branch lengths.
923
924     o The output of diversi.gof() was blurred by useless prints when a
925       user distribution was specified. This has been corrected, and
926       the help page of this function has been expanded.
927
928
929 OTHER CHANGES
930
931     o The internal structure of the class "phylo" has been changed:
932       see the document "Definition of Formats for Coding Phylogenetic
933       Trees in R" for the details. In addition, the code of most
934       functions has been improved.
935
936     o Several functions have been improved by replacing some R codes
937       by C codes: pic, plot.phylo, and reorder.phylo.
938
939     o There is now a citation information: see citation("ape") in R.
940
941     o write.tree() now does not add extra 0's to branch lengths so
942       that 1.23 is printed "1.23" by default, not "1.2300000000".
943
944     o The syntax of bind.tree() has been simplified. This function now
945       accepts trees with no branch lengths, and handles correctly node
946       labels.
947
948     o The option `as.numeric' of mrca() has been removed.
949
950     o The unused options `format' and `rooted' of read.tree() have
951       been removed.
952
953     o The unused option `format' of write.tree() has been removed.
954
955     o The use of node.depth() has been simplified.
956
957
958
959                 CHANGES IN APE VERSION 1.8-5
960
961
962 NEW FEATURES
963
964     o Two new functions read.nexus.data() and write.nexus.data(),
965       contributed by Johan Nylander, allow to read and write molecular
966       sequences in NEXUS files.
967
968     o The new function reorder.phylo() reorders the internal structure
969       of a tree of class "phylo". It is used as the generic, e.g.,
970       reorder(tr).
971
972     o read.tree() and read.nexus() can now read trees with a single
973       edge.
974
975     o The new data set `cynipids' supplies a set of protein sequences
976       in NEXUS format.
977
978
979 BUG FIXES
980
981     o The code of all.equal.phylo() has been completely rewritten
982       (thanks to Benoît Durand) which fixes several bugs.
983
984     o read.tree() and read.nexus() now checks the labels of the tree
985       to remove or substitute any characters that are illegal in the
986       Newick format (parentheses, etc.)
987
988     o A negative P-value could be returned by mantel.test(): this is
989       now fixed.
990
991
992
993                 CHANGES IN APE VERSION 1.8-4
994
995
996 NEW FEATURES
997
998     o The new function sh.test() computes the Shimodaira-
999       Hasegawa test.
1000
1001     o The new function collapse.singles() removes the nodes with a
1002       single descendant from a tree.
1003
1004     o plot.phylo() has a new argument `tip.color' to specify the
1005       colours of the tips.
1006
1007     o mlphylo() has now an option `quiet' to control the display of
1008       the progress of the analysis (the default is FALSE).
1009
1010
1011 BUG FIXES
1012
1013     o read.dna() did not read correctly sequences in sequential format
1014       with leading alignment gaps "-": this is fixed.
1015
1016     o ace() returned a list with no class so that the generic
1017       functions (anova, logLik, ...) could not be used directly. This
1018       is fixed as ace() now returns an object of class "ace".
1019
1020     o anova.ace() had a small bug when computing the number of degrees
1021       of freedom: this is fixed.
1022
1023     o mlphylo() did not work when the sequences were in a matrix or
1024       a data frame: this is fixed.
1025
1026     o rtree() did not work correctly when trying to simulate an
1027       unrooted tree with two tips: an error message is now issued.
1028
1029
1030 OTHER CHANGES
1031
1032     o The algorithm of rtree() has been changed: it is now about 40,
1033       100, and 130 times faster for 10, 100, and 1000 tips,
1034       respectively.
1035
1036
1037
1038                 CHANGES IN APE VERSION 1.8-3
1039
1040
1041 NEW FEATURES
1042
1043     o There are four new `method' functions to be used with the
1044       results of ace(): logLik(), deviance(), AIC(), and anova().
1045
1046     o The plot method of phymltest has two new arguments: `main' to
1047       change the title, and `col' to control the colour of the
1048       segments showing the AIC values.
1049
1050     o ace() has a new argument `ip' that gives the initial values used
1051       in the ML estimation with discrete characters (see the examples
1052       in ?ace). This function now returns a matrix giving the indices
1053       of the estimated rates when analysing discrete characters.
1054
1055     o nodelabels() and tiplabels() have a new argument `pie' to
1056       represent proportions, with any number of categories, as
1057       piecharts. The use of the option `thermo' has been improved:
1058       there is now no limitation on the number of categories.
1059
1060
1061 BUG FIXES
1062
1063     o mlphylo() did not work with more than two partitions: this is
1064       fixed.
1065
1066     o root() failed if the proposed outgroup was already an outgroup
1067       in the tree: this is fixed.
1068
1069     o The `col' argument in nodelabels() and tiplabels() was not
1070       correctly passed when `text' was used: this is fixed.
1071
1072     o Two bugs were fixed in mlphylo(): parameters were not always
1073       correctly output, and the estimation failed in some cases.
1074
1075     o plot.phylo() was stuck when given a tree with a single tip: this
1076       is fixed and a message error is now returned.
1077
1078     o An error was corrected in the help page of gammaStat regarding
1079       the calculation of P-values.
1080
1081     o Using gls() could crash R when the number of species in the tree
1082       and in the variables were different: this is fixed.
1083
1084
1085
1086                 CHANGES IN APE VERSION 1.8-2
1087
1088
1089 NEW FEATURES
1090
1091     o The new function mlphylo() fits a phylogenetic tree by maximum
1092       likelihood from DNA sequences. Its companion function DNAmodel()
1093       is used to define the substitution model which may include
1094       partitioning. There are methods for logLik(), deviance(), and
1095       AIC(), and the summary() method has been extended to display in
1096       a friendly way the results of this model fitting. Currently, the
1097       functionality is limited to estimating the substitution and
1098       associated parameters and computing the likelihood.
1099
1100     o The new function drop1.compar.gee (used as, e.g., drop1(m))
1101       tests for single effects in GEE-based comparative method. A
1102       warning message is printed if there is not enough degrees of
1103       freedom.
1104
1105
1106 BUG FIXES
1107
1108     o An error message was sometimes issued by plot.multi.tree(),
1109       though with no consequence.
1110
1111
1112
1113                 CHANGES IN APE VERSION 1.8-1
1114
1115
1116 NEW FEATURES
1117
1118     o There is a new plot method for lists of trees (objects of class
1119       "multi.tree"): it calls plot.phylo() internally and is
1120       documented on the same help page.
1121
1122
1123 BUG FIXES
1124
1125     o A bug was fixed in the C code that analyzes bipartitions: this
1126       has impact on several functions like prop.part, prop.clades,
1127       boot.phylo, or consensus.
1128
1129     o root() did not work correctly when the specified outgroup had
1130       more than one element: this is fixed.
1131
1132     o dist.dna() sometimes returned a warning inappropriately: this
1133       has been corrected.
1134
1135     o If the distance object given to nj() had no rownames, nj()
1136       returned a tree with no tip labels: it now returns tips labelled
1137       "1", "2", ..., corresponding to the row numbers.
1138
1139
1140 OTHER CHANGES
1141
1142     o nj() has been slightly changed so that tips with a zero distance
1143       are first aggregated with zero-lengthed branches; the usual NJ
1144       procedure is then performed on a distance matrix without 0's.
1145
1146
1147
1148                 CHANGES IN APE VERSION 1.8
1149
1150
1151 NEW FEATURES
1152
1153     o The new function chronopl() estimates dates using the penalized
1154       likelihood method by Sanderson (2002; Mol. Biol. Evol., 19:101).
1155
1156     o The new function consensus() calculates the consensus tree of a
1157       list of trees.
1158
1159     o The new function evolve.phylo() simulates the evolution of
1160       continuous characters along a phylogeny under a Brownian model.
1161
1162     o The new plot method for objects of class "ancestral" displays a
1163       tree together with ancestral values, as returned by the above
1164       function.
1165
1166     o The new function as.phylo.formula() returns a phylogeny from a
1167       set of nested taxonomic variables given as a formula.
1168
1169     o The new function read.caic() reads trees in CAIC format.
1170
1171     o The new function tiplabels() allows to add labels to the tips
1172       of a tree using text or plotting symbols in a flexible way.
1173
1174     o The new function unroot() unroots a phylogeny.
1175
1176     o multi2di() has a new option, `random', which specifies whether
1177       to resolve the multichotomies randomly (the default) or not.
1178
1179     o prop.part() now returns an object of class "prop.part" for which
1180       there are print (to display a partition in a more friendly way)
1181       and summary (to extract the numbers) methods.
1182
1183     o plot.phylo() has a new option, `show.tip.label', specifying
1184       whether to print the labels of the tips. The default is TRUE.
1185
1186     o The code of nj() has been replaced by a faster C code: it is now
1187       about 10, 25, and 40 times faster for 50, 100, and 200 taxa,
1188       respectively.
1189
1190     o write.nexus() now writes whether a tree is rooted or not.
1191
1192
1193 BUG FIXES
1194
1195     o Two bugs have been fixed in root(): unrooted trees are now
1196       handled corretly, and node labels are now output normally.
1197
1198     o A bug was fixed in phymltest(): the executable couldn't be found
1199       in some cases.
1200
1201     o Three bug have been fixed in ace(): computing the likelihood of
1202       ancestral states of discrete characters failed, custom models
1203       did not work, and the function failed with a null gradient (a
1204       warning message is now returned; this latter bug was also
1205       present in yule.cov() as well and is now fixed).
1206
1207     o pic() hanged out when missing data were present: a message error
1208       is now returned.
1209
1210     o A small bug was fixed in dist.dna() where the gamma correction
1211       was not always correctly dispatched.
1212
1213     o plot.phylo() plotted correctly the root edge only when the tree
1214       was plotted rightwards: this works now for all directions.
1215
1216
1217 OTHER CHANGES
1218
1219     o dist.taxo() has been renamed as weight.taxo().
1220
1221     o Various error and warning messages have been improved.
1222
1223
1224
1225                 CHANGES IN APE VERSION 1.7
1226 NEW FEATURES
1227
1228     o The new function ace() estimates ancestral character states for
1229       continuous characters (with ML, GLS, and contrasts methods), and
1230       discrete characters (with ML only) for any number of states.
1231
1232     o The new function compar.ou() fits the Ornstein-Uhlenbeck model
1233       of directional evolution for continuous characters. The user
1234       specifies the node(s) of the tree where the character optimum
1235       changes.
1236
1237     o The new function is.rooted() tests whether a tree (of class
1238       "phylo") is rooted.
1239
1240     o The new function rcoal() generates random ultrametric trees with
1241       the possibility to specify the function that generates the
1242       inter-nodes distances.
1243
1244     o The new function mrca() gives for all pairs of tips in a tree
1245       (and optionally nodes too) the most recent common ancestor.
1246
1247     o nodelabels() has a new option `thermo' to plot proportions (up
1248       to three classes) on the nodes of a tree.
1249
1250     o rtree() has been improved: it can now generate rooted or
1251       unrooted trees, and the mathematical function that generates the
1252       branch lengths may be specified by the user. The tip labels may
1253       be given directly in the call to rtree. The limit cases (n = 2,
1254       3) are now handled correctly.
1255
1256     o dist.topo() has a new argument `method' with two choices: "PH85"
1257       for Penny and Henny's method (already available before and now
1258       the default), and "BHV01" for the geometric distance by Billera
1259       et al. (2001, Adv. Appl. Math. 27:733).
1260
1261     o write.tree() has a new option, `digits', which specifies the
1262       number of digits to be printed in the Newick tree. By default
1263       digits = 10. The numbers are now always printed in decimal form
1264       (i.e., 1.0e-1 is now avoided).
1265
1266     o dist.dna() can now compute the raw distances between pairs of
1267       DNA sequences by specifying model = "raw".
1268
1269     o dist.phylo() has a new option `full' to possibly compute the
1270       distances among all tips and nodes of the tree. The default if
1271       `full = FALSE'.
1272
1273
1274 BUG FIXES
1275
1276     o Several bugs were fixed in all.equal.phylo().
1277
1278     o dist.dna() did not handle correctly gaps ("-") in alignments:
1279       they are now considered as missing data.
1280
1281     o rotate() did not work if the tips were not ordered: this is
1282       fixed.
1283
1284     o mantel.test() returned NA in some special cases: this is fixed
1285       and the function has been improved and is now faster.
1286
1287     o A bug was fixed in diversi.gof() where the calculation of A² was
1288       incorrect.
1289
1290     o cherry() did not work correctly under some OSs (mainly Linux):
1291       this is fixed.
1292
1293     o is.binary.tree() has been modified so that it works with both
1294       rooted and unrooted trees.
1295
1296     o The documentation of theta.s() was not correct: this has been
1297       fixed.
1298
1299     o plot.mst() did not work correctly: this is fixed.
1300
1301
1302
1303                 CHANGES IN APE VERSION 1.6
1304
1305
1306 NEW FEATURES
1307
1308     o The new function dist.topo() computes the topological distances
1309       between two trees.
1310
1311     o The new function boot.phylo() performs a bootstrap analysis on
1312       phylogeny estimation.
1313
1314     o The new functions prop.part() and prop.clades() analyse
1315       bipartitions from a series of trees.
1316
1317
1318 OTHER CHANGES
1319
1320     o read.GenBank() now uses the EFetch utility of NCBI instead of
1321       the usual Web interface: it is now much faster (e.g., 12 times
1322       faster to retrieve 8 sequences, 37 times for 60 sequences).
1323
1324
1325 BUG FIXES
1326
1327     o Several bugs were fixed in read.dna().
1328
1329     o Several bugs were fixed in diversi.time().
1330
1331     o is.binary.tree() did not work correctly if the tree has no edge
1332       lengths: this is fixed.
1333
1334     o drop.tip() did not correctly propagated the `node.label' of a
1335       tree: this is fixed.
1336
1337
1338
1339                 CHANGES IN APE VERSION 1.5
1340
1341
1342 NEW FEATURES
1343
1344     o Two new functions, as.matching.phylo() and as.phylo.matching(),
1345       convert objects between the classes "phylo" and "matching". The
1346       latter implements the representation of binary trees introduced by
1347       Diaconis and Holmes (1998; PNAS 95:14600). The generic function
1348       as.matching() has been introduced as well.
1349
1350     o Two new functions, multi2di() and di2multi(), allow to resolve
1351       and collapse multichotomies with branches of length zero.
1352
1353     o The new function nuc.div() computes the nucleotide diversity
1354       from a sample a DNA sequences.
1355
1356     o dist.dna() has been completely rewritten with a much faster
1357       (particularly for large data sets) C code. Eight models are
1358       available: JC69, K80, F81, K81, F84, T92, TN93, and GG95 (the
1359       option `method' has been renamed `model'). Computation of variance
1360       is available for all models. A gamma-correction is possible for
1361       JC69, K80, F81, and TN93. There is a new option, pairwise.deletion,
1362       to remove sites with missing data on a pairwise basis. The option
1363       `GCcontent' has been removed.
1364
1365     o read.GenBank() has a new option (species.names) which specifies
1366       whether to return the species names of the organisms in addition
1367       to the accession numbers of the sequences (this is the default
1368       behaviour).
1369
1370     o write.nexus() can now write several trees in the same NEXUS file.
1371
1372     o drop.tip() has a new option `root.edge' that allows to specify the
1373       new root edge if internal branches are trimmed.
1374
1375
1376 BUG FIXES
1377
1378     o as.phylo.hclust() failed if some labels had parentheses: this
1379       is fixed.
1380
1381     o Several bugs were fixed in all.equal.phylo(). This function now
1382       returns the logical TRUE if the trees are identical but with
1383       different representations (a report was printed previously).
1384
1385     o read.GenBank() did not correctly handle ambiguous base codes:
1386       this is fixed.
1387
1388
1389 OTHER CHANGES
1390
1391     o birthdeath() now returns an object of class "birthdeath" for
1392       which there is a print method.
1393
1394
1395
1396                 CHANGES IN APE VERSION 1.4
1397
1398
1399 NEW FEATURES
1400
1401     o The new function nj() performs phylogeny estimation with the
1402       neighbor-joining method of Saitou and Nei (1987; Mol. Biol.
1403       Evol., 4:406).
1404
1405     o The new function which.edge() identifies the edges of a tree
1406       that belong to a group specified as a set of tips.
1407
1408     o The new function as.phylo.phylog() converts an object of class
1409       "phylog" (from the package ade4) into an object of class
1410       "phylo".
1411
1412     o The new function axisPhylo() draws axes on the side of a
1413       phylogeny plot.
1414
1415     o The new function howmanytrees() calculates the number of trees
1416       in different cases and giving a number of tips.
1417
1418     o write.tree() has a new option `multi.line' (TRUE by default) to
1419       write a Newick tree on several lines rather than on a single
1420       line.
1421
1422     o The functionalities of zoom() have been extended. Several
1423       subtrees can be visualized at the same time, and they are marked
1424       on the main tree with colors. The context of the subtrees can be
1425       marked with the option `subtree' (see below).
1426
1427     o drop.tip() has a new option `subtree' (FALSE by default) which
1428       specifies whether to output in the tree how many tips have been
1429       deleted and where.
1430
1431     o The arguments of add.scale.bar() have been redefined and have
1432       now default values (see ?add.scale.bar for details). This
1433       function now works even if the plotted tree has no edge length.
1434
1435     o plot.phylo() can now plot radial trees, but this does not take
1436       edge lengths into account.
1437
1438     o In plot.phylo() with `type = "phylogram"', if the values of
1439       `edge.color' and `edge.width' are identical for sister-branches,
1440       they are propagated to the vertical line that link them.
1441
1442
1443 BUG FIXES
1444
1445     o Repeated calls to as.phylo.hclust() or as.hclust.phylo() made R
1446       crashing. This is fixed.
1447
1448     o In plot.phylo(), the options `edge.color' and `edge.width' are
1449       now properly recycled; their default values are now "black" and
1450       1, respectively.
1451
1452     o A bug has been fixed in write.nexus().
1453
1454
1455 OTHER CHANGES
1456
1457     o The function node.depth.edgelength() has been removed and
1458       replaced by a C code.
1459
1460
1461
1462                 CHANGES IN APE VERSION 1.3-1
1463
1464
1465 NEW FEATURES
1466
1467     o The new function nodelabels() allows to add labels to the nodes
1468       of a tree using text or plotting symbols in a flexible way.
1469
1470     o In plot.phylo() the arguments `x.lim' and `y.lim' can now be two
1471       numeric values specifying the lower and upper limits on the x-
1472       and y-axes. This allows to leave some space on any side of the
1473       tree. If a single value is given, this is taken as the upper
1474       limit (as before).
1475
1476
1477
1478                 CHANGES IN APE VERSION 1.3
1479
1480
1481 NEW FEATURES
1482
1483     o The new function phymltest() calls the software PHYML and fits
1484       28 models of DNA sequence evolution. There are a print method to
1485       display likelihood and AIC values, a summary method to compute
1486       the hierarchical likelihood ratio tests, and a plot method to
1487       display graphically the AIC values of each model.
1488
1489     o The new function yule.cov() fits the Yule model with covariates,
1490       a model where the speciation rate is affected by several species
1491       traits through a generalized linear model. The parameters are
1492       estimated by maximum likelihood.
1493
1494     o Three new functions, corBrownian(), corGrafen(), and
1495       corMartins(), compute the expected correlation structures among
1496       species given a phylogeny under different models of evolution.
1497       These can be used for GLS comparative phylogenetic methods (see
1498       the examples). There are coef() and corMatrix() methods and an
1499       Initialize.corPhyl() function associated.
1500
1501     o The new function compar.cheverud() implements Cheverud et al.'s
1502       (1985; Evolution 39:1335) phylogenetic comparative method.
1503
1504     o The new function varcomp() estimates variance components; it has
1505       a plot method.
1506
1507     o Two new functions, panel.superpose.correlogram() and
1508       plot.correlogramList(), allow to plot several phylogenetic
1509       correlograms.
1510
1511     o The new function node.leafnumber() computes the number of leaves
1512       of a subtree defined by a particular node.
1513
1514     o The new function node.sons() gets all tags of son nodes from a
1515       given parent node.
1516
1517     o The new function compute.brlen() computes the branch lengths of
1518       a tree according to a specified method.
1519
1520     o plot.phylo() has three new options: "cex" controls the size of
1521       the (tip and node) labels (thus it is no more needed to change
1522       the global graphical parameter), "direction" which allows to
1523       plot the tree rightwards, leftwards, upwards, or downwards, and
1524       "y.lim" which sets the upper limit on the y-axis.
1525
1526
1527 BUG FIXES
1528
1529     o Some functions which try to match tip labels and names of
1530       additional data (e.g. vector) are likely to fail if there are
1531       typing or syntax errors. If both series of names do not perfectly
1532       match, they are ignored and a warning message is now issued.
1533       These functions are bd.ext, compar.gee, pic. Their help pages
1534       have been clarified on this point.
1535
1536
1537
1538                 CHANGES IN APE VERSION 1.2-7
1539
1540
1541 NEW FEATURES
1542
1543     o The new function root() reroots a phylogenetic tree with respect
1544       to a specified outgroup.
1545
1546     o The new function rotate() rotates an internal branch of a tree.
1547
1548     o In plot.phylo(), the new argument "lab4ut" (labels for unrooted
1549       trees) controls the display of the tip labels in unrooted trees.
1550       This display has been greatly improved: the tip labels are now not
1551       expected to overlap with the tree (particularly if lab4ut =
1552       "axial"). In all cases, combining appropriate values of "lab4ut"
1553       and the font size (via "par(cex = )") should result in readable
1554       unrooted trees. See ?plot.phylo for some examples.
1555
1556     o In drop.tip(), the argument `tip' can now be numeric or character.
1557
1558
1559 BUG FIXES
1560
1561     o drop.tip() did not work correctly with trees with no branch
1562       lengths: this is fixed.
1563
1564     o A bug in plot.phylo(..., type = "unrooted") made some trees being
1565       plotted with some line crossings: this is now fixed.
1566
1567
1568
1569                 CHANGES IN APE VERSION 1.2-6
1570
1571
1572 NEW FEATURES
1573
1574     o Six new functions (Moran.I, correlogram.formula, discrete.dist,
1575       correlogram.phylo, dist.taxo, plot.correlogram) have been added
1576       to implement comparative methods with an autocorrelation approach.
1577
1578     o A new data set describing some life history traits of Carnivores
1579       has been included.
1580
1581
1582 BUG FIXES
1583
1584     o A fix was made on mcmc.popsize() to conform to R 2.0.0.
1585
1586
1587 OTHER CHANGES
1588
1589     o When plotting a tree with plot.phylo(), the new default of the
1590       option `label.offset' is now 0, so the labels are always visible.
1591
1592
1593
1594                 CHANGES IN APE VERSION 1.2-5
1595
1596
1597 NEW FEATURES
1598
1599     o The new function bd.ext() fits a birth-death model with combined
1600       phylogenetic and taxonomic data, and estimates the corresponding
1601       speciation and extinction rates.
1602
1603
1604 OTHER CHANGES
1605
1606     o The package gee is no more required by ape but only suggested
1607       since only the function compar.gee() calls gee.
1608
1609
1610
1611                 CHANGES IN APE VERSION 1.2-4
1612
1613
1614 NEW FEATURES
1615
1616     o Four new functions (mcmc.popsize, extract.popsize, plot.popsize,
1617       and lines.popsize) implementing a new approach for inferring the
1618       demographic history from genealogies using a reversible jump
1619       MCMC have been introduced.
1620
1621     o The unit of time in the skyline plot and in the new plots can
1622       now be chosen to be actual years, rather than substitutions.
1623
1624
1625
1626                 CHANGES IN APE VERSION 1.2-3
1627
1628
1629 NEW FEATURES
1630
1631     o The new function rtree() generates a random binary tree with or
1632       without branch lengths.
1633
1634     o Two new functions for drawing lineages-through-time (LTT) plots
1635       are provided: ltt.lines() adds a LTT curve to an existing plot,
1636       and mltt.plot() does a multiple LTT plot giving several trees as
1637       arguments (see `?ltt.plot' for details).
1638
1639
1640 BUG FIXES
1641
1642     o Some taxon names made R crashing when calling as.phylo.hclust():
1643       this is fixed.
1644
1645     o dist.dna() returned an error with two identical DNA sequences
1646       (only using the Jukes-Cantor method returned 0): this is fixed.
1647
1648
1649 OTHER CHANGES
1650
1651     o The function dist.phylo() has been re-written using a different
1652       algorithm: it is now about four times faster.
1653
1654     o The code of branching.times() has been improved: it is now about
1655       twice faster.
1656
1657
1658
1659                 CHANGES IN APE VERSION 1.2-2
1660
1661
1662 NEW FEATURES
1663
1664     o The new function seg.sites() finds the segregating sites in a
1665       sample of DNA sequences.
1666
1667
1668 BUG FIXES
1669
1670     o A bug introduced in read.tree() and in read.nexus() with version
1671       1.2-1 was fixed.
1672
1673     o A few errors were corrected and a few examples were added in the
1674       help pages.
1675
1676
1677
1678                 CHANGES IN APE VERSION 1.2-1
1679
1680
1681 NEW FEATURES
1682
1683     o plot.phylo() can now draw the edge of the root of a tree if it
1684       has one (see the new option `root.edge', its default is FALSE).
1685
1686
1687 BUG FIXES
1688
1689     o A bug was fixed in read.nexus(): files with semicolons inside
1690       comment blocks were not read correctly.
1691
1692     o The behaviour of read.tree() and read.nexus() was corrected so
1693       that tree files with badly represented root edges (e.g., with
1694       an extra pair of parentheses, see the help pages for details)
1695       are now correctly represented in the object of class "phylo";
1696       a warning message is now issued.
1697
1698
1699
1700                 CHANGES IN APE VERSION 1.2
1701
1702
1703 NEW FEATURES
1704
1705     o plot.phylo() has been completely re-written and offers several
1706       new functionalities. Three types of trees can now be drawn:
1707       phylogram (as previously), cladogram, and unrooted tree; in
1708       all three types the branch lengths can be drawn using the edge
1709       lengths of the phylogeny or not (e.g., if the latter is absent).
1710       The vertical position of the nodes can be adjusted with two
1711       choices (see option `node.pos'). The code has been re-structured,
1712       and two new functions (potentially useful for developpers) are
1713       documented separately: node.depth.edgelength() and node.depth();
1714       see the respective help pages for details.
1715
1716     o The new function zoom() allows to explore very large trees by
1717       focusing on a small portion of it.
1718
1719     o The new function yule() fits by maximum likelihood the Yule model
1720       (birth-only process) to a phylogenetic tree.
1721
1722     o Support for writing DNA sequences in FASTA format has been
1723       introduced in write.dna() (support for reading sequences in
1724       this format was introduced in read.dna() in version 1.1-2).
1725       The function has been completely re-written, fixing some bugs
1726       (see below); the default behaviour is no more to display the
1727       sequences on the standard output. Several options have been
1728       introduced to control the sequence printing in a flexible
1729       way. The help page has been extended.
1730
1731     o A new data set is included: a supertree of bats in NEXUS format.
1732
1733
1734 BUG FIXES
1735
1736     o In theta.s(), the default of the option `variance' has
1737       been changed to `FALSE' (as was indicated in the help page).
1738
1739     o Several bugs were fixed in the code of all.equal.phylo().
1740
1741     o Several bugs were fixed in write.dna(), particularly this
1742       function did not work with `format = "interleaved"'.
1743
1744     o Various errors were corrected in the help pages.
1745
1746
1747 OTHER CHANGES
1748
1749     o The argument names of as.hclust.phylo() have been changed
1750       from "(phy)" to "(x, ...)" to conform to the definition of
1751       the corresponding generic function.
1752
1753     o gamma.stat() has been renamed gammaStat() to avoid confusion
1754       since gamma() is a generic function.
1755
1756
1757
1758                 CHANGES IN APE VERSION 1.1-3
1759
1760
1761 BUG FIXES
1762
1763     o base.freq() previously did not return a value of 0 for
1764       bases absent in the data (e.g., a vector of length 3 was
1765       returned if one base was absent). This is now fixed (a
1766       vector of length 4 is always returned).
1767
1768     o Several bugs were fixed in read.nexus(), including that this
1769       function did not work in this absence of a "TRANSLATE"
1770       command in the NEXUS file, and that the commands were
1771       case-sensitive.
1772
1773
1774
1775                 CHANGES IN APE VERSION 1.1-2
1776
1777
1778 NEW FEATURES
1779
1780     o The Tamura and Nei (1993) model of DNA distance is now implemented
1781       in dist.dna(): five models are now available in this function.
1782
1783     o A new data set is included: a set of 15 sequences of the
1784       cytochrome b mitochondrial gene of the woodmouse (Apodemus
1785       sylvaticus).
1786
1787
1788 BUG FIXES
1789
1790     o A bug in read.nexus() was fixed.
1791
1792     o read.dna() previously did not work correctly in most cases.
1793       The function has been completely re-written and its help page
1794       has been considerably extended (see ?read.dna for details).
1795       Underscores (_) in taxon names are no more replaced with
1796       spaces (this behaviour was undocumented).
1797
1798     o A bug was fixed in write.dna().
1799
1800
1801
1802                 CHANGES IN APE VERSION 1.1-1
1803
1804
1805 BUG FIXES
1806
1807     o A bug in read.tree() introduced in APE 1.1 was fixed.
1808
1809     o A bug in compar.gee() resulted in an error when trying to fit
1810       a model with `family = "binomial"'. This is now fixed.
1811
1812
1813
1814                 CHANGES IN APE VERSION 1.1
1815
1816
1817 NEW FEATURES
1818
1819     o The Klastorin (1982) method as suggested by Misawa and Tajima
1820       (2000, Mol. Biol. Evol. 17:1879-1884) for classifying genes
1821       on the basis of phylogenetic trees has been implemented (see
1822       the function klastorin()).
1823
1824     o Functions have been added to convert APE's "phylo" objects in
1825       "hclust" cluster objects and vice versa (see the help page of
1826       as.phylo for details).
1827
1828     o Three new functions, ratogram(), chronogram() and NPRS.criterion(),
1829       are introduced for the estimation of absolute evolutionary rates
1830       (ratogram) and dated clock-like trees (chronogram) from
1831       phylogenetic trees using the non-parametric rate smoothing approach
1832       by MJ Sanderson (1997, Mol. Biol. Evol. 14:1218-1231).
1833
1834     o A summary method is now provided printing a summary information on a
1835       phylogenetic tree with, for instance, `summary(tree)'.
1836
1837     o The behaviour of read.tree() was changed so that all spaces and
1838       tabulations in tree files are now ignored. Consequently, spaces in tip
1839       labels are no more allowed. Another side effect is that read.nexus()
1840       now does not replace the underscores (_) in tip labels with spaces
1841       (this behaviour was undocumented).
1842
1843     o The function plot.phylo() has a new option (`underscore') which
1844       specifies whether the underscores in tip labels should be written on
1845       the plot as such or replaced with spaces (the default).
1846
1847     o The function birthdeath() now computes 95% confidence intervals of
1848       the estimated parameters using profile likelihood.
1849
1850     o Three new data sets are included: a gene tree estimated from 36
1851       landplant rbcL sequences, a gene tree estimated from 32 opsin
1852       sequences, and a gene tree for 50 BRCA1 mammalian sequences.
1853
1854
1855 BUG FIXES
1856
1857     o A bug was fixed in dist.gene() where nothing was returned.
1858
1859     o A bug in plot.mst() was fixed.
1860
1861     o A bug in vcv.phylo() resulted in false correlations when the
1862       option `cor = TRUE' was used (now fixed).
1863
1864
1865
1866                 CHANGES IN APE VERSION 1.0
1867
1868
1869 NEW FEATURES
1870
1871     o Two new functions, read.dna() and write.dna(), read/write in a file
1872       DNA sequences in interleaved or in sequential format.
1873
1874     o Two new functions, read.nexus() and write.nexus(), read/write trees
1875       in a NEXUS file.
1876
1877     o The new function bind.tree() allows to bind two trees together,
1878       possibly handling root edges to give internal branches.
1879
1880     o The new function drop.tip() removes the tips in a phylogenetic tree,
1881       and trims (or not) the corresponding internal branches.
1882
1883     o The new function is.ultrametric() tests if a tree is ultrametric.
1884
1885     o The function plot.phylo() has more functionalities such as drawing the
1886       branches with different colours and/or different widths, showing the
1887       node labels, controling the position and font of the labels, rotating
1888       the labels, and controling the space around the plot.
1889
1890     o The function read.tree() can now read trees with no branch length,
1891       such as "(a,b),c);". Consequently, the element `edge.length' in
1892       objects of class "phylo" is now optional.
1893
1894     o The function write.tree() has a new default behaviour: if the default
1895       for the option `file' is used (i.e. file = ""), then a variable of
1896       mode character containing the tree in Newick format is returned which
1897       can thus be assigned (e.g., tree <- write.tree(phy)).
1898
1899     o The function read.tree() has a new argument `text' which allows
1900       to read the tree in a variable of mode character.
1901
1902     o A new data set is included: the phylogenetic relationships among
1903       the orders of birds from Sibley and Ahlquist (1990).
1904
1905
1906
1907                 CHANGES IN APE VERSION 0.2-1
1908
1909
1910 BUG FIXES
1911
1912     o Several bugs were fixed in the help pages.
1913
1914
1915
1916                 CHANGES IN APE VERSION 0.2
1917
1918
1919 NEW FEATURES
1920
1921     o The function write.tree() writes phylogenetic trees (objects of class
1922       "phylo") in an ASCII file using the Newick parenthetic format.
1923
1924     o The function birthdeath() fits a birth-death model to branching times
1925       by maximum likelihood, and estimates the corresponding speciation and
1926       extinction rates.
1927
1928     o The function scale.bar() adds a scale bar to a plot of a phylogenetic
1929       tree.
1930
1931     o The function is.binary.tree() tests whether a phylogeny is binary.
1932
1933     o Two generic functions, coalescent.intervals() and collapsed.intervals(),
1934       as well as some methods are introduced.
1935
1936     o Several functions, including some generics and methods, for computing
1937       skyline plot estimates (classic and generalized) of effective
1938       population size through time are introduced and replace the function
1939       skyline.plot() in version 0.1.
1940
1941     o Two data sets are now included: the phylogenetic relationships among
1942       the families of birds from Sibley and Ahlquist (1990), and an
1943       estimated clock-like phylogeny of HIV sequences sampled in the
1944       Democratic Republic of Congo.
1945
1946
1947 DEPRECATED & DEFUNCT
1948
1949     o The function skyline.plot() in ape 0.1 has been deprecated and
1950       replaced by more elaborate functions (see above).
1951
1952
1953 BUG FIXES
1954
1955     o Two important bugs were fixed in plot.phylo(): phylogenies with
1956       multichotomies not at the root or not with only terminal branches,
1957       and phylogenies with a single node (i.e. only terminal branches)
1958       did not plot. These trees should be plotted correctly now.
1959
1960     o Several bugs were fixed in diversi.time() in the computation of
1961       AICs and LRTs.
1962
1963     o Various errors were corrected in the help pages.