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

Input data from file "261-aggressive-phylip/7254_hg18.phylip" (29 taxa, 710 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	-1617.9573	3376.8138	0.5694	0.5694
GTRIG	65	-1617.3070	3377.9369	0.3247	0.8941
HKYG	60	-1624.7114	3380.7018	0.0815	0.9756
HKYIG	61	-1624.7191	3383.1111	0.0244	1.0000
K2PG	57	-1645.9152	3415.9715	0.0000	1.0000
K2PIG	58	-1645.9230	3418.3590	0.0000	1.0000
SYMG	61	-1643.5365	3420.7458	0.0000	1.0000
SYMIG	62	-1642.9929	3422.0600	0.0000	1.0000
F81G	59	-1648.8893	3426.6708	0.0000	1.0000
F81IG	60	-1648.3910	3428.0610	0.0000	1.0000
GTR	63	-1654.7194	3447.9218	0.0000	1.0000
GTRI	64	-1654.1273	3449.1537	0.0000	1.0000
HKY	59	-1661.4507	3451.7937	0.0000	1.0000
HKYI	60	-1661.0009	3453.2808	0.0000	1.0000
JC69G	56	-1666.9240	3455.6244	0.0000	1.0000
JC69IG	57	-1666.9317	3458.0045	0.0000	1.0000
K2P	56	-1682.2477	3486.2717	0.0000	1.0000
K2PI	57	-1681.9379	3488.0170	0.0000	1.0000
SYM	60	-1679.1025	3489.4838	0.0000	1.0000
SYMI	61	-1679.1067	3491.8863	0.0000	1.0000
F81	58	-1684.1336	3494.7802	0.0000	1.0000
F81I	59	-1684.1380	3497.1683	0.0000	1.0000
JC69	55	-1702.6962	3524.8113	0.0000	1.0000
JC69I	56	-1702.7005	3527.1773	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1617.9573	3363.9145	0.5581	0.5581
GTRIG	65	-1617.3070	3364.6140	0.3934	0.9515
HKYG	60	-1624.7114	3369.4229	0.0355	0.9870
HKYIG	61	-1624.7191	3371.4383	0.0130	1.0000
K2PG	57	-1645.9152	3405.8304	0.0000	1.0000
K2PIG	58	-1645.9230	3407.8460	0.0000	1.0000
SYMG	61	-1643.5365	3409.0730	0.0000	1.0000
SYMIG	62	-1642.9929	3409.9858	0.0000	1.0000
F81G	59	-1648.8893	3415.7785	0.0000	1.0000
F81IG	60	-1648.3910	3416.7821	0.0000	1.0000
GTR	63	-1654.7194	3435.4389	0.0000	1.0000
GTRI	64	-1654.1273	3436.2545	0.0000	1.0000
HKY	59	-1661.4507	3440.9014	0.0000	1.0000
HKYI	60	-1661.0009	3442.0019	0.0000	1.0000
JC69G	56	-1666.9240	3445.8480	0.0000	1.0000
JC69IG	57	-1666.9317	3447.8634	0.0000	1.0000
K2P	56	-1682.2477	3476.4953	0.0000	1.0000
K2PI	57	-1681.9379	3477.8759	0.0000	1.0000
SYM	60	-1679.1025	3478.2049	0.0000	1.0000
SYMI	61	-1679.1067	3480.2135	0.0000	1.0000
F81	58	-1684.1336	3484.2672	0.0000	1.0000
F81I	59	-1684.1380	3486.2760	0.0000	1.0000
JC69	55	-1702.6962	3515.3923	0.0000	1.0000
JC69I	56	-1702.7005	3517.4009	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1624.7114	3643.3388	0.9624	0.9624
HKYIG	61	-1624.7191	3649.9194	0.0358	0.9982
GTRG	64	-1617.9573	3656.0915	0.0016	0.9999
GTRIG	65	-1617.3070	3661.3562	0.0001	1.0000
K2PG	57	-1645.9152	3666.0505	0.0000	1.0000
K2PIG	58	-1645.9230	3672.6313	0.0000	1.0000
F81G	59	-1648.8893	3685.1292	0.0000	1.0000
SYMG	61	-1643.5365	3687.5541	0.0000	1.0000
F81IG	60	-1648.3910	3690.6980	0.0000	1.0000
SYMIG	62	-1642.9929	3693.0323	0.0000	1.0000
JC69G	56	-1666.9240	3701.5029	0.0000	1.0000
JC69IG	57	-1666.9317	3708.0835	0.0000	1.0000
HKY	59	-1661.4507	3710.2520	0.0000	1.0000
HKYI	60	-1661.0009	3715.9178	0.0000	1.0000
GTR	63	-1654.7194	3723.0506	0.0000	1.0000
GTRI	64	-1654.1273	3728.4315	0.0000	1.0000
K2P	56	-1682.2477	3732.1502	0.0000	1.0000
K2PI	57	-1681.9379	3738.0960	0.0000	1.0000
F81	58	-1684.1336	3749.0525	0.0000	1.0000
SYM	60	-1679.1025	3752.1208	0.0000	1.0000
F81I	59	-1684.1380	3755.6266	0.0000	1.0000
SYMI	61	-1679.1067	3758.6946	0.0000	1.0000
JC69	55	-1702.6962	3766.4819	0.0000	1.0000
JC69I	56	-1702.7005	3773.0557	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/7254_hg18.phylip = 1 - 710;
 Partition Dummy = 1:261-aggressive-phylip/7254_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/7254_hg18.phylip = 1 - 710;
 Partition Dummy = 1:261-aggressive-phylip/7254_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/7254_hg18.phylip = 1 - 710;
 Partition Dummy = 1:261-aggressive-phylip/7254_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

