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

Input data from file "591-tomodel-phylip/chr2_18614.phylip" (7 taxa, 815 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	-2415.8298	4872.7176	0.5679	0.5679
GTRIG	21	-2415.4178	4874.0008	0.2990	0.8668
GTRI	20	-2417.3813	4875.8205	0.1204	0.9872
HKYG	16	-2424.3990	4881.4797	0.0071	0.9943
HKYIG	17	-2424.1408	4883.0495	0.0032	0.9975
HKYI	16	-2425.4496	4883.5809	0.0025	1.0000
SYMG	17	-2434.7284	4904.2247	0.0000	1.0000
K2PG	13	-2439.3414	4905.1372	0.0000	1.0000
SYMI	17	-2435.2946	4905.3571	0.0000	1.0000
SYMIG	18	-2434.5704	4906.0000	0.0000	1.0000
K2PI	13	-2439.9691	4906.3925	0.0000	1.0000
K2PIG	14	-2439.1789	4906.8829	0.0000	1.0000
F81G	15	-2455.9952	4942.5912	0.0000	1.0000
F81I	15	-2456.7435	4944.0877	0.0000	1.0000
F81IG	16	-2455.7680	4944.2176	0.0000	1.0000
JC69G	12	-2467.6440	4959.6769	0.0000	1.0000
JC69I	12	-2468.1219	4960.6328	0.0000	1.0000
JC69IG	13	-2467.4707	4961.3958	0.0000	1.0000
GTR	19	-2465.8440	4970.6440	0.0000	1.0000
HKY	15	-2472.9834	4976.5676	0.0000	1.0000
SYM	16	-2480.6842	4994.0502	0.0000	1.0000
K2P	12	-2485.7888	4995.9667	0.0000	1.0000
F81	14	-2501.9637	5032.4525	0.0000	1.0000
JC69	11	-2512.4101	5047.1490	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2415.8298	4871.6597	0.5598	0.5598
GTRIG	21	-2415.4178	4872.8356	0.3109	0.8708
GTRI	20	-2417.3813	4874.7626	0.1186	0.9894
HKYG	16	-2424.3990	4880.7980	0.0058	0.9952
HKYIG	17	-2424.1408	4882.2817	0.0028	0.9980
HKYI	16	-2425.4496	4882.8992	0.0020	1.0000
SYMG	17	-2434.7284	4903.4569	0.0000	1.0000
SYMI	17	-2435.2946	4904.5892	0.0000	1.0000
K2PG	13	-2439.3414	4904.6827	0.0000	1.0000
SYMIG	18	-2434.5704	4905.1407	0.0000	1.0000
K2PI	13	-2439.9691	4905.9381	0.0000	1.0000
K2PIG	14	-2439.1789	4906.3579	0.0000	1.0000
F81G	15	-2455.9952	4941.9904	0.0000	1.0000
F81I	15	-2456.7435	4943.4869	0.0000	1.0000
F81IG	16	-2455.7680	4943.5359	0.0000	1.0000
JC69G	12	-2467.6440	4959.2879	0.0000	1.0000
JC69I	12	-2468.1219	4960.2438	0.0000	1.0000
JC69IG	13	-2467.4707	4960.9414	0.0000	1.0000
GTR	19	-2465.8440	4969.6880	0.0000	1.0000
HKY	15	-2472.9834	4975.9668	0.0000	1.0000
SYM	16	-2480.6842	4993.3685	0.0000	1.0000
K2P	12	-2485.7888	4995.5777	0.0000	1.0000
F81	14	-2501.9637	5031.9275	0.0000	1.0000
JC69	11	-2512.4101	5046.8202	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2424.3990	4956.0490	0.7057	0.7057
HKYI	16	-2425.4496	4958.1502	0.2468	0.9525
HKYIG	17	-2424.1408	4962.2359	0.0320	0.9845
GTRG	20	-2415.8298	4965.7234	0.0056	0.9901
K2PG	13	-2439.3414	4965.8242	0.0053	0.9955
K2PI	13	-2439.9691	4967.0795	0.0028	0.9983
GTRI	20	-2417.3813	4968.8263	0.0012	0.9995
GTRIG	21	-2415.4178	4971.6026	0.0003	0.9998
K2PIG	14	-2439.1789	4972.2025	0.0002	1.0000
SYMG	17	-2434.7284	4983.4111	0.0000	1.0000
SYMI	17	-2435.2946	4984.5434	0.0000	1.0000
SYMIG	18	-2434.5704	4989.7981	0.0000	1.0000
F81G	15	-2455.9952	5012.5382	0.0000	1.0000
F81I	15	-2456.7435	5014.0347	0.0000	1.0000
JC69G	12	-2467.6440	5015.7262	0.0000	1.0000
JC69I	12	-2468.1219	5016.6820	0.0000	1.0000
F81IG	16	-2455.7680	5018.7869	0.0000	1.0000
JC69IG	13	-2467.4707	5022.0828	0.0000	1.0000
HKY	15	-2472.9834	5046.5146	0.0000	1.0000
K2P	12	-2485.7888	5052.0159	0.0000	1.0000
GTR	19	-2465.8440	5059.0486	0.0000	1.0000
SYM	16	-2480.6842	5068.6195	0.0000	1.0000
F81	14	-2501.9637	5097.7721	0.0000	1.0000
JC69	11	-2512.4101	5098.5553	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/chr2_18614.phylip = 1 - 815;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18614.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/chr2_18614.phylip = 1 - 815;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18614.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/chr2_18614.phylip = 1 - 815;
 Partition Dummy = 1:591-tomodel-phylip/chr2_18614.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

