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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-2346.5656	4836.3692	0.6512	0.6512
GTRIG	65	-2345.9524	4837.6479	0.3436	0.9949
GTRI	64	-2351.8292	4846.8966	0.0034	0.9982
HKYG	60	-2357.8923	4849.0936	0.0011	0.9994
HKYIG	61	-2357.2523	4850.2823	0.0006	1.0000
HKYI	60	-2362.6918	4858.6926	0.0000	1.0000
SYMG	61	-2391.0051	4917.7881	0.0000	1.0000
SYMIG	62	-2390.3189	4918.8933	0.0000	1.0000
F81G	59	-2395.6062	4922.0617	0.0000	1.0000
F81IG	60	-2395.0772	4923.4635	0.0000	1.0000
K2PG	57	-2402.0145	4929.9857	0.0000	1.0000
SYMI	61	-2397.3678	4930.5134	0.0000	1.0000
F81I	59	-2400.0262	4930.9018	0.0000	1.0000
K2PIG	58	-2401.2553	4930.9093	0.0000	1.0000
GTR	63	-2397.5042	4935.7507	0.0000	1.0000
K2PI	57	-2407.1581	4940.2728	0.0000	1.0000
HKY	59	-2410.6926	4952.2346	0.0000	1.0000
JC69G	56	-2431.4112	4986.3458	0.0000	1.0000
JC69IG	57	-2430.8266	4987.6099	0.0000	1.0000
JC69I	56	-2436.1356	4995.7947	0.0000	1.0000
SYM	60	-2439.8938	5013.0967	0.0000	1.0000
F81	58	-2443.7089	5015.8163	0.0000	1.0000
K2P	56	-2453.5155	5030.5544	0.0000	1.0000
JC69	55	-2479.6994	5080.4979	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2346.5656	4821.1311	0.5933	0.5933
GTRIG	65	-2345.9524	4821.9047	0.4030	0.9963
GTRI	64	-2351.8292	4831.6585	0.0031	0.9993
HKYG	60	-2357.8923	4835.7845	0.0004	0.9997
HKYIG	61	-2357.2523	4836.5046	0.0003	1.0000
HKYI	60	-2362.6918	4845.3835	0.0000	1.0000
SYMG	61	-2391.0051	4904.0103	0.0000	1.0000
SYMIG	62	-2390.3189	4904.6379	0.0000	1.0000
F81G	59	-2395.6062	4909.2123	0.0000	1.0000
F81IG	60	-2395.0772	4910.1544	0.0000	1.0000
SYMI	61	-2397.3678	4916.7356	0.0000	1.0000
K2PG	57	-2402.0145	4918.0291	0.0000	1.0000
F81I	59	-2400.0262	4918.0525	0.0000	1.0000
K2PIG	58	-2401.2553	4918.5107	0.0000	1.0000
GTR	63	-2397.5042	4921.0085	0.0000	1.0000
K2PI	57	-2407.1581	4928.3162	0.0000	1.0000
HKY	59	-2410.6926	4939.3852	0.0000	1.0000
JC69G	56	-2431.4112	4974.8224	0.0000	1.0000
JC69IG	57	-2430.8266	4975.6533	0.0000	1.0000
JC69I	56	-2436.1356	4984.2712	0.0000	1.0000
SYM	60	-2439.8938	4999.7876	0.0000	1.0000
F81	58	-2443.7089	5003.4178	0.0000	1.0000
K2P	56	-2453.5155	5019.0309	0.0000	1.0000
JC69	55	-2479.6994	5069.3988	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2357.8923	5100.6903	0.7547	0.7547
GTRG	64	-2346.5656	5103.6973	0.1678	0.9225
HKYIG	61	-2357.2523	5105.8255	0.0579	0.9804
GTRIG	65	-2345.9524	5108.8860	0.0125	0.9929
HKYI	60	-2362.6918	5110.2894	0.0062	0.9991
GTRI	64	-2351.8292	5114.2247	0.0009	1.0000
K2PG	57	-2402.0145	5169.6896	0.0000	1.0000
F81G	59	-2395.6062	5169.7031	0.0000	1.0000
SYMG	61	-2391.0051	5173.3312	0.0000	1.0000
K2PIG	58	-2401.2553	5174.5863	0.0000	1.0000
F81IG	60	-2395.0772	5175.0602	0.0000	1.0000
SYMIG	62	-2390.3189	5178.3739	0.0000	1.0000
F81I	59	-2400.0262	5178.5432	0.0000	1.0000
K2PI	57	-2407.1581	5179.9767	0.0000	1.0000
SYMI	61	-2397.3678	5186.0565	0.0000	1.0000
GTR	63	-2397.5042	5199.1596	0.0000	1.0000
HKY	59	-2410.6926	5199.8760	0.0000	1.0000
JC69G	56	-2431.4112	5222.0678	0.0000	1.0000
JC69IG	57	-2430.8266	5227.3138	0.0000	1.0000
JC69I	56	-2436.1356	5231.5167	0.0000	1.0000
F81	58	-2443.7089	5259.4934	0.0000	1.0000
SYM	60	-2439.8938	5264.6935	0.0000	1.0000
K2P	56	-2453.5155	5266.2764	0.0000	1.0000
JC69	55	-2479.6994	5312.2291	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/8993_hg18.phylip = 1 - 611;
 Partition Dummy = 1:261-aggressive-phylip/8993_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/8993_hg18.phylip = 1 - 611;
 Partition Dummy = 1:261-aggressive-phylip/8993_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 (HKYG)]

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


-------------------------------------------------------------
End of Output

