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

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

Minimum AIC  model: GTRG
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	-1262.5293	2664.3725	0.3375	0.3375
GTRG	64	-1257.5359	2665.2585	0.2167	0.5543
HKYI	60	-1263.0907	2665.4953	0.1925	0.7468
GTRI	64	-1258.2472	2666.6810	0.1064	0.8532
HKYIG	61	-1262.5090	2667.0285	0.0894	0.9427
GTRIG	65	-1257.4939	2667.9289	0.0570	0.9997
HKY	59	-1271.1176	2678.8668	0.0002	0.9999
GTR	63	-1266.8652	2681.1771	0.0001	1.0000
F81G	59	-1282.2080	2701.0475	0.0000	1.0000
K2PG	57	-1285.3233	2701.9555	0.0000	1.0000
F81I	59	-1282.7381	2702.1078	0.0000	1.0000
K2PI	57	-1285.8369	2702.9827	0.0000	1.0000
F81IG	60	-1282.1868	2703.6876	0.0000	1.0000
K2PIG	58	-1285.2985	2704.5602	0.0000	1.0000
SYMG	61	-1283.8581	2709.7267	0.0000	1.0000
SYMI	61	-1284.3812	2710.7730	0.0000	1.0000
SYMIG	62	-1283.8310	2712.3834	0.0000	1.0000
F81	58	-1290.7905	2715.5442	0.0000	1.0000
K2P	56	-1294.1738	2717.0160	0.0000	1.0000
SYM	60	-1292.7266	2724.7672	0.0000	1.0000
JC69G	56	-1302.0525	2732.7733	0.0000	1.0000
JC69I	56	-1302.5942	2733.8568	0.0000	1.0000
JC69IG	57	-1302.0245	2735.3579	0.0000	1.0000
JC69	55	-1310.8560	2747.7536	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1257.5359	2643.0718	0.3852	0.3852
GTRI	64	-1258.2472	2644.4944	0.1892	0.5744
GTRIG	65	-1257.4939	2644.9877	0.1478	0.7222
HKYG	60	-1262.5293	2645.0586	0.1427	0.8649
HKYI	60	-1263.0907	2646.1813	0.0814	0.9463
HKYIG	61	-1262.5090	2647.0179	0.0536	0.9998
GTR	63	-1266.8652	2659.7303	0.0001	0.9999
HKY	59	-1271.1176	2660.2352	0.0001	1.0000
F81G	59	-1282.2080	2682.4159	0.0000	1.0000
F81I	59	-1282.7381	2683.4762	0.0000	1.0000
F81IG	60	-1282.1868	2684.3737	0.0000	1.0000
K2PG	57	-1285.3233	2684.6466	0.0000	1.0000
K2PI	57	-1285.8369	2685.6738	0.0000	1.0000
K2PIG	58	-1285.2985	2686.5969	0.0000	1.0000
SYMG	61	-1283.8581	2689.7161	0.0000	1.0000
SYMI	61	-1284.3812	2690.7624	0.0000	1.0000
SYMIG	62	-1283.8310	2691.6619	0.0000	1.0000
F81	58	-1290.7905	2697.5810	0.0000	1.0000
K2P	56	-1294.1738	2700.3476	0.0000	1.0000
SYM	60	-1292.7266	2705.4532	0.0000	1.0000
JC69G	56	-1302.0525	2716.1049	0.0000	1.0000
JC69I	56	-1302.5942	2717.1884	0.0000	1.0000
JC69IG	57	-1302.0245	2718.0490	0.0000	1.0000
JC69	55	-1310.8560	2731.7119	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1262.5293	2890.2650	0.6157	0.6157
HKYI	60	-1263.0907	2891.3878	0.3512	0.9669
HKYIG	61	-1262.5090	2896.3112	0.0300	0.9969
HKY	59	-1271.1176	2901.3549	0.0024	0.9993
GTRG	64	-1257.5359	2904.6254	0.0005	0.9997
GTRI	64	-1258.2472	2906.0480	0.0002	1.0000
GTRIG	65	-1257.4939	2910.6281	0.0000	1.0000
GTR	63	-1266.8652	2917.1971	0.0000	1.0000
K2PG	57	-1285.3233	2917.5928	0.0000	1.0000
K2PI	57	-1285.8369	2918.6199	0.0000	1.0000
F81G	59	-1282.2080	2923.5356	0.0000	1.0000
K2PIG	58	-1285.2985	2923.6299	0.0000	1.0000
F81I	59	-1282.7381	2924.5959	0.0000	1.0000
K2P	56	-1294.1738	2929.2070	0.0000	1.0000
F81IG	60	-1282.1868	2929.5801	0.0000	1.0000
F81	58	-1290.7905	2934.6139	0.0000	1.0000
SYMG	61	-1283.8581	2939.0094	0.0000	1.0000
SYMI	61	-1284.3812	2940.0557	0.0000	1.0000
JC69G	56	-1302.0525	2944.9643	0.0000	1.0000
SYMIG	62	-1283.8310	2945.0420	0.0000	1.0000
JC69I	56	-1302.5942	2946.0477	0.0000	1.0000
SYM	60	-1292.7266	2950.6597	0.0000	1.0000
JC69IG	57	-1302.0245	2950.9952	0.0000	1.0000
JC69	55	-1310.8560	2956.4845	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/chr7_10849.phylip = 1 - 440;
 Partition Dummy = 1:183-phylip/chr7_10849.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 (HKYG)]

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

