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

Input data from file "683-tomodel-phylip/chr17_1772.phylip" (7 taxa, 713 characters)

Minimum AIC  model: SYMG
Minimum AICc model: SYMG
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
SYMG	17	-3153.0574	6340.9954	0.4169	0.4169
GTRG	20	-3150.4398	6342.0935	0.2408	0.6577
SYMIG	18	-3153.0732	6343.1320	0.1433	0.8010
GTRIG	21	-3150.4173	6344.1717	0.0852	0.8862
SYMI	17	-3154.9829	6344.8464	0.0608	0.9469
GTRI	20	-3152.0424	6345.2987	0.0485	0.9954
K2PG	13	-3162.5288	6351.5784	0.0021	0.9975
HKYG	16	-3160.2937	6353.3691	0.0009	0.9984
K2PIG	14	-3162.5208	6353.6433	0.0007	0.9991
K2PI	13	-3164.2696	6355.0600	0.0004	0.9995
HKYIG	17	-3160.2657	6355.4120	0.0003	0.9998
HKYI	16	-3161.8348	6356.4512	0.0002	1.0000
SYM	16	-3184.2716	6401.3249	0.0000	1.0000
GTR	19	-3181.6172	6402.3310	0.0000	1.0000
K2P	12	-3193.3680	6411.1818	0.0000	1.0000
HKY	15	-3190.9081	6412.5048	0.0000	1.0000
JC69I	12	-3275.8255	6576.0966	0.0000	1.0000
JC69G	12	-3275.9962	6576.4381	0.0000	1.0000
JC69IG	13	-3275.5185	6577.5578	0.0000	1.0000
F81I	15	-3273.5826	6577.8539	0.0000	1.0000
F81G	15	-3273.9394	6578.5675	0.0000	1.0000
F81IG	16	-3273.3509	6579.4834	0.0000	1.0000
JC69	11	-3298.9250	6620.2266	0.0000	1.0000
F81	14	-3296.9452	6622.4920	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMG	17	-3153.0574	6340.1149	0.3855	0.3855
GTRG	20	-3150.4398	6340.8796	0.2630	0.6485
SYMIG	18	-3153.0732	6342.1464	0.1396	0.7882
GTRIG	21	-3150.4173	6342.8345	0.0990	0.8871
SYMI	17	-3154.9829	6343.9658	0.0562	0.9433
GTRI	20	-3152.0424	6344.0848	0.0530	0.9963
K2PG	13	-3162.5288	6351.0577	0.0016	0.9979
HKYG	16	-3160.2937	6352.5875	0.0008	0.9987
K2PIG	14	-3162.5208	6353.0416	0.0006	0.9993
HKYIG	17	-3160.2657	6354.5314	0.0003	0.9996
K2PI	13	-3164.2696	6354.5392	0.0003	0.9998
HKYI	16	-3161.8348	6355.6696	0.0002	1.0000
SYM	16	-3184.2716	6400.5433	0.0000	1.0000
GTR	19	-3181.6172	6401.2344	0.0000	1.0000
K2P	12	-3193.3680	6410.7361	0.0000	1.0000
HKY	15	-3190.9081	6411.8162	0.0000	1.0000
JC69I	12	-3275.8255	6575.6509	0.0000	1.0000
JC69G	12	-3275.9962	6575.9924	0.0000	1.0000
JC69IG	13	-3275.5185	6577.0371	0.0000	1.0000
F81I	15	-3273.5826	6577.1653	0.0000	1.0000
F81G	15	-3273.9394	6577.8788	0.0000	1.0000
F81IG	16	-3273.3509	6578.7018	0.0000	1.0000
JC69	11	-3298.9250	6619.8500	0.0000	1.0000
F81	14	-3296.9452	6621.8903	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3162.5288	6410.4609	0.8039	0.8039
K2PI	13	-3164.2696	6413.9425	0.1410	0.9449
K2PIG	14	-3162.5208	6417.0143	0.0303	0.9752
SYMG	17	-3153.0574	6417.7960	0.0205	0.9957
SYMI	17	-3154.9829	6421.6470	0.0030	0.9987
SYMIG	18	-3153.0732	6424.3971	0.0008	0.9995
HKYG	16	-3160.2937	6425.6992	0.0004	0.9999
HKYI	16	-3161.8348	6428.7813	0.0001	1.0000
HKYIG	17	-3160.2657	6432.2126	0.0000	1.0000
GTRG	20	-3150.4398	6432.2692	0.0000	1.0000
GTRI	20	-3152.0424	6435.4744	0.0000	1.0000
GTRIG	21	-3150.4173	6438.7936	0.0000	1.0000
K2P	12	-3193.3680	6465.5699	0.0000	1.0000
SYM	16	-3184.2716	6473.6550	0.0000	1.0000
HKY	15	-3190.9081	6480.3584	0.0000	1.0000
GTR	19	-3181.6172	6488.0545	0.0000	1.0000
JC69I	12	-3275.8255	6630.4847	0.0000	1.0000
JC69G	12	-3275.9962	6630.8262	0.0000	1.0000
JC69IG	13	-3275.5185	6636.4403	0.0000	1.0000
F81I	15	-3273.5826	6645.7075	0.0000	1.0000
F81G	15	-3273.9394	6646.4210	0.0000	1.0000
F81IG	16	-3273.3509	6651.8135	0.0000	1.0000
JC69	11	-3298.9250	6670.1143	0.0000	1.0000
F81	14	-3296.9452	6685.8630	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (SYMG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr17_1772.phylip = 1 - 713;
 Partition Dummy = 1:683-tomodel-phylip/chr17_1772.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=Fixed(Equal);
END;


[Mrbayes block for the best AICc model (SYMG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr17_1772.phylip = 1 - 713;
 Partition Dummy = 1:683-tomodel-phylip/chr17_1772.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=Fixed(Equal);
END;


[Mrbayes block for the best BIC model (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr17_1772.phylip = 1 - 713;
 Partition Dummy = 1:683-tomodel-phylip/chr17_1772.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

