
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:41 04/04/2011

Input data from file "591-tomodel-phylip/chr6_4126.phylip" (7 taxa, 725 characters)

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYI


Note: number of parameters compared to sample size
is low (Nchar/Nparams < 40). Consider the use of AICc.


(Output sorted by AICc. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AICc		wAICc	cwAICc
HKYIG	17	-1594.3327	3223.5310	0.5153	0.5153
GTRIG	21	-1590.9154	3225.1452	0.2299	0.7452
HKYI	16	-1597.1023	3226.9731	0.0922	0.8374
HKYG	16	-1597.2734	3227.3152	0.0777	0.9151
GTRI	20	-1593.2131	3227.6194	0.0667	0.9818
GTRG	20	-1594.5472	3230.2875	0.0176	0.9994
F81IG	16	-1602.2908	3237.3499	0.0005	0.9999
F81I	15	-1605.0586	3240.7943	0.0001	1.0000
F81G	15	-1607.7055	3246.0880	0.0000	1.0000
K2PI	13	-1623.1699	3272.8517	0.0000	1.0000
K2PIG	14	-1622.7120	3274.0156	0.0000	1.0000
SYMIG	18	-1619.1124	3275.1936	0.0000	1.0000
SYMI	17	-1621.3547	3277.5751	0.0000	1.0000
K2PG	13	-1625.7245	3277.9609	0.0000	1.0000
SYMG	17	-1622.1979	3279.2614	0.0000	1.0000
JC69IG	13	-1632.9520	3292.4160	0.0000	1.0000
JC69I	12	-1635.9065	3296.2511	0.0000	1.0000
GTR	19	-1629.4005	3297.8790	0.0000	1.0000
JC69G	12	-1637.9831	3300.4044	0.0000	1.0000
HKY	15	-1636.3197	3303.3164	0.0000	1.0000
F81	14	-1645.6527	3319.8970	0.0000	1.0000
SYM	16	-1656.1561	3345.0806	0.0000	1.0000
K2P	12	-1663.6254	3351.6889	0.0000	1.0000
JC69	11	-1674.5312	3371.4326	0.0000	1.0000



(Output sorted by AIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	17	-1594.3327	3222.6654	0.4840	0.4840
GTRIG	21	-1590.9154	3223.8308	0.2702	0.7542
HKYI	16	-1597.1023	3226.2047	0.0825	0.8367
GTRI	20	-1593.2131	3226.4262	0.0738	0.9105
HKYG	16	-1597.2734	3226.5468	0.0695	0.9800
GTRG	20	-1594.5472	3229.0943	0.0194	0.9995
F81IG	16	-1602.2908	3236.5816	0.0005	0.9999
F81I	15	-1605.0586	3240.1173	0.0001	1.0000
F81G	15	-1607.7055	3245.4110	0.0000	1.0000
K2PI	13	-1623.1699	3272.3398	0.0000	1.0000
K2PIG	14	-1622.7120	3273.4240	0.0000	1.0000
SYMIG	18	-1619.1124	3274.2247	0.0000	1.0000
SYMI	17	-1621.3547	3276.7094	0.0000	1.0000
K2PG	13	-1625.7245	3277.4489	0.0000	1.0000
SYMG	17	-1622.1979	3278.3958	0.0000	1.0000
JC69IG	13	-1632.9520	3291.9040	0.0000	1.0000
JC69I	12	-1635.9065	3295.8129	0.0000	1.0000
GTR	19	-1629.4005	3296.8010	0.0000	1.0000
JC69G	12	-1637.9831	3299.9662	0.0000	1.0000
HKY	15	-1636.3197	3302.6394	0.0000	1.0000
F81	14	-1645.6527	3319.3055	0.0000	1.0000
SYM	16	-1656.1561	3344.3123	0.0000	1.0000
K2P	12	-1663.6254	3351.2507	0.0000	1.0000
JC69	11	-1674.5312	3371.0623	0.0000	1.0000



(Output sorted by BIC. w: Akaike weight, cw: cumulative w.)

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1597.1023	3299.5834	0.4080	0.4080
HKYG	16	-1597.2734	3299.9256	0.3438	0.7518
HKYIG	17	-1594.3327	3300.6303	0.2417	0.9935
F81I	15	-1605.0586	3308.9099	0.0038	0.9974
F81IG	16	-1602.2908	3309.9603	0.0023	0.9997
F81G	15	-1607.7055	3314.2036	0.0003	0.9999
GTRI	20	-1593.2131	3318.1496	0.0000	1.0000
GTRIG	21	-1590.9154	3320.1404	0.0000	1.0000
GTRG	20	-1594.5472	3320.8177	0.0000	1.0000
K2PI	13	-1623.1699	3331.9600	0.0000	1.0000
K2PG	13	-1625.7245	3337.0692	0.0000	1.0000
K2PIG	14	-1622.7120	3337.6304	0.0000	1.0000
JC69I	12	-1635.9065	3350.8470	0.0000	1.0000
JC69IG	13	-1632.9520	3351.5242	0.0000	1.0000
SYMI	17	-1621.3547	3354.6744	0.0000	1.0000
JC69G	12	-1637.9831	3355.0003	0.0000	1.0000
SYMG	17	-1622.1979	3356.3607	0.0000	1.0000
SYMIG	18	-1619.1124	3356.7758	0.0000	1.0000
HKY	15	-1636.3197	3371.4319	0.0000	1.0000
F81	14	-1645.6527	3383.5119	0.0000	1.0000
GTR	19	-1629.4005	3383.9383	0.0000	1.0000
K2P	12	-1663.6254	3406.2848	0.0000	1.0000
SYM	16	-1656.1561	3417.6910	0.0000	1.0000
JC69	11	-1674.5312	3421.5102	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (HKYIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_4126.phylip = 1 - 725;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4126.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (HKYIG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_4126.phylip = 1 - 725;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4126.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr6_4126.phylip = 1 - 725;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4126.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

