Migrate Analysis Diagnosis

127 views
Skip to first unread message

Dandára Yasmim Bonfim De Oliveira Silva

unread,
Jan 18, 2024, 8:28:46 PM1/18/24
to migrate...@googlegroups.com
Dear Professor Peter and all MIGRATE users,

I started my interaction with MIGRATE about 4 months ago, using version 3.7.2. My goal is to estimate the values of Theta and gene flow among six populations.

Data summary: I have six populations, 33 SNPs, and a variable number of individuals per population. However, my main concern is determining which parameter in the output file I should rely on to identify the right moment to stop – that is, the number of iterations necessary to obtain reliable data. Is there a specific "Stop Diagnostic" value?
image.png

For instance, we might reach a point where further iterations don't significantly improve the parameters. Identifying this "stopping point" as the "ideal number of iterations" is crucial (see the graph below).

Attached is the result of a recent analysis I did.
Thank you
Best regards,
Dandara Yasmim B. O. Silva
Doutoranda em Ciência Florestal (Unesp/FCA)
Mestre em Ciência Florestal (Unesp/FCA)
Engenharia Florestal  (UFPI/CPCE)
Contato: (89) 999156915
   
Resultados_002.pdf

Peter Beerli

unread,
Jan 18, 2024, 8:41:33 PM1/18/24
to migrate...@googlegroups.com
Dandara,

Your graph suggests that you want to run more than 200,000 iterations without heating,
but I suggest that you repeat your experiment with heating
You also want to consider whether the n-island model is really the model that you want to run, also may want to consider to use 
use the 5.0 version of migrate, which is a Bayesian inference instead of maximum likelihood (using 3.7 is fine, but lack capability such as model selection among different evolutionary models).

Peter


On Jan 18, 2024, at 8:26 PM, Dandára Yasmim Bonfim De Oliveira Silva <dandara...@unesp.br> wrote:

Dear Professor Peter and all MIGRATE users,

I started my interaction with MIGRATE about 4 months ago, using version 3.7.2. My goal is to estimate the values of Theta and gene flow among six populations.

Data summary: I have six populations, 33 SNPs, and a variable number of individuals per population. However, my main concern is determining which parameter in the output file I should rely on to identify the right moment to stop – that is, the number of iterations necessary to obtain reliable data. Is there a specific "Stop Diagnostic" value?
<image.png>

For instance, we might reach a point where further iterations don't significantly improve the parameters. Identifying this "stopping point" as the "ideal number of iterations" is crucial (see the graph below).

Attached is the result of a recent analysis I did.
Thank you
Best regards,
Dandara Yasmim B. O. Silva
Doutoranda em Ciência Florestal (Unesp/FCA)
Mestre em Ciência Florestal (Unesp/FCA)
Engenharia Florestal  (UFPI/CPCE)
Contato: (89) 999156915
   

--
You received this message because you are subscribed to the Google Groups "migrate-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email to migrate-suppo...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/migrate-support/CAG7Z15rjc%3DAktZfgT88NK4g1%3DLrrDUessU3BULK4sUZhaXj_9w%40mail.gmail.com.
<Resultados_002.pdf>

Reply all
Reply to author
Forward
0 new messages