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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1477.8186	3103.7826	0.5611	0.5611
HKYG	60	-1484.0841	3105.7221	0.2128	0.7739
HKYI	60	-1484.5558	3106.6656	0.1327	0.9066
HKYIG	61	-1484.0966	3108.3760	0.0564	0.9630
GTRI	64	-1481.0479	3110.2411	0.0222	0.9853
GTRIG	65	-1480.1179	3111.0611	0.0147	1.0000
HKY	59	-1506.1449	3147.2276	0.0000	1.0000
GTR	63	-1502.9465	3151.3713	0.0000	1.0000
F81G	59	-1514.6755	3164.2889	0.0000	1.0000
F81I	59	-1515.1416	3165.2210	0.0000	1.0000
F81IG	60	-1514.6041	3166.7623	0.0000	1.0000
SYMG	61	-1520.9575	3182.0977	0.0000	1.0000
SYMI	61	-1523.1892	3186.5611	0.0000	1.0000
K2PG	57	-1529.2875	3188.3178	0.0000	1.0000
K2PI	57	-1529.3572	3188.4573	0.0000	1.0000
SYMIG	62	-1522.8279	3188.4800	0.0000	1.0000
K2PIG	58	-1529.0727	3190.4796	0.0000	1.0000
F81	58	-1540.3810	3213.0962	0.0000	1.0000
SYM	60	-1544.1395	3225.8330	0.0000	1.0000
K2P	56	-1551.8231	3230.8101	0.0000	1.0000
JC69G	56	-1556.2574	3239.6786	0.0000	1.0000
JC69I	56	-1556.4017	3239.9673	0.0000	1.0000
JC69IG	57	-1556.0907	3241.9242	0.0000	1.0000
JC69	55	-1578.4759	3281.5489	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1477.8186	3083.6373	0.7796	0.7796
HKYG	60	-1484.0841	3088.1682	0.0809	0.8605
HKYI	60	-1484.5558	3089.1116	0.0505	0.9110
GTRI	64	-1481.0479	3090.0958	0.0309	0.9418
HKYIG	61	-1484.0966	3090.1933	0.0294	0.9712
GTRIG	65	-1480.1179	3090.2359	0.0288	1.0000
HKY	59	-1506.1449	3130.2898	0.0000	1.0000
GTR	63	-1502.9465	3131.8930	0.0000	1.0000
F81G	59	-1514.6755	3147.3511	0.0000	1.0000
F81I	59	-1515.1416	3148.2832	0.0000	1.0000
F81IG	60	-1514.6041	3149.2083	0.0000	1.0000
SYMG	61	-1520.9575	3163.9150	0.0000	1.0000
SYMI	61	-1523.1892	3168.3784	0.0000	1.0000
SYMIG	62	-1522.8279	3169.6559	0.0000	1.0000
K2PG	57	-1529.2875	3172.5749	0.0000	1.0000
K2PI	57	-1529.3572	3172.7144	0.0000	1.0000
K2PIG	58	-1529.0727	3174.1455	0.0000	1.0000
F81	58	-1540.3810	3196.7621	0.0000	1.0000
SYM	60	-1544.1395	3208.2790	0.0000	1.0000
K2P	56	-1551.8231	3215.6462	0.0000	1.0000
JC69G	56	-1556.2574	3224.5147	0.0000	1.0000
JC69I	56	-1556.4017	3224.8034	0.0000	1.0000
JC69IG	57	-1556.0907	3226.1814	0.0000	1.0000
JC69	55	-1578.4759	3266.9517	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1484.0841	3338.3448	0.5983	0.5983
HKYI	60	-1484.5558	3339.2883	0.3733	0.9715
HKYIG	61	-1484.0966	3344.5396	0.0270	0.9986
GTRG	64	-1477.8186	3350.4924	0.0014	0.9999
GTRI	64	-1481.0479	3356.9509	0.0001	1.0000
GTRIG	65	-1480.1179	3361.2606	0.0000	1.0000
HKY	59	-1506.1449	3376.2968	0.0000	1.0000
F81G	59	-1514.6755	3393.3581	0.0000	1.0000
F81I	59	-1515.1416	3394.2902	0.0000	1.0000
GTR	63	-1502.9465	3394.5785	0.0000	1.0000
F81IG	60	-1514.6041	3399.3849	0.0000	1.0000
K2PG	57	-1529.2875	3410.2428	0.0000	1.0000
K2PI	57	-1529.3572	3410.3822	0.0000	1.0000
K2PIG	58	-1529.0727	3415.9829	0.0000	1.0000
SYMG	61	-1520.9575	3418.2613	0.0000	1.0000
SYMI	61	-1523.1892	3422.7246	0.0000	1.0000
SYMIG	62	-1522.8279	3428.1717	0.0000	1.0000
F81	58	-1540.3810	3438.5995	0.0000	1.0000
K2P	56	-1551.8231	3449.1444	0.0000	1.0000
JC69G	56	-1556.2574	3458.0129	0.0000	1.0000
JC69I	56	-1556.4017	3458.3016	0.0000	1.0000
SYM	60	-1544.1395	3458.4557	0.0000	1.0000
JC69IG	57	-1556.0907	3463.8492	0.0000	1.0000
JC69	55	-1578.4759	3496.2803	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/chrZ_532.phylip = 1 - 478;
 Partition Dummy = 1:183-phylip/chrZ_532.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 (GTRG)]

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_532.phylip = 1 - 478;
 Partition Dummy = 1:183-phylip/chrZ_532.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 BIC model (HKYG)]

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

