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

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

Minimum AIC  model: F81
Minimum AICc model: JC69
Minimum BIC  model: JC69


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
JC69	55	-454.1657	1061.7117	0.3259	0.3259
F81	58	-448.3585	1061.9544	0.2887	0.6146
K2P	56	-453.5573	1064.3912	0.0854	0.6999
HKY	59	-447.7273	1064.7589	0.0710	0.7710
JC69I	56	-454.1408	1065.5581	0.0476	0.8186
JC69G	56	-454.1454	1065.5673	0.0474	0.8660
F81I	59	-448.3491	1066.0025	0.0381	0.9042
F81G	59	-448.3513	1066.0070	0.0381	0.9422
K2PI	57	-453.5297	1068.2880	0.0122	0.9544
K2PG	57	-453.5350	1068.2985	0.0121	0.9665
HKYI	60	-447.7176	1068.8658	0.0091	0.9756
HKYG	60	-447.7196	1068.8698	0.0091	0.9847
JC69IG	57	-454.1386	1069.5058	0.0066	0.9913
F81IG	60	-448.3488	1070.1282	0.0048	0.9961
K2PIG	58	-453.5284	1072.2941	0.0016	0.9978
HKYIG	61	-447.7172	1073.0521	0.0011	0.9989
GTR	63	-444.1223	1074.4237	0.0006	0.9995
SYM	60	-451.0879	1075.6065	0.0003	0.9998
GTRG	64	-444.1227	1078.8018	0.0001	0.9998
GTRI	64	-444.1238	1078.8040	0.0001	0.9999
SYMG	61	-451.0880	1079.7937	0.0000	0.9999
SYMI	61	-451.0891	1079.7958	0.0000	1.0000
GTRIG	65	-444.1240	1083.2480	0.0000	1.0000
SYMIG	62	-451.0892	1084.0450	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
F81	58	-448.3585	1012.7170	0.2349	0.2349
HKY	59	-447.7273	1013.4546	0.1624	0.3973
GTR	63	-444.1223	1014.2446	0.1094	0.5068
F81I	59	-448.3491	1014.6981	0.0872	0.5940
F81G	59	-448.3513	1014.7027	0.0870	0.6810
HKYI	60	-447.7176	1015.4351	0.0603	0.7414
HKYG	60	-447.7196	1015.4392	0.0602	0.8016
GTRG	64	-444.1227	1016.2454	0.0402	0.8418
GTRI	64	-444.1238	1016.2476	0.0402	0.8820
F81IG	60	-448.3488	1016.6975	0.0321	0.9141
HKYIG	61	-447.7172	1017.4344	0.0222	0.9364
GTRIG	65	-444.1240	1018.2480	0.0148	0.9511
JC69	55	-454.1657	1018.3314	0.0142	0.9653
K2P	56	-453.5573	1019.1146	0.0096	0.9749
JC69I	56	-454.1408	1020.2815	0.0053	0.9803
JC69G	56	-454.1454	1020.2907	0.0053	0.9856
K2PI	57	-453.5297	1021.0594	0.0036	0.9892
K2PG	57	-453.5350	1021.0699	0.0036	0.9928
SYM	60	-451.0879	1022.1758	0.0021	0.9949
JC69IG	57	-454.1386	1022.2772	0.0020	0.9969
K2PIG	58	-453.5284	1023.0567	0.0013	0.9982
SYMG	61	-451.0880	1024.1761	0.0008	0.9990
SYMI	61	-451.0891	1024.1781	0.0008	0.9997
SYMIG	62	-451.0892	1026.1784	0.0003	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
JC69	55	-454.1657	1199.1861	0.6866	0.6866
K2P	56	-453.5573	1203.2575	0.0897	0.7763
F81	58	-448.3585	1203.4365	0.0820	0.8583
JC69I	56	-454.1408	1204.4245	0.0500	0.9083
JC69G	56	-454.1454	1204.4337	0.0498	0.9581
HKY	59	-447.7273	1207.4623	0.0110	0.9691
K2PI	57	-453.5297	1208.4906	0.0066	0.9756
K2PG	57	-453.5350	1208.5012	0.0065	0.9821
F81I	59	-448.3491	1208.7059	0.0059	0.9880
F81G	59	-448.3513	1208.7104	0.0059	0.9939
JC69IG	57	-454.1386	1209.7084	0.0036	0.9974
HKYI	60	-447.7176	1212.7312	0.0008	0.9982
HKYG	60	-447.7196	1212.7352	0.0008	0.9990
K2PIG	58	-453.5284	1213.7762	0.0005	0.9995
F81IG	60	-448.3488	1213.9936	0.0004	0.9999
HKYIG	61	-447.7172	1218.0187	0.0001	1.0000
SYM	60	-451.0879	1219.4718	0.0000	1.0000
GTR	63	-444.1223	1221.4054	0.0000	1.0000
SYMG	61	-451.0880	1224.7604	0.0000	1.0000
SYMI	61	-451.0891	1224.7624	0.0000	1.0000
GTRG	64	-444.1227	1226.6945	0.0000	1.0000
GTRI	64	-444.1238	1226.6967	0.0000	1.0000
SYMIG	62	-451.0892	1230.0509	0.0000	1.0000
GTRIG	65	-444.1240	1231.9853	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_8059.phylip = 1 - 198;
 Partition Dummy = 1:183-phylip/chrZ_8059.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_8059.phylip = 1 - 198;
 Partition Dummy = 1:183-phylip/chrZ_8059.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_8059.phylip = 1 - 198;
 Partition Dummy = 1:183-phylip/chrZ_8059.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=1 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal);
END;


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

