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

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

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
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
HKYI	60	-442.8893	1062.0863	0.3966	0.3966
HKYG	60	-443.1766	1062.6609	0.2976	0.6942
K2PI	57	-450.6327	1064.9796	0.0933	0.7875
HKYIG	61	-442.4787	1065.5931	0.0687	0.8562
K2PG	57	-450.9538	1065.6219	0.0677	0.9239
F81I	59	-448.0140	1068.0739	0.0199	0.9438
K2PIG	58	-450.1346	1068.1176	0.0194	0.9633
F81G	59	-448.2515	1068.5487	0.0157	0.9789
JC69I	56	-455.0594	1069.7606	0.0085	0.9875
JC69G	56	-455.3333	1070.3083	0.0065	0.9940
F81IG	60	-447.5815	1071.4708	0.0036	0.9976
JC69IG	57	-454.5619	1072.8381	0.0018	0.9994
GTRI	64	-441.7479	1077.5276	0.0002	0.9996
GTRG	64	-441.9360	1077.9037	0.0001	0.9998
HKY	59	-453.2678	1078.5814	0.0001	0.9999
SYMI	61	-449.9615	1080.5587	0.0000	0.9999
SYMG	61	-450.2297	1081.0951	0.0000	0.9999
K2P	56	-461.0350	1081.7119	0.0000	1.0000
GTRIG	65	-441.5623	1081.7647	0.0000	1.0000
F81	58	-457.6268	1083.1021	0.0000	1.0000
SYMIG	62	-449.6740	1084.3793	0.0000	1.0000
JC69	55	-464.8077	1085.2451	0.0000	1.0000
GTR	63	-452.0548	1093.6056	0.0000	1.0000
SYM	60	-460.3128	1096.9334	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	60	-442.8893	1005.7786	0.4008	0.4008
HKYG	60	-443.1766	1006.3532	0.3007	0.7015
HKYIG	61	-442.4787	1006.9575	0.2223	0.9238
GTRI	64	-441.7479	1011.4958	0.0230	0.9468
GTRG	64	-441.9360	1011.8720	0.0190	0.9658
GTRIG	65	-441.5623	1013.1247	0.0102	0.9760
F81I	59	-448.0140	1014.0281	0.0065	0.9825
F81G	59	-448.2515	1014.5029	0.0051	0.9876
F81IG	60	-447.5815	1015.1631	0.0037	0.9913
K2PI	57	-450.6327	1015.2653	0.0035	0.9948
K2PG	57	-450.9538	1015.9076	0.0025	0.9973
K2PIG	58	-450.1346	1016.2691	0.0021	0.9994
SYMI	61	-449.9615	1021.9230	0.0001	0.9995
JC69I	56	-455.0594	1022.1188	0.0001	0.9997
SYMG	61	-450.2297	1022.4594	0.0001	0.9997
JC69G	56	-455.3333	1022.6665	0.0001	0.9998
JC69IG	57	-454.5619	1023.1239	0.0001	0.9999
SYMIG	62	-449.6740	1023.3481	0.0001	1.0000
HKY	59	-453.2678	1024.5356	0.0000	1.0000
GTR	63	-452.0548	1030.1095	0.0000	1.0000
F81	58	-457.6268	1031.2536	0.0000	1.0000
K2P	56	-461.0350	1034.0701	0.0000	1.0000
JC69	55	-464.8077	1039.6155	0.0000	1.0000
SYM	60	-460.3128	1040.6257	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	57	-450.6327	1200.6449	0.2553	0.2553
HKYI	60	-442.8893	1200.9150	0.2230	0.4783
K2PG	57	-450.9538	1201.2872	0.1852	0.6635
HKYG	60	-443.1766	1201.4896	0.1673	0.8309
JC69I	56	-455.0594	1204.2461	0.0422	0.8730
JC69G	56	-455.3333	1204.7938	0.0321	0.9051
K2PIG	58	-450.1346	1204.9010	0.0304	0.9355
HKYIG	61	-442.4787	1205.3462	0.0243	0.9598
F81I	59	-448.0140	1205.9122	0.0183	0.9782
F81G	59	-448.2515	1206.3871	0.0145	0.9926
JC69IG	57	-454.5619	1208.5034	0.0050	0.9976
F81IG	60	-447.5815	1210.2995	0.0020	0.9997
K2P	56	-461.0350	1216.1974	0.0001	0.9998
HKY	59	-453.2678	1216.4198	0.0001	0.9999
JC69	55	-464.8077	1218.4905	0.0000	0.9999
GTRI	64	-441.7479	1219.6413	0.0000	0.9999
F81	58	-457.6268	1219.8855	0.0000	1.0000
GTRG	64	-441.9360	1220.0175	0.0000	1.0000
SYMI	61	-449.9615	1220.3117	0.0000	1.0000
SYMG	61	-450.2297	1220.8481	0.0000	1.0000
GTRIG	65	-441.5623	1224.5225	0.0000	1.0000
SYMIG	62	-449.6740	1224.9890	0.0000	1.0000
GTR	63	-452.0548	1235.0028	0.0000	1.0000
SYM	60	-460.3128	1235.7621	0.0000	1.0000



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

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

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


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

BEGIN MRBAYES;
 Charset 183-phylip/chr19_365.phylip = 1 - 191;
 Partition Dummy = 1:183-phylip/chr19_365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) 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/chr19_365.phylip = 1 - 191;
 Partition Dummy = 1:183-phylip/chr19_365.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

