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

Input data from file "591-tomodel-phylip/chr5_1689.phylip" (7 taxa, 733 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	-1744.6507	3530.4811	0.5087	0.5087
GTRI	20	-1745.2170	3531.6139	0.2887	0.7975
GTRIG	21	-1744.5332	3532.3659	0.1982	0.9957
SYMG	17	-1753.6150	3542.0860	0.0015	0.9972
SYMI	17	-1754.1157	3543.0873	0.0009	0.9982
SYMIG	18	-1753.5312	3544.0204	0.0006	0.9987
HKYG	16	-1755.8185	3544.3967	0.0005	0.9992
HKYI	16	-1756.2736	3545.3070	0.0003	0.9995
GTR	19	-1753.2636	3545.5931	0.0003	0.9998
HKYIG	17	-1755.7180	3546.2920	0.0002	1.0000
K2PG	13	-1763.8554	3554.2171	0.0000	1.0000
K2PI	13	-1764.2708	3555.0479	0.0000	1.0000
SYM	16	-1761.4700	3555.6998	0.0000	1.0000
K2PIG	14	-1763.7782	3556.1414	0.0000	1.0000
HKY	15	-1764.2742	3559.2178	0.0000	1.0000
K2P	12	-1771.7403	3567.9140	0.0000	1.0000
F81G	15	-1769.4834	3569.6362	0.0000	1.0000
F81I	15	-1769.8725	3570.4145	0.0000	1.0000
F81IG	16	-1769.4147	3571.5891	0.0000	1.0000
JC69G	12	-1776.6407	3577.7148	0.0000	1.0000
JC69I	12	-1776.9966	3578.4265	0.0000	1.0000
JC69IG	13	-1776.5845	3579.6753	0.0000	1.0000
F81	14	-1777.1469	3582.8788	0.0000	1.0000
JC69	11	-1783.9982	3590.3626	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1744.6507	3529.3014	0.5029	0.5029
GTRI	20	-1745.2170	3530.4341	0.2854	0.7883
GTRIG	21	-1744.5332	3531.0663	0.2081	0.9964
SYMG	17	-1753.6150	3541.2300	0.0013	0.9977
SYMI	17	-1754.1157	3542.2314	0.0008	0.9984
SYMIG	18	-1753.5312	3543.0624	0.0005	0.9990
HKYG	16	-1755.8185	3543.6370	0.0004	0.9993
GTR	19	-1753.2636	3544.5272	0.0002	0.9996
HKYI	16	-1756.2736	3544.5472	0.0002	0.9998
HKYIG	17	-1755.7180	3545.4361	0.0002	1.0000
K2PG	13	-1763.8554	3553.7108	0.0000	1.0000
K2PI	13	-1764.2708	3554.5417	0.0000	1.0000
SYM	16	-1761.4700	3554.9400	0.0000	1.0000
K2PIG	14	-1763.7782	3555.5564	0.0000	1.0000
HKY	15	-1764.2742	3558.5483	0.0000	1.0000
K2P	12	-1771.7403	3567.4807	0.0000	1.0000
F81G	15	-1769.4834	3568.9667	0.0000	1.0000
F81I	15	-1769.8725	3569.7450	0.0000	1.0000
F81IG	16	-1769.4147	3570.8293	0.0000	1.0000
JC69G	12	-1776.6407	3577.2814	0.0000	1.0000
JC69I	12	-1776.9966	3577.9931	0.0000	1.0000
JC69IG	13	-1776.5845	3579.1690	0.0000	1.0000
F81	14	-1777.1469	3582.2938	0.0000	1.0000
JC69	11	-1783.9982	3589.9964	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1763.8554	3613.4737	0.4781	0.4781
K2PI	13	-1764.2708	3614.3046	0.3155	0.7936
HKYG	16	-1755.8185	3617.1913	0.0745	0.8681
HKYI	16	-1756.2736	3618.1016	0.0473	0.9154
SYMG	17	-1753.6150	3619.3815	0.0249	0.9403
K2PIG	14	-1763.7782	3619.9165	0.0191	0.9594
SYMI	17	-1754.1157	3620.3828	0.0151	0.9745
GTRG	20	-1744.6507	3621.2443	0.0098	0.9843
GTRI	20	-1745.2170	3622.3770	0.0056	0.9899
K2P	12	-1771.7403	3622.6464	0.0049	0.9947
HKYIG	17	-1755.7180	3623.5875	0.0030	0.9978
SYMIG	18	-1753.5312	3625.8110	0.0010	0.9988
HKY	15	-1764.2742	3627.5055	0.0004	0.9992
GTRIG	21	-1744.5332	3627.6064	0.0004	0.9996
SYM	16	-1761.4700	3628.4943	0.0003	0.9999
GTR	19	-1753.2636	3631.8729	0.0000	0.9999
JC69G	12	-1776.6407	3632.4472	0.0000	1.0000
JC69I	12	-1776.9966	3633.1589	0.0000	1.0000
F81G	15	-1769.4834	3637.9239	0.0000	1.0000
F81I	15	-1769.8725	3638.7022	0.0000	1.0000
JC69IG	13	-1776.5845	3638.9319	0.0000	1.0000
JC69	11	-1783.9982	3640.5650	0.0000	1.0000
F81IG	16	-1769.4147	3644.3837	0.0000	1.0000
F81	14	-1777.1469	3646.6539	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/chr5_1689.phylip = 1 - 733;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1689.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/chr5_1689.phylip = 1 - 733;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1689.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/chr5_1689.phylip = 1 - 733;
 Partition Dummy = 1:591-tomodel-phylip/chr5_1689.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

