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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: K2PI


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
GTRI	20	-2645.7833	5332.8549	0.3645	0.3645
GTRG	20	-2646.2372	5333.7627	0.2315	0.5960
GTRIG	21	-2645.3702	5334.1597	0.1898	0.7859
SYMI	17	-2650.1879	5335.3102	0.1068	0.8927
SYMG	17	-2650.8521	5336.6385	0.0550	0.9477
SYMIG	18	-2649.9278	5336.9015	0.0482	0.9959
K2PI	13	-2658.4719	5343.4961	0.0018	0.9976
K2PG	13	-2659.2048	5344.9620	0.0009	0.9985
K2PIG	14	-2658.2015	5345.0412	0.0008	0.9993
HKYI	16	-2656.9390	5346.7072	0.0004	0.9997
HKYG	16	-2657.7289	5348.2870	0.0002	0.9998
HKYIG	17	-2656.6821	5348.2986	0.0002	1.0000
GTR	19	-2691.4456	5422.0551	0.0000	1.0000
SYM	16	-2697.9780	5428.7852	0.0000	1.0000
K2P	12	-2709.2722	5443.0172	0.0000	1.0000
HKY	15	-2706.5337	5443.7980	0.0000	1.0000
JC69I	12	-2714.9385	5454.3498	0.0000	1.0000
JC69IG	13	-2714.8241	5456.2005	0.0000	1.0000
JC69G	12	-2716.0557	5456.5841	0.0000	1.0000
F81I	15	-2713.2726	5457.2759	0.0000	1.0000
F81IG	16	-2713.1881	5459.2055	0.0000	1.0000
F81G	15	-2714.5185	5459.7675	0.0000	1.0000
JC69	11	-2763.1660	5548.7313	0.0000	1.0000
F81	14	-2760.8652	5550.3686	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2645.7833	5331.5666	0.3720	0.3720
GTRG	20	-2646.2372	5332.4743	0.2363	0.6083
GTRIG	21	-2645.3702	5332.7404	0.2068	0.8151
SYMI	17	-2650.1879	5334.3758	0.0913	0.9064
SYMG	17	-2650.8521	5335.7041	0.0470	0.9534
SYMIG	18	-2649.9278	5335.8556	0.0436	0.9970
K2PI	13	-2658.4719	5342.9437	0.0013	0.9982
K2PIG	14	-2658.2015	5344.4029	0.0006	0.9988
K2PG	13	-2659.2048	5344.4096	0.0006	0.9994
HKYI	16	-2656.9390	5345.8780	0.0003	0.9997
HKYIG	17	-2656.6821	5347.3642	0.0001	0.9999
HKYG	16	-2657.7289	5347.4578	0.0001	1.0000
GTR	19	-2691.4456	5420.8913	0.0000	1.0000
SYM	16	-2697.9780	5427.9559	0.0000	1.0000
K2P	12	-2709.2722	5442.5445	0.0000	1.0000
HKY	15	-2706.5337	5443.0674	0.0000	1.0000
JC69I	12	-2714.9385	5453.8771	0.0000	1.0000
JC69IG	13	-2714.8241	5455.6482	0.0000	1.0000
JC69G	12	-2716.0557	5456.1114	0.0000	1.0000
F81I	15	-2713.2726	5456.5453	0.0000	1.0000
F81IG	16	-2713.1881	5458.3762	0.0000	1.0000
F81G	15	-2714.5185	5459.0369	0.0000	1.0000
JC69	11	-2763.1660	5548.3319	0.0000	1.0000
F81	14	-2760.8652	5549.7303	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-2658.4719	5401.5964	0.6472	0.6472
K2PG	13	-2659.2048	5403.0623	0.3110	0.9582
K2PIG	14	-2658.2015	5407.5673	0.0327	0.9908
SYMI	17	-2650.1879	5411.0755	0.0057	0.9965
SYMG	17	-2650.8521	5412.4038	0.0029	0.9994
SYMIG	18	-2649.9278	5417.0670	0.0003	0.9997
HKYI	16	-2656.9390	5418.0659	0.0002	0.9999
HKYG	16	-2657.7289	5419.6457	0.0001	0.9999
GTRI	20	-2645.7833	5421.8015	0.0000	1.0000
GTRG	20	-2646.2372	5422.7092	0.0000	1.0000
HKYIG	17	-2656.6821	5424.0639	0.0000	1.0000
GTRIG	21	-2645.3702	5427.4870	0.0000	1.0000
K2P	12	-2709.2722	5496.6854	0.0000	1.0000
SYM	16	-2697.9780	5500.1439	0.0000	1.0000
GTR	19	-2691.4456	5506.6144	0.0000	1.0000
JC69I	12	-2714.9385	5508.0180	0.0000	1.0000
JC69G	12	-2716.0557	5510.2523	0.0000	1.0000
HKY	15	-2706.5337	5510.7436	0.0000	1.0000
JC69IG	13	-2714.8241	5514.3009	0.0000	1.0000
F81I	15	-2713.2726	5524.2215	0.0000	1.0000
F81G	15	-2714.5185	5526.7131	0.0000	1.0000
F81IG	16	-2713.1881	5530.5641	0.0000	1.0000
JC69	11	-2763.1660	5597.9611	0.0000	1.0000
F81	14	-2760.8652	5612.8948	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3129.phylip = 1 - 673;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3129.phylip = 1 - 673;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3129.phylip = 1 - 673;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

