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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-734.4015	1654.1824	0.8750	0.8750
GTRI	64	-737.2896	1659.9586	0.0487	0.9237
HKYG	60	-746.2473	1661.6221	0.0212	0.9449
GTRIG	65	-736.1636	1661.9106	0.0184	0.9633
K2PG	57	-752.5481	1662.5963	0.0130	0.9763
SYMG	61	-745.0246	1663.1573	0.0098	0.9862
GTR	63	-741.8749	1664.9826	0.0040	0.9901
SYMI	61	-746.1524	1665.4129	0.0032	0.9933
K2PI	57	-754.3039	1666.1079	0.0023	0.9956
K2PIG	58	-752.4400	1666.2045	0.0021	0.9977
SYMIG	62	-745.0077	1667.1583	0.0013	0.9990
HKYI	60	-750.3450	1669.8174	0.0004	0.9994
HKYIG	61	-748.5317	1670.1715	0.0003	0.9997
K2P	56	-758.9991	1671.7238	0.0001	0.9998
SYM	60	-751.3003	1671.7282	0.0001	1.0000
HKY	59	-754.3777	1673.9554	0.0000	1.0000
JC69G	56	-772.2063	1698.1382	0.0000	1.0000
F81G	59	-767.8210	1700.8419	0.0000	1.0000
JC69I	56	-773.6482	1701.0219	0.0000	1.0000
JC69IG	57	-772.1298	1701.7596	0.0000	1.0000
F81I	59	-769.5154	1704.2307	0.0000	1.0000
F81IG	60	-768.1160	1705.3596	0.0000	1.0000
JC69	55	-778.1579	1706.3157	0.0000	1.0000
F81	58	-773.1606	1707.6457	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-734.4015	1596.8030	0.8915	0.8915
GTRIG	65	-736.1636	1602.3273	0.0563	0.9479
GTRI	64	-737.2896	1602.5793	0.0496	0.9975
GTR	63	-741.8749	1609.7497	0.0014	0.9989
SYMG	61	-745.0246	1612.0492	0.0004	0.9993
HKYG	60	-746.2473	1612.4946	0.0003	0.9997
SYMIG	62	-745.0077	1614.0155	0.0002	0.9998
SYMI	61	-746.1524	1614.3048	0.0001	1.0000
HKYIG	61	-748.5317	1619.0634	0.0000	1.0000
K2PG	57	-752.5481	1619.0963	0.0000	1.0000
HKYI	60	-750.3450	1620.6899	0.0000	1.0000
K2PIG	58	-752.4400	1620.8800	0.0000	1.0000
SYM	60	-751.3003	1622.6007	0.0000	1.0000
K2PI	57	-754.3039	1622.6079	0.0000	1.0000
HKY	59	-754.3777	1626.7554	0.0000	1.0000
K2P	56	-758.9991	1629.9983	0.0000	1.0000
F81G	59	-767.8210	1653.6419	0.0000	1.0000
F81IG	60	-768.1160	1656.2321	0.0000	1.0000
JC69G	56	-772.2063	1656.4127	0.0000	1.0000
F81I	59	-769.5154	1657.0307	0.0000	1.0000
JC69IG	57	-772.1298	1658.2596	0.0000	1.0000
JC69I	56	-773.6482	1659.2964	0.0000	1.0000
F81	58	-773.1606	1662.3212	0.0000	1.0000
JC69	55	-778.1579	1666.3157	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-752.5481	1809.8814	0.4704	0.4704
GTRG	64	-734.4015	1811.0179	0.2665	0.7369
HKYG	60	-746.2473	1813.3210	0.0843	0.8212
K2PI	57	-754.3039	1813.3930	0.0813	0.9025
K2PIG	58	-752.4400	1815.0122	0.0362	0.9386
SYMG	61	-745.0246	1816.2228	0.0197	0.9584
GTRI	64	-737.2896	1816.7941	0.0148	0.9732
K2P	56	-758.9991	1817.4363	0.0108	0.9840
SYMI	61	-746.1524	1818.4784	0.0064	0.9904
GTRIG	65	-736.1636	1819.8892	0.0032	0.9935
GTR	63	-741.8749	1820.6175	0.0022	0.9957
HKYI	60	-750.3450	1821.5164	0.0014	0.9971
SYMIG	62	-745.0077	1821.5361	0.0014	0.9985
HKYIG	61	-748.5317	1823.2370	0.0006	0.9991
SYM	60	-751.3003	1823.4271	0.0005	0.9996
HKY	59	-754.3777	1824.2348	0.0004	1.0000
JC69G	56	-772.2063	1843.8507	0.0000	1.0000
JC69I	56	-773.6482	1846.7344	0.0000	1.0000
JC69IG	57	-772.1298	1849.0447	0.0000	1.0000
JC69	55	-778.1579	1850.4066	0.0000	1.0000
F81G	59	-767.8210	1851.1212	0.0000	1.0000
F81I	59	-769.5154	1854.5101	0.0000	1.0000
F81	58	-773.1606	1856.4535	0.0000	1.0000
F81IG	60	-768.1160	1857.0586	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/11139_hg18.phylip = 1 - 210;
 Partition Dummy = 1:261-aggressive-phylip/11139_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11139_hg18.phylip = 1 - 210;
 Partition Dummy = 1:261-aggressive-phylip/11139_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 (K2PG)]

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

