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

Input data from file "591-tomodel-phylip/chrZ_7819.phylip" (7 taxa, 803 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	20	-2356.2311	4753.5363	0.4995	0.4995
GTRIG	21	-2356.1138	4755.4107	0.1957	0.6952
GTRI	20	-2357.5830	4756.2401	0.1292	0.8245
HKYG	16	-2361.9684	4756.6290	0.1064	0.9309
HKYIG	17	-2361.8740	4758.5277	0.0412	0.9721
HKYI	16	-2363.3055	4759.3031	0.0279	1.0000
SYMG	17	-2373.1013	4780.9823	0.0000	1.0000
SYMIG	18	-2373.0096	4782.8917	0.0000	1.0000
SYMI	17	-2374.3302	4783.4400	0.0000	1.0000
K2PG	13	-2379.6389	4785.7391	0.0000	1.0000
K2PIG	14	-2379.5442	4787.6214	0.0000	1.0000
K2PI	13	-2380.9605	4788.3822	0.0000	1.0000
GTR	19	-2381.6856	4802.3418	0.0000	1.0000
HKY	15	-2387.9091	4806.4280	0.0000	1.0000
SYM	16	-2398.9305	4830.5531	0.0000	1.0000
K2P	12	-2405.9399	4836.2747	0.0000	1.0000
F81G	15	-2407.6901	4845.9900	0.0000	1.0000
F81IG	16	-2407.5905	4847.8731	0.0000	1.0000
F81I	15	-2409.0463	4848.7025	0.0000	1.0000
JC69G	12	-2424.5702	4873.5354	0.0000	1.0000
JC69I	12	-2425.8075	4876.0100	0.0000	1.0000
JC69IG	13	-2427.4794	4881.4201	0.0000	1.0000
F81	14	-2432.8355	4894.2040	0.0000	1.0000
JC69	11	-2449.6152	4921.5642	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2356.2311	4752.4621	0.5087	0.5087
GTRIG	21	-2356.1138	4754.2276	0.2104	0.7191
GTRI	20	-2357.5830	4755.1660	0.1316	0.8508
HKYG	16	-2361.9684	4755.9369	0.0895	0.9403
HKYIG	17	-2361.8740	4757.7481	0.0362	0.9765
HKYI	16	-2363.3055	4758.6109	0.0235	1.0000
SYMG	17	-2373.1013	4780.2027	0.0000	1.0000
SYMIG	18	-2373.0096	4782.0193	0.0000	1.0000
SYMI	17	-2374.3302	4782.6604	0.0000	1.0000
K2PG	13	-2379.6389	4785.2777	0.0000	1.0000
K2PIG	14	-2379.5442	4787.0884	0.0000	1.0000
K2PI	13	-2380.9605	4787.9209	0.0000	1.0000
GTR	19	-2381.6856	4801.3712	0.0000	1.0000
HKY	15	-2387.9091	4805.8181	0.0000	1.0000
SYM	16	-2398.9305	4829.8610	0.0000	1.0000
K2P	12	-2405.9399	4835.8798	0.0000	1.0000
F81G	15	-2407.6901	4845.3801	0.0000	1.0000
F81IG	16	-2407.5905	4847.1810	0.0000	1.0000
F81I	15	-2409.0463	4848.0926	0.0000	1.0000
JC69G	12	-2424.5702	4873.1405	0.0000	1.0000
JC69I	12	-2425.8075	4875.6150	0.0000	1.0000
JC69IG	13	-2427.4794	4880.9587	0.0000	1.0000
F81	14	-2432.8355	4893.6710	0.0000	1.0000
JC69	11	-2449.6152	4921.2304	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2361.9684	4830.9505	0.7677	0.7677
HKYI	16	-2363.3055	4833.6246	0.2016	0.9693
HKYIG	17	-2361.8740	4837.4501	0.0298	0.9990
K2PG	13	-2379.6389	4846.2264	0.0004	0.9994
GTRG	20	-2356.2311	4846.2292	0.0004	0.9998
K2PI	13	-2380.9605	4848.8695	0.0001	0.9999
GTRI	20	-2357.5830	4848.9331	0.0001	1.0000
GTRIG	21	-2356.1138	4852.6831	0.0000	1.0000
K2PIG	14	-2379.5442	4852.7254	0.0000	1.0000
SYMG	17	-2373.1013	4859.9047	0.0000	1.0000
SYMI	17	-2374.3302	4862.3625	0.0000	1.0000
SYMIG	18	-2373.0096	4866.4096	0.0000	1.0000
HKY	15	-2387.9091	4876.1434	0.0000	1.0000
GTR	19	-2381.6856	4890.4499	0.0000	1.0000
K2P	12	-2405.9399	4892.1400	0.0000	1.0000
SYM	16	-2398.9305	4904.8747	0.0000	1.0000
F81G	15	-2407.6901	4915.7054	0.0000	1.0000
F81I	15	-2409.0463	4918.4179	0.0000	1.0000
F81IG	16	-2407.5905	4922.1947	0.0000	1.0000
JC69G	12	-2424.5702	4929.4007	0.0000	1.0000
JC69I	12	-2425.8075	4931.8753	0.0000	1.0000
JC69IG	13	-2427.4794	4941.9074	0.0000	1.0000
F81	14	-2432.8355	4959.3080	0.0000	1.0000
JC69	11	-2449.6152	4972.8023	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 591-tomodel-phylip/chrZ_7819.phylip = 1 - 803;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7819.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 591-tomodel-phylip/chrZ_7819.phylip = 1 - 803;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7819.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 591-tomodel-phylip/chrZ_7819.phylip = 1 - 803;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_7819.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

