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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2861.4053	5762.8106	0.4065	0.4065
GTRI	20	-2861.7078	5763.4156	0.3004	0.7069
GTRIG	21	-2860.7374	5763.4747	0.2917	0.9986
HKYG	16	-2871.8143	5775.6286	0.0007	0.9993
HKYI	16	-2872.1083	5776.2166	0.0005	0.9998
HKYIG	17	-2871.8731	5777.7463	0.0002	1.0000
GTR	19	-2889.4094	5816.8188	0.0000	1.0000
HKY	15	-2899.8825	5829.7650	0.0000	1.0000
SYMIG	18	-2904.1200	5844.2399	0.0000	1.0000
SYMI	17	-2905.5713	5845.1427	0.0000	1.0000
SYMG	17	-2905.7429	5845.4858	0.0000	1.0000
K2PG	13	-2912.8178	5851.6356	0.0000	1.0000
K2PI	13	-2912.9994	5851.9988	0.0000	1.0000
K2PIG	14	-2914.2308	5856.4615	0.0000	1.0000
SYM	16	-2933.9091	5899.8181	0.0000	1.0000
K2P	12	-2941.4546	5906.9093	0.0000	1.0000
F81G	15	-2944.2914	5918.5827	0.0000	1.0000
F81I	15	-2944.4324	5918.8648	0.0000	1.0000
F81IG	16	-2944.3546	5920.7091	0.0000	1.0000
F81	14	-2970.1802	5968.3603	0.0000	1.0000
JC69G	12	-2977.0460	5978.0919	0.0000	1.0000
JC69I	12	-2977.1544	5978.3089	0.0000	1.0000
JC69IG	13	-2977.0460	5980.0920	0.0000	1.0000
JC69	11	-3003.2922	6028.5845	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2871.8143	5852.9571	0.5418	0.5418
HKYI	16	-2872.1083	5853.5451	0.4038	0.9456
GTRG	20	-2861.4053	5859.4712	0.0209	0.9665
HKYIG	17	-2871.8731	5859.9078	0.0168	0.9832
GTRI	20	-2861.7078	5860.0762	0.0154	0.9987
GTRIG	21	-2860.7374	5864.9684	0.0013	1.0000
HKY	15	-2899.8825	5902.2605	0.0000	1.0000
GTR	19	-2889.4094	5908.6464	0.0000	1.0000
K2PG	13	-2912.8178	5914.4651	0.0000	1.0000
K2PI	13	-2912.9994	5914.8282	0.0000	1.0000
K2PIG	14	-2914.2308	5924.1239	0.0000	1.0000
SYMI	17	-2905.5713	5927.3042	0.0000	1.0000
SYMG	17	-2905.7429	5927.6473	0.0000	1.0000
SYMIG	18	-2904.1200	5931.2345	0.0000	1.0000
K2P	12	-2941.4546	5964.9057	0.0000	1.0000
SYM	16	-2933.9091	5977.1466	0.0000	1.0000
F81G	15	-2944.2914	5991.0782	0.0000	1.0000
F81I	15	-2944.4324	5991.3603	0.0000	1.0000
F81IG	16	-2944.3546	5998.0376	0.0000	1.0000
F81	14	-2970.1802	6036.0227	0.0000	1.0000
JC69G	12	-2977.0460	6036.0883	0.0000	1.0000
JC69I	12	-2977.1544	6036.3053	0.0000	1.0000
JC69IG	13	-2977.0460	6042.9214	0.0000	1.0000
JC69	11	-3003.2922	6081.7478	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_10445.phylip = 1 - 928;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10445.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_10445.phylip = 1 - 928;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10445.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_10445.phylip = 1 - 928;
 Partition Dummy = 1:591-tomodel-phylip/chr7_10445.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

