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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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
GTRIG	65	-4983.0032	10107.9729	0.7232	0.7232
GTRG	64	-4985.1534	10109.8946	0.2767	0.9999
HKYIG	61	-4996.8020	10126.0950	0.0001	1.0000
HKYG	60	-4999.1143	10128.3672	0.0000	1.0000
GTRI	64	-5019.4317	10178.4511	0.0000	1.0000
SYMIG	62	-5022.8580	10180.5661	0.0000	1.0000
SYMG	61	-5026.3724	10185.2358	0.0000	1.0000
K2PIG	58	-5032.4726	10190.3982	0.0000	1.0000
HKYI	60	-5030.5406	10191.2196	0.0000	1.0000
K2PG	57	-5035.4398	10193.9997	0.0000	1.0000
SYMI	61	-5050.2335	10232.9579	0.0000	1.0000
K2PI	57	-5064.5802	10252.2804	0.0000	1.0000
F81IG	60	-5132.9481	10396.0348	0.0000	1.0000
F81G	59	-5135.8659	10399.5244	0.0000	1.0000
JC69G	56	-5148.0850	10416.9635	0.0000	1.0000
JC69IG	57	-5148.4880	10420.0959	0.0000	1.0000
F81I	59	-5164.6487	10457.0899	0.0000	1.0000
JC69I	56	-5176.4124	10473.6181	0.0000	1.0000
GTR	63	-5202.6514	10542.5184	0.0000	1.0000
HKY	59	-5216.4315	10560.6556	0.0000	1.0000
SYM	60	-5243.5172	10617.1729	0.0000	1.0000
K2P	56	-5257.5472	10635.8877	0.0000	1.0000
F81	58	-5350.8688	10827.1906	0.0000	1.0000
JC69	55	-5363.9937	10846.4606	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4983.0032	10096.0063	0.7595	0.7595
GTRG	64	-4985.1534	10098.3069	0.2404	0.9999
HKYIG	61	-4996.8020	10115.6040	0.0000	1.0000
HKYG	60	-4999.1143	10118.2287	0.0000	1.0000
GTRI	64	-5019.4317	10166.8634	0.0000	1.0000
SYMIG	62	-5022.8580	10169.7161	0.0000	1.0000
SYMG	61	-5026.3724	10174.7448	0.0000	1.0000
K2PIG	58	-5032.4726	10180.9451	0.0000	1.0000
HKYI	60	-5030.5406	10181.0811	0.0000	1.0000
K2PG	57	-5035.4398	10184.8797	0.0000	1.0000
SYMI	61	-5050.2335	10222.4669	0.0000	1.0000
K2PI	57	-5064.5802	10243.1604	0.0000	1.0000
F81IG	60	-5132.9481	10385.8963	0.0000	1.0000
F81G	59	-5135.8659	10389.7318	0.0000	1.0000
JC69G	56	-5148.0850	10408.1701	0.0000	1.0000
JC69IG	57	-5148.4880	10410.9759	0.0000	1.0000
F81I	59	-5164.6487	10447.2973	0.0000	1.0000
JC69I	56	-5176.4124	10464.8248	0.0000	1.0000
GTR	63	-5202.6514	10531.3029	0.0000	1.0000
HKY	59	-5216.4315	10550.8631	0.0000	1.0000
SYM	60	-5243.5172	10607.0344	0.0000	1.0000
K2P	56	-5257.5472	10627.0944	0.0000	1.0000
F81	58	-5350.8688	10817.7376	0.0000	1.0000
JC69	55	-5363.9937	10837.9874	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4985.1534	10396.7474	0.4930	0.4930
HKYG	60	-4999.1143	10398.0167	0.2614	0.7544
GTRIG	65	-4983.0032	10399.1100	0.1513	0.9057
HKYIG	61	-4996.8020	10400.0551	0.0943	1.0000
K2PG	57	-5035.4398	10450.6783	0.0000	1.0000
K2PIG	58	-5032.4726	10451.4068	0.0000	1.0000
SYMIG	62	-5022.8580	10458.8303	0.0000	1.0000
SYMG	61	-5026.3724	10459.1959	0.0000	1.0000
HKYI	60	-5030.5406	10460.8691	0.0000	1.0000
GTRI	64	-5019.4317	10465.3039	0.0000	1.0000
SYMI	61	-5050.2335	10506.9180	0.0000	1.0000
K2PI	57	-5064.5802	10508.9590	0.0000	1.0000
F81G	59	-5135.8659	10664.8566	0.0000	1.0000
F81IG	60	-5132.9481	10665.6842	0.0000	1.0000
JC69G	56	-5148.0850	10669.3055	0.0000	1.0000
JC69IG	57	-5148.4880	10676.7745	0.0000	1.0000
F81I	59	-5164.6487	10722.4221	0.0000	1.0000
JC69I	56	-5176.4124	10725.9602	0.0000	1.0000
GTR	63	-5202.6514	10825.0802	0.0000	1.0000
HKY	59	-5216.4315	10825.9879	0.0000	1.0000
SYM	60	-5243.5172	10886.8223	0.0000	1.0000
K2P	56	-5257.5472	10888.2298	0.0000	1.0000
F81	58	-5350.8688	11088.1993	0.0000	1.0000
JC69	55	-5363.9937	11094.4597	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_3515.phylip = 1 - 783;
 Partition Dummy = 1:183-phylip/chrZ_3515.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chrZ_3515.phylip = 1 - 783;
 Partition Dummy = 1:183-phylip/chrZ_3515.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


-------------------------------------------------------------
End of Output

