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

Input data from file "591-tomodel-phylip/chr1_24433.phylip" (7 taxa, 404 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	20	-1306.2237	2654.6405	0.4426	0.4426
GTRIG	21	-1305.7744	2655.9677	0.2280	0.6706
GTRI	20	-1307.0258	2656.2447	0.1985	0.8691
SYMG	17	-1312.1355	2659.8566	0.0326	0.9017
SYMI	17	-1312.6019	2660.7892	0.0205	0.9222
SYMIG	18	-1311.7152	2661.2069	0.0166	0.9388
K2PG	13	-1317.1587	2661.2508	0.0162	0.9550
HKYG	16	-1313.9387	2661.2831	0.0160	0.9710
K2PIG	14	-1316.8190	2662.7177	0.0078	0.9788
HKYIG	17	-1313.5826	2662.7508	0.0077	0.9865
K2PI	13	-1317.9654	2662.8641	0.0072	0.9937
HKYI	16	-1314.8710	2663.1477	0.0063	1.0000
GTR	19	-1342.3143	2724.6077	0.0000	1.0000
SYM	16	-1345.8924	2725.1905	0.0000	1.0000
JC69G	12	-1350.5331	2725.8642	0.0000	1.0000
JC69I	12	-1351.0234	2726.8447	0.0000	1.0000
JC69IG	13	-1350.1781	2727.2895	0.0000	1.0000
F81G	15	-1348.0417	2727.3205	0.0000	1.0000
F81I	15	-1348.5218	2728.2808	0.0000	1.0000
F81IG	16	-1347.6780	2728.7616	0.0000	1.0000
K2P	12	-1352.0701	2728.9382	0.0000	1.0000
HKY	15	-1349.7936	2730.8243	0.0000	1.0000
JC69	11	-1380.9314	2784.5362	0.0000	1.0000
F81	14	-1378.5982	2786.2761	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1306.2237	2652.4473	0.4486	0.4486
GTRIG	21	-1305.7744	2653.5489	0.2586	0.7073
GTRI	20	-1307.0258	2654.0515	0.2012	0.9085
SYMG	17	-1312.1355	2658.2711	0.0244	0.9329
SYMI	17	-1312.6019	2659.2037	0.0153	0.9482
SYMIG	18	-1311.7152	2659.4303	0.0137	0.9618
HKYG	16	-1313.9387	2659.8774	0.0109	0.9727
K2PG	13	-1317.1587	2660.3175	0.0088	0.9815
HKYIG	17	-1313.5826	2661.1653	0.0057	0.9873
K2PIG	14	-1316.8190	2661.6380	0.0045	0.9918
HKYI	16	-1314.8710	2661.7420	0.0043	0.9961
K2PI	13	-1317.9654	2661.9307	0.0039	1.0000
GTR	19	-1342.3143	2722.6285	0.0000	1.0000
SYM	16	-1345.8924	2723.7848	0.0000	1.0000
JC69G	12	-1350.5331	2725.0662	0.0000	1.0000
JC69I	12	-1351.0234	2726.0468	0.0000	1.0000
F81G	15	-1348.0417	2726.0834	0.0000	1.0000
JC69IG	13	-1350.1781	2726.3562	0.0000	1.0000
F81I	15	-1348.5218	2727.0437	0.0000	1.0000
F81IG	16	-1347.6780	2727.3559	0.0000	1.0000
K2P	12	-1352.0701	2728.1402	0.0000	1.0000
HKY	15	-1349.7936	2729.5871	0.0000	1.0000
JC69	11	-1380.9314	2783.8628	0.0000	1.0000
F81	14	-1378.5982	2785.1964	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1317.1587	2712.3359	0.6569	0.6569
K2PI	13	-1317.9654	2713.9491	0.2932	0.9501
K2PIG	14	-1316.8190	2717.6578	0.0459	0.9960
HKYG	16	-1313.9387	2723.9000	0.0020	0.9980
HKYI	16	-1314.8710	2725.7647	0.0008	0.9988
SYMG	17	-1312.1355	2726.2951	0.0006	0.9994
SYMI	17	-1312.6019	2727.2278	0.0004	0.9998
HKYIG	17	-1313.5826	2729.1893	0.0001	0.9999
SYMIG	18	-1311.7152	2731.4558	0.0000	1.0000
GTRG	20	-1306.2237	2732.4756	0.0000	1.0000
GTRI	20	-1307.0258	2734.0798	0.0000	1.0000
GTRIG	21	-1305.7744	2737.5786	0.0000	1.0000
JC69G	12	-1350.5331	2773.0832	0.0000	1.0000
JC69I	12	-1351.0234	2774.0637	0.0000	1.0000
K2P	12	-1352.0701	2776.1572	0.0000	1.0000
JC69IG	13	-1350.1781	2778.3746	0.0000	1.0000
F81G	15	-1348.0417	2786.1046	0.0000	1.0000
F81I	15	-1348.5218	2787.0649	0.0000	1.0000
SYM	16	-1345.8924	2787.8074	0.0000	1.0000
HKY	15	-1349.7936	2789.6084	0.0000	1.0000
F81IG	16	-1347.6780	2791.3785	0.0000	1.0000
GTR	19	-1342.3143	2798.6554	0.0000	1.0000
JC69	11	-1380.9314	2827.8783	0.0000	1.0000
F81	14	-1378.5982	2841.2162	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/chr1_24433.phylip = 1 - 404;
 Partition Dummy = 1:591-tomodel-phylip/chr1_24433.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/chr1_24433.phylip = 1 - 404;
 Partition Dummy = 1:591-tomodel-phylip/chr1_24433.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 591-tomodel-phylip/chr1_24433.phylip = 1 - 404;
 Partition Dummy = 1:591-tomodel-phylip/chr1_24433.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

