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

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

Minimum AIC  model: K2PIG
Minimum AICc model: K2PIG
Minimum BIC  model: K2PI


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
K2PIG	58	-1224.9475	2585.9067	0.5637	0.5637
K2PI	57	-1227.3748	2588.0267	0.1953	0.7590
K2PG	57	-1227.4297	2588.1364	0.1849	0.9438
HKYIG	61	-1224.0460	2592.4046	0.0219	0.9657
SYMIG	62	-1223.3592	2593.8309	0.0107	0.9764
HKYI	60	-1226.6205	2594.7704	0.0067	0.9831
HKYG	60	-1226.7199	2594.9691	0.0061	0.9892
SYMI	61	-1225.6313	2595.5753	0.0045	0.9937
SYMG	61	-1225.6520	2595.6167	0.0044	0.9981
GTRIG	65	-1221.8992	2599.4103	0.0007	0.9987
JC69IG	57	-1233.3820	2600.0410	0.0005	0.9992
GTRG	64	-1224.3320	2601.4259	0.0002	0.9995
GTRI	64	-1224.3439	2601.4497	0.0002	0.9997
JC69G	56	-1236.0075	2602.5731	0.0001	0.9998
JC69I	56	-1236.0179	2602.5939	0.0001	1.0000
F81IG	60	-1232.5786	2606.6867	0.0000	1.0000
F81G	59	-1235.2373	2609.2371	0.0000	1.0000
F81I	59	-1235.2437	2609.2499	0.0000	1.0000
K2P	56	-1268.0695	2666.6971	0.0000	1.0000
HKY	59	-1267.4536	2673.6697	0.0000	1.0000
SYM	60	-1266.2035	2673.9365	0.0000	1.0000
GTR	63	-1265.3371	2680.6030	0.0000	1.0000
JC69	55	-1276.4969	2680.8488	0.0000	1.0000
F81	58	-1275.8123	2687.6362	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	58	-1224.9475	2565.8950	0.5535	0.5535
K2PI	57	-1227.3748	2568.7497	0.1328	0.6863
K2PG	57	-1227.4297	2568.8594	0.1257	0.8120
HKYIG	61	-1224.0460	2570.0919	0.0679	0.8799
SYMIG	62	-1223.3592	2570.7185	0.0496	0.9295
HKYI	60	-1226.6205	2573.2410	0.0141	0.9436
SYMI	61	-1225.6313	2573.2626	0.0139	0.9575
SYMG	61	-1225.6520	2573.3040	0.0136	0.9711
HKYG	60	-1226.7199	2573.4397	0.0127	0.9838
GTRIG	65	-1221.8992	2573.7984	0.0106	0.9944
GTRG	64	-1224.3320	2576.6640	0.0025	0.9970
GTRI	64	-1224.3439	2576.6878	0.0025	0.9995
JC69IG	57	-1233.3820	2580.7641	0.0003	0.9998
JC69G	56	-1236.0075	2584.0150	0.0001	0.9999
JC69I	56	-1236.0179	2584.0358	0.0001	0.9999
F81IG	60	-1232.5786	2585.1572	0.0000	1.0000
F81G	59	-1235.2373	2588.4747	0.0000	1.0000
F81I	59	-1235.2437	2588.4875	0.0000	1.0000
K2P	56	-1268.0695	2648.1390	0.0000	1.0000
SYM	60	-1266.2035	2652.4071	0.0000	1.0000
HKY	59	-1267.4536	2652.9072	0.0000	1.0000
GTR	63	-1265.3371	2656.6742	0.0000	1.0000
JC69	55	-1276.4969	2662.9937	0.0000	1.0000
F81	58	-1275.8123	2667.6245	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-1227.3748	2796.4055	0.3964	0.3964
K2PG	57	-1227.4297	2796.5152	0.3752	0.7716
K2PIG	58	-1224.9475	2797.5448	0.2243	0.9959
JC69G	56	-1236.0075	2807.6768	0.0014	0.9973
JC69I	56	-1236.0179	2807.6976	0.0014	0.9987
JC69IG	57	-1233.3820	2808.4199	0.0010	0.9997
HKYI	60	-1226.6205	2812.8787	0.0001	0.9998
HKYG	60	-1226.7199	2813.0774	0.0001	0.9999
HKYIG	61	-1224.0460	2813.7236	0.0001	1.0000
SYMI	61	-1225.6313	2816.8943	0.0000	1.0000
SYMG	61	-1225.6520	2816.9357	0.0000	1.0000
SYMIG	62	-1223.3592	2818.3441	0.0000	1.0000
F81G	59	-1235.2373	2824.1184	0.0000	1.0000
F81I	59	-1235.2437	2824.1312	0.0000	1.0000
F81IG	60	-1232.5786	2824.7949	0.0000	1.0000
GTRG	64	-1224.3320	2832.2776	0.0000	1.0000
GTRI	64	-1224.3439	2832.3013	0.0000	1.0000
GTRIG	65	-1221.8992	2833.4059	0.0000	1.0000
K2P	56	-1268.0695	2871.8008	0.0000	1.0000
JC69	55	-1276.4969	2882.6616	0.0000	1.0000
HKY	59	-1267.4536	2888.5509	0.0000	1.0000
SYM	60	-1266.2035	2892.0448	0.0000	1.0000
F81	58	-1275.8123	2899.2743	0.0000	1.0000
GTR	63	-1265.3371	2908.2938	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr3_3073.phylip = 1 - 401;
 Partition Dummy = 1:183-phylip/chr3_3073.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) 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 (K2PIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_3073.phylip = 1 - 401;
 Partition Dummy = 1:183-phylip/chr3_3073.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) 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 (K2PI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_3073.phylip = 1 - 401;
 Partition Dummy = 1:183-phylip/chr3_3073.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

