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

Input data from file "183-phylip/chr9_3633.phylip" (29 taxa, 631 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	-3450.9855	7044.6707	0.7005	0.7005
GTRIG	65	-3450.7657	7046.7173	0.2518	0.9523
HKYG	60	-3458.9268	7050.6958	0.0344	0.9867
HKYIG	61	-3458.7256	7052.7447	0.0124	0.9991
SYMG	61	-3461.7998	7058.8930	0.0006	0.9996
SYMIG	62	-3461.6419	7061.0373	0.0002	0.9998
K2PG	57	-3468.1793	7061.8978	0.0001	1.0000
K2PIG	58	-3468.0276	7064.0203	0.0000	1.0000
GTRI	64	-3464.2724	7071.2444	0.0000	1.0000
HKYI	60	-3471.8047	7076.4514	0.0000	1.0000
SYMI	61	-3473.2804	7081.8544	0.0000	1.0000
K2PI	57	-3479.8254	7085.1901	0.0000	1.0000
GTR	63	-3520.8108	7181.8439	0.0000	1.0000
HKY	59	-3528.6893	7187.7779	0.0000	1.0000
SYM	60	-3529.9611	7192.7643	0.0000	1.0000
K2P	56	-3536.9515	7197.0250	0.0000	1.0000
F81G	59	-3544.2524	7218.9041	0.0000	1.0000
JC69G	56	-3548.7919	7220.7057	0.0000	1.0000
F81IG	60	-3544.1368	7221.1156	0.0000	1.0000
JC69IG	57	-3548.6851	7222.9095	0.0000	1.0000
F81I	59	-3557.3441	7245.0876	0.0000	1.0000
JC69I	56	-3561.7386	7246.5991	0.0000	1.0000
JC69	55	-3616.4135	7353.5401	0.0000	1.0000
F81	58	-3613.9383	7355.8416	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3450.9855	7029.9711	0.6724	0.6724
GTRIG	65	-3450.7657	7031.5314	0.3082	0.9806
HKYG	60	-3458.9268	7037.8537	0.0131	0.9937
HKYIG	61	-3458.7256	7039.4512	0.0059	0.9996
SYMG	61	-3461.7998	7045.5995	0.0003	0.9998
SYMIG	62	-3461.6419	7047.2838	0.0001	1.0000
K2PG	57	-3468.1793	7050.3585	0.0000	1.0000
K2PIG	58	-3468.0276	7052.0553	0.0000	1.0000
GTRI	64	-3464.2724	7056.5447	0.0000	1.0000
HKYI	60	-3471.8047	7063.6093	0.0000	1.0000
SYMI	61	-3473.2804	7068.5609	0.0000	1.0000
K2PI	57	-3479.8254	7073.6508	0.0000	1.0000
GTR	63	-3520.8108	7167.6217	0.0000	1.0000
HKY	59	-3528.6893	7175.3786	0.0000	1.0000
SYM	60	-3529.9611	7179.9222	0.0000	1.0000
K2P	56	-3536.9515	7185.9030	0.0000	1.0000
F81G	59	-3544.2524	7206.5048	0.0000	1.0000
F81IG	60	-3544.1368	7208.2735	0.0000	1.0000
JC69G	56	-3548.7919	7209.5838	0.0000	1.0000
JC69IG	57	-3548.6851	7211.3703	0.0000	1.0000
F81I	59	-3557.3441	7232.6883	0.0000	1.0000
JC69I	56	-3561.7386	7235.4771	0.0000	1.0000
JC69	55	-3616.4135	7342.8271	0.0000	1.0000
F81	58	-3613.9383	7343.8766	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3468.1793	7303.8550	0.5738	0.5738
HKYG	60	-3458.9268	7304.6921	0.3776	0.9513
K2PIG	58	-3468.0276	7309.9990	0.0266	0.9779
HKYIG	61	-3458.7256	7310.7368	0.0184	0.9963
GTRG	64	-3450.9855	7314.5986	0.0027	0.9990
SYMG	61	-3461.7998	7316.8852	0.0008	0.9998
GTRIG	65	-3450.7657	7320.6063	0.0001	1.0000
SYMIG	62	-3461.6419	7323.0167	0.0000	1.0000
K2PI	57	-3479.8254	7327.1473	0.0000	1.0000
HKYI	60	-3471.8047	7330.4477	0.0000	1.0000
SYMI	61	-3473.2804	7339.8465	0.0000	1.0000
GTRI	64	-3464.2724	7341.1723	0.0000	1.0000
K2P	56	-3536.9515	7434.9521	0.0000	1.0000
HKY	59	-3528.6893	7437.7696	0.0000	1.0000
SYM	60	-3529.9611	7446.7606	0.0000	1.0000
GTR	63	-3520.8108	7447.8019	0.0000	1.0000
JC69G	56	-3548.7919	7458.6329	0.0000	1.0000
JC69IG	57	-3548.6851	7464.8667	0.0000	1.0000
F81G	59	-3544.2524	7468.8959	0.0000	1.0000
F81IG	60	-3544.1368	7475.1119	0.0000	1.0000
JC69I	56	-3561.7386	7484.5263	0.0000	1.0000
F81I	59	-3557.3441	7495.0793	0.0000	1.0000
JC69	55	-3616.4135	7587.4289	0.0000	1.0000
F81	58	-3613.9383	7601.8203	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/chr9_3633.phylip = 1 - 631;
 Partition Dummy = 1:183-phylip/chr9_3633.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/chr9_3633.phylip = 1 - 631;
 Partition Dummy = 1:183-phylip/chr9_3633.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/chr9_3633.phylip = 1 - 631;
 Partition Dummy = 1:183-phylip/chr9_3633.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

