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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYG
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
HKYG	60	-947.4079	2039.6294	0.3891	0.3891
HKYI	60	-947.6137	2040.0410	0.3167	0.7058
HKYIG	61	-947.3074	2042.3428	0.1002	0.8060
GTRG	64	-943.0345	2042.6602	0.0855	0.8915
GTRI	64	-943.3981	2043.3874	0.0594	0.9510
GTRIG	65	-942.8814	2045.3491	0.0223	0.9732
K2PG	57	-955.2101	2046.6081	0.0119	0.9851
K2PI	57	-955.3940	2046.9758	0.0099	0.9950
K2PIG	58	-955.0792	2049.2022	0.0032	0.9983
SYMG	61	-952.0603	2051.8484	0.0009	0.9991
SYMI	61	-952.4115	2052.5509	0.0006	0.9997
SYMIG	62	-951.8652	2054.3925	0.0002	1.0000
HKY	59	-958.7151	2059.3491	0.0000	1.0000
GTR	63	-953.4783	2060.5730	0.0000	1.0000
K2P	56	-967.2695	2067.8901	0.0000	1.0000
SYM	60	-962.9623	2070.7381	0.0000	1.0000
F81G	59	-971.2142	2084.3472	0.0000	1.0000
F81I	59	-971.3947	2084.7083	0.0000	1.0000
F81IG	60	-971.1052	2087.0240	0.0000	1.0000
JC69G	56	-978.2956	2089.9423	0.0000	1.0000
JC69I	56	-978.4768	2090.3048	0.0000	1.0000
JC69IG	57	-978.1658	2092.5195	0.0000	1.0000
F81	58	-983.1491	2105.3420	0.0000	1.0000
JC69	55	-990.4953	2111.5239	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-943.0345	2014.0691	0.2724	0.2724
GTRI	64	-943.3981	2014.7963	0.1894	0.4618
HKYG	60	-947.4079	2014.8158	0.1875	0.6493
HKYI	60	-947.6137	2015.2274	0.1527	0.8020
GTRIG	65	-942.8814	2015.7629	0.1168	0.9188
HKYIG	61	-947.3074	2016.6149	0.0763	0.9951
K2PG	57	-955.2101	2024.4201	0.0015	0.9966
K2PI	57	-955.3940	2024.7879	0.0013	0.9979
SYMG	61	-952.0603	2026.1205	0.0007	0.9986
K2PIG	58	-955.0792	2026.1584	0.0006	0.9992
SYMI	61	-952.4115	2026.8230	0.0005	0.9997
SYMIG	62	-951.8652	2027.7304	0.0003	1.0000
GTR	63	-953.4783	2032.9565	0.0000	1.0000
HKY	59	-958.7151	2035.4302	0.0000	1.0000
SYM	60	-962.9623	2045.9245	0.0000	1.0000
K2P	56	-967.2695	2046.5390	0.0000	1.0000
F81G	59	-971.2142	2060.4283	0.0000	1.0000
F81I	59	-971.3947	2060.7894	0.0000	1.0000
F81IG	60	-971.1052	2062.2105	0.0000	1.0000
JC69G	56	-978.2956	2068.5912	0.0000	1.0000
JC69I	56	-978.4768	2068.9536	0.0000	1.0000
JC69IG	57	-978.1658	2070.3316	0.0000	1.0000
F81	58	-983.1491	2082.2982	0.0000	1.0000
JC69	55	-990.4953	2090.9906	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-955.2101	2245.2912	0.3883	0.3883
K2PI	57	-955.3940	2245.6590	0.3231	0.7114
HKYG	60	-947.4079	2247.3117	0.1414	0.8528
HKYI	60	-947.6137	2247.7233	0.1151	0.9679
K2PIG	58	-955.0792	2250.9044	0.0235	0.9914
HKYIG	61	-947.3074	2252.9856	0.0083	0.9996
GTRG	64	-943.0345	2262.0647	0.0001	0.9997
SYMG	61	-952.0603	2262.4913	0.0001	0.9998
GTRI	64	-943.3981	2262.7919	0.0001	0.9999
SYMI	61	-952.4115	2263.1938	0.0001	0.9999
K2P	56	-967.2695	2263.5351	0.0000	1.0000
HKY	59	-958.7151	2264.0511	0.0000	1.0000
GTRIG	65	-942.8814	2267.6334	0.0000	1.0000
SYMIG	62	-951.8652	2267.9761	0.0000	1.0000
GTR	63	-953.4783	2277.0772	0.0000	1.0000
SYM	60	-962.9623	2278.4204	0.0000	1.0000
JC69G	56	-978.2956	2285.5873	0.0000	1.0000
JC69I	56	-978.4768	2285.9497	0.0000	1.0000
F81G	59	-971.2142	2289.0492	0.0000	1.0000
F81I	59	-971.3947	2289.4103	0.0000	1.0000
JC69IG	57	-978.1658	2291.2027	0.0000	1.0000
F81IG	60	-971.1052	2294.7063	0.0000	1.0000
JC69	55	-990.4953	2304.1118	0.0000	1.0000
F81	58	-983.1491	2307.0442	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/3416_hg18.phylip = 1 - 356;
 Partition Dummy = 1:261-aggressive-phylip/3416_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 (HKYG)]

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

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

