
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:28 04/04/2011

Input data from file "591-tomodel-phylip/chr1_28786.phylip" (7 taxa, 827 characters)

Minimum AIC  model: HKYG
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	16	-2571.3604	5175.3924	0.6150	0.6150
HKYIG	17	-2571.3517	5177.4599	0.2187	0.8337
GTRG	20	-2569.1231	5179.2884	0.0877	0.9214
HKYI	16	-2574.0396	5180.7508	0.0422	0.9636
GTRIG	21	-2569.1145	5181.3769	0.0309	0.9944
GTRI	20	-2571.8760	5184.7942	0.0056	1.0000
HKY	15	-2597.7896	5226.1711	0.0000	1.0000
GTR	19	-2595.2782	5229.4982	0.0000	1.0000
F81G	15	-2627.0861	5284.7641	0.0000	1.0000
F81IG	16	-2627.0922	5286.8560	0.0000	1.0000
F81I	15	-2628.3037	5287.1992	0.0000	1.0000
F81	14	-2654.0099	5336.5371	0.0000	1.0000
SYMG	17	-2658.6949	5352.1462	0.0000	1.0000
SYMIG	18	-2658.6964	5354.2394	0.0000	1.0000
SYMI	17	-2660.5699	5355.8962	0.0000	1.0000
K2PG	13	-2664.7926	5356.0329	0.0000	1.0000
K2PIG	14	-2664.7926	5358.1025	0.0000	1.0000
K2PI	13	-2666.2953	5359.0382	0.0000	1.0000
SYM	16	-2689.4324	5411.5364	0.0000	1.0000
K2P	12	-2694.3402	5413.0637	0.0000	1.0000
JC69G	12	-2714.5252	5453.4336	0.0000	1.0000
JC69IG	13	-2714.5282	5455.5041	0.0000	1.0000
JC69I	12	-2715.7922	5455.9676	0.0000	1.0000
JC69	11	-2742.2098	5506.7436	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2571.3604	5174.7208	0.5932	0.5932
HKYIG	17	-2571.3517	5176.7034	0.2201	0.8133
GTRG	20	-2569.1231	5178.2462	0.1018	0.9151
HKYI	16	-2574.0396	5180.0792	0.0407	0.9558
GTRIG	21	-2569.1145	5180.2291	0.0378	0.9935
GTRI	20	-2571.8760	5183.7521	0.0065	1.0000
HKY	15	-2597.7896	5225.5792	0.0000	1.0000
GTR	19	-2595.2782	5228.5565	0.0000	1.0000
F81G	15	-2627.0861	5284.1722	0.0000	1.0000
F81IG	16	-2627.0922	5286.1844	0.0000	1.0000
F81I	15	-2628.3037	5286.6073	0.0000	1.0000
F81	14	-2654.0099	5336.0199	0.0000	1.0000
SYMG	17	-2658.6949	5351.3897	0.0000	1.0000
SYMIG	18	-2658.6964	5353.3929	0.0000	1.0000
SYMI	17	-2660.5699	5355.1397	0.0000	1.0000
K2PG	13	-2664.7926	5355.5852	0.0000	1.0000
K2PIG	14	-2664.7926	5357.5853	0.0000	1.0000
K2PI	13	-2666.2953	5358.5905	0.0000	1.0000
SYM	16	-2689.4324	5410.8648	0.0000	1.0000
K2P	12	-2694.3402	5412.6804	0.0000	1.0000
JC69G	12	-2714.5252	5453.0503	0.0000	1.0000
JC69IG	13	-2714.5282	5455.0564	0.0000	1.0000
JC69I	12	-2715.7922	5455.5843	0.0000	1.0000
JC69	11	-2742.2098	5506.4196	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2571.3604	5250.2056	0.9060	0.9060
HKYI	16	-2574.0396	5255.5641	0.0622	0.9682
HKYIG	17	-2571.3517	5256.9061	0.0318	1.0000
GTRG	20	-2569.1231	5272.6023	0.0000	1.0000
GTRI	20	-2571.8760	5278.1082	0.0000	1.0000
GTRIG	21	-2569.1145	5279.3030	0.0000	1.0000
HKY	15	-2597.7896	5296.3463	0.0000	1.0000
GTR	19	-2595.2782	5318.1948	0.0000	1.0000
F81G	15	-2627.0861	5354.9393	0.0000	1.0000
F81I	15	-2628.3037	5357.3744	0.0000	1.0000
F81IG	16	-2627.0922	5361.6693	0.0000	1.0000
F81	14	-2654.0099	5402.0691	0.0000	1.0000
K2PG	13	-2664.7926	5416.9166	0.0000	1.0000
K2PI	13	-2666.2953	5419.9220	0.0000	1.0000
K2PIG	14	-2664.7926	5423.6345	0.0000	1.0000
SYMG	17	-2658.6949	5431.5924	0.0000	1.0000
SYMI	17	-2660.5699	5435.3424	0.0000	1.0000
SYMIG	18	-2658.6964	5438.3134	0.0000	1.0000
K2P	12	-2694.3402	5469.2941	0.0000	1.0000
SYM	16	-2689.4324	5486.3496	0.0000	1.0000
JC69G	12	-2714.5252	5509.6640	0.0000	1.0000
JC69I	12	-2715.7922	5512.1980	0.0000	1.0000
JC69IG	13	-2714.5282	5516.3879	0.0000	1.0000
JC69	11	-2742.2098	5558.3155	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 591-tomodel-phylip/chr1_28786.phylip = 1 - 827;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28786.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 591-tomodel-phylip/chr1_28786.phylip = 1 - 827;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28786.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 591-tomodel-phylip/chr1_28786.phylip = 1 - 827;
 Partition Dummy = 1:591-tomodel-phylip/chr1_28786.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

