IMa Island Model Load Mode difficulty

58 views
Skip to first unread message

Clifford Cunningham

unread,
May 11, 2015, 11:28:06 AM5/11/15
to isolation-wi...@googlegroups.com
Dear Jody,

I run IMa in Island model (-j4) .

In the Load trees mode model testing, it correctly gives T as very large, but for the alternate models, it seems to be estimating t, and the values of M are very strange.

This means I can't compare the log values either, correct?

Is there another instruction I should add to the load trees mode to do the model comparison correctly?  Thank you!

-i /home/abyss/cliff/IMaLo/LeFruBiEx.b1.s.txt -o /home/abyss/cliff/IMaLo/LeFruBiEx.b1.s.95.Eq.PST.out -r LeFruBiEx.b1.s.95.Eq.out -q1 100 -q2 100 -qA 100 -m1 100 -m2 100 -j4 -n 10 -c02 -b 20000 -L 100000 -p1234


Multidimensional Peak Location and Posterior Probability
--------------------------------------------------------
 Set   t       log(P)      q1         q2         m1         m2     
   0  1000000.0000  -1.1661     1.8538     1.1616     2.1456   0.000119
   1  1000000.0000  -1.2908     2.8910     1.2721     2.4471   0.000107
 ALL  1000000.0000  -1.3919     2.5014     1.1509     2.3627   0.000100


Nested Models Peak Location and Posterior Probability
------------------------------------------------------
Model               t          log(P)      q1          q2          m1          m2         df       2LLR
 AB_DD    0.009566    -4.6496     2.5014      1.1509  [   0.0001] [     -nan] [     -nan]    1        6.5153
 AA_DE    0.125329    -2.0768     2.5014  [   2.5014] [   0.0001] [   1.1509] [     -nan]    1        1.3698
 AA_DD    0.004498    -5.4042     2.5014  [   2.5014] [   0.0001] [   1.1509] [   1.1509]    2        8.0246

Reply all
Reply to author
Forward
0 new messages