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

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

Minimum AIC  model: GTRIG
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	-1417.2815	2979.5426	0.4017	0.4017
GTRIG	65	-1416.0285	2979.6030	0.3898	0.7915
GTRI	64	-1418.0303	2981.0401	0.1900	0.9815
HKYG	60	-1426.3481	2987.5139	0.0075	0.9890
HKYIG	61	-1425.2702	2987.8833	0.0062	0.9952
HKYI	60	-1426.7937	2988.4052	0.0048	1.0000
SYMIG	62	-1441.0585	3021.9951	0.0000	1.0000
SYMG	61	-1442.6689	3022.6806	0.0000	1.0000
SYMI	61	-1443.1088	3023.5604	0.0000	1.0000
GTR	63	-1442.0311	3026.4858	0.0000	1.0000
HKY	59	-1450.4610	3033.2251	0.0000	1.0000
K2PIG	58	-1452.6695	3035.1373	0.0000	1.0000
K2PG	57	-1453.9764	3035.2567	0.0000	1.0000
K2PI	57	-1454.2955	3035.8948	0.0000	1.0000
F81IG	60	-1453.4610	3041.7399	0.0000	1.0000
F81G	59	-1454.7803	3041.8637	0.0000	1.0000
F81I	59	-1455.2315	3042.7660	0.0000	1.0000
SYM	60	-1469.1710	3073.1598	0.0000	1.0000
K2P	56	-1479.4527	3083.7247	0.0000	1.0000
F81	58	-1479.2622	3088.3227	0.0000	1.0000
JC69IG	57	-1482.1373	3091.5784	0.0000	1.0000
JC69G	56	-1483.4254	3091.6702	0.0000	1.0000
JC69I	56	-1483.7959	3092.4111	0.0000	1.0000
JC69	55	-1508.3714	3139.0875	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1416.0285	2962.0569	0.4636	0.4636
GTRG	64	-1417.2815	2962.5630	0.3600	0.8236
GTRI	64	-1418.0303	2964.0605	0.1702	0.9938
HKYIG	61	-1425.2702	2972.5405	0.0025	0.9963
HKYG	60	-1426.3481	2972.6961	0.0023	0.9985
HKYI	60	-1426.7937	2973.5874	0.0015	1.0000
SYMIG	62	-1441.0585	3006.1171	0.0000	1.0000
SYMG	61	-1442.6689	3007.3378	0.0000	1.0000
SYMI	61	-1443.1088	3008.2176	0.0000	1.0000
GTR	63	-1442.0311	3010.0621	0.0000	1.0000
HKY	59	-1450.4610	3018.9220	0.0000	1.0000
K2PIG	58	-1452.6695	3021.3389	0.0000	1.0000
K2PG	57	-1453.9764	3021.9529	0.0000	1.0000
K2PI	57	-1454.2955	3022.5910	0.0000	1.0000
F81IG	60	-1453.4610	3026.9220	0.0000	1.0000
F81G	59	-1454.7803	3027.5606	0.0000	1.0000
F81I	59	-1455.2315	3028.4630	0.0000	1.0000
SYM	60	-1469.1710	3058.3420	0.0000	1.0000
K2P	56	-1479.4527	3070.9054	0.0000	1.0000
F81	58	-1479.2622	3074.5243	0.0000	1.0000
JC69IG	57	-1482.1373	3078.2746	0.0000	1.0000
JC69G	56	-1483.4254	3078.8509	0.0000	1.0000
JC69I	56	-1483.7959	3079.5919	0.0000	1.0000
JC69	55	-1508.3714	3126.7428	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1426.3481	3231.8342	0.5523	0.5523
HKYI	60	-1426.7937	3232.7254	0.3537	0.9059
HKYIG	61	-1425.2702	3235.9975	0.0689	0.9748
GTRG	64	-1417.2815	3238.9769	0.0155	0.9903
GTRI	64	-1418.0303	3240.4745	0.0073	0.9977
GTRIG	65	-1416.0285	3242.7899	0.0023	1.0000
K2PG	57	-1453.9764	3268.1340	0.0000	1.0000
K2PI	57	-1454.2955	3268.7721	0.0000	1.0000
SYMG	61	-1442.6689	3270.7948	0.0000	1.0000
SYMI	61	-1443.1088	3271.6747	0.0000	1.0000
K2PIG	58	-1452.6695	3271.8391	0.0000	1.0000
HKY	59	-1450.4610	3273.7412	0.0000	1.0000
SYMIG	62	-1441.0585	3273.8931	0.0000	1.0000
GTR	63	-1442.0311	3282.1571	0.0000	1.0000
F81G	59	-1454.7803	3282.3797	0.0000	1.0000
F81I	59	-1455.2315	3283.2821	0.0000	1.0000
F81IG	60	-1453.4610	3286.0601	0.0000	1.0000
K2P	56	-1479.4527	3312.7676	0.0000	1.0000
SYM	60	-1469.1710	3317.4801	0.0000	1.0000
JC69G	56	-1483.4254	3320.7131	0.0000	1.0000
JC69I	56	-1483.7959	3321.4541	0.0000	1.0000
JC69IG	57	-1482.1373	3324.4558	0.0000	1.0000
F81	58	-1479.2622	3325.0245	0.0000	1.0000
JC69	55	-1508.3714	3364.2861	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/10972_hg18.phylip = 1 - 555;
 Partition Dummy = 1:261-aggressive-phylip/10972_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 (GTRG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/10972_hg18.phylip = 1 - 555;
 Partition Dummy = 1:261-aggressive-phylip/10972_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/10972_hg18.phylip = 1 - 555;
 Partition Dummy = 1:261-aggressive-phylip/10972_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

