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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1741.2799	3523.7844	0.3980	0.3980
GTRI	20	-1741.8169	3524.8583	0.2326	0.6307
GTRIG	21	-1741.2657	3525.8802	0.1396	0.7702
HKYG	16	-1747.0950	3526.9785	0.0806	0.8508
HKYI	16	-1747.6320	3528.0525	0.0471	0.8979
SYMG	17	-1747.0789	3529.0461	0.0287	0.9266
HKYIG	17	-1747.0845	3529.0572	0.0285	0.9551
SYMI	17	-1747.6083	3530.1048	0.0169	0.9720
SYMIG	18	-1747.0631	3531.1203	0.0102	0.9821
K2PG	13	-1752.6735	3531.8723	0.0070	0.9891
K2PI	13	-1753.2078	3532.9409	0.0041	0.9932
GTR	19	-1747.0756	3533.2575	0.0035	0.9967
K2PIG	14	-1752.6557	3533.9183	0.0025	0.9992
HKY	15	-1753.0655	3536.8256	0.0006	0.9998
SYM	16	-1753.2522	3539.2928	0.0002	1.0000
K2P	12	-1758.9568	3542.3632	0.0000	1.0000
F81G	15	-1780.7118	3592.1182	0.0000	1.0000
F81I	15	-1781.2489	3593.1924	0.0000	1.0000
F81IG	16	-1780.7086	3594.2056	0.0000	1.0000
JC69G	12	-1785.6769	3595.8033	0.0000	1.0000
JC69I	12	-1786.2139	3596.8773	0.0000	1.0000
JC69IG	13	-1785.6725	3597.8703	0.0000	1.0000
F81	14	-1786.4776	3601.5622	0.0000	1.0000
JC69	11	-1791.5690	3605.5178	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1741.2799	3522.5599	0.4116	0.4116
GTRI	20	-1741.8169	3523.6338	0.2406	0.6521
GTRIG	21	-1741.2657	3524.5313	0.1536	0.8057
HKYG	16	-1747.0950	3526.1901	0.0670	0.8727
HKYI	16	-1747.6320	3527.2641	0.0392	0.9119
SYMG	17	-1747.0789	3528.1578	0.0251	0.9370
HKYIG	17	-1747.0845	3528.1689	0.0249	0.9619
SYMI	17	-1747.6083	3529.2165	0.0148	0.9766
SYMIG	18	-1747.0631	3530.1261	0.0094	0.9860
K2PG	13	-1752.6735	3531.3470	0.0051	0.9911
GTR	19	-1747.0756	3532.1512	0.0034	0.9945
K2PI	13	-1753.2078	3532.4157	0.0030	0.9975
K2PIG	14	-1752.6557	3533.3114	0.0019	0.9994
HKY	15	-1753.0655	3536.1310	0.0005	0.9998
SYM	16	-1753.2522	3538.5044	0.0001	1.0000
K2P	12	-1758.9568	3541.9136	0.0000	1.0000
F81G	15	-1780.7118	3591.4235	0.0000	1.0000
F81I	15	-1781.2489	3592.4977	0.0000	1.0000
F81IG	16	-1780.7086	3593.4172	0.0000	1.0000
JC69G	12	-1785.6769	3595.3538	0.0000	1.0000
JC69I	12	-1786.2139	3596.4277	0.0000	1.0000
JC69IG	13	-1785.6725	3597.3450	0.0000	1.0000
F81	14	-1786.4776	3600.9553	0.0000	1.0000
JC69	11	-1791.5690	3605.1379	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1752.6735	3590.6404	0.5887	0.5887
K2PI	13	-1753.2078	3591.7091	0.3450	0.9337
K2P	12	-1758.9568	3596.6460	0.0292	0.9629
K2PIG	14	-1752.6557	3597.1658	0.0225	0.9854
HKYG	16	-1747.0950	3599.1666	0.0083	0.9937
HKYI	16	-1747.6320	3600.2406	0.0048	0.9986
HKY	15	-1753.0655	3604.5465	0.0006	0.9991
SYMG	17	-1747.0789	3605.6954	0.0003	0.9995
HKYIG	17	-1747.0845	3605.7064	0.0003	0.9998
SYMI	17	-1747.6083	3606.7540	0.0002	1.0000
SYM	16	-1753.2522	3611.4809	0.0000	1.0000
SYMIG	18	-1747.0631	3612.2247	0.0000	1.0000
GTRG	20	-1741.2799	3613.7805	0.0000	1.0000
GTRI	20	-1741.8169	3614.8544	0.0000	1.0000
GTR	19	-1747.0756	3618.8108	0.0000	1.0000
GTRIG	21	-1741.2657	3620.3130	0.0000	1.0000
JC69G	12	-1785.6769	3650.0861	0.0000	1.0000
JC69I	12	-1786.2139	3651.1601	0.0000	1.0000
JC69	11	-1791.5690	3655.3093	0.0000	1.0000
JC69IG	13	-1785.6725	3656.6384	0.0000	1.0000
F81G	15	-1780.7118	3659.8390	0.0000	1.0000
F81I	15	-1781.2489	3660.9132	0.0000	1.0000
F81	14	-1786.4776	3664.8097	0.0000	1.0000
F81IG	16	-1780.7086	3666.3937	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_4670.phylip = 1 - 707;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4670.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_4670.phylip = 1 - 707;
 Partition Dummy = 1:591-tomodel-phylip/chr7_4670.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 (K2PG)]

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

