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

Input data from file "183-phylip/chr2_23548.phylip" (29 taxa, 456 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	-1457.5316	3064.3421	0.3679	0.3679
HKYG	60	-1463.1276	3064.7869	0.2945	0.6624
GTRIG	65	-1457.2610	3066.5219	0.1237	0.7861
HKYIG	61	-1462.8441	3066.8862	0.1031	0.8892
GTRI	64	-1459.3294	3067.9375	0.0609	0.9502
HKYI	60	-1464.9042	3068.3400	0.0498	1.0000
K2PG	57	-1479.0010	3088.6151	0.0000	1.0000
K2PIG	58	-1478.6844	3090.6082	0.0000	1.0000
K2PI	57	-1480.7341	3092.0813	0.0000	1.0000
GTR	63	-1473.5935	3093.7584	0.0000	1.0000
HKY	59	-1479.6777	3095.2342	0.0000	1.0000
SYMG	61	-1477.0699	3095.3377	0.0000	1.0000
SYMIG	62	-1476.7741	3097.4261	0.0000	1.0000
SYMI	61	-1478.8526	3098.9032	0.0000	1.0000
K2P	56	-1496.0039	3120.0078	0.0000	1.0000
SYM	60	-1493.6504	3125.8325	0.0000	1.0000
F81G	59	-1508.3500	3152.5787	0.0000	1.0000
F81IG	60	-1508.0572	3154.6461	0.0000	1.0000
F81I	59	-1510.1404	3156.1596	0.0000	1.0000
JC69G	56	-1522.4977	3172.9953	0.0000	1.0000
JC69IG	57	-1522.1835	3174.9800	0.0000	1.0000
JC69I	56	-1524.2409	3176.4819	0.0000	1.0000
F81	58	-1524.7103	3182.6599	0.0000	1.0000
JC69	55	-1539.2380	3203.8759	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1457.5316	3043.0633	0.5040	0.5040
GTRIG	65	-1457.2610	3044.5219	0.2431	0.7471
HKYG	60	-1463.1276	3046.2552	0.1022	0.8493
GTRI	64	-1459.3294	3046.6588	0.0835	0.9328
HKYIG	61	-1462.8441	3047.6883	0.0499	0.9827
HKYI	60	-1464.9042	3049.8083	0.0173	1.0000
K2PG	57	-1479.0010	3072.0020	0.0000	1.0000
GTR	63	-1473.5935	3073.1870	0.0000	1.0000
K2PIG	58	-1478.6844	3073.3689	0.0000	1.0000
K2PI	57	-1480.7341	3075.4683	0.0000	1.0000
SYMG	61	-1477.0699	3076.1398	0.0000	1.0000
HKY	59	-1479.6777	3077.3555	0.0000	1.0000
SYMIG	62	-1476.7741	3077.5483	0.0000	1.0000
SYMI	61	-1478.8526	3079.7053	0.0000	1.0000
K2P	56	-1496.0039	3104.0078	0.0000	1.0000
SYM	60	-1493.6504	3107.3009	0.0000	1.0000
F81G	59	-1508.3500	3134.6999	0.0000	1.0000
F81IG	60	-1508.0572	3136.1145	0.0000	1.0000
F81I	59	-1510.1404	3138.2808	0.0000	1.0000
JC69G	56	-1522.4977	3156.9953	0.0000	1.0000
JC69IG	57	-1522.1835	3158.3670	0.0000	1.0000
JC69I	56	-1524.2409	3160.4819	0.0000	1.0000
F81	58	-1524.7103	3165.4206	0.0000	1.0000
JC69	55	-1539.2380	3188.4759	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1463.1276	3293.6048	0.8100	0.8100
HKYI	60	-1464.9042	3297.1579	0.1371	0.9471
HKYIG	61	-1462.8441	3299.1603	0.0504	0.9975
GTRG	64	-1457.5316	3306.9028	0.0010	0.9985
K2PG	57	-1479.0010	3306.9841	0.0010	0.9995
K2PI	57	-1480.7341	3310.4504	0.0002	0.9997
GTRI	64	-1459.3294	3310.4983	0.0002	0.9999
K2PIG	58	-1478.6844	3312.4734	0.0001	0.9999
GTRIG	65	-1457.2610	3312.4840	0.0001	1.0000
HKY	59	-1479.6777	3320.5825	0.0000	1.0000
SYMG	61	-1477.0699	3327.6118	0.0000	1.0000
SYMI	61	-1478.8526	3331.1773	0.0000	1.0000
GTR	63	-1473.5935	3332.9040	0.0000	1.0000
SYMIG	62	-1476.7741	3333.1428	0.0000	1.0000
K2P	56	-1496.0039	3334.8674	0.0000	1.0000
SYM	60	-1493.6504	3354.6504	0.0000	1.0000
F81G	59	-1508.3500	3377.9270	0.0000	1.0000
F81I	59	-1510.1404	3381.5079	0.0000	1.0000
F81IG	60	-1508.0572	3383.4640	0.0000	1.0000
JC69G	56	-1522.4977	3387.8549	0.0000	1.0000
JC69I	56	-1524.2409	3391.3415	0.0000	1.0000
JC69IG	57	-1522.1835	3393.3491	0.0000	1.0000
F81	58	-1524.7103	3404.5252	0.0000	1.0000
JC69	55	-1539.2380	3415.2130	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/chr2_23548.phylip = 1 - 456;
 Partition Dummy = 1:183-phylip/chr2_23548.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/chr2_23548.phylip = 1 - 456;
 Partition Dummy = 1:183-phylip/chr2_23548.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/chr2_23548.phylip = 1 - 456;
 Partition Dummy = 1:183-phylip/chr2_23548.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

