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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-4794.9227	9727.0796	0.9658	0.9658
GTRIG	65	-4797.1141	9733.7615	0.0342	1.0000
HKYG	60	-4822.0255	9772.1394	0.0000	1.0000
GTRI	64	-4824.6967	9786.6276	0.0000	1.0000
HKYIG	61	-4829.7155	9789.7982	0.0000	1.0000
SYMG	61	-4834.2175	9798.8022	0.0000	1.0000
K2PG	57	-4842.2894	9805.8607	0.0000	1.0000
SYMIG	62	-4841.2496	9815.1504	0.0000	1.0000
K2PIG	58	-4857.5211	9838.5880	0.0000	1.0000
HKYI	60	-4868.1820	9864.4524	0.0000	1.0000
SYMI	61	-4877.6769	9885.7211	0.0000	1.0000
K2PI	57	-4889.0799	9899.4418	0.0000	1.0000
F81IG	60	-4956.5673	10041.2230	0.0000	1.0000
F81G	59	-4961.0060	10047.8265	0.0000	1.0000
JC69IG	57	-4972.8284	10066.9388	0.0000	1.0000
JC69G	56	-4975.3871	10069.7974	0.0000	1.0000
GTR	63	-4979.1126	10093.1653	0.0000	1.0000
F81I	59	-5014.4690	10154.7526	0.0000	1.0000
HKY	59	-5017.1838	10160.1822	0.0000	1.0000
SYM	60	-5017.2479	10162.5842	0.0000	1.0000
JC69I	56	-5026.3308	10171.6847	0.0000	1.0000
K2P	56	-5057.0116	10233.0464	0.0000	1.0000
F81	58	-5151.8161	10427.1780	0.0000	1.0000
JC69	55	-5162.5859	10441.9411	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-4794.9227	9717.8454	0.9605	0.9605
GTRIG	65	-4797.1141	9724.2281	0.0395	1.0000
HKYG	60	-4822.0255	9764.0510	0.0000	1.0000
GTRI	64	-4824.6967	9777.3934	0.0000	1.0000
HKYIG	61	-4829.7155	9781.4310	0.0000	1.0000
SYMG	61	-4834.2175	9790.4350	0.0000	1.0000
K2PG	57	-4842.2894	9798.5788	0.0000	1.0000
SYMIG	62	-4841.2496	9806.4993	0.0000	1.0000
K2PIG	58	-4857.5211	9831.0423	0.0000	1.0000
HKYI	60	-4868.1820	9856.3640	0.0000	1.0000
SYMI	61	-4877.6769	9877.3538	0.0000	1.0000
K2PI	57	-4889.0799	9892.1598	0.0000	1.0000
F81IG	60	-4956.5673	10033.1346	0.0000	1.0000
F81G	59	-4961.0060	10040.0120	0.0000	1.0000
JC69IG	57	-4972.8284	10059.6569	0.0000	1.0000
JC69G	56	-4975.3871	10062.7743	0.0000	1.0000
GTR	63	-4979.1126	10084.2251	0.0000	1.0000
F81I	59	-5014.4690	10146.9380	0.0000	1.0000
HKY	59	-5017.1838	10152.3677	0.0000	1.0000
SYM	60	-5017.2479	10154.4958	0.0000	1.0000
JC69I	56	-5026.3308	10164.6616	0.0000	1.0000
K2P	56	-5057.0116	10226.0233	0.0000	1.0000
F81	58	-5151.8161	10419.6322	0.0000	1.0000
JC69	55	-5162.5859	10435.1719	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-4794.9227	10029.7279	0.9964	0.9964
GTRIG	65	-4797.1141	10040.9838	0.0036	1.0000
HKYG	60	-4822.0255	10056.4408	0.0000	1.0000
K2PG	57	-4842.2894	10076.3491	0.0000	1.0000
HKYIG	61	-4829.7155	10078.6940	0.0000	1.0000
SYMG	61	-4834.2175	10087.6980	0.0000	1.0000
GTRI	64	-4824.6967	10089.2759	0.0000	1.0000
SYMIG	62	-4841.2496	10108.6354	0.0000	1.0000
K2PIG	58	-4857.5211	10113.6858	0.0000	1.0000
HKYI	60	-4868.1820	10148.7539	0.0000	1.0000
K2PI	57	-4889.0799	10169.9302	0.0000	1.0000
SYMI	61	-4877.6769	10174.6168	0.0000	1.0000
F81IG	60	-4956.5673	10325.5244	0.0000	1.0000
F81G	59	-4961.0060	10327.5286	0.0000	1.0000
JC69G	56	-4975.3871	10335.6715	0.0000	1.0000
JC69IG	57	-4972.8284	10337.4272	0.0000	1.0000
GTR	63	-4979.1126	10391.2345	0.0000	1.0000
F81I	59	-5014.4690	10434.4547	0.0000	1.0000
JC69I	56	-5026.3308	10437.5588	0.0000	1.0000
HKY	59	-5017.1838	10439.8843	0.0000	1.0000
SYM	60	-5017.2479	10446.8857	0.0000	1.0000
K2P	56	-5057.0116	10498.9204	0.0000	1.0000
F81	58	-5151.8161	10702.2757	0.0000	1.0000
JC69	55	-5162.5859	10703.1959	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/chr5_3310.phylip = 1 - 966;
 Partition Dummy = 1:183-phylip/chr5_3310.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/chr5_3310.phylip = 1 - 966;
 Partition Dummy = 1:183-phylip/chr5_3310.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 (GTRG)]

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

