Hi Daniel,
Thank you for your report. It's not very clear how or why the penalized
likelihood becomes NA (maybe a log of a negative number is calculated
somewhere). chronos() is still under progress and I take good note of
your fix.
Best,
Emmanuel
Le 05/09/2016 à 10:27, Daniel Lang a
Hi,
I hope I did not overlook ape's bug tracker - which would be the most
appropriate place to post the following:
I ran into this error while dating a tree using apes' chronos function
(discrete model):
Error in if (new.ploglik - current.ploglik > 1e-06 && i <=
dual.iter.max) { :
missing
Hello Kirston,
try this:
x <- read.fas("example.fasta")
y <- x[grep("1011",rownames(x)),]
Greetings,
Andreas
Am 2016-09-05 08:07, schrieb Kirston Barton:
Hi,
I have my data in a fasta file and am importing it into R using
read.dna, which creates a DNAbin matrix object. I would like to subset