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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-1292.0026	2741.5051	0.4661	0.4661
GTRG	64	-1293.9719	2742.5253	0.2798	0.7459
GTRI	64	-1294.6320	2743.8455	0.1446	0.8906
SYMIG	62	-1298.4855	2745.7709	0.0552	0.9458
SYMG	61	-1300.4811	2746.8988	0.0314	0.9772
SYMI	61	-1300.8666	2747.6699	0.0214	0.9986
K2PIG	58	-1308.3181	2754.0907	0.0009	0.9994
K2PG	57	-1310.8375	2756.3375	0.0003	0.9997
K2PI	57	-1311.4425	2757.5475	0.0002	0.9999
HKYIG	61	-1306.4448	2758.8263	0.0001	0.9999
HKYG	60	-1308.8060	2760.7034	0.0000	1.0000
HKYI	60	-1309.3090	2761.7094	0.0000	1.0000
JC69G	56	-1337.9248	2807.7375	0.0000	1.0000
F81IG	60	-1332.9825	2809.0565	0.0000	1.0000
F81G	59	-1335.2795	2810.8232	0.0000	1.0000
JC69I	56	-1339.5226	2810.9331	0.0000	1.0000
JC69IG	57	-1338.4219	2811.5062	0.0000	1.0000
F81I	59	-1337.0278	2814.3198	0.0000	1.0000
GTR	63	-1356.5296	2864.7408	0.0000	1.0000
SYM	60	-1362.2073	2867.5062	0.0000	1.0000
K2P	56	-1371.8934	2875.6747	0.0000	1.0000
HKY	59	-1369.5343	2879.3327	0.0000	1.0000
JC69	55	-1399.2714	2927.6732	0.0000	1.0000
F81	58	-1396.7763	2931.0071	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1292.0026	2714.0051	0.6153	0.6153
GTRG	64	-1293.9719	2715.9439	0.2334	0.8487
GTRI	64	-1294.6320	2717.2640	0.1206	0.9693
SYMIG	62	-1298.4855	2720.9709	0.0189	0.9882
SYMG	61	-1300.4811	2722.9621	0.0070	0.9952
SYMI	61	-1300.8666	2723.7332	0.0047	0.9999
K2PIG	58	-1308.3181	2732.6362	0.0001	1.0000
HKYIG	61	-1306.4448	2734.8896	0.0000	1.0000
K2PG	57	-1310.8375	2735.6750	0.0000	1.0000
K2PI	57	-1311.4425	2736.8850	0.0000	1.0000
HKYG	60	-1308.8060	2737.6119	0.0000	1.0000
HKYI	60	-1309.3090	2738.6179	0.0000	1.0000
F81IG	60	-1332.9825	2785.9651	0.0000	1.0000
JC69G	56	-1337.9248	2787.8496	0.0000	1.0000
F81G	59	-1335.2795	2788.5590	0.0000	1.0000
JC69IG	57	-1338.4219	2790.8437	0.0000	1.0000
JC69I	56	-1339.5226	2791.0452	0.0000	1.0000
F81I	59	-1337.0278	2792.0557	0.0000	1.0000
GTR	63	-1356.5296	2839.0592	0.0000	1.0000
SYM	60	-1362.2073	2844.4147	0.0000	1.0000
K2P	56	-1371.8934	2855.7868	0.0000	1.0000
HKY	59	-1369.5343	2857.0685	0.0000	1.0000
JC69	55	-1399.2714	2908.5428	0.0000	1.0000
F81	58	-1396.7763	2909.5525	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1310.8375	2959.9640	0.3734	0.3734
K2PIG	58	-1308.3181	2960.8601	0.2385	0.6119
K2PI	57	-1311.4425	2961.1740	0.2039	0.8158
SYMG	61	-1300.4811	2962.9907	0.0822	0.8980
SYMI	61	-1300.8666	2963.7617	0.0559	0.9539
SYMIG	62	-1298.4855	2964.9344	0.0311	0.9850
GTRG	64	-1293.9719	2967.7771	0.0075	0.9925
GTRI	64	-1294.6320	2969.0973	0.0039	0.9964
GTRIG	65	-1292.0026	2969.7732	0.0028	0.9992
HKYG	60	-1308.8060	2973.7056	0.0004	0.9996
HKYI	60	-1309.3090	2974.7116	0.0002	0.9998
HKYIG	61	-1306.4448	2974.9181	0.0002	1.0000
JC69G	56	-1337.9248	3008.2037	0.0000	1.0000
JC69I	56	-1339.5226	3011.3993	0.0000	1.0000
JC69IG	57	-1338.4219	3015.1327	0.0000	1.0000
F81G	59	-1335.2795	3020.7178	0.0000	1.0000
F81IG	60	-1332.9825	3022.0587	0.0000	1.0000
F81I	59	-1337.0278	3024.2144	0.0000	1.0000
K2P	56	-1371.8934	3076.1409	0.0000	1.0000
SYM	60	-1362.2073	3080.5083	0.0000	1.0000
GTR	63	-1356.5296	3086.9576	0.0000	1.0000
HKY	59	-1369.5343	3089.2273	0.0000	1.0000
JC69	55	-1399.2714	3124.9619	0.0000	1.0000
F81	58	-1396.7763	3137.7764	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_31767.phylip = 1 - 378;
 Partition Dummy = 1:183-phylip/chr1_31767.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr1_31767.phylip = 1 - 378;
 Partition Dummy = 1:183-phylip/chr1_31767.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

