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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRG
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
GTRG	64	-592.7977	1362.2505	0.2291	0.2291
HKYI	60	-600.4705	1362.7696	0.1767	0.4057
HKYG	60	-600.4774	1362.7834	0.1755	0.5812
HKYIG	61	-598.8318	1363.1349	0.1472	0.7284
GTRI	64	-593.2862	1363.2275	0.1405	0.8689
GTRIG	65	-591.4582	1363.3870	0.1298	0.9987
F81I	59	-608.2919	1374.8112	0.0004	0.9991
F81G	59	-608.4628	1375.1528	0.0004	0.9995
F81IG	60	-606.7661	1375.3609	0.0003	0.9998
K2PI	57	-613.6788	1378.5036	0.0001	0.9999
K2PIG	58	-612.1350	1378.9367	0.0001	0.9999
K2PG	57	-614.0286	1379.2033	0.0000	1.0000
GTR	63	-605.9600	1384.8037	0.0000	1.0000
SYMI	61	-610.3433	1386.1578	0.0000	1.0000
SYMG	61	-610.4005	1386.2723	0.0000	1.0000
HKY	59	-614.2723	1386.7719	0.0000	1.0000
SYMIG	62	-609.3125	1387.7810	0.0000	1.0000
JC69I	56	-621.2069	1390.0785	0.0000	1.0000
JC69IG	57	-619.7412	1390.6284	0.0000	1.0000
JC69G	56	-621.5296	1390.7240	0.0000	1.0000
F81	58	-623.0468	1400.7604	0.0000	1.0000
K2P	56	-629.5765	1406.8178	0.0000	1.0000
SYM	60	-624.1844	1410.1974	0.0000	1.0000
JC69	55	-636.9445	1418.1112	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-591.4582	1312.9164	0.4506	0.4506
GTRG	64	-592.7977	1313.5955	0.3209	0.7714
GTRI	64	-593.2862	1314.5725	0.1969	0.9683
HKYIG	61	-598.8318	1319.6636	0.0154	0.9837
HKYI	60	-600.4705	1320.9410	0.0082	0.9919
HKYG	60	-600.4774	1320.9548	0.0081	1.0000
F81IG	60	-606.7661	1333.5323	0.0000	1.0000
F81I	59	-608.2919	1334.5839	0.0000	1.0000
F81G	59	-608.4628	1334.9256	0.0000	1.0000
GTR	63	-605.9600	1337.9200	0.0000	1.0000
K2PIG	58	-612.1350	1340.2701	0.0000	1.0000
K2PI	57	-613.6788	1341.3575	0.0000	1.0000
K2PG	57	-614.0286	1342.0572	0.0000	1.0000
SYMIG	62	-609.3125	1342.6250	0.0000	1.0000
SYMI	61	-610.3433	1342.6865	0.0000	1.0000
SYMG	61	-610.4005	1342.8010	0.0000	1.0000
HKY	59	-614.2723	1346.5446	0.0000	1.0000
JC69IG	57	-619.7412	1353.4824	0.0000	1.0000
JC69I	56	-621.2069	1354.4137	0.0000	1.0000
JC69G	56	-621.5296	1355.0592	0.0000	1.0000
F81	58	-623.0468	1362.0937	0.0000	1.0000
SYM	60	-624.1844	1368.3688	0.0000	1.0000
K2P	56	-629.5765	1371.1530	0.0000	1.0000
JC69	55	-636.9445	1383.8890	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-600.4705	1528.7709	0.4121	0.4121
HKYG	60	-600.4774	1528.7847	0.4092	0.8213
HKYIG	61	-598.8318	1530.9574	0.1381	0.9594
GTRG	64	-592.7977	1535.2807	0.0159	0.9753
GTRI	64	-593.2862	1536.2577	0.0098	0.9850
GTRIG	65	-591.4582	1538.0655	0.0040	0.9890
K2PI	57	-613.6788	1538.7960	0.0027	0.9917
F81I	59	-608.2919	1538.9500	0.0025	0.9943
F81G	59	-608.4628	1539.2916	0.0021	0.9964
K2PG	57	-614.0286	1539.4956	0.0019	0.9983
K2PIG	58	-612.1350	1541.1723	0.0008	0.9992
F81IG	60	-606.7661	1541.3622	0.0008	0.9999
JC69I	56	-621.2069	1548.3883	0.0000	1.0000
JC69G	56	-621.5296	1549.0338	0.0000	1.0000
HKY	59	-614.2723	1550.9107	0.0000	1.0000
JC69IG	57	-619.7412	1550.9208	0.0000	1.0000
SYMI	61	-610.3433	1553.9802	0.0000	1.0000
SYMG	61	-610.4005	1554.0948	0.0000	1.0000
GTR	63	-605.9600	1556.1414	0.0000	1.0000
SYMIG	62	-609.3125	1557.3825	0.0000	1.0000
F81	58	-623.0468	1562.9959	0.0000	1.0000
K2P	56	-629.5765	1565.1276	0.0000	1.0000
JC69	55	-636.9445	1574.3997	0.0000	1.0000
SYM	60	-624.1844	1576.1987	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1623_hg18.phylip = 1 - 236;
 Partition Dummy = 1:261-aggressive-phylip/1623_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1623_hg18.phylip = 1 - 236;
 Partition Dummy = 1:261-aggressive-phylip/1623_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 BIC model (HKYI)]

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

