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

Input data from file "591-tomodel-phylip/chr23_707.phylip" (7 taxa, 755 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	-2328.4986	4698.1417	0.4607	0.4607
GTRI	20	-2328.8987	4698.9417	0.3088	0.7695
GTRIG	21	-2328.5039	4700.2684	0.1591	0.9286
HKYG	16	-2335.0681	4702.8733	0.0432	0.9718
HKYIG	17	-2335.0776	4704.9856	0.0150	0.9869
HKYI	16	-2336.2599	4705.2568	0.0131	1.0000
GTR	19	-2348.4506	4735.9352	0.0000	1.0000
HKY	15	-2354.8414	4740.3323	0.0000	1.0000
F81G	15	-2380.7061	4792.0617	0.0000	1.0000
F81I	15	-2381.6299	4793.9093	0.0000	1.0000
F81IG	16	-2380.7121	4794.1613	0.0000	1.0000
SYMG	17	-2384.7453	4804.3210	0.0000	1.0000
SYMI	17	-2385.5122	4805.8548	0.0000	1.0000
SYMIG	18	-2384.6933	4806.3159	0.0000	1.0000
F81	14	-2398.6728	4825.9132	0.0000	1.0000
K2PG	13	-2402.2351	4830.9615	0.0000	1.0000
K2PIG	14	-2402.2197	4833.0069	0.0000	1.0000
K2PI	13	-2403.3149	4833.1211	0.0000	1.0000
SYM	16	-2406.7922	4846.3215	0.0000	1.0000
K2P	12	-2422.1885	4868.7974	0.0000	1.0000
JC69G	12	-2438.2475	4900.9155	0.0000	1.0000
JC69IG	13	-2438.2433	4902.9777	0.0000	1.0000
JC69I	12	-2439.2964	4903.0133	0.0000	1.0000
JC69	11	-2456.9163	4936.1879	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2328.4986	4696.9972	0.4621	0.4621
GTRI	20	-2328.8987	4697.7973	0.3098	0.7719
GTRIG	21	-2328.5039	4699.0078	0.1691	0.9410
HKYG	16	-2335.0681	4702.1362	0.0354	0.9764
HKYIG	17	-2335.0776	4704.1552	0.0129	0.9893
HKYI	16	-2336.2599	4704.5197	0.0107	1.0000
GTR	19	-2348.4506	4734.9011	0.0000	1.0000
HKY	15	-2354.8414	4739.6828	0.0000	1.0000
F81G	15	-2380.7061	4791.4121	0.0000	1.0000
F81I	15	-2381.6299	4793.2598	0.0000	1.0000
F81IG	16	-2380.7121	4793.4242	0.0000	1.0000
SYMG	17	-2384.7453	4803.4906	0.0000	1.0000
SYMI	17	-2385.5122	4805.0244	0.0000	1.0000
SYMIG	18	-2384.6933	4805.3866	0.0000	1.0000
F81	14	-2398.6728	4825.3456	0.0000	1.0000
K2PG	13	-2402.2351	4830.4703	0.0000	1.0000
K2PIG	14	-2402.2197	4832.4393	0.0000	1.0000
K2PI	13	-2403.3149	4832.6299	0.0000	1.0000
SYM	16	-2406.7922	4845.5844	0.0000	1.0000
K2P	12	-2422.1885	4868.3770	0.0000	1.0000
JC69G	12	-2438.2475	4900.4950	0.0000	1.0000
JC69IG	13	-2438.2433	4902.4865	0.0000	1.0000
JC69I	12	-2439.2964	4902.5929	0.0000	1.0000
JC69	11	-2456.9163	4935.8326	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2335.0681	4776.1636	0.7452	0.7452
HKYI	16	-2336.2599	4778.5472	0.2263	0.9715
HKYIG	17	-2335.0776	4782.8094	0.0269	0.9984
GTRG	20	-2328.4986	4789.5316	0.0009	0.9993
GTRI	20	-2328.8987	4790.3317	0.0006	1.0000
GTRIG	21	-2328.5039	4796.1689	0.0000	1.0000
HKY	15	-2354.8414	4809.0836	0.0000	1.0000
GTR	19	-2348.4506	4822.8088	0.0000	1.0000
F81G	15	-2380.7061	4860.8129	0.0000	1.0000
F81I	15	-2381.6299	4862.6606	0.0000	1.0000
F81IG	16	-2380.7121	4867.4516	0.0000	1.0000
SYMG	17	-2384.7453	4882.1448	0.0000	1.0000
SYMI	17	-2385.5122	4883.6786	0.0000	1.0000
SYMIG	18	-2384.6933	4888.6675	0.0000	1.0000
F81	14	-2398.6728	4890.1197	0.0000	1.0000
K2PG	13	-2402.2351	4890.6176	0.0000	1.0000
K2PI	13	-2403.3149	4892.7772	0.0000	1.0000
K2PIG	14	-2402.2197	4897.2133	0.0000	1.0000
SYM	16	-2406.7922	4919.6119	0.0000	1.0000
K2P	12	-2422.1885	4923.8976	0.0000	1.0000
JC69G	12	-2438.2475	4956.0157	0.0000	1.0000
JC69I	12	-2439.2964	4958.1135	0.0000	1.0000
JC69IG	13	-2438.2433	4962.6339	0.0000	1.0000
JC69	11	-2456.9163	4986.7265	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/chr23_707.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr23_707.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/chr23_707.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr23_707.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/chr23_707.phylip = 1 - 755;
 Partition Dummy = 1:591-tomodel-phylip/chr23_707.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

