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

Input data from file "261-aggressive-phylip/1970_hg18.phylip" (29 taxa, 798 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	-1469.6807	3069.2936	0.6904	0.6904
HKYIG	61	-1469.7074	3071.6920	0.2081	0.8985
GTRG	64	-1467.1779	3073.7065	0.0760	0.9745
GTRIG	65	-1467.2047	3076.1307	0.0226	0.9971
HKY	59	-1476.6560	3080.9056	0.0021	0.9992
HKYI	60	-1476.6812	3083.2945	0.0006	0.9998
GTR	63	-1474.7631	3086.5125	0.0001	1.0000
GTRI	64	-1474.7881	3088.9268	0.0000	1.0000
F81G	59	-1494.3203	3116.2340	0.0000	1.0000
F81IG	60	-1494.3496	3118.6313	0.0000	1.0000
F81	58	-1500.3959	3126.0530	0.0000	1.0000
F81I	59	-1500.4226	3128.4387	0.0000	1.0000
K2PG	57	-1513.3014	3149.5380	0.0000	1.0000
K2PIG	58	-1513.3282	3151.9176	0.0000	1.0000
SYMG	61	-1512.6869	3157.6510	0.0000	1.0000
K2P	56	-1519.6190	3159.8534	0.0000	1.0000
SYMIG	62	-1512.7137	3160.0560	0.0000	1.0000
K2PI	57	-1519.6442	3162.2235	0.0000	1.0000
SYM	60	-1518.8275	3167.5872	0.0000	1.0000
SYMI	61	-1518.8528	3169.9828	0.0000	1.0000
JC69G	56	-1531.6956	3184.0066	0.0000	1.0000
JC69IG	57	-1531.7241	3186.3833	0.0000	1.0000
JC69	55	-1538.0541	3194.4101	0.0000	1.0000
JC69I	56	-1538.0812	3196.7778	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1469.6807	3059.3614	0.6003	0.6003
HKYIG	61	-1469.7074	3061.4148	0.2150	0.8153
GTRG	64	-1467.1779	3062.3559	0.1343	0.9496
GTRIG	65	-1467.2047	3064.4094	0.0481	0.9977
HKY	59	-1476.6560	3071.3121	0.0015	0.9992
HKYI	60	-1476.6812	3073.3623	0.0005	0.9997
GTR	63	-1474.7631	3075.5261	0.0002	0.9999
GTRI	64	-1474.7881	3077.5762	0.0001	1.0000
F81G	59	-1494.3203	3106.6405	0.0000	1.0000
F81IG	60	-1494.3496	3108.6992	0.0000	1.0000
F81	58	-1500.3959	3116.7918	0.0000	1.0000
F81I	59	-1500.4226	3118.8452	0.0000	1.0000
K2PG	57	-1513.3014	3140.6029	0.0000	1.0000
K2PIG	58	-1513.3282	3142.6564	0.0000	1.0000
SYMG	61	-1512.6869	3147.3738	0.0000	1.0000
SYMIG	62	-1512.7137	3149.4274	0.0000	1.0000
K2P	56	-1519.6190	3151.2380	0.0000	1.0000
K2PI	57	-1519.6442	3153.2884	0.0000	1.0000
SYM	60	-1518.8275	3157.6550	0.0000	1.0000
SYMI	61	-1518.8528	3159.7056	0.0000	1.0000
JC69G	56	-1531.6956	3175.3912	0.0000	1.0000
JC69IG	57	-1531.7241	3177.4481	0.0000	1.0000
JC69	55	-1538.0541	3186.1082	0.0000	1.0000
JC69I	56	-1538.0812	3188.1624	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1469.6807	3340.2880	0.9418	0.9418
HKYIG	61	-1469.7074	3347.0234	0.0325	0.9743
HKY	59	-1476.6560	3347.5565	0.0249	0.9991
HKYI	60	-1476.6812	3354.2888	0.0009	1.0000
GTRG	64	-1467.1779	3362.0108	0.0000	1.0000
GTRIG	65	-1467.2047	3368.7465	0.0000	1.0000
GTR	63	-1474.7631	3370.4989	0.0000	1.0000
GTRI	64	-1474.7881	3377.2312	0.0000	1.0000
F81G	59	-1494.3203	3382.8849	0.0000	1.0000
F81	58	-1500.3959	3388.3541	0.0000	1.0000
F81IG	60	-1494.3496	3389.6257	0.0000	1.0000
F81I	59	-1500.4226	3395.0896	0.0000	1.0000
K2PG	57	-1513.3014	3407.4831	0.0000	1.0000
K2P	56	-1519.6190	3413.4361	0.0000	1.0000
K2PIG	58	-1513.3282	3414.2187	0.0000	1.0000
K2PI	57	-1519.6442	3420.1686	0.0000	1.0000
SYMG	61	-1512.6869	3432.9824	0.0000	1.0000
JC69G	56	-1531.6956	3437.5893	0.0000	1.0000
SYM	60	-1518.8275	3438.5815	0.0000	1.0000
SYMIG	62	-1512.7137	3439.7181	0.0000	1.0000
JC69	55	-1538.0541	3443.6242	0.0000	1.0000
JC69IG	57	-1531.7241	3444.3283	0.0000	1.0000
SYMI	61	-1518.8528	3445.3142	0.0000	1.0000
JC69I	56	-1538.0812	3450.3605	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/1970_hg18.phylip = 1 - 798;
 Partition Dummy = 1:261-aggressive-phylip/1970_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/1970_hg18.phylip = 1 - 798;
 Partition Dummy = 1:261-aggressive-phylip/1970_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/1970_hg18.phylip = 1 - 798;
 Partition Dummy = 1:261-aggressive-phylip/1970_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

