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

Input data from file "591-tomodel-phylip/chr1_32365.phylip" (7 taxa, 963 characters)

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2536.6078	5113.2157	0.4581	0.4581
GTRG	20	-2536.7536	5113.5073	0.3960	0.8541
GTRIG	21	-2536.7536	5115.5073	0.1457	0.9998
HKYI	16	-2549.1453	5130.2907	0.0001	0.9999
HKYG	16	-2549.2732	5130.5465	0.0001	1.0000
HKYIG	17	-2549.2733	5132.5465	0.0000	1.0000
GTR	19	-2568.1823	5174.3647	0.0000	1.0000
HKY	15	-2580.8409	5191.6818	0.0000	1.0000
SYMI	17	-2579.8100	5193.6200	0.0000	1.0000
SYMG	17	-2579.9485	5193.8969	0.0000	1.0000
SYMIG	18	-2579.9484	5195.8969	0.0000	1.0000
K2PI	13	-2586.2922	5198.5845	0.0000	1.0000
K2PG	13	-2586.3472	5198.6944	0.0000	1.0000
K2PIG	14	-2586.3472	5200.6944	0.0000	1.0000
F81I	15	-2594.6559	5219.3118	0.0000	1.0000
F81G	15	-2594.9382	5219.8764	0.0000	1.0000
F81IG	16	-2594.9383	5221.8765	0.0000	1.0000
SYM	16	-2611.5020	5255.0040	0.0000	1.0000
K2P	12	-2618.7518	5261.5036	0.0000	1.0000
F81	14	-2624.0908	5276.1815	0.0000	1.0000
JC69I	12	-2627.1955	5278.3910	0.0000	1.0000
JC69G	12	-2627.4248	5278.8496	0.0000	1.0000
JC69IG	13	-2627.4249	5280.8498	0.0000	1.0000
JC69	11	-2657.4360	5336.8719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2549.1453	5208.2115	0.4038	0.4038
HKYG	16	-2549.2732	5208.4673	0.3553	0.7590
GTRI	20	-2536.6078	5210.6167	0.1213	0.8803
GTRG	20	-2536.7536	5210.9083	0.1048	0.9852
HKYIG	17	-2549.2733	5215.3374	0.0114	0.9966
GTRIG	21	-2536.7536	5217.7784	0.0034	1.0000
K2PI	13	-2586.2922	5261.8952	0.0000	1.0000
K2PG	13	-2586.3472	5262.0051	0.0000	1.0000
HKY	15	-2580.8409	5264.7326	0.0000	1.0000
GTR	19	-2568.1823	5266.8957	0.0000	1.0000
K2PIG	14	-2586.3472	5268.8752	0.0000	1.0000
SYMI	17	-2579.8100	5276.4109	0.0000	1.0000
SYMG	17	-2579.9485	5276.6878	0.0000	1.0000
SYMIG	18	-2579.9484	5283.5578	0.0000	1.0000
F81I	15	-2594.6559	5292.3626	0.0000	1.0000
F81G	15	-2594.9382	5292.9272	0.0000	1.0000
F81IG	16	-2594.9383	5299.7974	0.0000	1.0000
K2P	12	-2618.7518	5319.9443	0.0000	1.0000
SYM	16	-2611.5020	5332.9248	0.0000	1.0000
JC69I	12	-2627.1955	5336.8316	0.0000	1.0000
JC69G	12	-2627.4248	5337.2902	0.0000	1.0000
JC69IG	13	-2627.4249	5344.1605	0.0000	1.0000
F81	14	-2624.0908	5344.3623	0.0000	1.0000
JC69	11	-2657.4360	5390.4425	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32365.phylip = 1 - 963;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32365.phylip = 1 - 963;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_32365.phylip = 1 - 963;
 Partition Dummy = 1:591-tomodel-phylip/chr1_32365.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


-------------------------------------------------------------
End of Output

