
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:49 04/04/2011

Input data from file "591-tomodel-phylip/chr8_3277.phylip" (7 taxa, 473 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	16	-1428.4055	2890.0040	0.4580	0.4580
HKYI	16	-1428.7776	2890.7482	0.3157	0.7737
HKYIG	17	-1428.4061	2892.1572	0.1561	0.9297
GTRG	20	-1426.6716	2895.2016	0.0341	0.9638
GTRI	20	-1426.9916	2895.8416	0.0247	0.9885
GTRIG	21	-1426.6657	2897.3802	0.0115	1.0000
HKY	15	-1441.6713	2914.3929	0.0000	1.0000
GTR	19	-1440.0783	2919.8342	0.0000	1.0000
F81G	15	-1454.5177	2940.0857	0.0000	1.0000
F81I	15	-1454.5995	2940.2493	0.0000	1.0000
F81IG	16	-1454.5436	2942.2802	0.0000	1.0000
F81	14	-1467.0050	2962.9270	0.0000	1.0000
K2PG	13	-1485.7347	2998.2624	0.0000	1.0000
K2PI	13	-1485.9686	2998.7302	0.0000	1.0000
K2PIG	14	-1485.7400	3000.3970	0.0000	1.0000
SYMG	17	-1483.2227	3001.7905	0.0000	1.0000
SYMI	17	-1483.4483	3002.2416	0.0000	1.0000
SYMIG	18	-1483.2333	3003.9733	0.0000	1.0000
K2P	12	-1499.9785	3024.6352	0.0000	1.0000
SYM	16	-1498.0092	3029.2113	0.0000	1.0000
JC69G	12	-1503.8731	3032.4244	0.0000	1.0000
JC69I	12	-1503.9424	3032.5630	0.0000	1.0000
JC69IG	13	-1503.9003	3034.5937	0.0000	1.0000
JC69	11	-1517.1132	3056.7991	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1428.4055	2888.8110	0.4397	0.4397
HKYI	16	-1428.7776	2889.5552	0.3031	0.7427
HKYIG	17	-1428.4061	2890.8121	0.1617	0.9044
GTRG	20	-1426.6716	2893.3432	0.0456	0.9500
GTRI	20	-1426.9916	2893.9832	0.0331	0.9831
GTRIG	21	-1426.6657	2895.3314	0.0169	1.0000
HKY	15	-1441.6713	2913.3426	0.0000	1.0000
GTR	19	-1440.0783	2918.1565	0.0000	1.0000
F81G	15	-1454.5177	2939.0354	0.0000	1.0000
F81I	15	-1454.5995	2939.1990	0.0000	1.0000
F81IG	16	-1454.5436	2941.0872	0.0000	1.0000
F81	14	-1467.0050	2962.0100	0.0000	1.0000
K2PG	13	-1485.7347	2997.4694	0.0000	1.0000
K2PI	13	-1485.9686	2997.9372	0.0000	1.0000
K2PIG	14	-1485.7400	2999.4799	0.0000	1.0000
SYMG	17	-1483.2227	3000.4454	0.0000	1.0000
SYMI	17	-1483.4483	3000.8966	0.0000	1.0000
SYMIG	18	-1483.2333	3002.4667	0.0000	1.0000
K2P	12	-1499.9785	3023.9569	0.0000	1.0000
SYM	16	-1498.0092	3028.0183	0.0000	1.0000
JC69G	12	-1503.8731	3031.7461	0.0000	1.0000
JC69I	12	-1503.9424	3031.8847	0.0000	1.0000
JC69IG	13	-1503.9003	3033.8007	0.0000	1.0000
JC69	11	-1517.1132	3056.2265	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1428.4055	2955.3566	0.5763	0.5763
HKYI	16	-1428.7776	2956.1007	0.3972	0.9735
HKYIG	17	-1428.4061	2961.5168	0.0265	1.0000
HKY	15	-1441.6713	2975.7290	0.0000	1.0000
GTRG	20	-1426.6716	2976.5251	0.0000	1.0000
GTRI	20	-1426.9916	2977.1651	0.0000	1.0000
GTRIG	21	-1426.6657	2982.6724	0.0000	1.0000
GTR	19	-1440.0783	2997.1793	0.0000	1.0000
F81G	15	-1454.5177	3001.4219	0.0000	1.0000
F81I	15	-1454.5995	3001.5854	0.0000	1.0000
F81IG	16	-1454.5436	3007.6327	0.0000	1.0000
F81	14	-1467.0050	3020.2373	0.0000	1.0000
K2PG	13	-1485.7347	3051.5376	0.0000	1.0000
K2PI	13	-1485.9686	3052.0054	0.0000	1.0000
K2PIG	14	-1485.7400	3057.7073	0.0000	1.0000
SYMG	17	-1483.2227	3071.1500	0.0000	1.0000
SYMI	17	-1483.4483	3071.6012	0.0000	1.0000
K2P	12	-1499.9785	3073.8660	0.0000	1.0000
SYMIG	18	-1483.2333	3077.3304	0.0000	1.0000
JC69G	12	-1503.8731	3081.6553	0.0000	1.0000
JC69I	12	-1503.9424	3081.7939	0.0000	1.0000
JC69IG	13	-1503.9003	3087.8689	0.0000	1.0000
SYM	16	-1498.0092	3094.5638	0.0000	1.0000
JC69	11	-1517.1132	3101.9765	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 591-tomodel-phylip/chr8_3277.phylip = 1 - 473;
 Partition Dummy = 1:591-tomodel-phylip/chr8_3277.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 591-tomodel-phylip/chr8_3277.phylip = 1 - 473;
 Partition Dummy = 1:591-tomodel-phylip/chr8_3277.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 591-tomodel-phylip/chr8_3277.phylip = 1 - 473;
 Partition Dummy = 1:591-tomodel-phylip/chr8_3277.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

