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

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

Minimum AIC  model: SYMG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1538.7242	3200.1255	0.4300	0.4300
SYMG	61	-1534.6405	3201.2599	0.2439	0.6739
K2PIG	58	-1538.7596	3202.5127	0.1303	0.8042
HKYG	60	-1536.8512	3203.3467	0.0859	0.8901
SYMIG	62	-1534.6762	3203.6721	0.0730	0.9631
HKYIG	61	-1536.8869	3205.7526	0.0258	0.9889
GTRG	64	-1534.4704	3207.9607	0.0086	0.9975
GTRIG	65	-1534.5062	3210.3918	0.0025	1.0000
JC69G	56	-1584.6928	3289.7526	0.0000	1.0000
JC69IG	57	-1584.7285	3292.1341	0.0000	1.0000
F81G	59	-1583.3888	3294.0933	0.0000	1.0000
F81IG	60	-1583.4245	3296.4933	0.0000	1.0000
K2P	56	-1595.2727	3310.9123	0.0000	1.0000
K2PI	57	-1595.3165	3313.3102	0.0000	1.0000
HKY	59	-1593.1815	3313.6788	0.0000	1.0000
SYM	60	-1592.0657	3313.7757	0.0000	1.0000
GTR	63	-1589.2880	3315.2427	0.0000	1.0000
HKYI	60	-1593.2116	3316.0674	0.0000	1.0000
SYMI	61	-1592.1105	3316.1999	0.0000	1.0000
GTRI	64	-1589.3181	3317.6560	0.0000	1.0000
JC69	55	-1640.3605	3398.7838	0.0000	1.0000
JC69I	56	-1640.4043	3401.1756	0.0000	1.0000
F81	58	-1638.6840	3402.3615	0.0000	1.0000
F81I	59	-1638.7143	3404.7444	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	61	-1534.6405	3191.2810	0.3241	0.3241
K2PG	57	-1538.7242	3191.4483	0.2981	0.6223
SYMIG	62	-1534.6762	3193.3525	0.1151	0.7373
K2PIG	58	-1538.7596	3193.5193	0.1059	0.8432
HKYG	60	-1536.8512	3193.7024	0.0966	0.9398
HKYIG	61	-1536.8869	3195.7737	0.0343	0.9741
GTRG	64	-1534.4704	3196.9408	0.0191	0.9932
GTRIG	65	-1534.5062	3199.0125	0.0068	1.0000
JC69G	56	-1584.6928	3281.3856	0.0000	1.0000
JC69IG	57	-1584.7285	3283.4569	0.0000	1.0000
F81G	59	-1583.3888	3284.7776	0.0000	1.0000
F81IG	60	-1583.4245	3286.8490	0.0000	1.0000
K2P	56	-1595.2727	3302.5454	0.0000	1.0000
SYM	60	-1592.0657	3304.1314	0.0000	1.0000
HKY	59	-1593.1815	3304.3630	0.0000	1.0000
GTR	63	-1589.2880	3304.5760	0.0000	1.0000
K2PI	57	-1595.3165	3304.6331	0.0000	1.0000
SYMI	61	-1592.1105	3306.2210	0.0000	1.0000
HKYI	60	-1593.2116	3306.4231	0.0000	1.0000
GTRI	64	-1589.3181	3306.6362	0.0000	1.0000
JC69	55	-1640.3605	3390.7210	0.0000	1.0000
JC69I	56	-1640.4043	3392.8087	0.0000	1.0000
F81	58	-1638.6840	3393.3680	0.0000	1.0000
F81I	59	-1638.7143	3395.4286	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1538.7242	3459.8786	0.9670	0.9670
K2PIG	58	-1538.7596	3466.6589	0.0326	0.9996
HKYG	60	-1536.8512	3476.2607	0.0003	0.9999
SYMG	61	-1534.6405	3478.5486	0.0001	1.0000
HKYIG	61	-1536.8869	3483.0413	0.0000	1.0000
SYMIG	62	-1534.6762	3485.3293	0.0000	1.0000
GTRG	64	-1534.4704	3498.3363	0.0000	1.0000
GTRIG	65	-1534.5062	3505.1173	0.0000	1.0000
JC69G	56	-1584.6928	3545.1066	0.0000	1.0000
JC69IG	57	-1584.7285	3551.8873	0.0000	1.0000
F81G	59	-1583.3888	3562.6265	0.0000	1.0000
K2P	56	-1595.2727	3566.2664	0.0000	1.0000
F81IG	60	-1583.4245	3569.4073	0.0000	1.0000
K2PI	57	-1595.3165	3573.0634	0.0000	1.0000
HKY	59	-1593.1815	3582.2120	0.0000	1.0000
SYM	60	-1592.0657	3586.6897	0.0000	1.0000
HKYI	60	-1593.2116	3588.9814	0.0000	1.0000
SYMI	61	-1592.1105	3593.4885	0.0000	1.0000
GTR	63	-1589.2880	3601.2622	0.0000	1.0000
GTRI	64	-1589.3181	3608.0316	0.0000	1.0000
JC69	55	-1640.3605	3649.7327	0.0000	1.0000
JC69I	56	-1640.4043	3656.5297	0.0000	1.0000
F81	58	-1638.6840	3666.5077	0.0000	1.0000
F81I	59	-1638.7143	3673.2775	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/4077_hg18.phylip = 1 - 820;
 Partition Dummy = 1:261-aggressive-phylip/4077_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=Fixed(Equal);
END;


[Mrbayes block for the best AICc model (K2PG)]

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


[Mrbayes block for the best BIC model (K2PG)]

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


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

