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

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

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


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	-2739.2831	5619.5055	0.7033	0.7033
GTRIG	65	-2738.9440	5621.2525	0.2936	0.9969
GTRI	64	-2745.1519	5631.2432	0.0020	0.9989
SYMG	61	-2749.6796	5633.0682	0.0008	0.9997
SYMIG	62	-2749.3413	5634.7942	0.0003	1.0000
SYMI	61	-2755.2223	5644.1537	0.0000	1.0000
HKYG	60	-2759.4911	5650.2960	0.0000	1.0000
HKYIG	61	-2758.9485	5651.6059	0.0000	1.0000
K2PG	57	-2765.0388	5654.2499	0.0000	1.0000
K2PIG	58	-2764.1192	5654.7838	0.0000	1.0000
HKYI	60	-2765.3323	5661.9783	0.0000	1.0000
K2PI	57	-2770.3893	5664.9508	0.0000	1.0000
GTR	63	-2791.4064	5721.3345	0.0000	1.0000
SYM	60	-2802.7578	5736.8293	0.0000	1.0000
HKY	59	-2816.7319	5762.3896	0.0000	1.0000
K2P	56	-2822.4129	5766.6322	0.0000	1.0000
F81G	59	-2835.6613	5800.2485	0.0000	1.0000
JC69G	56	-2839.6652	5801.1368	0.0000	1.0000
F81IG	60	-2835.3564	5802.0266	0.0000	1.0000
JC69IG	57	-2839.3568	5802.8859	0.0000	1.0000
F81I	59	-2840.2427	5809.4114	0.0000	1.0000
JC69I	56	-2844.2255	5810.2575	0.0000	1.0000
F81	58	-2890.1387	5906.8229	0.0000	1.0000
JC69	55	-2895.1158	5909.6795	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2739.2831	5606.5661	0.6578	0.6578
GTRIG	65	-2738.9440	5607.8881	0.3397	0.9975
GTRI	64	-2745.1519	5618.3039	0.0019	0.9994
SYMG	61	-2749.6796	5621.3592	0.0004	0.9998
SYMIG	62	-2749.3413	5622.6826	0.0002	1.0000
SYMI	61	-2755.2223	5632.4447	0.0000	1.0000
HKYG	60	-2759.4911	5638.9823	0.0000	1.0000
HKYIG	61	-2758.9485	5639.8970	0.0000	1.0000
K2PG	57	-2765.0388	5644.0776	0.0000	1.0000
K2PIG	58	-2764.1192	5644.2383	0.0000	1.0000
HKYI	60	-2765.3323	5650.6645	0.0000	1.0000
K2PI	57	-2770.3893	5654.7785	0.0000	1.0000
GTR	63	-2791.4064	5708.8127	0.0000	1.0000
SYM	60	-2802.7578	5725.5156	0.0000	1.0000
HKY	59	-2816.7319	5751.4637	0.0000	1.0000
K2P	56	-2822.4129	5756.8257	0.0000	1.0000
F81G	59	-2835.6613	5789.3225	0.0000	1.0000
F81IG	60	-2835.3564	5790.7128	0.0000	1.0000
JC69G	56	-2839.6652	5791.3303	0.0000	1.0000
JC69IG	57	-2839.3568	5792.7136	0.0000	1.0000
F81I	59	-2840.2427	5798.4855	0.0000	1.0000
JC69I	56	-2844.2255	5800.4510	0.0000	1.0000
F81	58	-2890.1387	5896.2774	0.0000	1.0000
JC69	55	-2895.1158	5900.2316	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2739.2831	5898.5626	0.5774	0.5774
SYMG	61	-2749.6796	5899.6683	0.3321	0.9095
K2PG	57	-2765.0388	5904.1369	0.0356	0.9451
GTRIG	65	-2738.9440	5904.4469	0.0305	0.9755
SYMIG	62	-2749.3413	5905.5542	0.0175	0.9930
K2PIG	58	-2764.1192	5908.8601	0.0034	0.9964
GTRI	64	-2745.1519	5910.3003	0.0016	0.9980
SYMI	61	-2755.2223	5910.7538	0.0013	0.9993
HKYG	60	-2759.4911	5912.7289	0.0005	0.9998
K2PI	57	-2770.3893	5914.8379	0.0002	1.0000
HKYIG	61	-2758.9485	5918.2060	0.0000	1.0000
HKYI	60	-2765.3323	5924.4112	0.0000	1.0000
GTR	63	-2791.4064	5996.2467	0.0000	1.0000
SYM	60	-2802.7578	5999.2622	0.0000	1.0000
K2P	56	-2822.4129	6012.3226	0.0000	1.0000
HKY	59	-2816.7319	6020.6479	0.0000	1.0000
JC69G	56	-2839.6652	6046.8272	0.0000	1.0000
JC69IG	57	-2839.3568	6052.7729	0.0000	1.0000
JC69I	56	-2844.2255	6055.9479	0.0000	1.0000
F81G	59	-2835.6613	6058.5067	0.0000	1.0000
F81IG	60	-2835.3564	6064.4594	0.0000	1.0000
F81I	59	-2840.2427	6067.6697	0.0000	1.0000
JC69	55	-2895.1158	6151.1661	0.0000	1.0000
F81	58	-2890.1387	6160.8992	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/7668_hg18.phylip = 1 - 708;
 Partition Dummy = 1:261-aggressive-phylip/7668_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/7668_hg18.phylip = 1 - 708;
 Partition Dummy = 1:261-aggressive-phylip/7668_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 (GTRG)]

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


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

