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

Input data from file "591-tomodel-phylip/chrZ_1330.phylip" (7 taxa, 767 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	-1815.3552	3671.8365	0.3005	0.3005
GTRI	20	-1815.6499	3672.4257	0.2238	0.5243
HKYG	16	-1820.1277	3672.9807	0.1696	0.6939
HKYI	16	-1820.3287	3673.3828	0.1387	0.8326
GTRIG	21	-1815.3341	3673.9085	0.1066	0.9393
HKYIG	17	-1820.1089	3675.0349	0.0607	1.0000
GTR	19	-1827.2986	3693.6147	0.0000	1.0000
HKY	15	-1832.6786	3695.9963	0.0000	1.0000
F81G	15	-1849.0272	3728.6936	0.0000	1.0000
F81I	15	-1849.3393	3729.3177	0.0000	1.0000
F81IG	16	-1848.9769	3730.6791	0.0000	1.0000
K2PG	13	-1858.5345	3743.5524	0.0000	1.0000
SYMG	17	-1854.4915	3743.8000	0.0000	1.0000
K2PI	13	-1858.8874	3744.2582	0.0000	1.0000
SYMI	17	-1854.8869	3744.5909	0.0000	1.0000
K2PIG	14	-1858.4162	3745.3910	0.0000	1.0000
SYMIG	18	-1854.3474	3745.6093	0.0000	1.0000
F81	14	-1861.0624	3750.6834	0.0000	1.0000
K2P	12	-1872.5030	3769.4198	0.0000	1.0000
SYM	16	-1868.6508	3770.0270	0.0000	1.0000
JC69G	12	-1884.6504	3793.7146	0.0000	1.0000
JC69I	12	-1885.0381	3794.4901	0.0000	1.0000
JC69IG	13	-1884.5455	3795.5743	0.0000	1.0000
JC69	11	-1897.6287	3817.6071	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1815.3552	3670.7105	0.3192	0.3192
GTRI	20	-1815.6499	3671.2997	0.2377	0.5569
HKYG	16	-1820.1277	3672.2554	0.1474	0.7043
HKYI	16	-1820.3287	3672.6575	0.1206	0.8248
GTRIG	21	-1815.3341	3672.6683	0.1199	0.9447
HKYIG	17	-1820.1089	3674.2178	0.0553	1.0000
GTR	19	-1827.2986	3692.5973	0.0000	1.0000
HKY	15	-1832.6786	3695.3572	0.0000	1.0000
F81G	15	-1849.0272	3728.0545	0.0000	1.0000
F81I	15	-1849.3393	3728.6785	0.0000	1.0000
F81IG	16	-1848.9769	3729.9537	0.0000	1.0000
SYMG	17	-1854.4915	3742.9829	0.0000	1.0000
K2PG	13	-1858.5345	3743.0690	0.0000	1.0000
SYMI	17	-1854.8869	3743.7738	0.0000	1.0000
K2PI	13	-1858.8874	3743.7748	0.0000	1.0000
SYMIG	18	-1854.3474	3744.6949	0.0000	1.0000
K2PIG	14	-1858.4162	3744.8325	0.0000	1.0000
F81	14	-1861.0624	3750.1249	0.0000	1.0000
K2P	12	-1872.5030	3769.0060	0.0000	1.0000
SYM	16	-1868.6508	3769.3016	0.0000	1.0000
JC69G	12	-1884.6504	3793.3008	0.0000	1.0000
JC69I	12	-1885.0381	3794.0763	0.0000	1.0000
JC69IG	13	-1884.5455	3795.0909	0.0000	1.0000
JC69	11	-1897.6287	3817.2574	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1820.1277	3746.5352	0.5390	0.5390
HKYI	16	-1820.3287	3746.9373	0.4409	0.9799
HKYIG	17	-1820.1089	3753.1401	0.0198	0.9998
GTRG	20	-1815.3552	3763.5602	0.0001	0.9999
GTRI	20	-1815.6499	3764.1495	0.0001	0.9999
HKY	15	-1832.6786	3764.9945	0.0001	1.0000
GTRIG	21	-1815.3341	3770.1605	0.0000	1.0000
GTR	19	-1827.2986	3780.8045	0.0000	1.0000
F81G	15	-1849.0272	3797.6918	0.0000	1.0000
F81I	15	-1849.3393	3798.3158	0.0000	1.0000
K2PG	13	-1858.5345	3803.4213	0.0000	1.0000
K2PI	13	-1858.8874	3804.1271	0.0000	1.0000
F81IG	16	-1848.9769	3804.2335	0.0000	1.0000
K2PIG	14	-1858.4162	3809.8273	0.0000	1.0000
F81	14	-1861.0624	3815.1197	0.0000	1.0000
SYMG	17	-1854.4915	3821.9052	0.0000	1.0000
SYMI	17	-1854.8869	3822.6961	0.0000	1.0000
K2P	12	-1872.5030	3824.7158	0.0000	1.0000
SYMIG	18	-1854.3474	3828.2596	0.0000	1.0000
SYM	16	-1868.6508	3843.5814	0.0000	1.0000
JC69G	12	-1884.6504	3849.0107	0.0000	1.0000
JC69I	12	-1885.0381	3849.7861	0.0000	1.0000
JC69IG	13	-1884.5455	3855.4432	0.0000	1.0000
JC69	11	-1897.6287	3868.3248	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/chrZ_1330.phylip = 1 - 767;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_1330.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/chrZ_1330.phylip = 1 - 767;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_1330.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/chrZ_1330.phylip = 1 - 767;
 Partition Dummy = 1:591-tomodel-phylip/chrZ_1330.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

