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

Input data from file "261-aggressive-phylip/10230_hg18.phylip" (29 taxa, 623 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	-1547.0798	3227.1846	0.7705	0.7705
HKYIG	61	-1547.0821	3229.6472	0.2249	0.9955
GTRG	64	-1547.5268	3237.9640	0.0035	0.9990
GTRIG	65	-1547.5292	3240.4623	0.0010	1.0000
HKY	59	-1562.7130	3256.0016	0.0000	1.0000
HKYI	60	-1562.7127	3258.4502	0.0000	1.0000
GTR	63	-1564.0878	3268.6014	0.0000	1.0000
GTRI	64	-1564.0874	3271.0853	0.0000	1.0000
K2PG	57	-1577.1885	3280.0796	0.0000	1.0000
K2PIG	58	-1577.1909	3282.5165	0.0000	1.0000
SYMG	61	-1577.0489	3289.5808	0.0000	1.0000
SYMIG	62	-1577.0513	3292.0525	0.0000	1.0000
F81G	59	-1586.4270	3303.4295	0.0000	1.0000
F81IG	60	-1586.4294	3305.8837	0.0000	1.0000
K2P	56	-1595.4416	3314.1623	0.0000	1.0000
K2PI	57	-1595.4411	3316.5849	0.0000	1.0000
SYM	60	-1594.8037	3322.6322	0.0000	1.0000
SYMI	61	-1594.8032	3325.0895	0.0000	1.0000
F81	58	-1603.6902	3335.5152	0.0000	1.0000
F81I	59	-1603.6898	3337.9550	0.0000	1.0000
JC69G	56	-1611.2067	3345.6925	0.0000	1.0000
JC69IG	57	-1611.2089	3348.1204	0.0000	1.0000
JC69	55	-1629.1083	3379.0809	0.0000	1.0000
JC69I	56	-1629.1077	3381.4945	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1547.0798	3214.1597	0.7230	0.7230
HKYIG	61	-1547.0821	3216.1642	0.2654	0.9884
GTRG	64	-1547.5268	3223.0536	0.0085	0.9969
GTRIG	65	-1547.5292	3225.0584	0.0031	1.0000
HKY	59	-1562.7130	3243.4261	0.0000	1.0000
HKYI	60	-1562.7127	3245.4253	0.0000	1.0000
GTR	63	-1564.0878	3254.1757	0.0000	1.0000
GTRI	64	-1564.0874	3256.1749	0.0000	1.0000
K2PG	57	-1577.1885	3268.3769	0.0000	1.0000
K2PIG	58	-1577.1909	3270.3817	0.0000	1.0000
SYMG	61	-1577.0489	3276.0978	0.0000	1.0000
SYMIG	62	-1577.0513	3278.1025	0.0000	1.0000
F81G	59	-1586.4270	3290.8540	0.0000	1.0000
F81IG	60	-1586.4294	3292.8588	0.0000	1.0000
K2P	56	-1595.4416	3302.8831	0.0000	1.0000
K2PI	57	-1595.4411	3304.8822	0.0000	1.0000
SYM	60	-1594.8037	3309.6073	0.0000	1.0000
SYMI	61	-1594.8032	3311.6064	0.0000	1.0000
F81	58	-1603.6902	3323.3804	0.0000	1.0000
F81I	59	-1603.6898	3325.3795	0.0000	1.0000
JC69G	56	-1611.2067	3334.4134	0.0000	1.0000
JC69IG	57	-1611.2089	3336.4178	0.0000	1.0000
JC69	55	-1629.1083	3368.2167	0.0000	1.0000
JC69I	56	-1629.1077	3370.2154	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1547.0798	3480.2325	0.9616	0.9616
HKYIG	61	-1547.0821	3486.6715	0.0384	1.0000
HKY	59	-1562.7130	3505.0643	0.0000	1.0000
GTRG	64	-1547.5268	3506.8646	0.0000	1.0000
HKYI	60	-1562.7127	3511.4981	0.0000	1.0000
GTRIG	65	-1547.5292	3513.3039	0.0000	1.0000
K2PG	57	-1577.1885	3521.1461	0.0000	1.0000
K2PIG	58	-1577.1909	3527.5854	0.0000	1.0000
GTR	63	-1564.0878	3533.5521	0.0000	1.0000
GTRI	64	-1564.0874	3539.9858	0.0000	1.0000
SYMG	61	-1577.0489	3546.6051	0.0000	1.0000
K2P	56	-1595.4416	3551.2177	0.0000	1.0000
F81G	59	-1586.4270	3552.4923	0.0000	1.0000
SYMIG	62	-1577.0513	3553.0444	0.0000	1.0000
K2PI	57	-1595.4411	3557.6514	0.0000	1.0000
F81IG	60	-1586.4294	3558.9316	0.0000	1.0000
SYM	60	-1594.8037	3575.6801	0.0000	1.0000
F81	58	-1603.6902	3580.5841	0.0000	1.0000
SYMI	61	-1594.8032	3582.1137	0.0000	1.0000
JC69G	56	-1611.2067	3582.7480	0.0000	1.0000
F81I	59	-1603.6898	3587.0178	0.0000	1.0000
JC69IG	57	-1611.2089	3589.1869	0.0000	1.0000
JC69	55	-1629.1083	3612.1167	0.0000	1.0000
JC69I	56	-1629.1077	3618.5500	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/10230_hg18.phylip = 1 - 623;
 Partition Dummy = 1:261-aggressive-phylip/10230_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/10230_hg18.phylip = 1 - 623;
 Partition Dummy = 1:261-aggressive-phylip/10230_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/10230_hg18.phylip = 1 - 623;
 Partition Dummy = 1:261-aggressive-phylip/10230_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

