Re: [R-sig-phylo] ape: chronos

2016-09-05 Thread Emmanuel Paradis
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

[R-sig-phylo] ape: chronos

2016-09-05 Thread Daniel Lang
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

Re: [R-sig-phylo] subset DNAbin

2016-09-05 Thread Andreas Kolter
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