
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:59 02/17/2011

Input data from file "183-phylip/chr8_5307.phylip" (29 taxa, 469 characters)

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-1590.8612	3319.6636	0.6415	0.6415
HKYIG	61	-1590.7749	3322.1346	0.1865	0.8280
GTRG	64	-1587.2340	3323.0620	0.1173	0.9453
GTRIG	65	-1587.1925	3325.6753	0.0318	0.9771
HKYI	60	-1594.3853	3326.7118	0.0189	0.9960
GTRI	64	-1590.6144	3329.8228	0.0040	1.0000
HKY	59	-1609.4459	3354.2024	0.0000	1.0000
GTR	63	-1605.1155	3356.1420	0.0000	1.0000
K2PG	57	-1626.8927	3383.8731	0.0000	1.0000
K2PIG	58	-1626.7180	3386.1286	0.0000	1.0000
SYMG	61	-1624.3185	3389.2219	0.0000	1.0000
K2PI	57	-1630.0950	3390.2777	0.0000	1.0000
SYMIG	62	-1624.1920	3391.6254	0.0000	1.0000
SYMI	61	-1627.1599	3394.9047	0.0000	1.0000
K2P	56	-1646.2528	3420.0007	0.0000	1.0000
SYM	60	-1643.0411	3424.0235	0.0000	1.0000
F81G	59	-1650.7045	3436.7195	0.0000	1.0000
F81IG	60	-1650.5908	3439.1227	0.0000	1.0000
F81I	59	-1653.8464	3443.0033	0.0000	1.0000
F81	58	-1668.6577	3470.0082	0.0000	1.0000
JC69G	56	-1678.8082	3485.1116	0.0000	1.0000
JC69IG	57	-1678.6567	3487.4010	0.0000	1.0000
JC69I	56	-1681.9058	3491.3068	0.0000	1.0000
JC69	55	-1697.4240	3519.7633	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1590.8612	3301.7225	0.4155	0.4155
GTRG	64	-1587.2340	3302.4679	0.2862	0.7017
HKYIG	61	-1590.7749	3303.5498	0.1666	0.8683
GTRIG	65	-1587.1925	3304.3850	0.1097	0.9780
HKYI	60	-1594.3853	3308.7706	0.0122	0.9903
GTRI	64	-1590.6144	3309.2287	0.0097	1.0000
GTR	63	-1605.1155	3336.2309	0.0000	1.0000
HKY	59	-1609.4459	3336.8918	0.0000	1.0000
K2PG	57	-1626.8927	3367.7855	0.0000	1.0000
K2PIG	58	-1626.7180	3369.4360	0.0000	1.0000
SYMG	61	-1624.3185	3370.6371	0.0000	1.0000
SYMIG	62	-1624.1920	3372.3840	0.0000	1.0000
K2PI	57	-1630.0950	3374.1901	0.0000	1.0000
SYMI	61	-1627.1599	3376.3199	0.0000	1.0000
K2P	56	-1646.2528	3404.5055	0.0000	1.0000
SYM	60	-1643.0411	3406.0823	0.0000	1.0000
F81G	59	-1650.7045	3419.4090	0.0000	1.0000
F81IG	60	-1650.5908	3421.1815	0.0000	1.0000
F81I	59	-1653.8464	3425.6928	0.0000	1.0000
F81	58	-1668.6577	3453.3155	0.0000	1.0000
JC69G	56	-1678.8082	3469.6165	0.0000	1.0000
JC69IG	57	-1678.6567	3471.3134	0.0000	1.0000
JC69I	56	-1681.9058	3475.8116	0.0000	1.0000
JC69	55	-1697.4240	3504.8480	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1590.8612	3550.7586	0.9259	0.9259
HKYIG	61	-1590.7749	3556.7366	0.0466	0.9725
HKYI	60	-1594.3853	3557.8068	0.0273	0.9998
GTRG	64	-1587.2340	3568.1065	0.0002	1.0000
GTRIG	65	-1587.1925	3574.1742	0.0000	1.0000
GTRI	64	-1590.6144	3574.8673	0.0000	1.0000
HKY	59	-1609.4459	3581.7774	0.0000	1.0000
GTR	63	-1605.1155	3597.7189	0.0000	1.0000
K2PG	57	-1626.8927	3604.3698	0.0000	1.0000
K2PIG	58	-1626.7180	3610.1709	0.0000	1.0000
K2PI	57	-1630.0950	3610.7744	0.0000	1.0000
SYMG	61	-1624.3185	3623.8239	0.0000	1.0000
SYMI	61	-1627.1599	3629.5067	0.0000	1.0000
SYMIG	62	-1624.1920	3629.7214	0.0000	1.0000
K2P	56	-1646.2528	3636.9393	0.0000	1.0000
SYM	60	-1643.0411	3655.1184	0.0000	1.0000
F81G	59	-1650.7045	3664.2945	0.0000	1.0000
F81IG	60	-1650.5908	3670.2177	0.0000	1.0000
F81I	59	-1653.8464	3670.5784	0.0000	1.0000
F81	58	-1668.6577	3694.0504	0.0000	1.0000
JC69G	56	-1678.8082	3702.0502	0.0000	1.0000
JC69IG	57	-1678.6567	3707.8978	0.0000	1.0000
JC69I	56	-1681.9058	3708.2454	0.0000	1.0000
JC69	55	-1697.4240	3733.1312	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_5307.phylip = 1 - 469;
 Partition Dummy = 1:183-phylip/chr8_5307.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_5307.phylip = 1 - 469;
 Partition Dummy = 1:183-phylip/chr8_5307.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr8_5307.phylip = 1 - 469;
 Partition Dummy = 1:183-phylip/chr8_5307.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

