
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:15 04/04/2011

Input data from file "683-tomodel-phylip/chr3_17735.phylip" (7 taxa, 626 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	-1824.6400	3690.6685	0.4436	0.4436
GTRI	20	-1824.8460	3691.0805	0.3610	0.8046
GTRIG	21	-1824.6078	3692.7454	0.1570	0.9616
HKYG	16	-1832.0257	3696.9446	0.0192	0.9809
HKYI	16	-1832.5402	3697.9736	0.0115	0.9924
HKYIG	17	-1832.0237	3699.0539	0.0067	0.9991
GTR	19	-1832.3234	3703.9009	0.0006	0.9997
SYMG	17	-1835.8794	3706.7653	0.0001	0.9998
SYMI	17	-1836.2323	3707.4712	0.0001	0.9999
SYMIG	18	-1835.8696	3708.8661	0.0000	1.0000
HKY	15	-1839.3264	3709.4396	0.0000	1.0000
K2PG	13	-1843.4581	3713.5109	0.0000	1.0000
K2PI	13	-1843.9693	3714.5334	0.0000	1.0000
K2PIG	14	-1843.4584	3715.6041	0.0000	1.0000
SYM	16	-1842.8328	3718.5589	0.0000	1.0000
K2P	12	-1850.2731	3725.0551	0.0000	1.0000
F81G	15	-1862.4704	3755.7276	0.0000	1.0000
F81I	15	-1862.7534	3756.2936	0.0000	1.0000
F81IG	16	-1862.4611	3757.8155	0.0000	1.0000
F81	14	-1869.0802	3766.8479	0.0000	1.0000
JC69G	12	-1872.5513	3769.6116	0.0000	1.0000
JC69I	12	-1872.8749	3770.2588	0.0000	1.0000
JC69IG	13	-1872.5476	3771.6901	0.0000	1.0000
JC69	11	-1878.9865	3780.4030	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1824.6400	3689.2801	0.4420	0.4420
GTRI	20	-1824.8460	3689.6921	0.3598	0.8018
GTRIG	21	-1824.6078	3691.2156	0.1679	0.9697
HKYG	16	-1832.0257	3696.0514	0.0150	0.9847
HKYI	16	-1832.5402	3697.0803	0.0089	0.9937
HKYIG	17	-1832.0237	3698.0473	0.0055	0.9992
GTR	19	-1832.3234	3702.6468	0.0006	0.9997
SYMG	17	-1835.8794	3705.7588	0.0001	0.9998
SYMI	17	-1836.2323	3706.4646	0.0001	0.9999
SYMIG	18	-1835.8696	3707.7392	0.0000	1.0000
HKY	15	-1839.3264	3708.6527	0.0000	1.0000
K2PG	13	-1843.4581	3712.9161	0.0000	1.0000
K2PI	13	-1843.9693	3713.9386	0.0000	1.0000
K2PIG	14	-1843.4584	3714.9167	0.0000	1.0000
SYM	16	-1842.8328	3717.6657	0.0000	1.0000
K2P	12	-1850.2731	3724.5462	0.0000	1.0000
F81G	15	-1862.4704	3754.9407	0.0000	1.0000
F81I	15	-1862.7534	3755.5068	0.0000	1.0000
F81IG	16	-1862.4611	3756.9222	0.0000	1.0000
F81	14	-1869.0802	3766.1605	0.0000	1.0000
JC69G	12	-1872.5513	3769.1027	0.0000	1.0000
JC69I	12	-1872.8749	3769.7499	0.0000	1.0000
JC69IG	13	-1872.5476	3771.0953	0.0000	1.0000
JC69	11	-1878.9865	3779.9731	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1832.0257	3767.0810	0.5136	0.5136
HKYI	16	-1832.5402	3768.1099	0.3071	0.8207
K2PG	13	-1843.4581	3770.6277	0.0872	0.9079
K2PI	13	-1843.9693	3771.6502	0.0523	0.9602
HKYIG	17	-1832.0237	3773.5163	0.0206	0.9807
HKY	15	-1839.3264	3775.2430	0.0087	0.9894
K2PIG	14	-1843.4584	3777.0676	0.0035	0.9929
K2P	12	-1850.2731	3777.8184	0.0024	0.9953
GTRG	20	-1824.6400	3778.0671	0.0021	0.9974
GTRI	20	-1824.8460	3778.4791	0.0017	0.9991
SYMG	17	-1835.8794	3781.2277	0.0004	0.9996
SYMI	17	-1836.2323	3781.9336	0.0003	0.9999
GTRIG	21	-1824.6078	3784.4420	0.0001	0.9999
GTR	19	-1832.3234	3786.9944	0.0000	1.0000
SYMIG	18	-1835.8696	3787.6475	0.0000	1.0000
SYM	16	-1842.8328	3788.6953	0.0000	1.0000
F81G	15	-1862.4704	3821.5310	0.0000	1.0000
F81I	15	-1862.7534	3822.0970	0.0000	1.0000
JC69G	12	-1872.5513	3822.3749	0.0000	1.0000
JC69I	12	-1872.8749	3823.0221	0.0000	1.0000
F81IG	16	-1862.4611	3827.9518	0.0000	1.0000
F81	14	-1869.0802	3828.3114	0.0000	1.0000
JC69	11	-1878.9865	3828.8059	0.0000	1.0000
JC69IG	13	-1872.5476	3828.8068	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 683-tomodel-phylip/chr3_17735.phylip = 1 - 626;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17735.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 683-tomodel-phylip/chr3_17735.phylip = 1 - 626;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17735.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 683-tomodel-phylip/chr3_17735.phylip = 1 - 626;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17735.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

