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

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

Minimum AIC  model: GTRG
Minimum AICc model: SYMG
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
SYMG	61	-1199.2659	2545.5781	0.4203	0.4203
GTRG	64	-1195.6784	2547.1829	0.1884	0.6087
SYMI	61	-1200.2585	2547.5634	0.1557	0.7644
SYMIG	62	-1199.1088	2548.1712	0.1149	0.8793
GTRIG	65	-1195.4455	2549.6830	0.0540	0.9333
GTRI	64	-1196.9682	2549.7626	0.0519	0.9852
K2PG	57	-1209.0203	2553.6485	0.0074	0.9926
K2PI	57	-1209.8771	2555.3620	0.0032	0.9958
K2PIG	58	-1208.8569	2556.1531	0.0021	0.9979
HKYG	60	-1206.5373	2557.2331	0.0012	0.9991
HKYI	60	-1207.3975	2558.9534	0.0005	0.9996
HKYIG	61	-1206.3291	2559.7046	0.0004	1.0000
SYM	60	-1217.2752	2578.7089	0.0000	1.0000
GTR	63	-1214.7963	2582.4726	0.0000	1.0000
K2P	56	-1226.6328	2586.0603	0.0000	1.0000
HKY	59	-1224.4144	2590.1182	0.0000	1.0000
JC69G	56	-1236.0577	2604.9101	0.0000	1.0000
JC69I	56	-1236.8744	2606.5435	0.0000	1.0000
JC69IG	57	-1235.9181	2607.4441	0.0000	1.0000
F81G	59	-1234.2661	2609.8217	0.0000	1.0000
F81I	59	-1235.0035	2611.2965	0.0000	1.0000
F81IG	60	-1234.0671	2612.2926	0.0000	1.0000
JC69	55	-1253.3480	2636.6959	0.0000	1.0000
F81	58	-1251.5091	2641.4576	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1195.6784	2519.3569	0.3642	0.3642
SYMG	61	-1199.2659	2520.5317	0.2024	0.5666
GTRIG	65	-1195.4455	2520.8911	0.1691	0.7357
GTRI	64	-1196.9682	2521.9365	0.1003	0.8360
SYMIG	62	-1199.1088	2522.2177	0.0871	0.9231
SYMI	61	-1200.2585	2522.5170	0.0750	0.9981
K2PG	57	-1209.0203	2532.0406	0.0006	0.9987
HKYG	60	-1206.5373	2533.0747	0.0004	0.9991
K2PIG	58	-1208.8569	2533.7137	0.0003	0.9994
K2PI	57	-1209.8771	2533.7542	0.0003	0.9997
HKYIG	61	-1206.3291	2534.6582	0.0002	0.9998
HKYI	60	-1207.3975	2534.7950	0.0002	1.0000
SYM	60	-1217.2752	2554.5505	0.0000	1.0000
GTR	63	-1214.7963	2555.5926	0.0000	1.0000
K2P	56	-1226.6328	2565.2656	0.0000	1.0000
HKY	59	-1224.4144	2566.8287	0.0000	1.0000
JC69G	56	-1236.0577	2584.1154	0.0000	1.0000
JC69I	56	-1236.8744	2585.7487	0.0000	1.0000
JC69IG	57	-1235.9181	2585.8363	0.0000	1.0000
F81G	59	-1234.2661	2586.5322	0.0000	1.0000
F81I	59	-1235.0035	2588.0071	0.0000	1.0000
F81IG	60	-1234.0671	2588.1342	0.0000	1.0000
JC69	55	-1253.3480	2616.6959	0.0000	1.0000
F81	58	-1251.5091	2619.0183	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1209.0203	2754.1784	0.5967	0.5967
K2PI	57	-1209.8771	2755.8920	0.2533	0.8500
SYMG	61	-1199.2659	2758.2581	0.0776	0.9276
K2PIG	58	-1208.8569	2759.7487	0.0368	0.9644
SYMI	61	-1200.2585	2760.2434	0.0288	0.9932
SYMIG	62	-1199.1088	2763.8412	0.0048	0.9979
HKYG	60	-1206.5373	2766.9039	0.0010	0.9990
HKYI	60	-1207.3975	2768.6242	0.0004	0.9994
GTRG	64	-1195.6784	2768.7747	0.0004	0.9998
GTRI	64	-1196.9682	2771.3543	0.0001	0.9999
HKYIG	61	-1206.3291	2772.3846	0.0001	1.0000
GTRIG	65	-1195.4455	2774.2061	0.0000	1.0000
K2P	56	-1226.6328	2783.5062	0.0000	1.0000
SYM	60	-1217.2752	2788.3797	0.0000	1.0000
HKY	59	-1224.4144	2796.7608	0.0000	1.0000
GTR	63	-1214.7963	2801.1133	0.0000	1.0000
JC69G	56	-1236.0577	2802.3560	0.0000	1.0000
JC69I	56	-1236.8744	2803.9894	0.0000	1.0000
JC69IG	57	-1235.9181	2807.9740	0.0000	1.0000
F81G	59	-1234.2661	2816.4643	0.0000	1.0000
F81I	59	-1235.0035	2817.9391	0.0000	1.0000
F81IG	60	-1234.0671	2821.9634	0.0000	1.0000
JC69	55	-1253.3480	2831.0394	0.0000	1.0000
F81	58	-1251.5091	2845.0532	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/chr3_14142.phylip = 1 - 364;
 Partition Dummy = 1:183-phylip/chr3_14142.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 (SYMG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_14142.phylip = 1 - 364;
 Partition Dummy = 1:183-phylip/chr3_14142.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=Fixed(Equal);
END;


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

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

