
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:19 08/01/2011

Input data from file "261-aggressive-phylip/10307_hg18.phylip" (29 taxa, 935 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-2767.2735	5674.4204	0.6878	0.6878
GTRG	64	-2769.6462	5676.8556	0.2036	0.8914
HKYIG	61	-2774.3256	5679.3156	0.0595	0.9508
GTRI	64	-2771.8017	5681.1666	0.0236	0.9744
HKYG	60	-2776.4743	5681.3240	0.0218	0.9962
HKYI	60	-2778.2276	5684.8306	0.0038	1.0000
K2PIG	58	-2815.0884	5753.9895	0.0000	1.0000
K2PG	57	-2817.3083	5756.1560	0.0000	1.0000
K2PI	57	-2818.8495	5759.2384	0.0000	1.0000
SYMIG	62	-2813.5994	5760.1574	0.0000	1.0000
SYMG	61	-2815.8401	5762.3445	0.0000	1.0000
SYMI	61	-2817.4555	5765.5755	0.0000	1.0000
F81IG	60	-2834.2615	5796.8983	0.0000	1.0000
F81G	59	-2836.3752	5798.8418	0.0000	1.0000
F81I	59	-2837.9811	5802.0537	0.0000	1.0000
GTR	63	-2835.5830	5806.4242	0.0000	1.0000
HKY	59	-2840.7180	5807.5274	0.0000	1.0000
JC69IG	57	-2868.1146	5857.7685	0.0000	1.0000
JC69G	56	-2870.3238	5859.9186	0.0000	1.0000
JC69I	56	-2871.9033	5863.0777	0.0000	1.0000
K2P	56	-2881.4156	5882.1023	0.0000	1.0000
SYM	60	-2879.8414	5888.0581	0.0000	1.0000
F81	58	-2899.2957	5922.4042	0.0000	1.0000
JC69	55	-2933.4376	5983.8832	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2767.2735	5664.5470	0.7420	0.7420
GTRG	64	-2769.6462	5667.2924	0.1881	0.9301
HKYIG	61	-2774.3256	5670.6512	0.0351	0.9652
GTRI	64	-2771.8017	5671.6034	0.0218	0.9870
HKYG	60	-2776.4743	5672.9487	0.0111	0.9981
HKYI	60	-2778.2276	5676.4553	0.0019	1.0000
K2PIG	58	-2815.0884	5746.1768	0.0000	1.0000
K2PG	57	-2817.3083	5748.6167	0.0000	1.0000
SYMIG	62	-2813.5994	5751.1987	0.0000	1.0000
K2PI	57	-2818.8495	5751.6991	0.0000	1.0000
SYMG	61	-2815.8401	5753.6801	0.0000	1.0000
SYMI	61	-2817.4555	5756.9111	0.0000	1.0000
F81IG	60	-2834.2615	5788.5230	0.0000	1.0000
F81G	59	-2836.3752	5790.7504	0.0000	1.0000
F81I	59	-2837.9811	5793.9622	0.0000	1.0000
GTR	63	-2835.5830	5797.1659	0.0000	1.0000
HKY	59	-2840.7180	5799.4359	0.0000	1.0000
JC69IG	57	-2868.1146	5850.2291	0.0000	1.0000
JC69G	56	-2870.3238	5852.6475	0.0000	1.0000
JC69I	56	-2871.9033	5855.8067	0.0000	1.0000
K2P	56	-2881.4156	5874.8312	0.0000	1.0000
SYM	60	-2879.8414	5879.6829	0.0000	1.0000
F81	58	-2899.2957	5914.5914	0.0000	1.0000
JC69	55	-2933.4376	5976.8753	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2776.4743	5963.3815	0.6872	0.6872
HKYIG	61	-2774.3256	5965.9245	0.1927	0.8799
HKYI	60	-2778.2276	5966.8881	0.1190	0.9989
GTRG	64	-2769.6462	5977.0874	0.0007	0.9997
GTRIG	65	-2767.2735	5979.1825	0.0003	0.9999
GTRI	64	-2771.8017	5981.3984	0.0001	1.0000
K2PG	57	-2817.3083	6024.5279	0.0000	1.0000
K2PIG	58	-2815.0884	6026.9285	0.0000	1.0000
K2PI	57	-2818.8495	6027.6103	0.0000	1.0000
SYMG	61	-2815.8401	6048.9535	0.0000	1.0000
SYMIG	62	-2813.5994	6051.3126	0.0000	1.0000
SYMI	61	-2817.4555	6052.1844	0.0000	1.0000
F81G	59	-2836.3752	6076.3427	0.0000	1.0000
F81IG	60	-2834.2615	6078.9558	0.0000	1.0000
F81I	59	-2837.9811	6079.5545	0.0000	1.0000
HKY	59	-2840.7180	6085.0282	0.0000	1.0000
GTR	63	-2835.5830	6102.1204	0.0000	1.0000
JC69G	56	-2870.3238	6123.7181	0.0000	1.0000
JC69IG	57	-2868.1146	6126.1403	0.0000	1.0000
JC69I	56	-2871.9033	6126.8773	0.0000	1.0000
K2P	56	-2881.4156	6145.9018	0.0000	1.0000
SYM	60	-2879.8414	6170.1157	0.0000	1.0000
F81	58	-2899.2957	6195.3431	0.0000	1.0000
JC69	55	-2933.4376	6243.1053	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10307_hg18.phylip = 1 - 935;
 Partition Dummy = 1:261-aggressive-phylip/10307_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10307_hg18.phylip = 1 - 935;
 Partition Dummy = 1:261-aggressive-phylip/10307_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10307_hg18.phylip = 1 - 935;
 Partition Dummy = 1:261-aggressive-phylip/10307_hg18.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

