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

Input data from file "591-tomodel-phylip/chr1_28694.phylip" (7 taxa, 806 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	-2450.0687	4941.2074	0.5906	0.5906
GTRIG	21	-2449.5690	4942.3165	0.3392	0.9298
GTRI	20	-2452.6681	4946.4063	0.0439	0.9737
HKYG	16	-2457.8327	4948.3548	0.0166	0.9902
HKYIG	17	-2457.4449	4949.6664	0.0086	0.9988
HKYI	16	-2460.4812	4953.6519	0.0012	1.0000
SYMG	17	-2467.3563	4969.4893	0.0000	1.0000
SYMIG	18	-2468.0111	4972.8914	0.0000	1.0000
K2PG	13	-2474.3286	4975.1167	0.0000	1.0000
SYMI	17	-2470.5789	4975.9345	0.0000	1.0000
K2PIG	14	-2473.9850	4976.5010	0.0000	1.0000
K2PI	13	-2476.7658	4979.9912	0.0000	1.0000
GTR	19	-2489.3582	5017.6832	0.0000	1.0000
HKY	15	-2495.3851	5021.3777	0.0000	1.0000
SYM	16	-2507.4707	5047.6309	0.0000	1.0000
K2P	12	-2511.9331	5048.2597	0.0000	1.0000
F81G	15	-2514.8877	5060.3831	0.0000	1.0000
F81IG	16	-2514.6541	5061.9977	0.0000	1.0000
F81I	15	-2516.9632	5064.5339	0.0000	1.0000
JC69G	12	-2528.2137	5080.8208	0.0000	1.0000
JC69IG	13	-2528.0000	5082.4596	0.0000	1.0000
JC69I	12	-2530.1901	5084.7736	0.0000	1.0000
F81	14	-2550.7356	5130.0022	0.0000	1.0000
JC69	11	-2564.0530	5150.4385	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2450.0687	4940.1374	0.5820	0.5820
GTRIG	21	-2449.5690	4941.1379	0.3529	0.9350
GTRI	20	-2452.6681	4945.3363	0.0433	0.9782
HKYG	16	-2457.8327	4947.6653	0.0135	0.9917
HKYIG	17	-2457.4449	4948.8898	0.0073	0.9990
HKYI	16	-2460.4812	4952.9624	0.0010	1.0000
SYMG	17	-2467.3563	4968.7127	0.0000	1.0000
SYMIG	18	-2468.0111	4972.0223	0.0000	1.0000
K2PG	13	-2474.3286	4974.6571	0.0000	1.0000
SYMI	17	-2470.5789	4975.1579	0.0000	1.0000
K2PIG	14	-2473.9850	4975.9700	0.0000	1.0000
K2PI	13	-2476.7658	4979.5316	0.0000	1.0000
GTR	19	-2489.3582	5016.7163	0.0000	1.0000
HKY	15	-2495.3851	5020.7701	0.0000	1.0000
SYM	16	-2507.4707	5046.9414	0.0000	1.0000
K2P	12	-2511.9331	5047.8662	0.0000	1.0000
F81G	15	-2514.8877	5059.7755	0.0000	1.0000
F81IG	16	-2514.6541	5061.3083	0.0000	1.0000
F81I	15	-2516.9632	5063.9263	0.0000	1.0000
JC69G	12	-2528.2137	5080.4273	0.0000	1.0000
JC69IG	13	-2528.0000	5082.0000	0.0000	1.0000
JC69I	12	-2530.1901	5084.3801	0.0000	1.0000
F81	14	-2550.7356	5129.4713	0.0000	1.0000
JC69	11	-2564.0530	5150.1060	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2457.8327	5022.7387	0.8861	0.8861
HKYI	16	-2460.4812	5028.0358	0.0627	0.9488
HKYIG	17	-2457.4449	5028.6552	0.0460	0.9948
GTRG	20	-2450.0687	5033.9791	0.0032	0.9980
K2PG	13	-2474.3286	5035.6542	0.0014	0.9994
GTRI	20	-2452.6681	5039.1779	0.0002	0.9996
GTRIG	21	-2449.5690	5039.6717	0.0002	0.9998
K2PI	13	-2476.7658	5040.5286	0.0001	0.9999
K2PIG	14	-2473.9850	5041.6592	0.0001	1.0000
SYMG	17	-2467.3563	5048.4781	0.0000	1.0000
SYMI	17	-2470.5789	5054.9233	0.0000	1.0000
SYMIG	18	-2468.0111	5056.4798	0.0000	1.0000
HKY	15	-2495.3851	5091.1514	0.0000	1.0000
K2P	12	-2511.9331	5104.1712	0.0000	1.0000
GTR	19	-2489.3582	5105.8659	0.0000	1.0000
SYM	16	-2507.4707	5122.0148	0.0000	1.0000
F81G	15	-2514.8877	5130.1567	0.0000	1.0000
F81I	15	-2516.9632	5134.3076	0.0000	1.0000
F81IG	16	-2514.6541	5136.3816	0.0000	1.0000
JC69G	12	-2528.2137	5136.7323	0.0000	1.0000
JC69I	12	-2530.1901	5140.6851	0.0000	1.0000
JC69IG	13	-2528.0000	5142.9971	0.0000	1.0000
F81	14	-2550.7356	5195.1604	0.0000	1.0000
JC69	11	-2564.0530	5201.7189	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/chr1_28694.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28694.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/chr1_28694.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28694.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/chr1_28694.phylip = 1 - 806;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28694.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

