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

Input data from file "591-tomodel-phylip/chr3_9696.phylip" (7 taxa, 620 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	20	-1840.6456	3722.6935	0.3726	0.3726
HKYG	16	-1845.3385	3723.5792	0.2393	0.6119
GTRIG	21	-1840.6463	3724.8377	0.1275	0.7394
GTRI	20	-1841.9608	3725.3240	0.1000	0.8394
HKYIG	17	-1845.3385	3725.6937	0.0831	0.9225
HKYI	16	-1846.4834	3725.8690	0.0762	0.9987
HKY	15	-1851.8887	3734.5722	0.0010	0.9997
GTR	19	-1848.7412	3736.7491	0.0003	1.0000
K2PG	13	-1863.5601	3753.7209	0.0000	1.0000
K2PIG	14	-1863.5601	3755.8145	0.0000	1.0000
K2PI	13	-1864.7004	3756.0015	0.0000	1.0000
SYMG	17	-1862.5572	3760.1310	0.0000	1.0000
SYMIG	18	-1862.5572	3762.2524	0.0000	1.0000
SYMI	17	-1863.7417	3762.5000	0.0000	1.0000
K2P	12	-1869.2248	3762.9636	0.0000	1.0000
SYM	16	-1868.4647	3769.8316	0.0000	1.0000
F81G	15	-1887.7874	3806.3695	0.0000	1.0000
F81I	15	-1888.7868	3808.3684	0.0000	1.0000
F81IG	16	-1887.7874	3808.4770	0.0000	1.0000
F81	14	-1892.9401	3814.5744	0.0000	1.0000
JC69G	12	-1901.3927	3827.2994	0.0000	1.0000
JC69I	12	-1902.3887	3829.2913	0.0000	1.0000
JC69IG	13	-1901.3927	3829.3861	0.0000	1.0000
JC69	11	-1906.2980	3835.0302	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1840.6456	3721.2912	0.4029	0.4029
HKYG	16	-1845.3385	3722.6770	0.2015	0.6044
GTRIG	21	-1840.6463	3723.2925	0.1481	0.7525
GTRI	20	-1841.9608	3723.9217	0.1081	0.8606
HKYIG	17	-1845.3385	3724.6771	0.0741	0.9348
HKYI	16	-1846.4834	3724.9669	0.0641	0.9989
HKY	15	-1851.8887	3733.7775	0.0008	0.9997
GTR	19	-1848.7412	3735.4825	0.0003	1.0000
K2PG	13	-1863.5601	3753.1202	0.0000	1.0000
K2PIG	14	-1863.5601	3755.1203	0.0000	1.0000
K2PI	13	-1864.7004	3755.4008	0.0000	1.0000
SYMG	17	-1862.5572	3759.1143	0.0000	1.0000
SYMIG	18	-1862.5572	3761.1143	0.0000	1.0000
SYMI	17	-1863.7417	3761.4834	0.0000	1.0000
K2P	12	-1869.2248	3762.4496	0.0000	1.0000
SYM	16	-1868.4647	3768.9294	0.0000	1.0000
F81G	15	-1887.7874	3805.5748	0.0000	1.0000
F81I	15	-1888.7868	3807.5737	0.0000	1.0000
F81IG	16	-1887.7874	3807.5749	0.0000	1.0000
F81	14	-1892.9401	3813.8801	0.0000	1.0000
JC69G	12	-1901.3927	3826.7854	0.0000	1.0000
JC69I	12	-1902.3887	3828.7773	0.0000	1.0000
JC69IG	13	-1901.3927	3828.7854	0.0000	1.0000
JC69	11	-1906.2980	3834.5960	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1845.3385	3793.5526	0.7170	0.7170
HKYI	16	-1846.4834	3795.8424	0.2282	0.9452
HKYIG	17	-1845.3385	3799.9823	0.0288	0.9740
HKY	15	-1851.8887	3800.2233	0.0255	0.9995
GTRG	20	-1840.6456	3809.8856	0.0002	0.9997
K2PG	13	-1863.5601	3810.7066	0.0001	0.9999
GTRI	20	-1841.9608	3812.5160	0.0001	0.9999
K2PI	13	-1864.7004	3812.9872	0.0000	1.0000
K2P	12	-1869.2248	3815.6063	0.0000	1.0000
GTRIG	21	-1840.6463	3816.3166	0.0000	1.0000
K2PIG	14	-1863.5601	3817.1363	0.0000	1.0000
GTR	19	-1848.7412	3819.6471	0.0000	1.0000
SYMG	17	-1862.5572	3834.4196	0.0000	1.0000
SYMI	17	-1863.7417	3836.7887	0.0000	1.0000
SYM	16	-1868.4647	3839.8049	0.0000	1.0000
SYMIG	18	-1862.5572	3840.8493	0.0000	1.0000
F81G	15	-1887.7874	3872.0206	0.0000	1.0000
F81I	15	-1888.7868	3874.0195	0.0000	1.0000
F81	14	-1892.9401	3875.8962	0.0000	1.0000
F81IG	16	-1887.7874	3878.4504	0.0000	1.0000
JC69G	12	-1901.3927	3879.9421	0.0000	1.0000
JC69I	12	-1902.3887	3881.9339	0.0000	1.0000
JC69	11	-1906.2980	3883.3229	0.0000	1.0000
JC69IG	13	-1901.3927	3886.3718	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 591-tomodel-phylip/chr3_9696.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9696.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 591-tomodel-phylip/chr3_9696.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9696.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 591-tomodel-phylip/chr3_9696.phylip = 1 - 620;
 Partition Dummy = 1:591-tomodel-phylip/chr3_9696.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

