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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
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
HKYIG	61	-1963.8204	4066.5625	0.3260	0.3260
GTRIG	65	-1958.8903	4067.1485	0.2432	0.5693
GTRG	64	-1960.2550	4067.2487	0.2314	0.8006
HKYG	60	-1965.6252	4067.5897	0.1951	0.9957
HKYI	60	-1969.9217	4076.1827	0.0027	0.9984
GTRI	64	-1965.2211	4077.1808	0.0016	1.0000
K2PIG	58	-1980.4058	4092.0205	0.0000	1.0000
SYMIG	62	-1975.5395	4092.5947	0.0000	1.0000
K2PG	57	-1982.0756	4092.8119	0.0000	1.0000
SYMG	61	-1977.0602	4093.0420	0.0000	1.0000
K2PI	57	-1985.9209	4100.5026	0.0000	1.0000
SYMI	61	-1981.1904	4101.3025	0.0000	1.0000
F81IG	60	-2016.9163	4170.1718	0.0000	1.0000
F81G	59	-2018.4576	4170.6835	0.0000	1.0000
F81I	59	-2022.3706	4178.5095	0.0000	1.0000
JC69IG	57	-2027.3423	4183.3454	0.0000	1.0000
JC69G	56	-2028.6634	4183.4506	0.0000	1.0000
JC69I	56	-2035.0002	4196.1243	0.0000	1.0000
GTR	63	-2029.1630	4202.4473	0.0000	1.0000
HKY	59	-2036.2189	4206.2061	0.0000	1.0000
SYM	60	-2042.9973	4222.3339	0.0000	1.0000
K2P	56	-2050.3395	4226.8028	0.0000	1.0000
F81	58	-2083.5047	4298.2183	0.0000	1.0000
JC69	55	-2094.8901	4313.3785	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1958.8903	4047.7806	0.4400	0.4400
GTRG	64	-1960.2550	4048.5100	0.3056	0.7456
HKYIG	61	-1963.8204	4049.6408	0.1736	0.9192
HKYG	60	-1965.6252	4051.2504	0.0776	0.9968
GTRI	64	-1965.2211	4058.4421	0.0021	0.9989
HKYI	60	-1969.9217	4059.8434	0.0011	1.0000
SYMIG	62	-1975.5395	4075.0790	0.0000	1.0000
SYMG	61	-1977.0602	4076.1203	0.0000	1.0000
K2PIG	58	-1980.4058	4076.8116	0.0000	1.0000
K2PG	57	-1982.0756	4078.1512	0.0000	1.0000
SYMI	61	-1981.1904	4084.3808	0.0000	1.0000
K2PI	57	-1985.9209	4085.8419	0.0000	1.0000
F81IG	60	-2016.9163	4153.8326	0.0000	1.0000
F81G	59	-2018.4576	4154.9151	0.0000	1.0000
F81I	59	-2022.3706	4162.7412	0.0000	1.0000
JC69IG	57	-2027.3423	4168.6846	0.0000	1.0000
JC69G	56	-2028.6634	4169.3267	0.0000	1.0000
JC69I	56	-2035.0002	4182.0004	0.0000	1.0000
GTR	63	-2029.1630	4184.3260	0.0000	1.0000
HKY	59	-2036.2189	4190.4377	0.0000	1.0000
SYM	60	-2042.9973	4205.9946	0.0000	1.0000
K2P	56	-2050.3395	4212.6789	0.0000	1.0000
F81	58	-2083.5047	4283.0094	0.0000	1.0000
JC69	55	-2094.8901	4299.7803	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1965.6252	4305.1973	0.7782	0.7782
HKYIG	61	-1963.8204	4307.8201	0.2097	0.9878
HKYI	60	-1969.9217	4313.7903	0.0106	0.9984
GTRG	64	-1960.2550	4319.3867	0.0006	0.9991
K2PG	57	-1982.0756	4319.4007	0.0006	0.9997
K2PIG	58	-1980.4058	4322.2936	0.0002	0.9999
GTRIG	65	-1958.8903	4322.8897	0.0001	1.0000
K2PI	57	-1985.9209	4327.0914	0.0000	1.0000
GTRI	64	-1965.2211	4329.3188	0.0000	1.0000
SYMG	61	-1977.0602	4334.2996	0.0000	1.0000
SYMIG	62	-1975.5395	4337.4908	0.0000	1.0000
SYMI	61	-1981.1904	4342.5602	0.0000	1.0000
F81G	59	-2018.4576	4404.6295	0.0000	1.0000
JC69G	56	-2028.6634	4406.3438	0.0000	1.0000
F81IG	60	-2016.9163	4407.7794	0.0000	1.0000
JC69IG	57	-2027.3423	4409.9342	0.0000	1.0000
F81I	59	-2022.3706	4412.4556	0.0000	1.0000
JC69I	56	-2035.0002	4419.0174	0.0000	1.0000
HKY	59	-2036.2189	4440.1522	0.0000	1.0000
K2P	56	-2050.3395	4449.6960	0.0000	1.0000
GTR	63	-2029.1630	4450.9702	0.0000	1.0000
SYM	60	-2042.9973	4459.9415	0.0000	1.0000
F81	58	-2083.5047	4528.4914	0.0000	1.0000
JC69	55	-2094.8901	4532.5649	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/8682_hg18.phylip = 1 - 509;
 Partition Dummy = 1:261-aggressive-phylip/8682_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 (HKYIG)]

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


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

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

