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

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

Minimum AIC  model: GTRIG
Minimum AICc model: HKYIG
Minimum BIC  model: HKYIG


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
HKYIG	61	-1512.5453	3179.0061	0.5245	0.5245
GTRIG	65	-1506.2686	3179.3612	0.4392	0.9638
HKYG	60	-1517.1095	3184.9752	0.0265	0.9903
GTRG	64	-1511.8018	3187.1592	0.0089	0.9992
K2PIG	58	-1523.9323	3192.3812	0.0007	0.9998
SYMIG	62	-1519.2552	3195.6122	0.0001	1.0000
K2PG	57	-1529.2051	3199.8458	0.0000	1.0000
SYMG	61	-1523.7622	3201.4401	0.0000	1.0000
F81IG	60	-1530.1342	3211.0247	0.0000	1.0000
F81G	59	-1533.6711	3214.9656	0.0000	1.0000
GTRI	64	-1530.5050	3224.5655	0.0000	1.0000
HKYI	60	-1537.0807	3224.9177	0.0000	1.0000
JC69IG	57	-1545.9780	3233.3917	0.0000	1.0000
JC69G	56	-1550.7299	3239.8399	0.0000	1.0000
K2PI	57	-1551.0171	3243.4700	0.0000	1.0000
SYMI	61	-1544.9309	3243.7773	0.0000	1.0000
F81I	59	-1555.0991	3257.8216	0.0000	1.0000
JC69I	56	-1570.1420	3278.6642	0.0000	1.0000
GTR	63	-1608.3621	3377.0390	0.0000	1.0000
HKY	59	-1618.1932	3384.0099	0.0000	1.0000
SYM	60	-1628.9446	3408.6456	0.0000	1.0000
K2P	56	-1635.9345	3410.2492	0.0000	1.0000
F81	58	-1636.5583	3417.6332	0.0000	1.0000
JC69	55	-1653.7813	3442.9125	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1506.2686	3142.5372	0.8958	0.8958
HKYIG	61	-1512.5453	3147.0905	0.0919	0.9877
GTRG	64	-1511.8018	3151.6037	0.0096	0.9973
HKYG	60	-1517.1095	3154.2189	0.0026	0.9999
SYMIG	62	-1519.2552	3162.5105	0.0000	1.0000
K2PIG	58	-1523.9323	3163.8646	0.0000	1.0000
SYMG	61	-1523.7622	3169.5245	0.0000	1.0000
K2PG	57	-1529.2051	3172.4102	0.0000	1.0000
F81IG	60	-1530.1342	3180.2684	0.0000	1.0000
F81G	59	-1533.6711	3185.3421	0.0000	1.0000
GTRI	64	-1530.5050	3189.0099	0.0000	1.0000
HKYI	60	-1537.0807	3194.1614	0.0000	1.0000
JC69IG	57	-1545.9780	3205.9560	0.0000	1.0000
SYMI	61	-1544.9309	3211.8617	0.0000	1.0000
JC69G	56	-1550.7299	3213.4597	0.0000	1.0000
K2PI	57	-1551.0171	3216.0343	0.0000	1.0000
F81I	59	-1555.0991	3228.1981	0.0000	1.0000
JC69I	56	-1570.1420	3252.2840	0.0000	1.0000
GTR	63	-1608.3621	3342.7241	0.0000	1.0000
HKY	59	-1618.1932	3354.3865	0.0000	1.0000
SYM	60	-1628.9446	3377.8893	0.0000	1.0000
K2P	56	-1635.9345	3383.8690	0.0000	1.0000
F81	58	-1636.5583	3389.1166	0.0000	1.0000
JC69	55	-1653.7813	3417.5627	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-1512.5453	3372.8176	0.7997	0.7997
HKYG	60	-1517.1095	3376.2455	0.1441	0.9438
K2PIG	58	-1523.9323	3378.4903	0.0469	0.9907
GTRIG	65	-1506.2686	3383.0660	0.0048	0.9954
K2PG	57	-1529.2051	3383.3354	0.0042	0.9996
GTRG	64	-1511.8018	3388.4320	0.0003	0.9999
SYMIG	62	-1519.2552	3391.9380	0.0001	1.0000
SYMG	61	-1523.7622	3395.2516	0.0000	1.0000
F81IG	60	-1530.1342	3402.2950	0.0000	1.0000
F81G	59	-1533.6711	3403.6683	0.0000	1.0000
HKYI	60	-1537.0807	3416.1880	0.0000	1.0000
JC69IG	57	-1545.9780	3416.8813	0.0000	1.0000
JC69G	56	-1550.7299	3420.6846	0.0000	1.0000
GTRI	64	-1530.5050	3425.8383	0.0000	1.0000
K2PI	57	-1551.0171	3426.9596	0.0000	1.0000
SYMI	61	-1544.9309	3437.5888	0.0000	1.0000
F81I	59	-1555.0991	3446.5243	0.0000	1.0000
JC69I	56	-1570.1420	3459.5089	0.0000	1.0000
HKY	59	-1618.1932	3572.7126	0.0000	1.0000
GTR	63	-1608.3621	3575.8521	0.0000	1.0000
K2P	56	-1635.9345	3591.0938	0.0000	1.0000
SYM	60	-1628.9446	3599.9159	0.0000	1.0000
F81	58	-1636.5583	3603.7423	0.0000	1.0000
JC69	55	-1653.7813	3621.0871	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9944_hg18.phylip = 1 - 299;
 Partition Dummy = 1:261-aggressive-phylip/9944_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9944_hg18.phylip = 1 - 299;
 Partition Dummy = 1:261-aggressive-phylip/9944_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/9944_hg18.phylip = 1 - 299;
 Partition Dummy = 1:261-aggressive-phylip/9944_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

