
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:53 04/04/2011

Input data from file "591-tomodel-phylip/chr7_2784.phylip" (7 taxa, 449 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1397.8937	2837.7501	0.5673	0.5673
GTRIG	21	-1397.8874	2839.9387	0.1899	0.7572
GTRI	20	-1399.5342	2841.0309	0.1100	0.8671
HKYG	16	-1404.1433	2841.5459	0.0850	0.9522
HKYIG	17	-1404.1373	2843.6945	0.0290	0.9812
HKYI	16	-1405.6529	2844.5650	0.0188	1.0000
GTR	19	-1410.4264	2860.6243	0.0000	1.0000
HKY	15	-1416.6672	2864.4429	0.0000	1.0000
F81G	15	-1425.8497	2882.8079	0.0000	1.0000
F81IG	16	-1425.8408	2884.9408	0.0000	1.0000
F81I	15	-1427.3476	2885.8037	0.0000	1.0000
K2PG	13	-1432.9200	2892.6769	0.0000	1.0000
SYMG	17	-1429.3072	2894.0343	0.0000	1.0000
K2PIG	14	-1432.9034	2894.7745	0.0000	1.0000
K2PI	13	-1434.2422	2895.3211	0.0000	1.0000
SYMIG	18	-1429.2777	2896.1461	0.0000	1.0000
SYMI	17	-1430.4545	2896.3290	0.0000	1.0000
F81	14	-1438.4240	2905.8156	0.0000	1.0000
K2P	12	-1446.6638	2918.0433	0.0000	1.0000
SYM	16	-1442.4261	2918.1115	0.0000	1.0000
JC69G	12	-1451.1896	2927.0948	0.0000	1.0000
JC69IG	13	-1451.1728	2929.1825	0.0000	1.0000
JC69I	12	-1452.6520	2930.0196	0.0000	1.0000
JC69	11	-1464.6309	2951.8659	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1397.8937	2835.7875	0.5774	0.5774
GTRIG	21	-1397.8874	2837.7748	0.2138	0.7912
GTRI	20	-1399.5342	2839.0683	0.1120	0.9031
HKYG	16	-1404.1433	2840.2867	0.0609	0.9640
HKYIG	17	-1404.1373	2842.2746	0.0225	0.9865
HKYI	16	-1405.6529	2843.3058	0.0135	1.0000
GTR	19	-1410.4264	2858.8528	0.0000	1.0000
HKY	15	-1416.6672	2863.3344	0.0000	1.0000
F81G	15	-1425.8497	2881.6993	0.0000	1.0000
F81IG	16	-1425.8408	2883.6815	0.0000	1.0000
F81I	15	-1427.3476	2884.6952	0.0000	1.0000
K2PG	13	-1432.9200	2891.8401	0.0000	1.0000
SYMG	17	-1429.3072	2892.6143	0.0000	1.0000
K2PIG	14	-1432.9034	2893.8068	0.0000	1.0000
K2PI	13	-1434.2422	2894.4843	0.0000	1.0000
SYMIG	18	-1429.2777	2894.5554	0.0000	1.0000
SYMI	17	-1430.4545	2894.9091	0.0000	1.0000
F81	14	-1438.4240	2904.8479	0.0000	1.0000
SYM	16	-1442.4261	2916.8522	0.0000	1.0000
K2P	12	-1446.6638	2917.3277	0.0000	1.0000
JC69G	12	-1451.1896	2926.3792	0.0000	1.0000
JC69IG	13	-1451.1728	2928.3457	0.0000	1.0000
JC69I	12	-1452.6520	2929.3040	0.0000	1.0000
JC69	11	-1464.6309	2951.2618	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1404.1433	2905.9990	0.7863	0.7863
HKYI	16	-1405.6529	2909.0181	0.1738	0.9601
HKYIG	17	-1404.1373	2912.0940	0.0373	0.9974
GTRG	20	-1397.8937	2917.9280	0.0020	0.9995
GTRI	20	-1399.5342	2921.2088	0.0004	0.9998
GTRIG	21	-1397.8874	2924.0223	0.0001	0.9999
HKY	15	-1416.6672	2924.9397	0.0001	1.0000
GTR	19	-1410.4264	2936.8862	0.0000	1.0000
F81G	15	-1425.8497	2943.3047	0.0000	1.0000
K2PG	13	-1432.9200	2945.2314	0.0000	1.0000
F81I	15	-1427.3476	2946.3005	0.0000	1.0000
K2PI	13	-1434.2422	2947.8756	0.0000	1.0000
F81IG	16	-1425.8408	2949.3939	0.0000	1.0000
K2PIG	14	-1432.9034	2951.3051	0.0000	1.0000
F81	14	-1438.4240	2962.3462	0.0000	1.0000
SYMG	17	-1429.3072	2962.4337	0.0000	1.0000
SYMI	17	-1430.4545	2964.7284	0.0000	1.0000
K2P	12	-1446.6638	2966.6119	0.0000	1.0000
SYMIG	18	-1429.2777	2968.4818	0.0000	1.0000
JC69G	12	-1451.1896	2975.6634	0.0000	1.0000
JC69I	12	-1452.6520	2978.5883	0.0000	1.0000
JC69IG	13	-1451.1728	2981.7370	0.0000	1.0000
SYM	16	-1442.4261	2982.5646	0.0000	1.0000
JC69	11	-1464.6309	2996.4390	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 591-tomodel-phylip/chr7_2784.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2784.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 591-tomodel-phylip/chr7_2784.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2784.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 (HKYG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_2784.phylip = 1 - 449;
 Partition Dummy = 1:591-tomodel-phylip/chr7_2784.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

