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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2929.8861	5891.7721	0.4077	0.4077
K2PG	13	-2933.8579	5893.7159	0.1543	0.5620
HKYIG	17	-2929.8853	5893.7705	0.1501	0.7121
GTRG	20	-2927.2904	5894.5808	0.1001	0.8122
K2PIG	14	-2933.8542	5895.7084	0.0570	0.8691
HKYI	16	-2932.0460	5896.0921	0.0470	0.9161
GTRIG	21	-2927.2971	5896.5943	0.0366	0.9527
SYMG	17	-2932.0898	5898.1796	0.0166	0.9693
K2PI	13	-2936.2540	5898.5080	0.0140	0.9833
GTRI	20	-2929.6567	5899.3134	0.0094	0.9927
SYMIG	18	-2932.0844	5900.1689	0.0061	0.9988
SYMI	17	-2934.7492	5903.4983	0.0012	1.0000
HKY	15	-2957.4900	5944.9801	0.0000	1.0000
GTR	19	-2954.3474	5946.6947	0.0000	1.0000
K2P	12	-2961.9418	5947.8835	0.0000	1.0000
SYM	16	-2959.6707	5951.3413	0.0000	1.0000
F81G	15	-2992.1004	6014.2009	0.0000	1.0000
JC69G	12	-2995.7775	6015.5550	0.0000	1.0000
F81IG	16	-2992.1052	6016.2104	0.0000	1.0000
F81I	15	-2993.6551	6017.3102	0.0000	1.0000
JC69IG	13	-2995.7815	6017.5631	0.0000	1.0000
JC69I	12	-2997.4426	6018.8852	0.0000	1.0000
F81	14	-3016.8874	6061.7747	0.0000	1.0000
JC69	11	-3020.7837	6063.5673	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2933.8579	5956.2047	0.8875	0.8875
K2PI	13	-2936.2540	5960.9968	0.0808	0.9684
K2PIG	14	-2933.8542	5963.0041	0.0296	0.9980
HKYG	16	-2929.8861	5968.6814	0.0017	0.9997
HKYI	16	-2932.0460	5973.0013	0.0002	0.9999
HKYIG	17	-2929.8853	5975.4866	0.0001	1.0000
SYMG	17	-2932.0898	5979.8957	0.0000	1.0000
SYMI	17	-2934.7492	5985.2144	0.0000	1.0000
SYMIG	18	-2932.0844	5986.6918	0.0000	1.0000
GTRG	20	-2927.2904	5990.7174	0.0000	1.0000
GTRI	20	-2929.6567	5995.4500	0.0000	1.0000
GTRIG	21	-2927.2971	5997.5377	0.0000	1.0000
K2P	12	-2961.9418	6005.5655	0.0000	1.0000
HKY	15	-2957.4900	6017.0825	0.0000	1.0000
SYM	16	-2959.6707	6028.2506	0.0000	1.0000
GTR	19	-2954.3474	6038.0245	0.0000	1.0000
JC69G	12	-2995.7775	6073.2369	0.0000	1.0000
JC69I	12	-2997.4426	6076.5671	0.0000	1.0000
JC69IG	13	-2995.7815	6080.0519	0.0000	1.0000
F81G	15	-2992.1004	6086.3033	0.0000	1.0000
F81I	15	-2993.6551	6089.4126	0.0000	1.0000
F81IG	16	-2992.1052	6093.1197	0.0000	1.0000
JC69	11	-3020.7837	6116.4425	0.0000	1.0000
F81	14	-3016.8874	6129.0704	0.0000	1.0000



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

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

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


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

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


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

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

