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

Input data from file "261-aggressive-phylip/1644_hg18.phylip" (29 taxa, 383 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	-1605.4767	3365.1169	0.5789	0.5789
HKYG	60	-1612.1711	3367.0751	0.2174	0.7963
GTRIG	65	-1605.4190	3367.9042	0.1437	0.9400
HKYIG	61	-1612.0443	3369.6525	0.0599	0.9999
GTRI	64	-1614.8210	3383.8055	0.0001	1.0000
HKYI	60	-1620.8331	3384.3992	0.0000	1.0000
GTR	63	-1632.5660	3416.4110	0.0000	1.0000
HKY	59	-1640.6567	3421.2328	0.0000	1.0000
K2PG	57	-1644.4616	3423.2678	0.0000	1.0000
K2PIG	58	-1644.2202	3425.5639	0.0000	1.0000
F81G	59	-1643.0790	3426.0775	0.0000	1.0000
F81IG	60	-1642.9381	3428.6091	0.0000	1.0000
SYMG	61	-1642.9183	3431.4005	0.0000	1.0000
SYMIG	62	-1642.7105	3433.8336	0.0000	1.0000
K2PI	57	-1651.2957	3436.9360	0.0000	1.0000
F81I	59	-1649.7476	3439.4147	0.0000	1.0000
SYMI	61	-1649.9937	3445.5513	0.0000	1.0000
JC69G	56	-1669.6301	3470.8431	0.0000	1.0000
JC69IG	57	-1669.4358	3473.2162	0.0000	1.0000
F81	58	-1670.3476	3477.8186	0.0000	1.0000
K2P	56	-1674.4659	3480.5146	0.0000	1.0000
JC69I	56	-1676.1285	3483.8398	0.0000	1.0000
SYM	60	-1672.0276	3486.7882	0.0000	1.0000
JC69	55	-1698.8189	3526.4757	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1605.4767	3338.9533	0.6731	0.6731
GTRIG	65	-1605.4190	3340.8379	0.2623	0.9354
HKYG	60	-1612.1711	3344.3422	0.0455	0.9809
HKYIG	61	-1612.0443	3346.0886	0.0190	0.9999
GTRI	64	-1614.8210	3357.6420	0.0001	1.0000
HKYI	60	-1620.8331	3361.6663	0.0000	1.0000
GTR	63	-1632.5660	3391.1320	0.0000	1.0000
HKY	59	-1640.6567	3399.3133	0.0000	1.0000
K2PG	57	-1644.4616	3402.9232	0.0000	1.0000
F81G	59	-1643.0790	3404.1580	0.0000	1.0000
K2PIG	58	-1644.2202	3404.4405	0.0000	1.0000
F81IG	60	-1642.9381	3405.8762	0.0000	1.0000
SYMG	61	-1642.9183	3407.8367	0.0000	1.0000
SYMIG	62	-1642.7105	3409.4211	0.0000	1.0000
K2PI	57	-1651.2957	3416.5914	0.0000	1.0000
F81I	59	-1649.7476	3417.4952	0.0000	1.0000
SYMI	61	-1649.9937	3421.9874	0.0000	1.0000
JC69G	56	-1669.6301	3451.2603	0.0000	1.0000
JC69IG	57	-1669.4358	3452.8716	0.0000	1.0000
F81	58	-1670.3476	3456.6951	0.0000	1.0000
K2P	56	-1674.4659	3460.9318	0.0000	1.0000
SYM	60	-1672.0276	3464.0553	0.0000	1.0000
JC69I	56	-1676.1285	3464.2570	0.0000	1.0000
JC69	55	-1698.8189	3507.6378	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1612.1711	3581.2243	0.9399	0.9399
HKYIG	61	-1612.0443	3586.9188	0.0545	0.9944
GTRG	64	-1605.4767	3591.6276	0.0052	0.9996
GTRIG	65	-1605.4190	3597.4602	0.0003	0.9998
HKYI	60	-1620.8331	3598.5484	0.0002	1.0000
GTRI	64	-1614.8210	3610.3162	0.0000	1.0000
K2PG	57	-1644.4616	3627.9612	0.0000	1.0000
HKY	59	-1640.6567	3632.2474	0.0000	1.0000
K2PIG	58	-1644.2202	3633.4265	0.0000	1.0000
F81G	59	-1643.0790	3637.0921	0.0000	1.0000
GTR	63	-1632.5660	3639.8582	0.0000	1.0000
K2PI	57	-1651.2957	3641.6294	0.0000	1.0000
F81IG	60	-1642.9381	3642.7583	0.0000	1.0000
SYMG	61	-1642.9183	3648.6668	0.0000	1.0000
F81I	59	-1649.7476	3650.4293	0.0000	1.0000
SYMIG	62	-1642.7105	3654.1992	0.0000	1.0000
SYMI	61	-1649.9937	3662.8176	0.0000	1.0000
JC69G	56	-1669.6301	3672.3503	0.0000	1.0000
JC69IG	57	-1669.4358	3677.9096	0.0000	1.0000
K2P	56	-1674.4659	3682.0218	0.0000	1.0000
JC69I	56	-1676.1285	3685.3470	0.0000	1.0000
F81	58	-1670.3476	3685.6812	0.0000	1.0000
SYM	60	-1672.0276	3700.9374	0.0000	1.0000
JC69	55	-1698.8189	3724.7797	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/1644_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/1644_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/1644_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/1644_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/1644_hg18.phylip = 1 - 383;
 Partition Dummy = 1:261-aggressive-phylip/1644_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

