Dear Sergei,

I was wondering if there was a bf to get substitution rates and divtimes from heterochronous sequences, similar to BEAST (Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login) - using a global clock or relaxing it.

Would your TipDate-type bf do the job? I read the two TipDate-related posts, but you say "dated contemporaneous sequences" - do you mean that your TipDate bf estimates MRCA dates from contemporaneous sequences or you just meant "non-contemporaneous" instead? If it's the latter, would it be possible to get substitution rates, brlengths and proceed to a RR test?

Congrats for such a great program (and server)!

Thanx for your time,
sergios

Dear Sergios,

I did mean to say dated NON-contemporanous sequences. That analysis is effectively an updated and more model-rich version of the global clock test a-la TipDate. We have some 'beta' Genetic Algorithm analyses which can do the test allowing a relaxed clock (i.e. multiple rates). Let me know if you are interested and I can e-mail you the files. Eventually (in the next month or so), I am actually going to add dating type analyses to datamonkey.org as well.

Cheers,
Sergei

Dear Sergei,

thank you for your reply. Yes, I'd be interested to check out the relaxed TipDate-type bf.

Is there a readme for the TipDate bf (input file, how to get substitution rates, MRCA divtimes)? I think I saw something on a post some time ago, but I can't seem to find it.

Thanks again,
sergios

Dear Sergious,

I updated the DatedTipsMolecularClock.bf file. It'll be rolled into the next build; in the meantime, you can download (Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login) and run it in the current version of HyPhy. File input are normal sequence data, with TipDate style node dates (or embedded branch lengths) - take a look at examples given in the prompts. The analysis outputs everything - rates, MRCA and internal node dates, confidence intervals, and dated clock p-values - to the console and graphically automatically.

Cheers,
Sergei