
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
10:59 08/01/2011

Input data from file "261-aggressive-phylip/10151_hg18.phylip" (29 taxa, 163 characters)

Minimum AIC  model: HKYG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-569.3195	1315.6103	0.4731	0.4731
K2PI	57	-569.3379	1315.6471	0.4645	0.9376
K2PIG	58	-569.3188	1320.4453	0.0422	0.9798
HKYG	60	-565.8940	1323.5526	0.0089	0.9887
K2P	56	-576.0542	1324.3348	0.0060	0.9947
HKYI	60	-566.5364	1324.8375	0.0047	0.9994
HKYIG	61	-566.4203	1329.7317	0.0004	0.9998
HKY	59	-572.8799	1332.4977	0.0001	0.9999
SYMI	61	-568.9147	1334.7205	0.0000	1.0000
SYMG	61	-569.0807	1335.0524	0.0000	1.0000
SYMIG	62	-569.0809	1340.2818	0.0000	1.0000
JC69G	56	-585.2083	1342.6431	0.0000	1.0000
GTRI	64	-565.5704	1344.0388	0.0000	1.0000
GTRG	64	-565.5761	1344.0502	0.0000	1.0000
JC69I	56	-586.0542	1344.3348	0.0000	1.0000
SYM	60	-576.7220	1345.2087	0.0000	1.0000
JC69IG	57	-586.0916	1349.1547	0.0000	1.0000
GTRIG	65	-565.5821	1349.6178	0.0000	1.0000
JC69	55	-591.6250	1350.8201	0.0000	1.0000
F81G	59	-582.0554	1350.8486	0.0000	1.0000
F81I	59	-583.1839	1353.1057	0.0000	1.0000
GTR	63	-573.0541	1353.5627	0.0000	1.0000
F81IG	60	-583.1802	1358.1250	0.0000	1.0000
F81	58	-588.3166	1358.4409	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-565.8940	1251.7879	0.2957	0.2957
K2PG	57	-569.3195	1252.6389	0.1932	0.4889
K2PI	57	-569.3379	1252.6757	0.1897	0.6786
HKYI	60	-566.5364	1253.0728	0.1555	0.8341
K2PIG	58	-569.3188	1254.6376	0.0711	0.9053
HKYIG	61	-566.4203	1254.8407	0.0643	0.9695
GTRI	64	-565.5704	1259.1408	0.0075	0.9770
GTRG	64	-565.5761	1259.1522	0.0074	0.9844
SYMI	61	-568.9147	1259.8294	0.0053	0.9897
SYMG	61	-569.0807	1260.1613	0.0045	0.9942
GTRIG	65	-565.5821	1261.1642	0.0027	0.9970
SYMIG	62	-569.0809	1262.1618	0.0017	0.9986
HKY	59	-572.8799	1263.7599	0.0007	0.9994
K2P	56	-576.0542	1264.1083	0.0006	1.0000
GTR	63	-573.0541	1272.1082	0.0000	1.0000
SYM	60	-576.7220	1273.4440	0.0000	1.0000
F81G	59	-582.0554	1282.1107	0.0000	1.0000
JC69G	56	-585.2083	1282.4166	0.0000	1.0000
JC69I	56	-586.0542	1284.1084	0.0000	1.0000
F81I	59	-583.1839	1284.3679	0.0000	1.0000
JC69IG	57	-586.0916	1286.1833	0.0000	1.0000
F81IG	60	-583.1802	1286.3603	0.0000	1.0000
F81	58	-588.3166	1292.6332	0.0000	1.0000
JC69	55	-591.6250	1293.2500	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-569.3195	1428.9827	0.4765	0.4765
K2PI	57	-569.3379	1429.0195	0.4678	0.9442
K2PIG	58	-569.3188	1434.0751	0.0373	0.9816
K2P	56	-576.0542	1437.3584	0.0072	0.9888
HKYG	60	-565.8940	1437.4129	0.0070	0.9958
HKYI	60	-566.5364	1438.6978	0.0037	0.9995
HKYIG	61	-566.4203	1443.5594	0.0003	0.9999
HKY	59	-572.8799	1446.2911	0.0001	0.9999
SYMI	61	-568.9147	1448.5481	0.0000	1.0000
SYMG	61	-569.0807	1448.8801	0.0000	1.0000
SYMIG	62	-569.0809	1453.9743	0.0000	1.0000
JC69G	56	-585.2083	1455.6667	0.0000	1.0000
GTRI	64	-565.5704	1457.1409	0.0000	1.0000
GTRG	64	-565.5761	1457.1522	0.0000	1.0000
JC69I	56	-586.0542	1457.3584	0.0000	1.0000
SYM	60	-576.7220	1459.0690	0.0000	1.0000
GTRIG	65	-565.5821	1462.2580	0.0000	1.0000
JC69IG	57	-586.0916	1462.5271	0.0000	1.0000
JC69	55	-591.6250	1463.4062	0.0000	1.0000
F81G	59	-582.0554	1464.6420	0.0000	1.0000
F81I	59	-583.1839	1466.8991	0.0000	1.0000
GTR	63	-573.0541	1467.0145	0.0000	1.0000
F81IG	60	-583.1802	1471.9853	0.0000	1.0000
F81	58	-588.3166	1472.0707	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10151_hg18.phylip = 1 - 163;
 Partition Dummy = 1:261-aggressive-phylip/10151_hg18.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 AICc model (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10151_hg18.phylip = 1 - 163;
 Partition Dummy = 1:261-aggressive-phylip/10151_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10151_hg18.phylip = 1 - 163;
 Partition Dummy = 1:261-aggressive-phylip/10151_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

