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

Input data from file "183-phylip/chr1_5478.phylip" (29 taxa, 470 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	-2251.3177	4651.1787	0.9797	0.9797
GTRIG	65	-2254.2141	4659.6659	0.0141	0.9937
SYMG	61	-2260.4117	4661.3626	0.0060	0.9997
HKYG	60	-2264.9689	4667.8351	0.0002	1.0000
K2PG	57	-2272.0371	4674.1228	0.0000	1.0000
SYMIG	62	-2266.0051	4675.2044	0.0000	1.0000
HKYIG	61	-2267.3603	4675.2598	0.0000	1.0000
GTRI	64	-2265.6880	4679.9192	0.0000	1.0000
K2PIG	58	-2278.4620	4689.5761	0.0000	1.0000
SYMI	61	-2277.1757	4694.8907	0.0000	1.0000
HKYI	60	-2281.1781	4700.2535	0.0000	1.0000
K2PI	57	-2292.3082	4714.6649	0.0000	1.0000
F81G	59	-2344.5819	4824.4321	0.0000	1.0000
F81IG	60	-2348.3929	4834.6832	0.0000	1.0000
JC69G	56	-2353.7618	4834.9813	0.0000	1.0000
GTR	63	-2345.2621	4836.3862	0.0000	1.0000
JC69I	56	-2356.7902	4841.0380	0.0000	1.0000
JC69IG	57	-2356.6261	4843.3008	0.0000	1.0000
SYM	60	-2353.5577	4845.0127	0.0000	1.0000
F81I	59	-2355.8924	4847.0531	0.0000	1.0000
K2P	56	-2371.0288	4869.5153	0.0000	1.0000
HKY	59	-2370.9482	4877.1647	0.0000	1.0000
F81	58	-2448.3636	5029.3793	0.0000	1.0000
JC69	55	-2454.0775	5033.0343	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2251.3177	4630.6355	0.9779	0.9779
GTRIG	65	-2254.2141	4638.4282	0.0199	0.9977
SYMG	61	-2260.4117	4642.8234	0.0022	0.9999
HKYG	60	-2264.9689	4649.9378	0.0001	1.0000
SYMIG	62	-2266.0051	4656.0103	0.0000	1.0000
HKYIG	61	-2267.3603	4656.7206	0.0000	1.0000
K2PG	57	-2272.0371	4658.0742	0.0000	1.0000
GTRI	64	-2265.6880	4659.3760	0.0000	1.0000
K2PIG	58	-2278.4620	4672.9240	0.0000	1.0000
SYMI	61	-2277.1757	4676.3515	0.0000	1.0000
HKYI	60	-2281.1781	4682.3562	0.0000	1.0000
K2PI	57	-2292.3082	4698.6163	0.0000	1.0000
F81G	59	-2344.5819	4807.1638	0.0000	1.0000
GTR	63	-2345.2621	4816.5242	0.0000	1.0000
F81IG	60	-2348.3929	4816.7859	0.0000	1.0000
JC69G	56	-2353.7618	4819.5236	0.0000	1.0000
JC69I	56	-2356.7902	4825.5803	0.0000	1.0000
SYM	60	-2353.5577	4827.1154	0.0000	1.0000
JC69IG	57	-2356.6261	4827.2523	0.0000	1.0000
F81I	59	-2355.8924	4829.7848	0.0000	1.0000
K2P	56	-2371.0288	4854.0577	0.0000	1.0000
HKY	59	-2370.9482	4859.8964	0.0000	1.0000
F81	58	-2448.3636	5012.7273	0.0000	1.0000
JC69	55	-2454.0775	5018.1551	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2272.0371	4894.7800	0.4840	0.4840
SYMG	61	-2260.4117	4896.1401	0.2452	0.7292
GTRG	64	-2251.3177	4896.4104	0.2142	0.9434
HKYG	60	-2264.9689	4899.1018	0.0558	0.9991
GTRIG	65	-2254.2141	4908.3559	0.0005	0.9997
HKYIG	61	-2267.3603	4910.0373	0.0002	0.9999
SYMIG	62	-2266.0051	4913.4797	0.0000	1.0000
K2PIG	58	-2278.4620	4913.7825	0.0000	1.0000
GTRI	64	-2265.6880	4925.1509	0.0000	1.0000
SYMI	61	-2277.1757	4929.6682	0.0000	1.0000
HKYI	60	-2281.1781	4931.5201	0.0000	1.0000
K2PI	57	-2292.3082	4935.3221	0.0000	1.0000
JC69G	56	-2353.7618	5052.0767	0.0000	1.0000
F81G	59	-2344.5819	5052.1751	0.0000	1.0000
JC69I	56	-2356.7902	5058.1334	0.0000	1.0000
JC69IG	57	-2356.6261	5063.9580	0.0000	1.0000
F81IG	60	-2348.3929	5065.9498	0.0000	1.0000
F81I	59	-2355.8924	5074.7960	0.0000	1.0000
SYM	60	-2353.5577	5076.2794	0.0000	1.0000
GTR	63	-2345.2621	5078.1463	0.0000	1.0000
K2P	56	-2371.0288	5086.6107	0.0000	1.0000
HKY	59	-2370.9482	5104.9076	0.0000	1.0000
JC69	55	-2454.0775	5246.5554	0.0000	1.0000
F81	58	-2448.3636	5253.5858	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/chr1_5478.phylip = 1 - 470;
 Partition Dummy = 1:183-phylip/chr1_5478.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/chr1_5478.phylip = 1 - 470;
 Partition Dummy = 1:183-phylip/chr1_5478.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/chr1_5478.phylip = 1 - 470;
 Partition Dummy = 1:183-phylip/chr1_5478.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

