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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1443.6803	3042.8946	0.2928	0.2928
HKYI	60	-1443.7981	3043.1302	0.2602	0.5530
K2PG	57	-1449.3011	3044.2385	0.1495	0.7026
K2PI	57	-1449.3772	3044.3907	0.1386	0.8411
HKYIG	61	-1443.1230	3045.1435	0.0951	0.9362
K2PIG	58	-1448.7042	3046.3123	0.0530	0.9893
SYMG	61	-1446.2800	3051.4576	0.0040	0.9933
SYMI	61	-1446.5648	3052.0272	0.0030	0.9963
SYMIG	62	-1445.7352	3053.7645	0.0013	0.9976
GTRG	64	-1442.3631	3053.9142	0.0012	0.9988
GTRI	64	-1442.6727	3054.5335	0.0009	0.9997
GTRIG	65	-1441.8998	3056.4861	0.0003	1.0000
HKY	59	-1464.9940	3082.1908	0.0000	1.0000
K2P	56	-1470.2287	3082.8574	0.0000	1.0000
SYM	60	-1468.8114	3093.1568	0.0000	1.0000
GTR	63	-1464.0715	3093.8671	0.0000	1.0000
F81I	59	-1497.1285	3146.4599	0.0000	1.0000
F81G	59	-1497.4672	3147.1373	0.0000	1.0000
F81IG	60	-1496.6279	3148.7898	0.0000	1.0000
JC69I	56	-1504.1462	3150.6925	0.0000	1.0000
JC69G	56	-1504.6601	3151.7201	0.0000	1.0000
JC69IG	57	-1503.7310	3153.0984	0.0000	1.0000
F81	58	-1516.5526	3182.0091	0.0000	1.0000
JC69	55	-1523.3489	3185.8921	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1443.6803	3007.3606	0.3397	0.3397
HKYI	60	-1443.7981	3007.5962	0.3020	0.6417
HKYIG	61	-1443.1230	3008.2460	0.2182	0.8599
K2PG	57	-1449.3011	3012.6021	0.0247	0.8846
GTRG	64	-1442.3631	3012.7261	0.0232	0.9078
K2PI	57	-1449.3772	3012.7543	0.0229	0.9307
GTRI	64	-1442.6727	3013.3454	0.0170	0.9477
K2PIG	58	-1448.7042	3013.4084	0.0165	0.9643
GTRIG	65	-1441.8998	3013.7995	0.0136	0.9778
SYMG	61	-1446.2800	3014.5601	0.0093	0.9871
SYMI	61	-1446.5648	3015.1296	0.0070	0.9941
SYMIG	62	-1445.7352	3015.4704	0.0059	1.0000
HKY	59	-1464.9940	3047.9879	0.0000	1.0000
K2P	56	-1470.2287	3052.4574	0.0000	1.0000
GTR	63	-1464.0715	3054.1429	0.0000	1.0000
SYM	60	-1468.8114	3057.6228	0.0000	1.0000
F81I	59	-1497.1285	3112.2570	0.0000	1.0000
F81G	59	-1497.4672	3112.9344	0.0000	1.0000
F81IG	60	-1496.6279	3113.2559	0.0000	1.0000
JC69I	56	-1504.1462	3120.2925	0.0000	1.0000
JC69G	56	-1504.6601	3121.3201	0.0000	1.0000
JC69IG	57	-1503.7310	3121.4620	0.0000	1.0000
F81	58	-1516.5526	3149.1053	0.0000	1.0000
JC69	55	-1523.3489	3156.6978	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1449.3011	3217.0753	0.4620	0.4620
K2PI	57	-1449.3772	3217.2275	0.4281	0.8901
K2PIG	58	-1448.7042	3221.4689	0.0514	0.9414
HKYG	60	-1443.6803	3222.5955	0.0292	0.9706
HKYI	60	-1443.7981	3222.8312	0.0260	0.9966
HKYIG	61	-1443.1230	3227.0681	0.0031	0.9998
SYMG	61	-1446.2800	3233.3822	0.0001	0.9999
SYMI	61	-1446.5648	3233.9518	0.0001	1.0000
SYMIG	62	-1445.7352	3237.8798	0.0000	1.0000
GTRG	64	-1442.3631	3242.3100	0.0000	1.0000
GTRI	64	-1442.6727	3242.9293	0.0000	1.0000
GTRIG	65	-1441.8998	3246.9707	0.0000	1.0000
K2P	56	-1470.2287	3253.3434	0.0000	1.0000
HKY	59	-1464.9940	3259.6356	0.0000	1.0000
SYM	60	-1468.8114	3272.8577	0.0000	1.0000
GTR	63	-1464.0715	3280.1396	0.0000	1.0000
JC69I	56	-1504.1462	3321.1784	0.0000	1.0000
JC69G	56	-1504.6601	3322.2061	0.0000	1.0000
F81I	59	-1497.1285	3323.9047	0.0000	1.0000
F81G	59	-1497.4672	3324.5821	0.0000	1.0000
JC69IG	57	-1503.7310	3325.9352	0.0000	1.0000
F81IG	60	-1496.6279	3328.4908	0.0000	1.0000
JC69	55	-1523.3489	3353.9965	0.0000	1.0000
F81	58	-1516.5526	3357.1657	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/chr3_24921.phylip = 1 - 267;
 Partition Dummy = 1:183-phylip/chr3_24921.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/chr3_24921.phylip = 1 - 267;
 Partition Dummy = 1:183-phylip/chr3_24921.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 (K2PG)]

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

