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

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

Minimum AIC  model: GTRG
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	-609.9166	1386.3093	0.6270	0.6270
K2PI	57	-611.0511	1388.5785	0.2016	0.8286
K2PIG	58	-609.6551	1390.0622	0.0960	0.9246
HKYG	60	-606.1953	1391.9029	0.0382	0.9628
HKYI	60	-607.0583	1393.6288	0.0161	0.9790
SYMG	61	-605.3713	1394.7426	0.0092	0.9882
HKYIG	61	-605.9507	1395.9015	0.0052	0.9934
SYMI	61	-606.4989	1396.9977	0.0030	0.9964
K2P	56	-618.1018	1398.4713	0.0014	0.9978
SYMIG	62	-605.1703	1398.9027	0.0012	0.9990
GTRG	64	-601.2323	1400.3805	0.0006	0.9995
GTRI	64	-602.1390	1402.1940	0.0002	0.9998
HKY	59	-614.0318	1403.1605	0.0001	0.9999
GTRIG	65	-600.9980	1404.7078	0.0001	1.0000
SYM	60	-613.4617	1406.4356	0.0000	1.0000
JC69G	56	-624.1460	1410.5597	0.0000	1.0000
GTR	63	-609.3828	1411.9656	0.0000	1.0000
JC69I	56	-625.1783	1412.6243	0.0000	1.0000
JC69IG	57	-623.9494	1414.3750	0.0000	1.0000
F81G	59	-620.2970	1415.6908	0.0000	1.0000
F81I	59	-621.1182	1417.3332	0.0000	1.0000
F81IG	60	-620.0883	1419.6888	0.0000	1.0000
JC69	55	-631.9423	1422.0096	0.0000	1.0000
F81	58	-627.9723	1426.6966	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-601.2323	1330.4645	0.2862	0.2862
GTRIG	65	-600.9980	1331.9959	0.1331	0.4193
GTRI	64	-602.1390	1332.2780	0.1156	0.5348
HKYG	60	-606.1953	1332.3907	0.1092	0.6441
SYMG	61	-605.3713	1332.7426	0.0916	0.7357
K2PG	57	-609.9166	1333.8331	0.0531	0.7888
HKYIG	61	-605.9507	1333.9015	0.0513	0.8401
HKYI	60	-607.0583	1334.1166	0.0461	0.8862
SYMIG	62	-605.1703	1334.3407	0.0412	0.9274
SYMI	61	-606.4989	1334.9977	0.0297	0.9571
K2PIG	58	-609.6551	1335.3102	0.0254	0.9825
K2PI	57	-611.0511	1336.1023	0.0171	0.9995
GTR	63	-609.3828	1344.7656	0.0002	0.9998
HKY	59	-614.0318	1346.0637	0.0001	0.9999
SYM	60	-613.4617	1346.9234	0.0001	1.0000
K2P	56	-618.1018	1348.2036	0.0000	1.0000
F81G	59	-620.2970	1358.5940	0.0000	1.0000
F81IG	60	-620.0883	1360.1766	0.0000	1.0000
F81I	59	-621.1182	1360.2364	0.0000	1.0000
JC69G	56	-624.1460	1360.2919	0.0000	1.0000
JC69IG	57	-623.9494	1361.8988	0.0000	1.0000
JC69I	56	-625.1783	1362.3565	0.0000	1.0000
F81	58	-627.9723	1371.9446	0.0000	1.0000
JC69	55	-631.9423	1373.8846	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-609.9166	1517.0845	0.6895	0.6895
K2PI	57	-611.0511	1519.3536	0.2217	0.9112
K2PIG	58	-609.6551	1521.7765	0.0660	0.9772
HKYG	60	-606.1953	1525.2868	0.0114	0.9886
HKYI	60	-607.0583	1527.0128	0.0048	0.9935
K2P	56	-618.1018	1528.2400	0.0026	0.9961
SYMG	61	-605.3713	1528.8537	0.0019	0.9980
HKYIG	61	-605.9507	1530.0125	0.0011	0.9991
SYMI	61	-606.4989	1531.1088	0.0006	0.9997
SYMIG	62	-605.1703	1533.6667	0.0002	0.9998
HKY	59	-614.0318	1535.7449	0.0001	0.9999
GTRG	64	-601.2323	1536.2204	0.0000	1.0000
GTRI	64	-602.1390	1538.0339	0.0000	1.0000
SYM	60	-613.4617	1539.8195	0.0000	1.0000
JC69G	56	-624.1460	1540.3283	0.0000	1.0000
GTRIG	65	-600.9980	1540.9667	0.0000	1.0000
JC69I	56	-625.1783	1542.3929	0.0000	1.0000
JC69IG	57	-623.9494	1545.1501	0.0000	1.0000
GTR	63	-609.3828	1547.3066	0.0000	1.0000
F81G	59	-620.2970	1548.2752	0.0000	1.0000
F81I	59	-621.1182	1549.9176	0.0000	1.0000
JC69	55	-631.9423	1550.7061	0.0000	1.0000
F81IG	60	-620.0883	1553.0727	0.0000	1.0000
F81	58	-627.9723	1558.4108	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 261-aggressive-phylip/7621_hg18.phylip = 1 - 184;
 Partition Dummy = 1:261-aggressive-phylip/7621_hg18.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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/7621_hg18.phylip = 1 - 184;
 Partition Dummy = 1:261-aggressive-phylip/7621_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/7621_hg18.phylip = 1 - 184;
 Partition Dummy = 1:261-aggressive-phylip/7621_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

