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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKY


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	-829.6913	1701.7757	0.2287	0.2287
GTRI	20	-829.7750	1701.9431	0.2104	0.4391
HKYG	16	-834.3622	1702.2569	0.1798	0.6189
HKYI	16	-834.7622	1703.0569	0.1205	0.7395
GTRIG	21	-829.6555	1703.9509	0.0771	0.8166
HKY	15	-836.4693	1704.2869	0.0652	0.8817
HKYIG	17	-834.3282	1704.3853	0.0620	0.9438
GTR	19	-832.2112	1704.5816	0.0562	1.0000
F81G	15	-850.5438	1732.4359	0.0000	1.0000
F81I	15	-850.9759	1733.3000	0.0000	1.0000
F81	14	-852.1912	1733.5589	0.0000	1.0000
F81IG	16	-850.5344	1734.6012	0.0000	1.0000
SYMI	17	-849.4457	1734.6201	0.0000	1.0000
SYMG	17	-849.4741	1734.6771	0.0000	1.0000
K2PG	13	-853.9182	1734.8532	0.0000	1.0000
K2PI	13	-854.3078	1735.6324	0.0000	1.0000
SYMIG	18	-849.4871	1736.9119	0.0000	1.0000
K2PIG	14	-853.8732	1736.9228	0.0000	1.0000
SYM	16	-851.7860	1737.1043	0.0000	1.0000
K2P	12	-856.1688	1737.2067	0.0000	1.0000
JC69G	12	-868.0461	1760.9613	0.0000	1.0000
JC69I	12	-868.1939	1761.2570	0.0000	1.0000
JC69	11	-869.5225	1761.7784	0.0000	1.0000
JC69IG	13	-868.0467	1763.1101	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-829.6913	1699.3826	0.2676	0.2676
GTRI	20	-829.7750	1699.5499	0.2461	0.5137
HKYG	16	-834.3622	1700.7245	0.1368	0.6505
GTRIG	21	-829.6555	1701.3109	0.1020	0.7525
HKYI	16	-834.7622	1701.5245	0.0917	0.8442
GTR	19	-832.2112	1702.4225	0.0585	0.9027
HKYIG	17	-834.3282	1702.6565	0.0521	0.9548
HKY	15	-836.4693	1702.9386	0.0452	1.0000
F81G	15	-850.5438	1731.0876	0.0000	1.0000
F81I	15	-850.9759	1731.9517	0.0000	1.0000
F81	14	-852.1912	1732.3824	0.0000	1.0000
SYMI	17	-849.4457	1732.8913	0.0000	1.0000
SYMG	17	-849.4741	1732.9482	0.0000	1.0000
F81IG	16	-850.5344	1733.0688	0.0000	1.0000
K2PG	13	-853.9182	1733.8364	0.0000	1.0000
K2PI	13	-854.3078	1734.6157	0.0000	1.0000
SYMIG	18	-849.4871	1734.9742	0.0000	1.0000
SYM	16	-851.7860	1735.5719	0.0000	1.0000
K2PIG	14	-853.8732	1735.7463	0.0000	1.0000
K2P	12	-856.1688	1736.3376	0.0000	1.0000
JC69G	12	-868.0461	1760.0923	0.0000	1.0000
JC69I	12	-868.1939	1760.3879	0.0000	1.0000
JC69	11	-869.5225	1761.0450	0.0000	1.0000
JC69IG	13	-868.0467	1762.0934	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	15	-836.4693	1761.7220	0.5759	0.5759
HKYG	16	-834.3622	1763.4268	0.2456	0.8215
HKYI	16	-834.7622	1764.2268	0.1646	0.9861
HKYIG	17	-834.3282	1769.2777	0.0132	0.9993
GTR	19	-832.2112	1776.8815	0.0003	0.9996
GTRG	20	-829.6913	1777.7605	0.0002	0.9998
GTRI	20	-829.7750	1777.9278	0.0002	1.0000
K2P	12	-856.1688	1783.3643	0.0000	1.0000
GTRIG	21	-829.6555	1783.6077	0.0000	1.0000
K2PG	13	-853.9182	1784.7820	0.0000	1.0000
K2PI	13	-854.3078	1785.5613	0.0000	1.0000
F81	14	-852.1912	1787.2470	0.0000	1.0000
F81G	15	-850.5438	1789.8710	0.0000	1.0000
K2PIG	14	-853.8732	1790.6108	0.0000	1.0000
F81I	15	-850.9759	1790.7351	0.0000	1.0000
F81IG	16	-850.5344	1795.7711	0.0000	1.0000
SYM	16	-851.7860	1798.2742	0.0000	1.0000
SYMI	17	-849.4457	1799.5125	0.0000	1.0000
SYMG	17	-849.4741	1799.5694	0.0000	1.0000
JC69	11	-869.5225	1804.1529	0.0000	1.0000
SYMIG	18	-849.4871	1805.5143	0.0000	1.0000
JC69G	12	-868.0461	1807.1190	0.0000	1.0000
JC69I	12	-868.1939	1807.4146	0.0000	1.0000
JC69IG	13	-868.0467	1813.0390	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/chr11_3432.phylip = 1 - 372;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3432.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/chr11_3432.phylip = 1 - 372;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3432.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 (HKY)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr11_3432.phylip = 1 - 372;
 Partition Dummy = 1:591-tomodel-phylip/chr11_3432.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

