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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: K2PIG


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
GTRIG	65	-1945.2301	4038.1510	0.9668	0.9668
GTRG	64	-1950.0311	4045.1815	0.0288	0.9956
HKYIG	61	-1955.8335	4049.1352	0.0040	0.9996
GTRI	64	-1954.5720	4054.2633	0.0003	0.9999
HKYG	60	-1961.3437	4057.6261	0.0001	0.9999
K2PIG	58	-1964.5784	4059.0673	0.0000	1.0000
SYMIG	62	-1959.6555	4059.3191	0.0000	1.0000
HKYI	60	-1964.7855	4064.5098	0.0000	1.0000
SYMG	61	-1963.6650	4064.7982	0.0000	1.0000
K2PG	57	-1968.9248	4065.2614	0.0000	1.0000
SYMI	61	-1966.9914	4071.4510	0.0000	1.0000
K2PI	57	-1972.3127	4072.0371	0.0000	1.0000
F81IG	60	-2001.6066	4138.1521	0.0000	1.0000
F81G	59	-2007.1798	4146.7791	0.0000	1.0000
JC69IG	57	-2011.1313	4149.6744	0.0000	1.0000
F81I	59	-2008.7834	4149.9863	0.0000	1.0000
JC69G	56	-2015.6120	4156.1471	0.0000	1.0000
JC69I	56	-2019.0632	4163.0494	0.0000	1.0000
GTR	63	-2030.7104	4203.9792	0.0000	1.0000
HKY	59	-2042.3353	4217.0902	0.0000	1.0000
SYM	60	-2044.1364	4223.2116	0.0000	1.0000
K2P	56	-2050.2926	4225.5082	0.0000	1.0000
F81	58	-2085.2438	4300.3982	0.0000	1.0000
JC69	55	-2096.7355	4315.9154	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1945.2301	4020.4603	0.9766	0.9766
GTRG	64	-1950.0311	4028.0621	0.0218	0.9984
HKYIG	61	-1955.8335	4033.6669	0.0013	0.9997
GTRI	64	-1954.5720	4037.1439	0.0002	1.0000
HKYG	60	-1961.3437	4042.6873	0.0000	1.0000
SYMIG	62	-1959.6555	4043.3109	0.0000	1.0000
K2PIG	58	-1964.5784	4045.1567	0.0000	1.0000
SYMG	61	-1963.6650	4049.3299	0.0000	1.0000
HKYI	60	-1964.7855	4049.5711	0.0000	1.0000
K2PG	57	-1968.9248	4051.8496	0.0000	1.0000
SYMI	61	-1966.9914	4055.9827	0.0000	1.0000
K2PI	57	-1972.3127	4058.6253	0.0000	1.0000
F81IG	60	-2001.6066	4123.2133	0.0000	1.0000
F81G	59	-2007.1798	4132.3595	0.0000	1.0000
F81I	59	-2008.7834	4135.5668	0.0000	1.0000
JC69IG	57	-2011.1313	4136.2626	0.0000	1.0000
JC69G	56	-2015.6120	4143.2240	0.0000	1.0000
JC69I	56	-2019.0632	4150.1263	0.0000	1.0000
GTR	63	-2030.7104	4187.4207	0.0000	1.0000
HKY	59	-2042.3353	4202.6706	0.0000	1.0000
SYM	60	-2044.1364	4208.2728	0.0000	1.0000
K2P	56	-2050.2926	4212.5851	0.0000	1.0000
F81	58	-2085.2438	4286.4876	0.0000	1.0000
JC69	55	-2096.7355	4303.4709	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-1964.5784	4295.2373	0.5195	0.5195
HKYIG	61	-1955.8335	4296.6827	0.2522	0.7717
K2PG	57	-1968.9248	4297.6185	0.1580	0.9296
GTRIG	65	-1945.2301	4300.7230	0.0334	0.9631
HKYG	60	-1961.3437	4301.3914	0.0239	0.9870
GTRG	64	-1950.0311	4304.0132	0.0065	0.9935
K2PI	57	-1972.3127	4304.3942	0.0053	0.9988
HKYI	60	-1964.7855	4308.2751	0.0008	0.9996
SYMIG	62	-1959.6555	4310.6385	0.0002	0.9998
SYMG	61	-1963.6650	4312.3457	0.0001	0.9999
GTRI	64	-1954.5720	4313.0949	0.0001	1.0000
SYMI	61	-1966.9914	4318.9985	0.0000	1.0000
F81IG	60	-2001.6066	4381.9174	0.0000	1.0000
JC69IG	57	-2011.1313	4382.0315	0.0000	1.0000
JC69G	56	-2015.6120	4384.6812	0.0000	1.0000
F81G	59	-2007.1798	4386.7519	0.0000	1.0000
F81I	59	-2008.7834	4389.9591	0.0000	1.0000
JC69I	56	-2019.0632	4391.5835	0.0000	1.0000
K2P	56	-2050.2926	4454.0423	0.0000	1.0000
HKY	59	-2042.3353	4457.0630	0.0000	1.0000
GTR	63	-2030.7104	4459.0600	0.0000	1.0000
SYM	60	-2044.1364	4466.9769	0.0000	1.0000
F81	58	-2085.2438	4536.5682	0.0000	1.0000
JC69	55	-2096.7355	4540.6163	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/11388_hg18.phylip = 1 - 551;
 Partition Dummy = 1:261-aggressive-phylip/11388_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 (GTRIG)]

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11388_hg18.phylip = 1 - 551;
 Partition Dummy = 1:261-aggressive-phylip/11388_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

