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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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
GTRG	64	-1432.9043	3019.2520	0.4424	0.4424
SYMG	61	-1438.1650	3021.2513	0.1628	0.6053
GTRI	64	-1434.1846	3021.8126	0.1230	0.7282
GTRIG	65	-1432.8853	3022.0896	0.1071	0.8353
SYMI	61	-1439.2872	3023.4957	0.0530	0.8883
SYMIG	62	-1438.1557	3024.0561	0.0401	0.9284
K2PG	57	-1445.1906	3024.1777	0.0377	0.9661
K2PI	57	-1446.1172	3026.0308	0.0149	0.9810
K2PIG	58	-1445.1824	3026.9173	0.0096	0.9906
HKYG	60	-1442.9030	3027.9209	0.0058	0.9964
HKYI	60	-1443.8858	3029.8865	0.0022	0.9986
HKYIG	61	-1442.8927	3030.7067	0.0014	1.0000
SYM	60	-1450.1080	3042.3307	0.0000	1.0000
GTR	63	-1446.3567	3043.2987	0.0000	1.0000
K2P	56	-1456.5050	3044.0666	0.0000	1.0000
HKY	59	-1454.5240	3048.3733	0.0000	1.0000
JC69G	56	-1475.6613	3082.3793	0.0000	1.0000
JC69I	56	-1476.6613	3084.3793	0.0000	1.0000
JC69IG	57	-1475.6471	3085.0905	0.0000	1.0000
F81G	59	-1474.4096	3088.1444	0.0000	1.0000
F81I	59	-1475.4459	3090.2170	0.0000	1.0000
F81IG	60	-1474.3959	3090.9067	0.0000	1.0000
JC69	55	-1486.8730	3102.0793	0.0000	1.0000
F81	58	-1485.6942	3107.9410	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1432.9043	2993.8086	0.5426	0.5426
GTRIG	65	-1432.8853	2995.7706	0.2034	0.7460
GTRI	64	-1434.1846	2996.3691	0.1508	0.8968
SYMG	61	-1438.1650	2998.3300	0.0566	0.9534
SYMIG	62	-1438.1557	3000.3114	0.0210	0.9744
SYMI	61	-1439.2872	3000.5745	0.0184	0.9928
K2PG	57	-1445.1906	3004.3813	0.0027	0.9955
HKYG	60	-1442.9030	3005.8061	0.0013	0.9969
K2PI	57	-1446.1172	3006.2344	0.0011	0.9980
K2PIG	58	-1445.1824	3006.3647	0.0010	0.9990
HKYI	60	-1443.8858	3007.7717	0.0005	0.9995
HKYIG	61	-1442.8927	3007.7855	0.0005	1.0000
GTR	63	-1446.3567	3018.7133	0.0000	1.0000
SYM	60	-1450.1080	3020.2159	0.0000	1.0000
K2P	56	-1456.5050	3025.0099	0.0000	1.0000
HKY	59	-1454.5240	3027.0480	0.0000	1.0000
JC69G	56	-1475.6613	3063.3226	0.0000	1.0000
JC69IG	57	-1475.6471	3065.2941	0.0000	1.0000
JC69I	56	-1476.6613	3065.3225	0.0000	1.0000
F81G	59	-1474.4096	3066.8191	0.0000	1.0000
F81IG	60	-1474.3959	3068.7919	0.0000	1.0000
F81I	59	-1475.4459	3068.8917	0.0000	1.0000
JC69	55	-1486.8730	3083.7460	0.0000	1.0000
F81	58	-1485.6942	3087.3884	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1445.1906	3230.7432	0.6853	0.6853
K2PI	57	-1446.1172	3232.5963	0.2713	0.9566
K2PIG	58	-1445.1824	3236.6979	0.0349	0.9915
SYMG	61	-1438.1650	3240.5770	0.0050	0.9965
SYMI	61	-1439.2872	3242.8215	0.0016	0.9982
HKYG	60	-1442.9030	3244.0818	0.0009	0.9990
HKYI	60	-1443.8858	3246.0474	0.0003	0.9994
SYMIG	62	-1438.1557	3246.5296	0.0003	0.9996
K2P	56	-1456.5050	3247.4006	0.0002	0.9998
GTRG	64	-1432.9043	3247.9694	0.0001	0.9999
HKYIG	61	-1442.8927	3250.0325	0.0000	1.0000
GTRI	64	-1434.1846	3250.5299	0.0000	1.0000
GTRIG	65	-1432.8853	3253.9026	0.0000	1.0000
SYM	60	-1450.1080	3258.4916	0.0000	1.0000
HKY	59	-1454.5240	3261.3524	0.0000	1.0000
GTR	63	-1446.3567	3268.9028	0.0000	1.0000
JC69G	56	-1475.6613	3285.7133	0.0000	1.0000
JC69I	56	-1476.6613	3287.7132	0.0000	1.0000
JC69IG	57	-1475.6471	3291.6561	0.0000	1.0000
F81G	59	-1474.4096	3301.1236	0.0000	1.0000
JC69	55	-1486.8730	3302.1654	0.0000	1.0000
F81I	59	-1475.4459	3303.1962	0.0000	1.0000
F81IG	60	-1474.3959	3307.0676	0.0000	1.0000
F81	58	-1485.6942	3317.7216	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_9749.phylip = 1 - 392;
 Partition Dummy = 1:183-phylip/chr3_9749.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_9749.phylip = 1 - 392;
 Partition Dummy = 1:183-phylip/chr3_9749.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_9749.phylip = 1 - 392;
 Partition Dummy = 1:183-phylip/chr3_9749.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

