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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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
GTRI	20	-2530.3203	5101.9190	0.6488	0.6488
GTRIG	21	-2530.3477	5104.1038	0.2176	0.8664
GTRG	20	-2531.9949	5105.2683	0.1216	0.9880
HKYI	16	-2538.9646	5110.7521	0.0078	0.9958
HKYIG	17	-2538.9933	5112.9138	0.0027	0.9985
HKYG	16	-2540.6087	5114.0405	0.0015	1.0000
GTR	19	-2544.2827	5127.7203	0.0000	1.0000
HKY	15	-2553.0235	5136.7720	0.0000	1.0000
SYMI	17	-2552.9195	5140.7662	0.0000	1.0000
SYMIG	18	-2552.9465	5142.9310	0.0000	1.0000
K2PI	13	-2558.4714	5143.4910	0.0000	1.0000
SYMG	17	-2554.4686	5143.8645	0.0000	1.0000
K2PIG	14	-2558.5040	5145.6415	0.0000	1.0000
K2PG	13	-2560.2536	5147.0553	0.0000	1.0000
SYM	16	-2566.1493	5165.1215	0.0000	1.0000
K2P	12	-2573.0176	5170.5045	0.0000	1.0000
F81I	15	-2590.2310	5211.1870	0.0000	1.0000
F81IG	16	-2590.2661	5213.3551	0.0000	1.0000
F81G	15	-2592.0629	5214.8509	0.0000	1.0000
F81	14	-2603.3571	5235.3477	0.0000	1.0000
JC69I	12	-2606.4525	5237.3742	0.0000	1.0000
JC69IG	13	-2606.4892	5239.5266	0.0000	1.0000
JC69G	12	-2608.3302	5241.1296	0.0000	1.0000
JC69	11	-2619.5369	5261.4702	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2530.3203	5100.6405	0.6409	0.6409
GTRIG	21	-2530.3477	5102.6953	0.2294	0.8703
GTRG	20	-2531.9949	5103.9898	0.1201	0.9904
HKYI	16	-2538.9646	5109.9291	0.0062	0.9966
HKYIG	17	-2538.9933	5111.9865	0.0022	0.9988
HKYG	16	-2540.6087	5113.2175	0.0012	1.0000
GTR	19	-2544.2827	5126.5653	0.0000	1.0000
HKY	15	-2553.0235	5136.0470	0.0000	1.0000
SYMI	17	-2552.9195	5139.8390	0.0000	1.0000
SYMIG	18	-2552.9465	5141.8930	0.0000	1.0000
SYMG	17	-2554.4686	5142.9372	0.0000	1.0000
K2PI	13	-2558.4714	5142.9428	0.0000	1.0000
K2PIG	14	-2558.5040	5145.0080	0.0000	1.0000
K2PG	13	-2560.2536	5146.5071	0.0000	1.0000
SYM	16	-2566.1493	5164.2985	0.0000	1.0000
K2P	12	-2573.0176	5170.0353	0.0000	1.0000
F81I	15	-2590.2310	5210.4619	0.0000	1.0000
F81IG	16	-2590.2661	5212.5321	0.0000	1.0000
F81G	15	-2592.0629	5214.1258	0.0000	1.0000
F81	14	-2603.3571	5234.7142	0.0000	1.0000
JC69I	12	-2606.4525	5236.9050	0.0000	1.0000
JC69IG	13	-2606.4892	5238.9784	0.0000	1.0000
JC69G	12	-2608.3302	5240.6605	0.0000	1.0000
JC69	11	-2619.5369	5261.0738	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2538.9646	5182.2355	0.8028	0.8028
HKYG	16	-2540.6087	5185.5238	0.1551	0.9578
HKYIG	17	-2538.9933	5188.8120	0.0300	0.9878
GTRI	20	-2530.3203	5191.0234	0.0099	0.9977
GTRG	20	-2531.9949	5194.3727	0.0019	0.9996
GTRIG	21	-2530.3477	5197.5974	0.0004	0.9999
K2PI	13	-2558.4714	5201.6917	0.0000	1.0000
HKY	15	-2553.0235	5203.8342	0.0000	1.0000
K2PG	13	-2560.2536	5205.2560	0.0000	1.0000
K2PIG	14	-2558.5040	5208.2760	0.0000	1.0000
GTR	19	-2544.2827	5212.4291	0.0000	1.0000
SYMI	17	-2552.9195	5216.6645	0.0000	1.0000
SYMG	17	-2554.4686	5219.7627	0.0000	1.0000
SYMIG	18	-2552.9465	5223.2377	0.0000	1.0000
K2P	12	-2573.0176	5224.2650	0.0000	1.0000
SYM	16	-2566.1493	5236.6049	0.0000	1.0000
F81I	15	-2590.2310	5278.2491	0.0000	1.0000
F81G	15	-2592.0629	5281.9130	0.0000	1.0000
F81IG	16	-2590.2661	5284.8385	0.0000	1.0000
JC69I	12	-2606.4525	5291.1348	0.0000	1.0000
JC69G	12	-2608.3302	5294.8902	0.0000	1.0000
JC69IG	13	-2606.4892	5297.7273	0.0000	1.0000
F81	14	-2603.3571	5297.9823	0.0000	1.0000
JC69	11	-2619.5369	5310.7844	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_11557.phylip = 1 - 678;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11557.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_11557.phylip = 1 - 678;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11557.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chrZ_11557.phylip = 1 - 678;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_11557.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

