Hi everyone,
I'm running different models in DIYABC using both microsats and mtDNA. The models run fine as well as the analysis, however I get some unexpected results when I estimate posterior distributions of parameters.
This is my scenario
N1 N4
0 sample 1
0 sample 2
ti varNe 2 Ni
t1 merge 1 2
t1 varNe 1 N3
And one of the conditions set is ti < t1. I point out that for both time parameters I'm using very relaxed and broad prior values with considerable overlap.
When estimating posterior distributions of parameters I would expect time parameters to be constrained by the conditions so that ti < t1. However ti is consistently higher than t1 (see picture below).
Am I missing something here?
Hi Francisco, have you tired changing the default settings eg mu rate of SSRs, mean mu value etc...? Do you know what the SSR mutation rate of your organism is?
Hi Samuel A. Logan,
I have some references for the SSR and I did not use the default (only for microsats). But I'm curious, how the settings might contribute to the violation of one of my conditions? It seems odd that I set ti < t1, and the model returns ti > t1.
Not completely sure but If your markers don't fit the model based on your priors, you simply can't make them fit. Have you ticked the box at bottom of the parameter set page to run until all parameters are met? There's two options at the end of that page..Also, how many times do you run it?
Perhaps email the authors and ask them. People are usually happy to reply to queries about their programs.
It might be something very simple.
- Badges
- Science topic
- Similar topics
- Computer Science and Engineering
- Bioinformatics
- Bioinformatic Tools