
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:39 02/17/2011

Input data from file "183-phylip/chr2_23124.phylip" (29 taxa, 365 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	64	-1548.7276	3253.1885	0.9516	0.9516
SYMG	61	-1556.6006	3260.1648	0.0291	0.9806
K2PG	57	-1563.7136	3262.9646	0.0072	0.9878
K2PIG	58	-1562.6085	3263.5830	0.0053	0.9931
HKYG	60	-1560.4395	3264.9579	0.0026	0.9957
HKYIG	61	-1559.2642	3265.4921	0.0020	0.9977
GTRIG	65	-1553.4957	3265.6871	0.0018	0.9996
SYMI	61	-1561.5262	3270.0161	0.0002	0.9998
SYMIG	62	-1560.1908	3270.2492	0.0002	1.0000
K2PI	57	-1569.7442	3275.0258	0.0000	1.0000
HKYI	60	-1567.0005	3278.0800	0.0000	1.0000
GTRI	64	-1563.0323	3281.7980	0.0000	1.0000
JC69IG	57	-1602.3764	3340.2902	0.0000	1.0000
JC69G	56	-1604.2123	3341.1518	0.0000	1.0000
F81G	59	-1601.7697	3344.7525	0.0000	1.0000
F81IG	60	-1600.5915	3345.2619	0.0000	1.0000
JC69I	56	-1610.0926	3352.9125	0.0000	1.0000
F81I	59	-1607.9344	3357.0819	0.0000	1.0000
GTR	63	-1603.4529	3359.6965	0.0000	1.0000
K2P	56	-1614.4934	3361.7142	0.0000	1.0000
HKY	59	-1610.3107	3361.8346	0.0000	1.0000
SYM	60	-1610.0558	3364.1905	0.0000	1.0000
JC69	55	-1653.9286	3437.7924	0.0000	1.0000
F81	58	-1651.1473	3440.6606	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1548.7276	3225.4552	0.9875	0.9875
SYMG	61	-1556.6006	3235.2011	0.0076	0.9951
GTRIG	65	-1553.4957	3236.9914	0.0031	0.9982
HKYIG	61	-1559.2642	3240.5284	0.0005	0.9987
HKYG	60	-1560.4395	3240.8790	0.0004	0.9992
K2PIG	58	-1562.6085	3241.2170	0.0004	0.9995
K2PG	57	-1563.7136	3241.4271	0.0003	0.9999
SYMIG	62	-1560.1908	3244.3817	0.0001	0.9999
SYMI	61	-1561.5262	3245.0524	0.0001	1.0000
K2PI	57	-1569.7442	3253.4884	0.0000	1.0000
HKYI	60	-1567.0005	3254.0010	0.0000	1.0000
GTRI	64	-1563.0323	3254.0647	0.0000	1.0000
JC69IG	57	-1602.3764	3318.7527	0.0000	1.0000
JC69G	56	-1604.2123	3320.4246	0.0000	1.0000
F81IG	60	-1600.5915	3321.1829	0.0000	1.0000
F81G	59	-1601.7697	3321.5394	0.0000	1.0000
JC69I	56	-1610.0926	3332.1852	0.0000	1.0000
GTR	63	-1603.4529	3332.9058	0.0000	1.0000
F81I	59	-1607.9344	3333.8687	0.0000	1.0000
HKY	59	-1610.3107	3338.6215	0.0000	1.0000
SYM	60	-1610.0558	3340.1115	0.0000	1.0000
K2P	56	-1614.4934	3340.9869	0.0000	1.0000
JC69	55	-1653.9286	3417.8571	0.0000	1.0000
F81	58	-1651.1473	3418.2946	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1563.7136	3463.7213	0.8491	0.8491
K2PIG	58	-1562.6085	3467.4110	0.1342	0.9834
SYMG	61	-1556.6006	3473.0949	0.0078	0.9912
HKYG	60	-1560.4395	3474.8728	0.0032	0.9944
GTRG	64	-1548.7276	3475.0486	0.0029	0.9973
K2PI	57	-1569.7442	3475.7825	0.0020	0.9994
HKYIG	61	-1559.2642	3478.4221	0.0005	0.9999
SYMI	61	-1561.5262	3482.9461	0.0001	1.0000
SYMIG	62	-1560.1908	3486.1753	0.0000	1.0000
HKYI	60	-1567.0005	3487.9949	0.0000	1.0000
GTRIG	65	-1553.4957	3490.4847	0.0000	1.0000
GTRI	64	-1563.0323	3503.6581	0.0000	1.0000
JC69G	56	-1604.2123	3538.8188	0.0000	1.0000
JC69IG	57	-1602.3764	3541.0469	0.0000	1.0000
JC69I	56	-1610.0926	3550.5794	0.0000	1.0000
F81G	59	-1601.7697	3551.6334	0.0000	1.0000
F81IG	60	-1600.5915	3555.1768	0.0000	1.0000
K2P	56	-1614.4934	3559.3811	0.0000	1.0000
F81I	59	-1607.9344	3563.9627	0.0000	1.0000
HKY	59	-1610.3107	3568.7154	0.0000	1.0000
SYM	60	-1610.0558	3574.1054	0.0000	1.0000
GTR	63	-1603.4529	3578.5994	0.0000	1.0000
JC69	55	-1653.9286	3632.3515	0.0000	1.0000
F81	58	-1651.1473	3644.4887	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 183-phylip/chr2_23124.phylip = 1 - 365;
 Partition Dummy = 1:183-phylip/chr2_23124.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 183-phylip/chr2_23124.phylip = 1 - 365;
 Partition Dummy = 1:183-phylip/chr2_23124.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 183-phylip/chr2_23124.phylip = 1 - 365;
 Partition Dummy = 1:183-phylip/chr2_23124.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

