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

Input data from file "261-aggressive-phylip/7721_hg18.phylip" (29 taxa, 741 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	-4513.2187	9166.7452	0.7688	0.7688
GTRIG	65	-4513.2217	9169.1544	0.2305	0.9993
HKYG	60	-4525.3228	9181.4102	0.0005	0.9998
HKYIG	61	-4525.3259	9183.7916	0.0002	1.0000
GTRI	64	-4571.2588	9282.8253	0.0000	1.0000
HKYI	60	-4576.7154	9284.1956	0.0000	1.0000
SYMG	61	-4589.8244	9312.7888	0.0000	1.0000
K2PG	57	-4595.6005	9314.8819	0.0000	1.0000
SYMIG	62	-4589.8275	9315.1772	0.0000	1.0000
K2PIG	58	-4595.6036	9317.2424	0.0000	1.0000
SYMI	61	-4643.5921	9420.3241	0.0000	1.0000
K2PI	57	-4652.8573	9429.3953	0.0000	1.0000
GTR	63	-4674.6509	9487.2132	0.0000	1.0000
F81G	59	-4688.0854	9504.5672	0.0000	1.0000
HKY	59	-4688.5109	9505.4182	0.0000	1.0000
F81IG	60	-4688.0709	9506.9065	0.0000	1.0000
JC69G	56	-4731.9146	9585.1626	0.0000	1.0000
JC69IG	57	-4731.9181	9587.5171	0.0000	1.0000
F81I	59	-4737.2957	9602.9879	0.0000	1.0000
SYM	60	-4756.8052	9644.3750	0.0000	1.0000
K2P	56	-4762.9972	9647.3277	0.0000	1.0000
JC69I	56	-4782.5693	9686.4719	0.0000	1.0000
F81	58	-4836.7841	9799.6035	0.0000	1.0000
JC69	55	-4882.3927	9883.7782	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4513.2187	9154.4375	0.7314	0.7314
GTRIG	65	-4513.2217	9156.4433	0.2683	0.9997
HKYG	60	-4525.3228	9170.6455	0.0002	0.9999
HKYIG	61	-4525.3259	9172.6517	0.0001	1.0000
GTRI	64	-4571.2588	9270.5176	0.0000	1.0000
HKYI	60	-4576.7154	9273.4309	0.0000	1.0000
SYMG	61	-4589.8244	9301.6489	0.0000	1.0000
SYMIG	62	-4589.8275	9303.6551	0.0000	1.0000
K2PG	57	-4595.6005	9305.2011	0.0000	1.0000
K2PIG	58	-4595.6036	9307.2072	0.0000	1.0000
SYMI	61	-4643.5921	9409.1842	0.0000	1.0000
K2PI	57	-4652.8573	9419.7145	0.0000	1.0000
GTR	63	-4674.6509	9475.3018	0.0000	1.0000
F81G	59	-4688.0854	9494.1707	0.0000	1.0000
HKY	59	-4688.5109	9495.0217	0.0000	1.0000
F81IG	60	-4688.0709	9496.1418	0.0000	1.0000
JC69G	56	-4731.9146	9575.8293	0.0000	1.0000
JC69IG	57	-4731.9181	9577.8363	0.0000	1.0000
F81I	59	-4737.2957	9592.5914	0.0000	1.0000
SYM	60	-4756.8052	9633.6103	0.0000	1.0000
K2P	56	-4762.9972	9637.9943	0.0000	1.0000
JC69I	56	-4782.5693	9677.1386	0.0000	1.0000
F81	58	-4836.7841	9789.5683	0.0000	1.0000
JC69	55	-4882.3927	9874.7855	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-4525.3228	9447.1256	0.7259	0.7259
GTRG	64	-4513.2187	9449.3495	0.2388	0.9647
HKYIG	61	-4525.3259	9453.7398	0.0266	0.9913
GTRIG	65	-4513.2217	9455.9634	0.0087	1.0000
HKYI	60	-4576.7154	9549.9109	0.0000	1.0000
GTRI	64	-4571.2588	9565.4296	0.0000	1.0000
K2PG	57	-4595.6005	9567.8571	0.0000	1.0000
K2PIG	58	-4595.6036	9574.4712	0.0000	1.0000
SYMG	61	-4589.8244	9582.7369	0.0000	1.0000
SYMIG	62	-4589.8275	9589.3511	0.0000	1.0000
K2PI	57	-4652.8573	9682.3705	0.0000	1.0000
SYMI	61	-4643.5921	9690.2722	0.0000	1.0000
GTR	63	-4674.6509	9765.6058	0.0000	1.0000
F81G	59	-4688.0854	9766.0427	0.0000	1.0000
HKY	59	-4688.5109	9766.8938	0.0000	1.0000
F81IG	60	-4688.0709	9772.6219	0.0000	1.0000
JC69G	56	-4731.9146	9833.8773	0.0000	1.0000
JC69IG	57	-4731.9181	9840.4923	0.0000	1.0000
F81I	59	-4737.2957	9864.4634	0.0000	1.0000
K2P	56	-4762.9972	9896.0424	0.0000	1.0000
SYM	60	-4756.8052	9910.0904	0.0000	1.0000
JC69I	56	-4782.5693	9935.1866	0.0000	1.0000
F81	58	-4836.7841	10056.8323	0.0000	1.0000
JC69	55	-4882.3927	10128.2255	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/7721_hg18.phylip = 1 - 741;
 Partition Dummy = 1:261-aggressive-phylip/7721_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/7721_hg18.phylip = 1 - 741;
 Partition Dummy = 1:261-aggressive-phylip/7721_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/7721_hg18.phylip = 1 - 741;
 Partition Dummy = 1:261-aggressive-phylip/7721_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

