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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-2908.8807	5945.5403	0.4818	0.4818
GTRG	64	-2905.1038	5947.0869	0.2223	0.7041
HKYIG	61	-2908.9178	5947.8825	0.1494	0.8535
SYMG	61	-2909.4415	5948.9299	0.0885	0.9420
GTRIG	65	-2906.1705	5951.5076	0.0244	0.9664
K2PG	57	-2915.6120	5952.2283	0.0170	0.9834
SYMIG	62	-2910.3614	5953.0422	0.0113	0.9947
K2PIG	58	-2915.6491	5954.5559	0.0053	1.0000
HKY	59	-2953.4423	6032.4005	0.0000	1.0000
GTR	63	-2948.9759	6032.5488	0.0000	1.0000
SYM	60	-2953.4863	6034.7516	0.0000	1.0000
GTRI	64	-2949.0014	6034.8821	0.0000	1.0000
SYMI	61	-2953.5080	6037.0627	0.0000	1.0000
K2P	56	-2959.2802	6037.3160	0.0000	1.0000
HKYI	60	-2957.3510	6042.4810	0.0000	1.0000
K2PI	57	-2963.1321	6047.2684	0.0000	1.0000
F81G	59	-2971.0037	6067.5233	0.0000	1.0000
F81IG	60	-2971.0408	6069.8606	0.0000	1.0000
JC69G	56	-2979.0835	6076.9226	0.0000	1.0000
JC69IG	57	-2979.1207	6079.2457	0.0000	1.0000
F81	58	-3015.9543	6155.1662	0.0000	1.0000
F81I	59	-3015.9798	6157.4756	0.0000	1.0000
JC69	55	-3023.6789	6163.8694	0.0000	1.0000
JC69I	56	-3023.7043	6166.1642	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-2908.8807	5937.7614	0.3954	0.3954
GTRG	64	-2905.1038	5938.2075	0.3164	0.7118
HKYIG	61	-2908.9178	5939.8357	0.1402	0.8519
SYMG	61	-2909.4415	5940.8831	0.0830	0.9349
GTRIG	65	-2906.1705	5942.3409	0.0401	0.9750
SYMIG	62	-2910.3614	5944.7227	0.0122	0.9872
K2PG	57	-2915.6120	5945.2240	0.0095	0.9966
K2PIG	58	-2915.6491	5947.2982	0.0034	1.0000
GTR	63	-2948.9759	6023.9518	0.0000	1.0000
HKY	59	-2953.4423	6024.8846	0.0000	1.0000
GTRI	64	-2949.0014	6026.0027	0.0000	1.0000
SYM	60	-2953.4863	6026.9726	0.0000	1.0000
SYMI	61	-2953.5080	6029.0159	0.0000	1.0000
K2P	56	-2959.2802	6030.5604	0.0000	1.0000
HKYI	60	-2957.3510	6034.7020	0.0000	1.0000
K2PI	57	-2963.1321	6040.2641	0.0000	1.0000
F81G	59	-2971.0037	6060.0074	0.0000	1.0000
F81IG	60	-2971.0408	6062.0816	0.0000	1.0000
JC69G	56	-2979.0835	6070.1670	0.0000	1.0000
JC69IG	57	-2979.1207	6072.2415	0.0000	1.0000
F81	58	-3015.9543	6147.9085	0.0000	1.0000
F81I	59	-3015.9798	6149.9596	0.0000	1.0000
JC69	55	-3023.6789	6157.3578	0.0000	1.0000
JC69I	56	-3023.7043	6159.4086	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2915.6120	6225.0800	0.9450	0.9450
K2PIG	58	-2915.6491	6232.0639	0.0288	0.9738
HKYG	60	-2908.8807	6232.3466	0.0250	0.9988
HKYIG	61	-2908.9178	6239.3306	0.0008	0.9995
SYMG	61	-2909.4415	6240.3780	0.0005	1.0000
SYMIG	62	-2910.3614	6249.1274	0.0000	1.0000
GTRG	64	-2905.1038	6252.4318	0.0000	1.0000
GTRIG	65	-2906.1705	6261.4749	0.0000	1.0000
K2P	56	-2959.2802	6305.5066	0.0000	1.0000
HKY	59	-2953.4423	6314.5600	0.0000	1.0000
K2PI	57	-2963.1321	6320.1201	0.0000	1.0000
SYM	60	-2953.4863	6321.5578	0.0000	1.0000
SYMI	61	-2953.5080	6328.5109	0.0000	1.0000
HKYI	60	-2957.3510	6329.2872	0.0000	1.0000
GTR	63	-2948.9759	6333.2662	0.0000	1.0000
GTRI	64	-2949.0014	6340.2270	0.0000	1.0000
JC69G	56	-2979.0835	6345.1132	0.0000	1.0000
F81G	59	-2971.0037	6349.6828	0.0000	1.0000
JC69IG	57	-2979.1207	6352.0974	0.0000	1.0000
F81IG	60	-2971.0408	6356.6668	0.0000	1.0000
JC69	55	-3023.6789	6427.3942	0.0000	1.0000
F81	58	-3015.9543	6432.6742	0.0000	1.0000
JC69I	56	-3023.7043	6434.3548	0.0000	1.0000
F81I	59	-3015.9798	6439.6351	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/6657_hg18.phylip = 1 - 1002;
 Partition Dummy = 1:261-aggressive-phylip/6657_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/6657_hg18.phylip = 1 - 1002;
 Partition Dummy = 1:261-aggressive-phylip/6657_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 (K2PG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/6657_hg18.phylip = 1 - 1002;
 Partition Dummy = 1:261-aggressive-phylip/6657_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

