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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	60	-484.7287	1137.9341	0.4164	0.4164
HKYI	60	-485.1101	1138.6970	0.2844	0.7008
HKY	59	-487.4221	1139.4232	0.1978	0.8986
HKYIG	61	-484.6279	1141.6825	0.0639	0.9625
GTRG	64	-480.1171	1144.8329	0.0132	0.9757
GTRI	64	-480.6176	1145.8339	0.0080	0.9837
GTR	63	-482.8750	1146.2366	0.0066	0.9903
F81G	59	-491.4783	1147.5356	0.0034	0.9937
F81I	59	-491.8160	1148.2109	0.0024	0.9961
GTRIG	65	-480.0130	1148.7931	0.0018	0.9980
F81	58	-494.2298	1149.1916	0.0015	0.9995
F81IG	60	-491.3635	1151.2039	0.0005	1.0000
K2PG	57	-505.4512	1167.8375	0.0000	1.0000
K2PI	57	-505.7490	1168.4331	0.0000	1.0000
K2P	56	-508.2697	1169.7265	0.0000	1.0000
K2PIG	58	-505.3406	1171.4132	0.0000	1.0000
JC69G	56	-510.3995	1173.9862	0.0000	1.0000
JC69I	56	-510.7007	1174.5884	0.0000	1.0000
JC69	55	-513.2208	1175.9288	0.0000	1.0000
JC69IG	57	-510.2794	1177.4938	0.0000	1.0000
SYMG	61	-503.5619	1179.5505	0.0000	1.0000
SYMI	61	-503.9054	1180.2375	0.0000	1.0000
SYM	60	-506.1080	1180.6928	0.0000	1.0000
SYMIG	62	-503.4735	1183.3766	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-480.1171	1088.2343	0.2920	0.2920
GTRI	64	-480.6176	1089.2352	0.1770	0.4691
HKYG	60	-484.7287	1089.4573	0.1584	0.6275
GTRIG	65	-480.0130	1090.0260	0.1192	0.7467
HKYI	60	-485.1101	1090.2201	0.1082	0.8549
HKYIG	61	-484.6279	1091.2559	0.0645	0.9194
GTR	63	-482.8750	1091.7501	0.0503	0.9697
HKY	59	-487.4221	1092.8443	0.0291	0.9988
F81G	59	-491.4783	1100.9567	0.0005	0.9993
F81I	59	-491.8160	1101.6320	0.0004	0.9997
F81IG	60	-491.3635	1102.7271	0.0002	0.9999
F81	58	-494.2298	1104.4596	0.0001	1.0000
K2PG	57	-505.4512	1124.9025	0.0000	1.0000
K2PI	57	-505.7490	1125.4981	0.0000	1.0000
K2PIG	58	-505.3406	1126.6812	0.0000	1.0000
K2P	56	-508.2697	1128.5394	0.0000	1.0000
SYMG	61	-503.5619	1129.1238	0.0000	1.0000
SYMI	61	-503.9054	1129.8109	0.0000	1.0000
SYMIG	62	-503.4735	1130.9471	0.0000	1.0000
SYM	60	-506.1080	1132.2160	0.0000	1.0000
JC69G	56	-510.3995	1132.7991	0.0000	1.0000
JC69I	56	-510.7007	1133.4013	0.0000	1.0000
JC69IG	57	-510.2794	1134.5588	0.0000	1.0000
JC69	55	-513.2208	1136.4416	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-484.7287	1290.8525	0.3577	0.3577
HKY	59	-487.4221	1290.8829	0.3523	0.7099
HKYI	60	-485.1101	1291.6153	0.2442	0.9542
HKYIG	61	-484.6279	1296.0076	0.0272	0.9813
F81G	59	-491.4783	1298.9953	0.0061	0.9874
F81	58	-494.2298	1299.1416	0.0057	0.9931
F81I	59	-491.8160	1299.6706	0.0044	0.9974
GTRG	64	-480.1171	1303.0558	0.0008	0.9982
GTR	63	-482.8750	1303.2150	0.0007	0.9990
GTRI	64	-480.6176	1304.0568	0.0005	0.9995
F81IG	60	-491.3635	1304.1223	0.0005	0.9999
GTRIG	65	-480.0130	1308.2041	0.0001	1.0000
K2PG	57	-505.4512	1316.2279	0.0000	1.0000
K2P	56	-508.2697	1316.5083	0.0000	1.0000
K2PI	57	-505.7490	1316.8235	0.0000	1.0000
JC69G	56	-510.3995	1320.7679	0.0000	1.0000
JC69	55	-513.2208	1321.0539	0.0000	1.0000
K2PIG	58	-505.3406	1321.3632	0.0000	1.0000
JC69I	56	-510.7007	1321.3702	0.0000	1.0000
JC69IG	57	-510.2794	1325.8842	0.0000	1.0000
SYM	60	-506.1080	1333.6112	0.0000	1.0000
SYMG	61	-503.5619	1333.8756	0.0000	1.0000
SYMI	61	-503.9054	1334.5626	0.0000	1.0000
SYMIG	62	-503.4735	1339.0554	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/chr7_10353.phylip = 1 - 212;
 Partition Dummy = 1:183-phylip/chr7_10353.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_10353.phylip = 1 - 212;
 Partition Dummy = 1:183-phylip/chr7_10353.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr7_10353.phylip = 1 - 212;
 Partition Dummy = 1:183-phylip/chr7_10353.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

