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

Input data from file "591-tomodel-phylip/chr2_24880.phylip" (7 taxa, 680 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	-2283.6895	4608.6536	0.3013	0.3013
HKYG	16	-2288.0789	4608.9782	0.2562	0.5575
GTRI	20	-2284.5249	4610.3244	0.1307	0.6882
GTRIG	21	-2283.5914	4610.5872	0.1146	0.8028
HKYI	16	-2288.9948	4610.8102	0.1025	0.9053
HKYIG	17	-2288.0221	4610.9686	0.0947	1.0000
HKY	15	-2307.9112	4646.5454	0.0000	1.0000
GTR	19	-2304.3381	4647.8278	0.0000	1.0000
SYMG	17	-2320.2164	4675.3573	0.0000	1.0000
K2PG	13	-2325.1814	4676.9094	0.0000	1.0000
SYMI	17	-2321.0950	4677.1144	0.0000	1.0000
SYMIG	18	-2320.2195	4677.4737	0.0000	1.0000
K2PI	13	-2326.0612	4678.6690	0.0000	1.0000
K2PIG	14	-2325.1723	4678.9762	0.0000	1.0000
SYM	16	-2339.9489	4712.7183	0.0000	1.0000
K2P	12	-2344.3149	4713.0976	0.0000	1.0000
F81G	15	-2341.3328	4713.3886	0.0000	1.0000
F81I	15	-2342.0556	4714.8340	0.0000	1.0000
F81IG	16	-2341.3026	4715.4257	0.0000	1.0000
F81	14	-2358.7304	4746.0923	0.0000	1.0000
JC69G	12	-2371.7901	4768.0479	0.0000	1.0000
JC69I	12	-2372.6478	4769.7634	0.0000	1.0000
JC69IG	13	-2371.8030	4770.1525	0.0000	1.0000
JC69	11	-2388.9534	4800.3020	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2283.6895	4607.3789	0.3277	0.3277
HKYG	16	-2288.0789	4608.1577	0.2220	0.5497
GTRI	20	-2284.5249	4609.0498	0.1421	0.6918
GTRIG	21	-2283.5914	4609.1829	0.1330	0.8247
HKYI	16	-2288.9948	4609.9896	0.0888	0.9136
HKYIG	17	-2288.0221	4610.0442	0.0864	1.0000
HKY	15	-2307.9112	4645.8225	0.0000	1.0000
GTR	19	-2304.3381	4646.6763	0.0000	1.0000
SYMG	17	-2320.2164	4674.4328	0.0000	1.0000
SYMI	17	-2321.0950	4676.1899	0.0000	1.0000
K2PG	13	-2325.1814	4676.3629	0.0000	1.0000
SYMIG	18	-2320.2195	4676.4389	0.0000	1.0000
K2PI	13	-2326.0612	4678.1224	0.0000	1.0000
K2PIG	14	-2325.1723	4678.3446	0.0000	1.0000
SYM	16	-2339.9489	4711.8978	0.0000	1.0000
K2P	12	-2344.3149	4712.6298	0.0000	1.0000
F81G	15	-2341.3328	4712.6657	0.0000	1.0000
F81I	15	-2342.0556	4714.1111	0.0000	1.0000
F81IG	16	-2341.3026	4714.6051	0.0000	1.0000
F81	14	-2358.7304	4745.4608	0.0000	1.0000
JC69G	12	-2371.7901	4767.5802	0.0000	1.0000
JC69I	12	-2372.6478	4769.2956	0.0000	1.0000
JC69IG	13	-2371.8030	4769.6059	0.0000	1.0000
JC69	11	-2388.9534	4799.9068	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2288.0789	4680.5112	0.6940	0.6940
HKYI	16	-2288.9948	4682.3431	0.2777	0.9717
HKYIG	17	-2288.0221	4686.9197	0.0282	0.9998
GTRG	20	-2283.6895	4697.8208	0.0001	0.9999
GTRI	20	-2284.5249	4699.4916	0.0001	1.0000
GTRIG	21	-2283.5914	4704.1468	0.0000	1.0000
HKY	15	-2307.9112	4713.6539	0.0000	1.0000
GTR	19	-2304.3381	4732.5960	0.0000	1.0000
K2PG	13	-2325.1814	4735.1501	0.0000	1.0000
K2PI	13	-2326.0612	4736.9096	0.0000	1.0000
K2PIG	14	-2325.1723	4741.6539	0.0000	1.0000
SYMG	17	-2320.2164	4751.3084	0.0000	1.0000
SYMI	17	-2321.0950	4753.0655	0.0000	1.0000
SYMIG	18	-2320.2195	4757.8366	0.0000	1.0000
K2P	12	-2344.3149	4766.8949	0.0000	1.0000
F81G	15	-2341.3328	4780.4971	0.0000	1.0000
F81I	15	-2342.0556	4781.9425	0.0000	1.0000
SYM	16	-2339.9489	4784.2512	0.0000	1.0000
F81IG	16	-2341.3026	4786.9586	0.0000	1.0000
F81	14	-2358.7304	4808.7701	0.0000	1.0000
JC69G	12	-2371.7901	4821.8453	0.0000	1.0000
JC69I	12	-2372.6478	4823.5607	0.0000	1.0000
JC69IG	13	-2371.8030	4828.3931	0.0000	1.0000
JC69	11	-2388.9534	4849.6498	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_24880.phylip = 1 - 680;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24880.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_24880.phylip = 1 - 680;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24880.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_24880.phylip = 1 - 680;
 Partition Dummy = 1:591-tomodel-phylip/chr2_24880.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

