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

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

Minimum AIC  model: K2PIG
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	57	-329.2093	902.0657	0.5030	0.5030
K2PG	57	-329.2689	902.1848	0.4739	0.9769
K2PIG	58	-327.7568	908.3937	0.0213	0.9982
JC69I	56	-339.9582	914.6856	0.0009	0.9991
JC69G	56	-340.0018	914.7728	0.0009	1.0000
JC69IG	57	-338.6015	920.8501	0.0000	1.0000
K2P	56	-345.3147	925.3986	0.0000	1.0000
HKYI	60	-328.9232	930.3463	0.0000	1.0000
HKYG	60	-328.9830	930.4660	0.0000	1.0000
SYMG	61	-325.9860	934.9081	0.0000	1.0000
SYMI	61	-326.0260	934.9882	0.0000	1.0000
JC69	55	-355.8083	937.8430	0.0000	1.0000
HKYIG	61	-327.5536	938.0434	0.0000	1.0000
F81I	59	-339.4268	941.3434	0.0000	1.0000
F81G	59	-339.4970	941.4838	0.0000	1.0000
SYMIG	62	-324.7804	943.3868	0.0000	1.0000
F81IG	60	-338.1500	948.8000	0.0000	1.0000
HKY	59	-344.9226	952.3350	0.0000	1.0000
SYM	60	-341.0380	954.5761	0.0000	1.0000
F81	58	-355.4772	963.8344	0.0000	1.0000
GTRG	64	-325.1950	967.4810	0.0000	1.0000
GTRI	64	-325.2570	967.6050	0.0000	1.0000
GTRIG	65	-324.1620	977.8589	0.0000	1.0000
GTR	63	-340.5051	986.2102	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	58	-327.7568	771.5137	0.2935	0.2935
K2PI	57	-329.2093	772.4187	0.1867	0.4802
K2PG	57	-329.2689	772.5377	0.1759	0.6561
SYMIG	62	-324.7804	773.5608	0.1055	0.7616
SYMG	61	-325.9860	773.9719	0.0859	0.8475
SYMI	61	-326.0260	774.0520	0.0825	0.9300
HKYIG	61	-327.5536	777.1072	0.0179	0.9479
HKYI	60	-328.9232	777.8463	0.0124	0.9603
HKYG	60	-328.9830	777.9660	0.0117	0.9719
GTRIG	65	-324.1620	778.3240	0.0097	0.9817
GTRG	64	-325.1950	778.3901	0.0094	0.9911
GTRI	64	-325.2570	778.5141	0.0089	1.0000
JC69IG	57	-338.6015	791.2030	0.0000	1.0000
JC69I	56	-339.9582	791.9163	0.0000	1.0000
JC69G	56	-340.0018	792.0035	0.0000	1.0000
F81IG	60	-338.1500	796.3000	0.0000	1.0000
F81I	59	-339.4268	796.8536	0.0000	1.0000
F81G	59	-339.4970	796.9940	0.0000	1.0000
SYM	60	-341.0380	802.0761	0.0000	1.0000
K2P	56	-345.3147	802.6293	0.0000	1.0000
GTR	63	-340.5051	807.0102	0.0000	1.0000
HKY	59	-344.9226	807.8452	0.0000	1.0000
JC69	55	-355.8083	821.6166	0.0000	1.0000
F81	58	-355.4772	826.9544	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-329.2093	925.8255	0.4238	0.4238
K2PG	57	-329.2689	925.9446	0.3993	0.8231
K2PIG	58	-327.7568	927.6119	0.1735	0.9966
SYMG	61	-325.9860	938.1441	0.0009	0.9975
SYMI	61	-326.0260	938.2242	0.0009	0.9983
HKYI	60	-328.9232	939.3272	0.0005	0.9988
HKYG	60	-328.9830	939.4468	0.0005	0.9993
SYMIG	62	-324.7804	940.4243	0.0003	0.9996
HKYIG	61	-327.5536	941.2795	0.0002	0.9998
JC69I	56	-339.9582	942.6318	0.0001	0.9999
JC69G	56	-340.0018	942.7190	0.0001	1.0000
JC69IG	57	-338.6015	944.6099	0.0000	1.0000
GTRG	64	-325.1950	950.6363	0.0000	1.0000
GTRI	64	-325.2570	950.7603	0.0000	1.0000
GTRIG	65	-324.1620	953.2616	0.0000	1.0000
K2P	56	-345.3147	953.3448	0.0000	1.0000
F81I	59	-339.4268	955.6432	0.0000	1.0000
F81G	59	-339.4970	955.7835	0.0000	1.0000
F81IG	60	-338.1500	957.7808	0.0000	1.0000
SYM	60	-341.0380	963.5569	0.0000	1.0000
HKY	59	-344.9226	966.6347	0.0000	1.0000
JC69	55	-355.8083	969.6407	0.0000	1.0000
GTR	63	-340.5051	976.5651	0.0000	1.0000
F81	58	-355.4772	983.0525	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_5495.phylip = 1 - 109;
 Partition Dummy = 1:183-phylip/chrZ_5495.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_5495.phylip = 1 - 109;
 Partition Dummy = 1:183-phylip/chrZ_5495.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_5495.phylip = 1 - 109;
 Partition Dummy = 1:183-phylip/chrZ_5495.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

