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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	64	-979.3182	2120.5955	0.6257	0.6257
GTRG	64	-980.2748	2122.5088	0.2404	0.8661
GTRIG	65	-979.3307	2123.8253	0.1245	0.9906
HKYI	60	-990.2142	2129.8261	0.0062	0.9968
HKYG	60	-991.3746	2132.1467	0.0019	0.9987
HKYIG	61	-990.2440	2132.9880	0.0013	1.0000
GTR	63	-992.5940	2143.9684	0.0000	1.0000
HKY	59	-1000.5646	2147.4492	0.0000	1.0000
K2PI	57	-1007.6862	2155.6106	0.0000	1.0000
K2PG	57	-1008.7462	2157.7305	0.0000	1.0000
K2PIG	58	-1007.7155	2158.6980	0.0000	1.0000
SYMI	61	-1004.0308	2160.5616	0.0000	1.0000
SYMG	61	-1005.1182	2162.7363	0.0000	1.0000
SYMIG	62	-1004.0601	2163.7477	0.0000	1.0000
F81I	59	-1012.0522	2170.4245	0.0000	1.0000
K2P	56	-1017.6702	2172.5737	0.0000	1.0000
F81G	59	-1013.1682	2172.6565	0.0000	1.0000
F81IG	60	-1012.0820	2173.5617	0.0000	1.0000
SYM	60	-1014.6436	2178.6848	0.0000	1.0000
F81	58	-1021.7984	2186.8637	0.0000	1.0000
JC69I	56	-1025.5102	2188.2535	0.0000	1.0000
JC69G	56	-1026.5764	2190.3860	0.0000	1.0000
JC69IG	57	-1025.5392	2191.3165	0.0000	1.0000
JC69	55	-1035.2366	2204.7251	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	64	-979.3182	2086.6363	0.5717	0.5717
GTRG	64	-980.2748	2088.5496	0.2196	0.7913
GTRIG	65	-979.3307	2088.6614	0.2077	0.9990
HKYI	60	-990.2142	2100.4285	0.0006	0.9996
HKYIG	61	-990.2440	2102.4880	0.0002	0.9998
HKYG	60	-991.3746	2102.7491	0.0002	1.0000
GTR	63	-992.5940	2111.1879	0.0000	1.0000
HKY	59	-1000.5646	2119.1292	0.0000	1.0000
K2PI	57	-1007.6862	2129.3725	0.0000	1.0000
SYMI	61	-1004.0308	2130.0616	0.0000	1.0000
K2PIG	58	-1007.7155	2131.4311	0.0000	1.0000
K2PG	57	-1008.7462	2131.4924	0.0000	1.0000
SYMIG	62	-1004.0601	2132.1202	0.0000	1.0000
SYMG	61	-1005.1182	2132.2363	0.0000	1.0000
F81I	59	-1012.0522	2142.1045	0.0000	1.0000
F81IG	60	-1012.0820	2144.1641	0.0000	1.0000
F81G	59	-1013.1682	2144.3365	0.0000	1.0000
K2P	56	-1017.6702	2147.3405	0.0000	1.0000
SYM	60	-1014.6436	2149.2872	0.0000	1.0000
F81	58	-1021.7984	2159.5968	0.0000	1.0000
JC69I	56	-1025.5102	2163.0203	0.0000	1.0000
JC69IG	57	-1025.5392	2165.0784	0.0000	1.0000
JC69G	56	-1026.5764	2165.1528	0.0000	1.0000
JC69	55	-1035.2366	2180.4731	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-990.2142	2324.6228	0.4591	0.4591
GTRI	64	-979.3182	2325.7770	0.2578	0.7169
HKYG	60	-991.3746	2326.9435	0.1439	0.8608
GTRG	64	-980.2748	2327.6902	0.0991	0.9599
HKYIG	61	-990.2440	2330.4189	0.0253	0.9852
GTRIG	65	-979.3307	2331.5386	0.0145	0.9996
HKY	59	-1000.5646	2339.5869	0.0003	0.9999
K2PI	57	-1007.6862	2342.3571	0.0001	1.0000
K2PG	57	-1008.7462	2344.4770	0.0000	1.0000
GTR	63	-992.5940	2346.5920	0.0000	1.0000
K2PIG	58	-1007.7155	2348.1523	0.0000	1.0000
K2P	56	-1017.6702	2356.5885	0.0000	1.0000
SYMI	61	-1004.0308	2357.9926	0.0000	1.0000
SYMG	61	-1005.1182	2360.1672	0.0000	1.0000
F81I	59	-1012.0522	2362.5623	0.0000	1.0000
SYMIG	62	-1004.0601	2363.7876	0.0000	1.0000
F81G	59	-1013.1682	2364.7942	0.0000	1.0000
F81IG	60	-1012.0820	2368.3584	0.0000	1.0000
JC69I	56	-1025.5102	2372.2683	0.0000	1.0000
SYM	60	-1014.6436	2373.4816	0.0000	1.0000
JC69G	56	-1026.5764	2374.4009	0.0000	1.0000
F81	58	-1021.7984	2376.3180	0.0000	1.0000
JC69IG	57	-1025.5392	2378.0630	0.0000	1.0000
JC69	55	-1035.2366	2385.9846	0.0000	1.0000



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

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

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


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

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


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

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

