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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	64	-1157.6414	2460.1933	0.4942	0.4942
GTRIG	65	-1157.4266	2462.3278	0.1700	0.6642
GTRG	64	-1158.7408	2462.3921	0.1646	0.8288
HKYI	60	-1164.3031	2463.3642	0.1012	0.9300
HKYG	60	-1165.2050	2465.1681	0.0411	0.9711
HKYIG	61	-1164.2953	2465.8714	0.0289	1.0000
F81G	59	-1178.0835	2488.4125	0.0000	1.0000
F81I	59	-1178.5139	2489.2734	0.0000	1.0000
GTR	63	-1173.6321	2489.6213	0.0000	1.0000
F81IG	60	-1178.4130	2491.5841	0.0000	1.0000
HKY	59	-1180.8653	2493.9761	0.0000	1.0000
SYMI	61	-1178.7060	2494.6928	0.0000	1.0000
K2PI	57	-1184.2158	2495.6821	0.0000	1.0000
SYMG	61	-1179.5003	2496.2813	0.0000	1.0000
SYMIG	62	-1178.3388	2496.4914	0.0000	1.0000
K2PG	57	-1185.2200	2497.6905	0.0000	1.0000
K2PIG	58	-1184.1540	2498.0510	0.0000	1.0000
F81	58	-1192.7455	2515.2339	0.0000	1.0000
JC69G	56	-1196.9563	2518.6807	0.0000	1.0000
JC69I	56	-1197.2268	2519.2215	0.0000	1.0000
JC69IG	57	-1196.9859	2521.2223	0.0000	1.0000
SYM	60	-1194.1900	2523.1380	0.0000	1.0000
K2P	56	-1200.9107	2526.5895	0.0000	1.0000
JC69	55	-1211.6381	2545.5717	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-1157.6414	2443.2827	0.5227	0.5227
GTRIG	65	-1157.4266	2444.8532	0.2384	0.7611
GTRG	64	-1158.7408	2445.4816	0.1741	0.9352
HKYI	60	-1164.3031	2448.6062	0.0365	0.9717
HKYG	60	-1165.2050	2450.4100	0.0148	0.9865
HKYIG	61	-1164.2953	2450.5906	0.0135	1.0000
GTR	63	-1173.6321	2473.2643	0.0000	1.0000
F81G	59	-1178.0835	2474.1670	0.0000	1.0000
F81I	59	-1178.5139	2475.0279	0.0000	1.0000
F81IG	60	-1178.4130	2476.8261	0.0000	1.0000
SYMI	61	-1178.7060	2479.4120	0.0000	1.0000
HKY	59	-1180.8653	2479.7306	0.0000	1.0000
SYMIG	62	-1178.3388	2480.6777	0.0000	1.0000
SYMG	61	-1179.5003	2481.0005	0.0000	1.0000
K2PI	57	-1184.2158	2482.4316	0.0000	1.0000
K2PIG	58	-1184.1540	2484.3080	0.0000	1.0000
K2PG	57	-1185.2200	2484.4400	0.0000	1.0000
F81	58	-1192.7455	2501.4910	0.0000	1.0000
JC69G	56	-1196.9563	2505.9127	0.0000	1.0000
JC69I	56	-1197.2268	2506.4535	0.0000	1.0000
JC69IG	57	-1196.9859	2507.9718	0.0000	1.0000
SYM	60	-1194.1900	2508.3800	0.0000	1.0000
K2P	56	-1200.9107	2513.8215	0.0000	1.0000
JC69	55	-1211.6381	2533.2763	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1164.3031	2707.9601	0.6887	0.6887
HKYG	60	-1165.2050	2709.7639	0.2795	0.9681
HKYIG	61	-1164.2953	2714.2671	0.0294	0.9975
GTRI	64	-1157.6414	2719.9269	0.0017	0.9993
GTRG	64	-1158.7408	2722.1258	0.0006	0.9999
GTRIG	65	-1157.4266	2725.8200	0.0001	0.9999
K2PI	57	-1184.2158	2728.8178	0.0000	1.0000
F81G	59	-1178.0835	2729.1984	0.0000	1.0000
F81I	59	-1178.5139	2730.0592	0.0000	1.0000
K2PG	57	-1185.2200	2730.8262	0.0000	1.0000
HKY	59	-1180.8653	2734.7620	0.0000	1.0000
K2PIG	58	-1184.1540	2735.0168	0.0000	1.0000
F81IG	60	-1178.4130	2736.1800	0.0000	1.0000
SYMI	61	-1178.7060	2743.0884	0.0000	1.0000
SYMG	61	-1179.5003	2744.6770	0.0000	1.0000
GTR	63	-1173.6321	2745.5859	0.0000	1.0000
JC69G	56	-1196.9563	2747.9763	0.0000	1.0000
JC69I	56	-1197.2268	2748.5172	0.0000	1.0000
SYMIG	62	-1178.3388	2748.6767	0.0000	1.0000
F81	58	-1192.7455	2752.1997	0.0000	1.0000
JC69IG	57	-1196.9859	2754.3581	0.0000	1.0000
K2P	56	-1200.9107	2755.8851	0.0000	1.0000
SYM	60	-1194.1900	2767.7339	0.0000	1.0000
JC69	55	-1211.6381	2771.0174	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_23668.phylip = 1 - 557;
 Partition Dummy = 1:183-phylip/chr2_23668.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_23668.phylip = 1 - 557;
 Partition Dummy = 1:183-phylip/chr2_23668.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

