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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-2308.3843	4747.9784	0.5689	0.5689
HKYIG	61	-2308.0819	4749.7651	0.2328	0.8017
GTRG	64	-2305.3494	4751.5185	0.0969	0.8986
HKYI	60	-2310.7405	4752.6907	0.0539	0.9525
GTRIG	65	-2305.0396	4753.3199	0.0394	0.9919
GTRI	64	-2307.8307	4756.4812	0.0081	1.0000
K2PG	57	-2325.7946	4775.6685	0.0000	1.0000
K2PIG	58	-2324.7640	4775.9769	0.0000	1.0000
K2PI	57	-2327.7423	4779.5638	0.0000	1.0000
SYMIG	62	-2323.7571	4783.5141	0.0000	1.0000
SYMG	61	-2324.9565	4783.5142	0.0000	1.0000
SYMI	61	-2326.6516	4786.9043	0.0000	1.0000
HKY	59	-2337.9456	4804.7169	0.0000	1.0000
GTR	63	-2335.1754	4808.7570	0.0000	1.0000
K2P	56	-2355.1757	4832.0684	0.0000	1.0000
SYM	60	-2354.0769	4839.3635	0.0000	1.0000
F81G	59	-2381.5304	4891.8864	0.0000	1.0000
F81IG	60	-2381.2317	4893.6732	0.0000	1.0000
F81I	59	-2383.7075	4896.2406	0.0000	1.0000
JC69G	56	-2394.9573	4911.6315	0.0000	1.0000
JC69IG	57	-2394.6974	4913.4741	0.0000	1.0000
JC69I	56	-2397.7625	4917.2420	0.0000	1.0000
F81	58	-2410.4415	4947.3318	0.0000	1.0000
JC69	55	-2423.4375	4966.2366	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2308.3843	4736.7686	0.4553	0.4553
HKYIG	61	-2308.0819	4738.1639	0.2266	0.6819
GTRG	64	-2305.3494	4738.6988	0.1734	0.8554
GTRIG	65	-2305.0396	4740.0791	0.0870	0.9423
HKYI	60	-2310.7405	4741.4809	0.0432	0.9855
GTRI	64	-2307.8307	4743.6614	0.0145	1.0000
K2PIG	58	-2324.7640	4765.5280	0.0000	1.0000
K2PG	57	-2325.7946	4765.5893	0.0000	1.0000
K2PI	57	-2327.7423	4769.4845	0.0000	1.0000
SYMIG	62	-2323.7571	4771.5141	0.0000	1.0000
SYMG	61	-2324.9565	4771.9129	0.0000	1.0000
SYMI	61	-2326.6516	4775.3031	0.0000	1.0000
HKY	59	-2337.9456	4793.8913	0.0000	1.0000
GTR	63	-2335.1754	4796.3508	0.0000	1.0000
K2P	56	-2355.1757	4822.3515	0.0000	1.0000
SYM	60	-2354.0769	4828.1537	0.0000	1.0000
F81G	59	-2381.5304	4881.0607	0.0000	1.0000
F81IG	60	-2381.2317	4882.4634	0.0000	1.0000
F81I	59	-2383.7075	4885.4149	0.0000	1.0000
JC69G	56	-2394.9573	4901.9146	0.0000	1.0000
JC69IG	57	-2394.6974	4903.3948	0.0000	1.0000
JC69I	56	-2397.7625	4907.5251	0.0000	1.0000
F81	58	-2410.4415	4936.8830	0.0000	1.0000
JC69	55	-2423.4375	4956.8749	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2308.3843	5011.0216	0.8725	0.8725
HKYI	60	-2310.7405	5015.7339	0.0827	0.9552
HKYIG	61	-2308.0819	5016.9878	0.0442	0.9994
K2PG	57	-2325.7946	5026.1296	0.0005	0.9998
K2PI	57	-2327.7423	5030.0248	0.0001	0.9999
K2PIG	58	-2324.7640	5030.6392	0.0000	1.0000
GTRG	64	-2305.3494	5031.2353	0.0000	1.0000
GTRI	64	-2307.8307	5036.1979	0.0000	1.0000
GTRIG	65	-2305.0396	5037.1865	0.0000	1.0000
SYMG	61	-2324.9565	5050.7368	0.0000	1.0000
SYMI	61	-2326.6516	5054.1270	0.0000	1.0000
SYMIG	62	-2323.7571	5054.9089	0.0000	1.0000
HKY	59	-2337.9456	5063.5734	0.0000	1.0000
K2P	56	-2355.1757	5078.3209	0.0000	1.0000
GTR	63	-2335.1754	5084.3164	0.0000	1.0000
SYM	60	-2354.0769	5102.4067	0.0000	1.0000
F81G	59	-2381.5304	5150.7428	0.0000	1.0000
F81I	59	-2383.7075	5155.0970	0.0000	1.0000
F81IG	60	-2381.2317	5156.7164	0.0000	1.0000
JC69G	56	-2394.9573	5157.8841	0.0000	1.0000
JC69I	56	-2397.7625	5163.4945	0.0000	1.0000
JC69IG	57	-2394.6974	5163.9352	0.0000	1.0000
F81	58	-2410.4415	5201.9942	0.0000	1.0000
JC69	55	-2423.4375	5208.2735	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12726_hg18.phylip = 1 - 714;
 Partition Dummy = 1:261-aggressive-phylip/12726_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 AICc model (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/12726_hg18.phylip = 1 - 714;
 Partition Dummy = 1:261-aggressive-phylip/12726_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 (HKYG)]

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

