
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:51 02/17/2011

Input data from file "183-phylip/chr5_10047.phylip" (29 taxa, 605 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-2362.1716	4867.7506	0.7077	0.7077
GTRIG	65	-2362.0185	4869.9554	0.2350	0.9427
SYMG	61	-2369.0647	4874.0594	0.0302	0.9729
GTRI	64	-2366.1945	4875.7964	0.0127	0.9856
SYMIG	62	-2368.9112	4876.2358	0.0102	0.9957
K2PG	57	-2376.7641	4879.6159	0.0019	0.9976
HKYG	60	-2373.5004	4880.4566	0.0012	0.9988
SYMI	61	-2372.3755	4880.6810	0.0011	0.9999
K2PIG	58	-2379.6454	4887.8255	0.0000	1.0000
HKYIG	61	-2376.3265	4888.5831	0.0000	1.0000
K2PI	57	-2382.6505	4891.3887	0.0000	1.0000
HKYI	60	-2379.4329	4892.3217	0.0000	1.0000
GTR	63	-2395.2007	4931.3071	0.0000	1.0000
SYM	60	-2400.9681	4935.3921	0.0000	1.0000
K2P	56	-2411.3669	4946.3834	0.0000	1.0000
HKY	59	-2408.1805	4947.3518	0.0000	1.0000
JC69G	56	-2426.3414	4976.3324	0.0000	1.0000
F81G	59	-2426.3398	4983.6704	0.0000	1.0000
F81IG	60	-2426.1950	4985.8459	0.0000	1.0000
JC69IG	57	-2429.9312	4985.9501	0.0000	1.0000
F81I	59	-2428.8041	4988.5990	0.0000	1.0000
JC69I	56	-2432.5532	4988.7561	0.0000	1.0000
F81	58	-2457.2816	5043.0981	0.0000	1.0000
JC69	55	-2461.0793	5043.3791	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2362.1716	4852.3432	0.6767	0.6767
GTRIG	65	-2362.0185	4854.0370	0.2901	0.9669
SYMG	61	-2369.0647	4860.1294	0.0138	0.9807
GTRI	64	-2366.1945	4860.3890	0.0121	0.9928
SYMIG	62	-2368.9112	4861.8225	0.0059	0.9987
SYMI	61	-2372.3755	4866.7510	0.0005	0.9992
HKYG	60	-2373.5004	4867.0007	0.0004	0.9996
K2PG	57	-2376.7641	4867.5281	0.0003	1.0000
HKYIG	61	-2376.3265	4874.6531	0.0000	1.0000
K2PIG	58	-2379.6454	4875.2907	0.0000	1.0000
HKYI	60	-2379.4329	4878.8658	0.0000	1.0000
K2PI	57	-2382.6505	4879.3009	0.0000	1.0000
GTR	63	-2395.2007	4916.4014	0.0000	1.0000
SYM	60	-2400.9681	4921.9363	0.0000	1.0000
HKY	59	-2408.1805	4934.3610	0.0000	1.0000
K2P	56	-2411.3669	4934.7337	0.0000	1.0000
JC69G	56	-2426.3414	4964.6828	0.0000	1.0000
F81G	59	-2426.3398	4970.6796	0.0000	1.0000
F81IG	60	-2426.1950	4972.3900	0.0000	1.0000
JC69IG	57	-2429.9312	4973.8624	0.0000	1.0000
F81I	59	-2428.8041	4975.6082	0.0000	1.0000
JC69I	56	-2432.5532	4977.1065	0.0000	1.0000
F81	58	-2457.2816	5030.5633	0.0000	1.0000
JC69	55	-2461.0793	5032.1587	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2376.7641	5118.6262	0.9864	0.9864
SYMG	61	-2369.0647	5128.8483	0.0059	0.9924
K2PI	57	-2382.6505	5130.3989	0.0027	0.9951
K2PIG	58	-2379.6454	5130.7940	0.0022	0.9973
HKYG	60	-2373.5004	5131.3144	0.0017	0.9991
GTRG	64	-2362.1716	5134.2778	0.0004	0.9995
SYMIG	62	-2368.9112	5134.9466	0.0003	0.9997
SYMI	61	-2372.3755	5135.4699	0.0002	1.0000
GTRIG	65	-2362.0185	5140.3769	0.0000	1.0000
GTRI	64	-2366.1945	5142.3236	0.0000	1.0000
HKYI	60	-2379.4329	5143.1795	0.0000	1.0000
HKYIG	61	-2376.3265	5143.3720	0.0000	1.0000
K2P	56	-2411.3669	5181.4265	0.0000	1.0000
SYM	60	-2400.9681	5186.2500	0.0000	1.0000
GTR	63	-2395.2007	5193.9308	0.0000	1.0000
HKY	59	-2408.1805	5194.2694	0.0000	1.0000
JC69G	56	-2426.3414	5211.3756	0.0000	1.0000
JC69I	56	-2432.5532	5223.7993	0.0000	1.0000
JC69IG	57	-2429.9312	5224.9604	0.0000	1.0000
F81G	59	-2426.3398	5230.5881	0.0000	1.0000
F81I	59	-2428.8041	5235.5167	0.0000	1.0000
F81IG	60	-2426.1950	5236.7037	0.0000	1.0000
JC69	55	-2461.0793	5274.4462	0.0000	1.0000
F81	58	-2457.2816	5286.0665	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 183-phylip/chr5_10047.phylip = 1 - 605;
 Partition Dummy = 1:183-phylip/chr5_10047.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 183-phylip/chr5_10047.phylip = 1 - 605;
 Partition Dummy = 1:183-phylip/chr5_10047.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 (K2PG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr5_10047.phylip = 1 - 605;
 Partition Dummy = 1:183-phylip/chr5_10047.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

