Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-29 Thread Roland Sookias
Hi Having another problem with the same function now. What you suggested worked for one tree, but for the attached tree and file it didn't... I get: ttree-date.phylo(tree, ages, rlen=10, method=equal) Error: subscript out of bounds traceback() 1: date.phylo(tree, ages, rlen = 10, method =

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Emmanuel Paradis
Hi Roland, Just after the error, type: traceback() This will show you, hopefully, where the error happened. Best, Emmanuel -Original Message- From: Roland Sookias r.sook...@gmail.com Sender: r-sig-phylo-boun...@r-project.org Date: Thu, 21 Jul 2011 12:44:30 To:

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Klaus Schliep
Dear Roland, it would be good if you add the datasets, that one can reproduce your results (archotreeresolved, ages). I would guess from the error message that your ages list may has the wrong format. Also traceback() sometimes tells you where the error happens exactly. Kind regards Klaus On

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Roland Sookias
Hi Thanks guys. Attached are the two files (will these work via the list?) When I type trackback() all I get is 1: date.phylo(archotreeresolved, ages, rlen = 1, method = equal). Roland On Thu, Jul 21, 2011 at 1:24 PM, Klaus Schliep klaus.schl...@gmail.comwrote: Dear Roland, it would be

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Klaus Schliep
Hi Roland, as I suspected, your ages list was not in the right format (has nothing to do with ape). Try this: tmp = read.csv(agescut.csv, header=FALSE) ages = matrix(tmp[, 2], ncol=1) rownames(ages) = tmp[,1] tree = read.tree(archotreeresolved6.tre) ttree-date.phylo(tree, ages, rlen=1,

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Roland Sookias
Great! It worked. What format should the file be in then? I thought I'd followed the format specified.. Thanks very much indeed :) On Thu, Jul 21, 2011 at 4:10 PM, Klaus Schliep klaus.schl...@gmail.comwrote: Hi Roland, as I suspected, your ages list was not in the right format (has nothing

Re: [R-sig-phylo] Dating trees using date.phylo - ape problem?

2011-07-21 Thread Klaus Schliep
Hi Roland! ages should be a matrix, with row names and the ages in the first column. On 7/21/11, Roland Sookias r.sook...@gmail.com wrote: Great! It worked. What format should the file be in then? I thought I'd followed the format specified.. Thanks very much indeed :) On Thu, Jul 21, 2011 at