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

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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	57	-1352.2495	2842.9879	0.6577	0.6577
K2PIG	58	-1352.0940	2845.6305	0.1755	0.8332
K2PI	57	-1354.4085	2847.3059	0.0759	0.9092
HKYG	60	-1350.3330	2848.0817	0.0515	0.9607
HKYIG	61	-1350.1963	2850.8286	0.0130	0.9737
SYMG	61	-1350.2011	2850.8383	0.0130	0.9867
HKYI	60	-1352.8308	2853.0773	0.0042	0.9909
GTRG	64	-1346.8468	2853.3286	0.0037	0.9947
SYMIG	62	-1350.0391	2853.5574	0.0033	0.9980
SYMI	61	-1352.8008	2856.0378	0.0010	0.9990
GTRIG	65	-1346.7079	2856.1640	0.0009	0.9999
GTRI	64	-1350.2388	2860.1127	0.0001	1.0000
K2P	56	-1381.1469	2897.8511	0.0000	1.0000
JC69G	56	-1383.2611	2902.0794	0.0000	1.0000
HKY	59	-1379.2064	2902.8307	0.0000	1.0000
SYM	60	-1377.7097	2902.8351	0.0000	1.0000
JC69IG	57	-1383.2219	2904.9326	0.0000	1.0000
GTR	63	-1374.3961	2905.3376	0.0000	1.0000
JC69I	56	-1385.5828	2906.7228	0.0000	1.0000
F81G	59	-1382.1856	2908.7891	0.0000	1.0000
F81IG	60	-1382.1389	2911.6935	0.0000	1.0000
F81I	59	-1385.3414	2915.1007	0.0000	1.0000
JC69	55	-1410.6583	2953.9636	0.0000	1.0000
F81	58	-1409.6373	2960.7169	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-1352.2495	2818.4990	0.3902	0.3902
K2PIG	58	-1352.0940	2820.1881	0.1677	0.5578
HKYG	60	-1350.3330	2820.6660	0.1320	0.6899
GTRG	64	-1346.8468	2821.6936	0.0790	0.7689
HKYIG	61	-1350.1963	2822.3925	0.0557	0.8245
SYMG	61	-1350.2011	2822.4022	0.0554	0.8800
K2PI	57	-1354.4085	2822.8170	0.0450	0.9250
GTRIG	65	-1346.7079	2823.4159	0.0334	0.9584
SYMIG	62	-1350.0391	2824.0782	0.0240	0.9824
HKYI	60	-1352.8308	2825.6615	0.0109	0.9932
SYMI	61	-1352.8008	2827.6017	0.0041	0.9973
GTRI	64	-1350.2388	2828.4777	0.0027	1.0000
K2P	56	-1381.1469	2874.2939	0.0000	1.0000
GTR	63	-1374.3961	2874.7922	0.0000	1.0000
SYM	60	-1377.7097	2875.4193	0.0000	1.0000
HKY	59	-1379.2064	2876.4128	0.0000	1.0000
JC69G	56	-1383.2611	2878.5222	0.0000	1.0000
JC69IG	57	-1383.2219	2880.4437	0.0000	1.0000
F81G	59	-1382.1856	2882.3712	0.0000	1.0000
JC69I	56	-1385.5828	2883.1656	0.0000	1.0000
F81IG	60	-1382.1389	2884.2778	0.0000	1.0000
F81I	59	-1385.3414	2888.6828	0.0000	1.0000
JC69	55	-1410.6583	2931.3165	0.0000	1.0000
F81	58	-1409.6373	2935.2745	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1352.2495	3034.7008	0.8465	0.8465
K2PI	57	-1354.4085	3039.0188	0.0977	0.9442
K2PIG	58	-1352.0940	3040.1829	0.0546	0.9988
HKYG	60	-1350.3330	3048.2468	0.0010	0.9998
HKYI	60	-1352.8308	3053.2423	0.0001	0.9999
HKYIG	61	-1350.1963	3053.7664	0.0001	0.9999
SYMG	61	-1350.2011	3053.7760	0.0001	1.0000
SYMI	61	-1352.8008	3058.9755	0.0000	1.0000
SYMIG	62	-1350.0391	3059.2450	0.0000	1.0000
GTRG	64	-1346.8468	3064.4465	0.0000	1.0000
GTRIG	65	-1346.7079	3069.9618	0.0000	1.0000
GTRI	64	-1350.2388	3071.2306	0.0000	1.0000
K2P	56	-1381.1469	3086.7026	0.0000	1.0000
JC69G	56	-1383.2611	3090.9310	0.0000	1.0000
JC69I	56	-1385.5828	3095.5743	0.0000	1.0000
JC69IG	57	-1383.2219	3096.6455	0.0000	1.0000
HKY	59	-1379.2064	3100.2006	0.0000	1.0000
SYM	60	-1377.7097	3103.0002	0.0000	1.0000
F81G	59	-1382.1856	3106.1590	0.0000	1.0000
F81IG	60	-1382.1389	3111.8586	0.0000	1.0000
F81I	59	-1385.3414	3112.4706	0.0000	1.0000
GTR	63	-1374.3961	3113.7520	0.0000	1.0000
JC69	55	-1410.6583	3139.9323	0.0000	1.0000
F81	58	-1409.6373	3155.2693	0.0000	1.0000



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

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

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


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

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


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

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

