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

Input data from file "183-phylip/chr2_27968.phylip" (29 taxa, 400 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	-1268.2661	2678.1251	0.4523	0.4523
HKYI	60	-1268.8298	2679.2524	0.2574	0.7097
HKYIG	61	-1267.8980	2680.1747	0.1623	0.8720
GTRG	64	-1264.5871	2682.0100	0.0648	0.9368
GTRI	64	-1265.0678	2682.9714	0.0401	0.9769
GTRIG	65	-1264.2471	2684.1828	0.0219	0.9988
K2PG	57	-1279.0391	2691.4116	0.0006	0.9993
K2PI	57	-1279.6354	2692.6042	0.0003	0.9997
K2PIG	58	-1278.7027	2693.4757	0.0002	0.9999
SYMG	61	-1275.7167	2695.8120	0.0001	0.9999
SYMI	61	-1276.4146	2697.2079	0.0000	1.0000
SYMIG	62	-1275.4150	2698.0109	0.0000	1.0000
HKY	59	-1289.4219	2717.6673	0.0000	1.0000
GTR	63	-1285.3945	2720.7890	0.0000	1.0000
F81G	59	-1291.0287	2720.8810	0.0000	1.0000
F81I	59	-1291.4939	2721.8113	0.0000	1.0000
F81IG	60	-1290.7259	2723.0448	0.0000	1.0000
K2P	56	-1298.6194	2727.8511	0.0000	1.0000
JC69G	56	-1299.2659	2729.1441	0.0000	1.0000
JC69I	56	-1299.8075	2730.2272	0.0000	1.0000
SYM	60	-1294.6302	2730.8533	0.0000	1.0000
JC69IG	57	-1298.9725	2731.2783	0.0000	1.0000
F81	58	-1311.0359	2758.1422	0.0000	1.0000
JC69	55	-1318.7145	2765.3360	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1268.2661	2656.5322	0.2739	0.2739
GTRG	64	-1264.5871	2657.1742	0.1987	0.4726
HKYI	60	-1268.8298	2657.6595	0.1559	0.6285
HKYIG	61	-1267.8980	2657.7960	0.1456	0.7741
GTRI	64	-1265.0678	2658.1356	0.1229	0.8969
GTRIG	65	-1264.2471	2658.4942	0.1027	0.9996
K2PG	57	-1279.0391	2672.0783	0.0001	0.9998
K2PI	57	-1279.6354	2673.2708	0.0001	0.9998
K2PIG	58	-1278.7027	2673.4053	0.0001	0.9999
SYMG	61	-1275.7167	2673.4333	0.0001	0.9999
SYMI	61	-1276.4146	2674.8292	0.0000	1.0000
SYMIG	62	-1275.4150	2674.8299	0.0000	1.0000
GTR	63	-1285.3945	2696.7890	0.0000	1.0000
HKY	59	-1289.4219	2696.8438	0.0000	1.0000
F81G	59	-1291.0287	2700.0575	0.0000	1.0000
F81I	59	-1291.4939	2700.9878	0.0000	1.0000
F81IG	60	-1290.7259	2701.4519	0.0000	1.0000
K2P	56	-1298.6194	2709.2389	0.0000	1.0000
SYM	60	-1294.6302	2709.2604	0.0000	1.0000
JC69G	56	-1299.2659	2710.5318	0.0000	1.0000
JC69I	56	-1299.8075	2711.6149	0.0000	1.0000
JC69IG	57	-1298.9725	2711.9450	0.0000	1.0000
F81	58	-1311.0359	2738.0718	0.0000	1.0000
JC69	55	-1318.7145	2747.4291	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1268.2661	2896.0201	0.5226	0.5226
HKYI	60	-1268.8298	2897.1474	0.2974	0.8200
K2PG	57	-1279.0391	2899.5917	0.0876	0.9076
K2PI	57	-1279.6354	2900.7843	0.0483	0.9559
HKYIG	61	-1267.8980	2901.2754	0.0378	0.9936
K2PIG	58	-1278.7027	2904.9103	0.0061	0.9998
GTRG	64	-1264.5871	2912.6279	0.0001	0.9999
GTRI	64	-1265.0678	2913.5894	0.0001	1.0000
SYMG	61	-1275.7167	2916.9127	0.0000	1.0000
GTRIG	65	-1264.2471	2917.9394	0.0000	1.0000
SYMI	61	-1276.4146	2918.3085	0.0000	1.0000
SYMIG	62	-1275.4150	2922.3007	0.0000	1.0000
HKY	59	-1289.4219	2932.3402	0.0000	1.0000
K2P	56	-1298.6194	2932.7609	0.0000	1.0000
JC69G	56	-1299.2659	2934.0539	0.0000	1.0000
JC69I	56	-1299.8075	2935.1370	0.0000	1.0000
F81G	59	-1291.0287	2935.5539	0.0000	1.0000
F81I	59	-1291.4939	2936.4842	0.0000	1.0000
JC69IG	57	-1298.9725	2939.4585	0.0000	1.0000
F81IG	60	-1290.7259	2940.9398	0.0000	1.0000
GTR	63	-1285.3945	2948.2513	0.0000	1.0000
SYM	60	-1294.6302	2948.7482	0.0000	1.0000
JC69	55	-1318.7145	2966.9596	0.0000	1.0000
F81	58	-1311.0359	2969.5767	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/chr2_27968.phylip = 1 - 400;
 Partition Dummy = 1:183-phylip/chr2_27968.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/chr2_27968.phylip = 1 - 400;
 Partition Dummy = 1:183-phylip/chr2_27968.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/chr2_27968.phylip = 1 - 400;
 Partition Dummy = 1:183-phylip/chr2_27968.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

