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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-3162.5601	6458.5022	0.7041	0.7041
HKYIG	61	-3162.2056	6460.2647	0.2917	0.9958
HKYI	60	-3167.6874	6468.7568	0.0042	1.0000
GTRG	64	-3167.5775	6478.4772	0.0000	1.0000
GTRIG	65	-3167.2142	6480.2586	0.0000	1.0000
GTRI	64	-3174.4979	6492.3180	0.0000	1.0000
HKY	59	-3207.3666	6545.6528	0.0000	1.0000
GTR	63	-3213.5678	6567.9591	0.0000	1.0000
K2PG	57	-3226.7114	6579.4447	0.0000	1.0000
K2PIG	58	-3226.4515	6581.3693	0.0000	1.0000
SYMG	61	-3225.4614	6586.7764	0.0000	1.0000
SYMIG	62	-3225.2155	6588.7649	0.0000	1.0000
K2PI	57	-3233.3116	6592.6451	0.0000	1.0000
SYMI	61	-3232.2485	6600.3505	0.0000	1.0000
K2P	56	-3275.1497	6673.8856	0.0000	1.0000
SYM	60	-3273.5850	6680.5521	0.0000	1.0000
F81G	59	-3282.0178	6694.9554	0.0000	1.0000
F81IG	60	-3281.7245	6696.8311	0.0000	1.0000
F81I	59	-3287.5742	6706.0681	0.0000	1.0000
JC69G	56	-3317.2383	6758.0628	0.0000	1.0000
JC69IG	57	-3316.9460	6759.9139	0.0000	1.0000
JC69I	56	-3323.1414	6769.8691	0.0000	1.0000
F81	58	-3326.1866	6780.8395	0.0000	1.0000
JC69	55	-3362.8930	6846.9454	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-3162.5601	6445.1201	0.6534	0.6534
HKYIG	61	-3162.2056	6446.4112	0.3426	0.9960
HKYI	60	-3167.6874	6455.3748	0.0039	0.9999
GTRG	64	-3167.5775	6463.1550	0.0001	1.0000
GTRIG	65	-3167.2142	6464.4284	0.0000	1.0000
GTRI	64	-3174.4979	6476.9958	0.0000	1.0000
HKY	59	-3207.3666	6532.7331	0.0000	1.0000
GTR	63	-3213.5678	6553.1356	0.0000	1.0000
K2PG	57	-3226.7114	6567.4229	0.0000	1.0000
K2PIG	58	-3226.4515	6568.9030	0.0000	1.0000
SYMG	61	-3225.4614	6572.9229	0.0000	1.0000
SYMIG	62	-3225.2155	6574.4309	0.0000	1.0000
K2PI	57	-3233.3116	6580.6233	0.0000	1.0000
SYMI	61	-3232.2485	6586.4971	0.0000	1.0000
K2P	56	-3275.1497	6662.2994	0.0000	1.0000
SYM	60	-3273.5850	6667.1700	0.0000	1.0000
F81G	59	-3282.0178	6682.0357	0.0000	1.0000
F81IG	60	-3281.7245	6683.4490	0.0000	1.0000
F81I	59	-3287.5742	6693.1484	0.0000	1.0000
JC69G	56	-3317.2383	6746.4766	0.0000	1.0000
JC69IG	57	-3316.9460	6747.8921	0.0000	1.0000
JC69I	56	-3323.1414	6758.2829	0.0000	1.0000
F81	58	-3326.1866	6768.3732	0.0000	1.0000
JC69	55	-3362.8930	6835.7860	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3162.5601	6709.7306	0.9401	0.9401
HKYIG	61	-3162.2056	6715.4319	0.0543	0.9944
HKYI	60	-3167.6874	6719.9853	0.0056	1.0000
GTRG	64	-3167.5775	6745.4062	0.0000	1.0000
GTRIG	65	-3167.2142	6751.0897	0.0000	1.0000
GTRI	64	-3174.4979	6759.2470	0.0000	1.0000
HKY	59	-3207.3666	6792.9334	0.0000	1.0000
K2PG	57	-3226.7114	6818.8028	0.0000	1.0000
K2PIG	58	-3226.4515	6824.6931	0.0000	1.0000
GTR	63	-3213.5678	6830.9766	0.0000	1.0000
K2PI	57	-3233.3116	6832.0032	0.0000	1.0000
SYMG	61	-3225.4614	6841.9436	0.0000	1.0000
SYMIG	62	-3225.2155	6847.8618	0.0000	1.0000
SYMI	61	-3232.2485	6855.5177	0.0000	1.0000
K2P	56	-3275.1497	6909.2692	0.0000	1.0000
SYM	60	-3273.5850	6931.7805	0.0000	1.0000
F81G	59	-3282.0178	6942.2360	0.0000	1.0000
F81IG	60	-3281.7245	6948.0595	0.0000	1.0000
F81I	59	-3287.5742	6953.3488	0.0000	1.0000
JC69G	56	-3317.2383	6993.4464	0.0000	1.0000
JC69IG	57	-3316.9460	6999.2720	0.0000	1.0000
JC69I	56	-3323.1414	7005.2527	0.0000	1.0000
F81	58	-3326.1866	7024.1634	0.0000	1.0000
JC69	55	-3362.8930	7078.3456	0.0000	1.0000



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

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

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


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

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


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

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

