Hi All,
I have tried using multi-species coalescent model with a goal of reconstructing a gene for a single locus (as a test) from 3 species (two ingroup and one outgroup). I fixed prior for TMRCA between ingroup and outgroup. But I can't get MRCA for ingroup. Why is it so?
I have mentioned the parameter and problem below. Please let me know if anyone have any idea on it.
***********Parameters****************
1. site model - gamma category count 4, shape checked for estimate, substitution model GTR
2. gene clock model - clock rate unchecked for estimate
3. Species clock model - species tree relaxed clock, stdev 0.1 & checked for estimate and clock rate 1.0 & checked for estimate
4. Priors - Yule model
5. Add prior i) MRCA prior ---> Tree.t:Species
a. internal node prior - uniform species (A & B) as monophyletic;
b. root node prior - lognormal MRCA prior for species A, B, C (M = 1, S = 5.5, offset = 42.0, and checked for mean in real space)
ii) MRCA prior ---> Tree.t.filesuffix (for gene tree and MRCA)
a. internal node prior - uniform alleles of a locus from all individuals species A & B as monophyletic
b. root node prior - lognormal MRCA prior for alleles of a locus from all individuals across species A, B, C (M = 1, S = 5.5, offset = 42.0, and checked for mean in real space)
6. MCMC - chain length 10 millions
*********** Warnings and Errors*******************
A. Warning : initial value <= 0, which is not appropriate for a LogTransform
B. Error message: Start likelihood: -Infinity after 10 initialisation attempts
P(posterior) = -Infinity (was -Infinity)
- java.lang.RuntimeException: Could not find a proper state to initialise. Perhaps try another seed.
When I re-ran without step 5 ii) from Add prior parameter; the run initialized and completed with good ESS value for species clade. However, I want MRCA for internal node.
When I ran loganalyzer I could not find MRCA for internal node (i.e. MRCA for alleles of a locus from all individuals of species A and B. Also, ESS for gene tree total height was very low.
Is possible to increase ESS for gene tree total height as well as extract MRCA from gene tree or loganalyzer output? Please let me know.
Regards,
Kazee