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

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

Minimum AIC  model: GTRG
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	-2969.8459	6072.0149	0.4463	0.4463
GTRG	64	-2965.3486	6072.7990	0.3016	0.7479
HKYIG	61	-2969.6593	6074.0741	0.1594	0.9073
GTRIG	65	-2965.2981	6075.1633	0.0925	0.9998
HKYI	60	-2978.1334	6088.5899	0.0001	0.9999
GTRI	64	-2973.3714	6088.8445	0.0001	1.0000
K2PG	57	-3005.1644	6135.4043	0.0000	1.0000
SYMG	61	-3000.3361	6135.4276	0.0000	1.0000
K2PIG	58	-3004.9661	6137.4155	0.0000	1.0000
SYMIG	62	-3000.1966	6137.5892	0.0000	1.0000
SYMI	61	-3008.2598	6151.2750	0.0000	1.0000
K2PI	57	-3013.1000	6151.2754	0.0000	1.0000
GTR	63	-3036.3106	6212.2659	0.0000	1.0000
HKY	59	-3042.9171	6215.7333	0.0000	1.0000
F81G	59	-3060.6369	6251.1729	0.0000	1.0000
F81IG	60	-3060.5444	6253.4121	0.0000	1.0000
F81I	59	-3067.4839	6264.8670	0.0000	1.0000
K2P	56	-3083.1572	6288.9900	0.0000	1.0000
SYM	60	-3078.6433	6289.6098	0.0000	1.0000
JC69G	56	-3090.8228	6304.3211	0.0000	1.0000
JC69IG	57	-3090.6849	6306.4451	0.0000	1.0000
JC69I	56	-3097.9068	6318.4891	0.0000	1.0000
F81	58	-3128.5156	6384.5144	0.0000	1.0000
JC69	55	-3161.3263	6442.9365	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2965.3486	6058.6973	0.4414	0.4414
HKYG	60	-2969.8459	6059.6917	0.2685	0.7099
GTRIG	65	-2965.2981	6060.5962	0.1708	0.8808
HKYIG	61	-2969.6593	6061.3186	0.1190	0.9998
GTRI	64	-2973.3714	6074.7428	0.0001	0.9999
HKYI	60	-2978.1334	6076.2667	0.0001	1.0000
SYMG	61	-3000.3361	6122.6722	0.0000	1.0000
K2PG	57	-3005.1644	6124.3289	0.0000	1.0000
SYMIG	62	-3000.1966	6124.3932	0.0000	1.0000
K2PIG	58	-3004.9661	6125.9323	0.0000	1.0000
SYMI	61	-3008.2598	6138.5196	0.0000	1.0000
K2PI	57	-3013.1000	6140.2000	0.0000	1.0000
GTR	63	-3036.3106	6198.6212	0.0000	1.0000
HKY	59	-3042.9171	6203.8342	0.0000	1.0000
F81G	59	-3060.6369	6239.2737	0.0000	1.0000
F81IG	60	-3060.5444	6241.0889	0.0000	1.0000
F81I	59	-3067.4839	6252.9679	0.0000	1.0000
SYM	60	-3078.6433	6277.2865	0.0000	1.0000
K2P	56	-3083.1572	6278.3144	0.0000	1.0000
JC69G	56	-3090.8228	6293.6455	0.0000	1.0000
JC69IG	57	-3090.6849	6295.3697	0.0000	1.0000
JC69I	56	-3097.9068	6307.8135	0.0000	1.0000
F81	58	-3128.5156	6373.0312	0.0000	1.0000
JC69	55	-3161.3263	6432.6527	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-2969.8459	6328.7698	0.9546	0.9546
HKYIG	61	-2969.6593	6334.8814	0.0449	0.9996
HKYI	60	-2978.1334	6345.3448	0.0002	0.9998
GTRG	64	-2965.3486	6345.7139	0.0002	1.0000
GTRIG	65	-2965.2981	6352.0975	0.0000	1.0000
GTRI	64	-2973.3714	6361.7595	0.0000	1.0000
K2PG	57	-3005.1644	6379.9531	0.0000	1.0000
K2PIG	58	-3004.9661	6386.0411	0.0000	1.0000
K2PI	57	-3013.1000	6395.8242	0.0000	1.0000
SYMG	61	-3000.3361	6396.2349	0.0000	1.0000
SYMIG	62	-3000.1966	6402.4406	0.0000	1.0000
SYMI	61	-3008.2598	6412.0823	0.0000	1.0000
HKY	59	-3042.9171	6468.4276	0.0000	1.0000
GTR	63	-3036.3106	6481.1532	0.0000	1.0000
F81G	59	-3060.6369	6503.8672	0.0000	1.0000
F81IG	60	-3060.5444	6510.1670	0.0000	1.0000
F81I	59	-3067.4839	6517.5613	0.0000	1.0000
K2P	56	-3083.1572	6529.4540	0.0000	1.0000
JC69G	56	-3090.8228	6544.7851	0.0000	1.0000
SYM	60	-3078.6433	6546.3647	0.0000	1.0000
JC69IG	57	-3090.6849	6550.9939	0.0000	1.0000
JC69I	56	-3097.9068	6558.9531	0.0000	1.0000
F81	58	-3128.5156	6633.1400	0.0000	1.0000
JC69	55	-3161.3263	6679.3076	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/chr1_5522.phylip = 1 - 655;
 Partition Dummy = 1:183-phylip/chr1_5522.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 (HKYG)]

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

