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

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

Minimum AIC  model: GTRG
Minimum AICc model: SYMG
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
SYMG	61	-1223.1572	2592.7933	0.3285	0.3285
GTRG	64	-1219.0545	2593.2984	0.2552	0.5836
K2PG	57	-1229.9903	2595.1053	0.1034	0.6870
SYMIG	62	-1223.0817	2595.5270	0.0837	0.7708
GTRIG	65	-1218.9536	2596.0384	0.0648	0.8356
SYMI	61	-1224.8598	2596.1986	0.0598	0.8954
GTRI	64	-1220.7539	2596.6974	0.0466	0.9421
K2PIG	58	-1229.9267	2597.7893	0.0270	0.9691
K2PI	57	-1231.7424	2598.6094	0.0179	0.9870
HKYG	60	-1228.1851	2599.9831	0.0090	0.9961
HKYIG	61	-1228.1115	2602.7020	0.0023	0.9984
HKYI	60	-1229.8965	2603.4059	0.0016	1.0000
SYM	60	-1239.9602	2623.5333	0.0000	1.0000
K2P	56	-1246.0346	2624.4005	0.0000	1.0000
GTR	63	-1237.7307	2627.7285	0.0000	1.0000
HKY	59	-1244.8256	2630.4165	0.0000	1.0000
JC69G	56	-1252.1246	2636.5804	0.0000	1.0000
JC69IG	57	-1252.0827	2639.2901	0.0000	1.0000
JC69I	56	-1253.8079	2639.9471	0.0000	1.0000
F81G	59	-1251.3594	2643.4841	0.0000	1.0000
F81IG	60	-1251.3161	2646.2450	0.0000	1.0000
F81I	59	-1253.0224	2646.8101	0.0000	1.0000
JC69	55	-1267.6933	2664.9422	0.0000	1.0000
F81	58	-1267.1110	2672.1578	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1219.0545	2566.1089	0.4643	0.4643
GTRIG	65	-1218.9536	2567.9073	0.1889	0.6532
SYMG	61	-1223.1572	2568.3143	0.1541	0.8074
GTRI	64	-1220.7539	2569.5079	0.0849	0.8922
SYMIG	62	-1223.0817	2570.1634	0.0611	0.9534
SYMI	61	-1224.8598	2571.7196	0.0281	0.9815
K2PG	57	-1229.9903	2573.9807	0.0091	0.9905
K2PIG	58	-1229.9267	2575.8534	0.0036	0.9941
HKYG	60	-1228.1851	2576.3702	0.0027	0.9968
K2PI	57	-1231.7424	2577.4848	0.0016	0.9984
HKYIG	61	-1228.1115	2578.2230	0.0011	0.9995
HKYI	60	-1229.8965	2579.7930	0.0005	1.0000
SYM	60	-1239.9602	2599.9204	0.0000	1.0000
GTR	63	-1237.7307	2601.4614	0.0000	1.0000
K2P	56	-1246.0346	2604.0693	0.0000	1.0000
HKY	59	-1244.8256	2607.6512	0.0000	1.0000
JC69G	56	-1252.1246	2616.2492	0.0000	1.0000
JC69IG	57	-1252.0827	2618.1655	0.0000	1.0000
JC69I	56	-1253.8079	2619.6159	0.0000	1.0000
F81G	59	-1251.3594	2620.7188	0.0000	1.0000
F81IG	60	-1251.3161	2622.6321	0.0000	1.0000
F81I	59	-1253.0224	2624.0449	0.0000	1.0000
JC69	55	-1267.6933	2645.3866	0.0000	1.0000
F81	58	-1267.1110	2650.2219	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1229.9903	2797.2042	0.8076	0.8076
K2PI	57	-1231.7424	2800.7083	0.1401	0.9477
K2PIG	58	-1229.9267	2802.9931	0.0447	0.9923
SYMG	61	-1223.1572	2807.2026	0.0054	0.9978
SYMI	61	-1224.8598	2810.6080	0.0010	0.9988
HKYG	60	-1228.1851	2811.3423	0.0007	0.9995
SYMIG	62	-1223.0817	2812.9679	0.0003	0.9998
HKYI	60	-1229.8965	2814.7651	0.0001	0.9999
GTRG	64	-1219.0545	2816.7458	0.0000	0.9999
HKYIG	61	-1228.1115	2817.1113	0.0000	1.0000
GTRI	64	-1220.7539	2820.1448	0.0000	1.0000
GTRIG	65	-1218.9536	2822.4604	0.0000	1.0000
K2P	56	-1246.0346	2823.3766	0.0000	1.0000
SYM	60	-1239.9602	2834.8925	0.0000	1.0000
JC69G	56	-1252.1246	2835.5565	0.0000	1.0000
HKY	59	-1244.8256	2838.7071	0.0000	1.0000
JC69I	56	-1253.8079	2838.9232	0.0000	1.0000
JC69IG	57	-1252.0827	2841.3890	0.0000	1.0000
GTR	63	-1237.7307	2848.1821	0.0000	1.0000
F81G	59	-1251.3594	2851.7747	0.0000	1.0000
F81I	59	-1253.0224	2855.1008	0.0000	1.0000
F81IG	60	-1251.3161	2857.6043	0.0000	1.0000
JC69	55	-1267.6933	2860.7778	0.0000	1.0000
F81	58	-1267.1110	2877.3616	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/chr2_18608.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr2_18608.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 (SYMG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_18608.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr2_18608.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_18608.phylip = 1 - 371;
 Partition Dummy = 1:183-phylip/chr2_18608.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

