
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:33 04/04/2011

Input data from file "591-tomodel-phylip/chr13_3033.phylip" (7 taxa, 264 characters)

Minimum AIC  model: HKYI
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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	16	-576.9209	1188.0443	0.4009	0.4009
HKYG	16	-576.9272	1188.0568	0.3984	0.7994
HKYIG	17	-576.9325	1190.3529	0.1264	0.9258
GTRG	20	-574.8552	1193.1671	0.0310	0.9567
GTRI	20	-574.8580	1193.1728	0.0309	0.9876
GTRIG	21	-574.8533	1195.5249	0.0095	0.9971
HKY	15	-584.3105	1200.5565	0.0008	0.9979
K2PI	13	-586.6592	1200.7744	0.0007	0.9986
K2PG	13	-586.7428	1200.9416	0.0006	0.9992
K2PIG	14	-586.6594	1203.0054	0.0002	0.9994
F81I	15	-585.9590	1203.8534	0.0001	0.9996
F81G	15	-585.9871	1203.9098	0.0001	0.9997
GTR	19	-581.4828	1204.0805	0.0001	0.9998
F81IG	16	-585.9592	1206.1208	0.0000	0.9999
SYMI	17	-584.8643	1206.2163	0.0000	0.9999
SYMG	17	-584.8645	1206.2167	0.0000	1.0000
SYMIG	18	-584.8644	1208.5207	0.0000	1.0000
K2P	12	-594.8480	1214.9390	0.0000	1.0000
JC69I	12	-595.3403	1215.9237	0.0000	1.0000
JC69G	12	-595.3949	1216.0328	0.0000	1.0000
F81	14	-593.6125	1216.9118	0.0000	1.0000
SYM	16	-591.9641	1218.1306	0.0000	1.0000
JC69IG	13	-595.3403	1218.1366	0.0000	1.0000
JC69	11	-603.2649	1229.5773	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYI	16	-576.9209	1185.8419	0.3696	0.3696
HKYG	16	-576.9272	1185.8544	0.3673	0.7370
HKYIG	17	-576.9325	1187.8651	0.1344	0.8714
GTRG	20	-574.8552	1189.7103	0.0534	0.9248
GTRI	20	-574.8580	1189.7160	0.0533	0.9781
GTRIG	21	-574.8533	1191.7067	0.0197	0.9978
HKY	15	-584.3105	1198.6211	0.0006	0.9984
K2PI	13	-586.6592	1199.3184	0.0004	0.9988
K2PG	13	-586.7428	1199.4856	0.0004	0.9993
GTR	19	-581.4828	1200.9657	0.0002	0.9994
K2PIG	14	-586.6594	1201.3187	0.0002	0.9996
F81I	15	-585.9590	1201.9179	0.0001	0.9997
F81G	15	-585.9871	1201.9743	0.0001	0.9998
SYMI	17	-584.8643	1203.7285	0.0000	0.9999
SYMG	17	-584.8645	1203.7289	0.0000	0.9999
F81IG	16	-585.9592	1203.9184	0.0000	1.0000
SYMIG	18	-584.8644	1205.7289	0.0000	1.0000
K2P	12	-594.8480	1213.6960	0.0000	1.0000
JC69I	12	-595.3403	1214.6807	0.0000	1.0000
JC69G	12	-595.3949	1214.7898	0.0000	1.0000
F81	14	-593.6125	1215.2251	0.0000	1.0000
SYM	16	-591.9641	1215.9282	0.0000	1.0000
JC69IG	13	-595.3403	1216.6806	0.0000	1.0000
JC69	11	-603.2649	1228.5297	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-576.9209	1243.0571	0.3887	0.3887
HKYG	16	-576.9272	1243.0696	0.3863	0.7750
K2PI	13	-586.6592	1245.8058	0.0983	0.8733
K2PG	13	-586.7428	1245.9730	0.0905	0.9638
HKYIG	17	-576.9325	1248.6562	0.0236	0.9874
K2PIG	14	-586.6594	1251.3820	0.0061	0.9935
HKY	15	-584.3105	1252.2603	0.0039	0.9974
F81I	15	-585.9590	1255.5572	0.0008	0.9981
F81G	15	-585.9871	1255.6135	0.0007	0.9988
K2P	12	-594.8480	1256.6074	0.0004	0.9993
JC69I	12	-595.3403	1257.5920	0.0003	0.9996
JC69G	12	-595.3949	1257.7011	0.0003	0.9998
F81IG	16	-585.9592	1261.1336	0.0000	0.9999
GTRG	20	-574.8552	1261.2293	0.0000	0.9999
GTRI	20	-574.8580	1261.2350	0.0000	1.0000
JC69IG	13	-595.3403	1263.1679	0.0000	1.0000
SYMI	17	-584.8643	1264.5197	0.0000	1.0000
SYMG	17	-584.8645	1264.5201	0.0000	1.0000
F81	14	-593.6125	1265.2883	0.0000	1.0000
GTRIG	21	-574.8533	1266.8016	0.0000	1.0000
JC69	11	-603.2649	1267.8652	0.0000	1.0000
GTR	19	-581.4828	1268.9087	0.0000	1.0000
SYMIG	18	-584.8644	1270.0960	0.0000	1.0000
SYM	16	-591.9641	1273.1433	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 591-tomodel-phylip/chr13_3033.phylip = 1 - 264;
 Partition Dummy = 1:591-tomodel-phylip/chr13_3033.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 591-tomodel-phylip/chr13_3033.phylip = 1 - 264;
 Partition Dummy = 1:591-tomodel-phylip/chr13_3033.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr13_3033.phylip = 1 - 264;
 Partition Dummy = 1:591-tomodel-phylip/chr13_3033.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;


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

