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

Input data from file "261-aggressive-phylip/6224_hg18.phylip" (29 taxa, 1211 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	-3341.1451	6817.5502	0.3421	0.3421
GTRIG	65	-3340.3409	6818.1752	0.2503	0.5924
GTRI	64	-3341.7843	6818.8286	0.1805	0.7730
HKYG	60	-3346.8995	6820.1642	0.0926	0.8656
HKYI	60	-3347.2093	6820.7838	0.0679	0.9335
HKYIG	61	-3346.1212	6820.8255	0.0665	1.0000
K2PG	57	-3375.1488	6870.0322	0.0000	1.0000
K2PI	57	-3375.1813	6870.0972	0.0000	1.0000
K2PIG	58	-3374.2974	6870.5357	0.0000	1.0000
SYMG	61	-3372.5988	6873.7806	0.0000	1.0000
SYMI	61	-3373.1028	6874.7887	0.0000	1.0000
SYMIG	62	-3372.2247	6875.2543	0.0000	1.0000
GTR	63	-3398.9626	6930.9558	0.0000	1.0000
HKY	59	-3403.4329	6931.0171	0.0000	1.0000
K2P	56	-3427.8108	6973.1537	0.0000	1.0000
SYM	60	-3424.9413	6976.2479	0.0000	1.0000
F81G	59	-3431.1112	6986.3737	0.0000	1.0000
F81I	59	-3431.3310	6986.8132	0.0000	1.0000
F81IG	60	-3431.0808	6988.5268	0.0000	1.0000
JC69IG	57	-3450.5846	7020.9038	0.0000	1.0000
JC69G	56	-3451.8071	7021.1462	0.0000	1.0000
JC69I	56	-3451.8966	7021.3253	0.0000	1.0000
F81	58	-3481.6931	7085.3272	0.0000	1.0000
JC69	55	-3501.5187	7118.3708	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3341.1451	6810.2901	0.3586	0.3586
GTRIG	65	-3340.3409	6810.6817	0.2949	0.6535
GTRI	64	-3341.7843	6811.5685	0.1893	0.8427
HKYG	60	-3346.8995	6813.7990	0.0620	0.9048
HKYIG	61	-3346.1212	6814.2424	0.0497	0.9545
HKYI	60	-3347.2093	6814.4185	0.0455	1.0000
K2PG	57	-3375.1488	6864.2976	0.0000	1.0000
K2PI	57	-3375.1813	6864.3626	0.0000	1.0000
K2PIG	58	-3374.2974	6864.5948	0.0000	1.0000
SYMG	61	-3372.5988	6867.1975	0.0000	1.0000
SYMI	61	-3373.1028	6868.2055	0.0000	1.0000
SYMIG	62	-3372.2247	6868.4495	0.0000	1.0000
GTR	63	-3398.9626	6923.9253	0.0000	1.0000
HKY	59	-3403.4329	6924.8659	0.0000	1.0000
K2P	56	-3427.8108	6967.6216	0.0000	1.0000
SYM	60	-3424.9413	6969.8827	0.0000	1.0000
F81G	59	-3431.1112	6980.2225	0.0000	1.0000
F81I	59	-3431.3310	6980.6621	0.0000	1.0000
F81IG	60	-3431.0808	6982.1616	0.0000	1.0000
JC69IG	57	-3450.5846	7015.1692	0.0000	1.0000
JC69G	56	-3451.8071	7015.6142	0.0000	1.0000
JC69I	56	-3451.8966	7015.7932	0.0000	1.0000
F81	58	-3481.6931	7079.3862	0.0000	1.0000
JC69	55	-3501.5187	7113.0374	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3346.8995	7119.7511	0.5566	0.5566
HKYI	60	-3347.2093	7120.3706	0.4083	0.9650
HKYIG	61	-3346.1212	7125.2937	0.0348	0.9998
GTRG	64	-3341.1451	7136.6390	0.0001	0.9999
GTRI	64	-3341.7843	7137.9174	0.0001	1.0000
GTRIG	65	-3340.3409	7142.1298	0.0000	1.0000
K2PG	57	-3375.1488	7154.9521	0.0000	1.0000
K2PI	57	-3375.1813	7155.0171	0.0000	1.0000
K2PIG	58	-3374.2974	7160.3485	0.0000	1.0000
SYMG	61	-3372.5988	7178.2488	0.0000	1.0000
SYMI	61	-3373.1028	7179.2568	0.0000	1.0000
SYMIG	62	-3372.2247	7184.6000	0.0000	1.0000
HKY	59	-3403.4329	7225.7188	0.0000	1.0000
GTR	63	-3398.9626	7245.1750	0.0000	1.0000
K2P	56	-3427.8108	7253.1769	0.0000	1.0000
SYM	60	-3424.9413	7275.8348	0.0000	1.0000
F81G	59	-3431.1112	7281.0754	0.0000	1.0000
F81I	59	-3431.3310	7281.5150	0.0000	1.0000
F81IG	60	-3431.0808	7288.1137	0.0000	1.0000
JC69G	56	-3451.8071	7301.1695	0.0000	1.0000
JC69I	56	-3451.8966	7301.3485	0.0000	1.0000
JC69IG	57	-3450.5846	7305.8237	0.0000	1.0000
F81	58	-3481.6931	7375.1399	0.0000	1.0000
JC69	55	-3501.5187	7393.4935	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/6224_hg18.phylip = 1 - 1211;
 Partition Dummy = 1:261-aggressive-phylip/6224_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/6224_hg18.phylip = 1 - 1211;
 Partition Dummy = 1:261-aggressive-phylip/6224_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/6224_hg18.phylip = 1 - 1211;
 Partition Dummy = 1:261-aggressive-phylip/6224_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

