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

Input data from file "183-phylip/chr1_32247.phylip" (29 taxa, 119 characters)

Minimum AIC  model: GTR
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	56	-307.5573	830.0824	0.8188	0.8188
K2PI	57	-306.2185	834.8305	0.0762	0.8950
K2PG	57	-306.2744	834.9422	0.0721	0.9671
JC69	55	-314.6792	837.1361	0.0241	0.9912
HKY	59	-301.7874	841.5747	0.0026	0.9938
JC69I	56	-313.4795	841.9267	0.0022	0.9960
JC69G	56	-313.5784	842.1245	0.0020	0.9980
K2PIG	58	-306.2838	842.6344	0.0015	0.9995
HKYI	60	-300.4975	847.2020	0.0002	0.9997
HKYG	60	-300.5863	847.3794	0.0001	0.9998
F81	58	-308.9140	847.8948	0.0001	0.9999
JC69IG	57	-313.5876	849.5687	0.0000	1.0000
SYM	60	-302.2871	850.7810	0.0000	1.0000
F81I	59	-307.8389	853.6777	0.0000	1.0000
F81G	59	-307.9740	853.9481	0.0000	1.0000
HKYIG	61	-300.5975	855.8967	0.0000	1.0000
SYMI	61	-301.1943	857.0903	0.0000	1.0000
SYMG	61	-301.3065	857.3147	0.0000	1.0000
F81IG	60	-307.8433	861.8935	0.0000	1.0000
GTR	63	-295.4185	863.4553	0.0000	1.0000
SYMIG	62	-301.3119	866.1237	0.0000	1.0000
GTRI	64	-294.5955	871.2652	0.0000	1.0000
GTRG	64	-294.7740	871.6222	0.0000	1.0000
GTRIG	65	-294.6007	881.0882	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-295.4185	716.8371	0.3015	0.3015
GTRI	64	-294.5955	717.1911	0.2526	0.5540
GTRG	64	-294.7740	717.5481	0.2113	0.7653
GTRIG	65	-294.6007	719.2014	0.0924	0.8577
HKYI	60	-300.4975	720.9951	0.0377	0.8954
HKYG	60	-300.5863	721.1725	0.0345	0.9299
HKY	59	-301.7874	721.5747	0.0282	0.9581
HKYIG	61	-300.5975	723.1949	0.0125	0.9707
SYMI	61	-301.1943	724.3886	0.0069	0.9776
SYM	60	-302.2871	724.5741	0.0063	0.9839
SYMG	61	-301.3065	724.6130	0.0062	0.9901
K2PI	57	-306.2185	726.4370	0.0025	0.9925
K2PG	57	-306.2744	726.5488	0.0023	0.9949
SYMIG	62	-301.3119	726.6237	0.0023	0.9972
K2P	56	-307.5573	727.1146	0.0018	0.9989
K2PIG	58	-306.2838	728.5677	0.0009	0.9998
F81I	59	-307.8389	733.6777	0.0001	0.9998
F81	58	-308.9140	733.8281	0.0001	0.9999
F81G	59	-307.9740	733.9481	0.0001	1.0000
F81IG	60	-307.8433	735.6866	0.0000	1.0000
JC69I	56	-313.4795	738.9589	0.0000	1.0000
JC69G	56	-313.5784	739.1568	0.0000	1.0000
JC69	55	-314.6792	739.3584	0.0000	1.0000
JC69IG	57	-313.5876	741.1752	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-307.5573	882.7456	0.4599	0.4599
K2PI	57	-306.2185	884.8471	0.1608	0.6207
K2PG	57	-306.2744	884.9588	0.1521	0.7728
HKY	59	-301.7874	885.5430	0.1136	0.8863
HKYI	60	-300.4975	887.7425	0.0378	0.9241
HKYG	60	-300.5863	887.9199	0.0346	0.9587
K2PIG	58	-306.2838	889.7569	0.0138	0.9726
SYM	60	-302.2871	891.3215	0.0063	0.9789
GTR	63	-295.4185	891.9219	0.0047	0.9836
JC69	55	-314.6792	892.2102	0.0040	0.9876
HKYIG	61	-300.5975	892.7215	0.0031	0.9907
SYMI	61	-301.1943	893.9151	0.0017	0.9925
SYMG	61	-301.3065	894.1395	0.0015	0.9940
JC69I	56	-313.4795	894.5898	0.0012	0.9952
JC69G	56	-313.5784	894.7877	0.0011	0.9964
F81	58	-308.9140	895.0173	0.0010	0.9974
GTRI	64	-294.5955	895.0550	0.0010	0.9983
GTRG	64	-294.7740	895.4120	0.0008	0.9991
F81I	59	-307.8389	897.6460	0.0003	0.9994
F81G	59	-307.9740	897.9164	0.0002	0.9996
SYMIG	62	-301.3119	898.9294	0.0001	0.9998
JC69IG	57	-313.5876	899.5853	0.0001	0.9999
GTRIG	65	-294.6007	899.8444	0.0001	1.0000
F81IG	60	-307.8433	902.4340	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32247.phylip = 1 - 119;
 Partition Dummy = 1:183-phylip/chr1_32247.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32247.phylip = 1 - 119;
 Partition Dummy = 1:183-phylip/chr1_32247.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_32247.phylip = 1 - 119;
 Partition Dummy = 1:183-phylip/chr1_32247.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

