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

Input data from file "591-tomodel-phylip/chrZ_4766.phylip" (7 taxa, 612 characters)

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
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
SYMG	17	-1446.4784	2927.9870	0.1756	0.1756
SYMI	17	-1446.6760	2928.3822	0.1441	0.3197
K2PG	13	-1451.0368	2928.6824	0.1240	0.4438
HKYG	16	-1448.0430	2929.0002	0.1058	0.5496
K2PI	13	-1451.2390	2929.0867	0.1013	0.6509
HKYI	16	-1448.2681	2929.4505	0.0845	0.7354
SYMIG	18	-1446.4035	2929.9605	0.0655	0.8009
GTRG	20	-1444.5342	2930.4897	0.0502	0.8511
K2PIG	14	-1450.9374	2930.5783	0.0481	0.8992
GTRI	20	-1444.7361	2930.8935	0.0411	0.9402
HKYIG	17	-1447.9350	2930.9004	0.0409	0.9811
GTRIG	21	-1444.4501	2932.4663	0.0187	0.9998
SYM	16	-1455.5186	2943.9515	0.0001	0.9999
K2P	12	-1460.1617	2944.8442	0.0000	0.9999
HKY	15	-1457.0520	2944.9093	0.0000	1.0000
GTR	19	-1453.3947	2946.0732	0.0000	1.0000
JC69G	12	-1465.2159	2954.9527	0.0000	1.0000
F81G	15	-1462.0875	2954.9804	0.0000	1.0000
JC69I	12	-1465.4767	2955.4742	0.0000	1.0000
F81I	15	-1462.3493	2955.5040	0.0000	1.0000
JC69IG	13	-1465.1439	2956.8965	0.0000	1.0000
F81IG	16	-1462.0002	2956.9146	0.0000	1.0000
JC69	11	-1473.2183	2968.8767	0.0000	1.0000
F81	14	-1470.1516	2969.0067	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-1446.4784	2926.9567	0.1813	0.1813
SYMI	17	-1446.6760	2927.3519	0.1488	0.3300
K2PG	13	-1451.0368	2928.0737	0.1037	0.4337
HKYG	16	-1448.0430	2928.0859	0.1031	0.5368
K2PI	13	-1451.2390	2928.4780	0.0847	0.6215
HKYI	16	-1448.2681	2928.5362	0.0823	0.7038
SYMIG	18	-1446.4035	2928.8071	0.0719	0.7757
GTRG	20	-1444.5342	2929.0684	0.0631	0.8387
GTRI	20	-1444.7361	2929.4722	0.0515	0.8902
HKYIG	17	-1447.9350	2929.8701	0.0422	0.9325
K2PIG	14	-1450.9374	2929.8748	0.0421	0.9746
GTRIG	21	-1444.4501	2930.9002	0.0252	0.9999
SYM	16	-1455.5186	2943.0372	0.0001	0.9999
HKY	15	-1457.0520	2944.1039	0.0000	0.9999
K2P	12	-1460.1617	2944.3234	0.0000	1.0000
GTR	19	-1453.3947	2944.7894	0.0000	1.0000
F81G	15	-1462.0875	2954.1750	0.0000	1.0000
JC69G	12	-1465.2159	2954.4318	0.0000	1.0000
F81I	15	-1462.3493	2954.6987	0.0000	1.0000
JC69I	12	-1465.4767	2954.9534	0.0000	1.0000
F81IG	16	-1462.0002	2956.0004	0.0000	1.0000
JC69IG	13	-1465.1439	2956.2878	0.0000	1.0000
F81	14	-1470.1516	2968.3032	0.0000	1.0000
JC69	11	-1473.2183	2968.4367	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1451.0368	2985.4912	0.5356	0.5356
K2PI	13	-1451.2390	2985.8956	0.4375	0.9731
K2PIG	14	-1450.9374	2991.7090	0.0239	0.9970
K2P	12	-1460.1617	2997.3241	0.0014	0.9984
HKYG	16	-1448.0430	2998.7536	0.0007	0.9991
HKYI	16	-1448.2681	2999.2039	0.0006	0.9997
SYMG	17	-1446.4784	3002.0412	0.0001	0.9998
SYMI	17	-1446.6760	3002.4363	0.0001	0.9999
HKYIG	17	-1447.9350	3004.9545	0.0000	1.0000
JC69G	12	-1465.2159	3007.4326	0.0000	1.0000
JC69I	12	-1465.4767	3007.9541	0.0000	1.0000
SYMIG	18	-1446.4035	3008.3083	0.0000	1.0000
HKY	15	-1457.0520	3010.3549	0.0000	1.0000
SYM	16	-1455.5186	3013.7049	0.0000	1.0000
JC69IG	13	-1465.1439	3013.7053	0.0000	1.0000
JC69	11	-1473.2183	3017.0207	0.0000	1.0000
GTRG	20	-1444.5342	3017.4030	0.0000	1.0000
GTRI	20	-1444.7361	3017.8068	0.0000	1.0000
F81G	15	-1462.0875	3020.4260	0.0000	1.0000
F81I	15	-1462.3493	3020.9497	0.0000	1.0000
GTRIG	21	-1444.4501	3023.6516	0.0000	1.0000
F81IG	16	-1462.0002	3026.6681	0.0000	1.0000
GTR	19	-1453.3947	3028.7073	0.0000	1.0000
F81	14	-1470.1516	3030.1374	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4766.phylip = 1 - 612;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4766.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4766.phylip = 1 - 612;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4766.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=Fixed(Equal);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_4766.phylip = 1 - 612;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_4766.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

