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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: K2PG


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	-1839.0468	3711.0598	0.3383	0.3383
GTRG	20	-1835.2989	3712.1004	0.2010	0.5393
HKYI	16	-1839.7063	3712.3789	0.1749	0.7142
HKYIG	17	-1839.0467	3713.1824	0.1170	0.8313
GTRI	20	-1836.1378	3713.7782	0.0869	0.9182
GTRIG	21	-1835.3000	3714.2560	0.0684	0.9866
K2PG	13	-1846.2402	3719.1234	0.0060	0.9926
K2PI	13	-1846.6094	3719.8619	0.0041	0.9968
K2PIG	14	-1846.2401	3721.2235	0.0021	0.9989
SYMG	17	-1844.4133	3723.9155	0.0005	0.9994
SYMI	17	-1844.6992	3724.4874	0.0004	0.9998
SYMIG	18	-1844.4487	3726.1166	0.0002	1.0000
HKY	15	-1870.6675	3772.1861	0.0000	1.0000
GTR	19	-1867.5767	3774.5105	0.0000	1.0000
K2P	12	-1877.3842	3779.3187	0.0000	1.0000
SYM	16	-1875.6668	3784.2998	0.0000	1.0000
F81G	15	-1877.9116	3786.6743	0.0000	1.0000
F81I	15	-1877.9959	3786.8429	0.0000	1.0000
F81IG	16	-1877.9117	3788.7896	0.0000	1.0000
JC69G	12	-1882.9090	3790.3682	0.0000	1.0000
JC69I	12	-1882.9595	3790.4693	0.0000	1.0000
JC69IG	13	-1882.9090	3792.4611	0.0000	1.0000
F81	14	-1906.1188	3840.9809	0.0000	1.0000
JC69	11	-1911.1907	3844.8461	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1839.0468	3710.0935	0.3014	0.3014
GTRG	20	-1835.2989	3710.5978	0.2342	0.5356
HKYI	16	-1839.7063	3711.4126	0.1558	0.6914
HKYIG	17	-1839.0467	3712.0935	0.1109	0.8023
GTRI	20	-1836.1378	3712.2755	0.1012	0.9035
GTRIG	21	-1835.3000	3712.6000	0.0861	0.9895
K2PG	13	-1846.2402	3718.4803	0.0045	0.9941
K2PI	13	-1846.6094	3719.2188	0.0031	0.9972
K2PIG	14	-1846.2401	3720.4802	0.0017	0.9989
SYMG	17	-1844.4133	3722.8265	0.0005	0.9994
SYMI	17	-1844.6992	3723.3984	0.0004	0.9998
SYMIG	18	-1844.4487	3724.8974	0.0002	1.0000
HKY	15	-1870.6675	3771.3350	0.0000	1.0000
GTR	19	-1867.5767	3773.1533	0.0000	1.0000
K2P	12	-1877.3842	3778.7684	0.0000	1.0000
SYM	16	-1875.6668	3783.3336	0.0000	1.0000
F81G	15	-1877.9116	3785.8232	0.0000	1.0000
F81I	15	-1877.9959	3785.9918	0.0000	1.0000
F81IG	16	-1877.9117	3787.8233	0.0000	1.0000
JC69G	12	-1882.9090	3789.8180	0.0000	1.0000
JC69I	12	-1882.9595	3789.9191	0.0000	1.0000
JC69IG	13	-1882.9090	3791.8180	0.0000	1.0000
F81	14	-1906.1188	3840.2376	0.0000	1.0000
JC69	11	-1911.1907	3844.3814	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1846.2402	3775.1997	0.5315	0.5315
K2PI	13	-1846.6094	3775.9382	0.3674	0.8990
HKYG	16	-1839.0468	3779.9020	0.0506	0.9496
HKYI	16	-1839.7063	3781.2211	0.0262	0.9758
K2PIG	14	-1846.2401	3781.5626	0.0221	0.9979
HKYIG	17	-1839.0467	3786.2650	0.0021	1.0000
SYMG	17	-1844.4133	3796.9980	0.0000	1.0000
SYMI	17	-1844.6992	3797.5699	0.0000	1.0000
GTRG	20	-1835.2989	3797.8583	0.0000	1.0000
GTRI	20	-1836.1378	3799.5361	0.0000	1.0000
SYMIG	18	-1844.4487	3803.4319	0.0000	1.0000
GTRIG	21	-1835.3000	3804.2236	0.0000	1.0000
K2P	12	-1877.3842	3831.1248	0.0000	1.0000
HKY	15	-1870.6675	3836.7805	0.0000	1.0000
JC69G	12	-1882.9090	3842.1743	0.0000	1.0000
JC69I	12	-1882.9595	3842.2754	0.0000	1.0000
JC69IG	13	-1882.9090	3848.5374	0.0000	1.0000
F81G	15	-1877.9116	3851.2686	0.0000	1.0000
F81I	15	-1877.9959	3851.4372	0.0000	1.0000
SYM	16	-1875.6668	3853.1420	0.0000	1.0000
GTR	19	-1867.5767	3856.0509	0.0000	1.0000
F81IG	16	-1877.9117	3857.6318	0.0000	1.0000
JC69	11	-1911.1907	3892.3747	0.0000	1.0000
F81	14	-1906.1188	3901.3200	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/chr2_11740.phylip = 1 - 580;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11740.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/chr2_11740.phylip = 1 - 580;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11740.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_11740.phylip = 1 - 580;
 Partition Dummy = 1:591-tomodel-phylip/chr2_11740.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

