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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-2774.3145	5589.6976	0.3016	0.3016
GTRI	20	-2775.1079	5591.2845	0.1364	0.4380
GTRIG	21	-2774.2670	5591.7111	0.1102	0.5482
K2PG	13	-2782.6756	5591.8101	0.1049	0.6531
K2PI	13	-2783.0662	5592.5914	0.0710	0.7240
SYMG	17	-2779.0651	5592.9059	0.0606	0.7847
HKYG	16	-2780.1992	5593.0869	0.0554	0.8400
K2PIG	14	-2782.5644	5593.6591	0.0416	0.8816
SYMI	17	-2779.4427	5593.6611	0.0416	0.9232
HKYI	16	-2780.7681	5594.2247	0.0314	0.9546
SYMIG	18	-2778.9345	5594.7371	0.0243	0.9788
HKYIG	17	-2780.1181	5595.0118	0.0212	1.0000
GTR	19	-2797.7008	5634.3673	0.0000	1.0000
K2P	12	-2806.6417	5637.6764	0.0000	1.0000
SYM	16	-2802.5767	5637.8420	0.0000	1.0000
HKY	15	-2803.7481	5638.1030	0.0000	1.0000
JC69G	12	-2834.6665	5693.7259	0.0000	1.0000
JC69I	12	-2834.8903	5694.1735	0.0000	1.0000
JC69IG	13	-2834.5248	5695.5087	0.0000	1.0000
F81G	15	-2832.5709	5695.7486	0.0000	1.0000
F81I	15	-2832.8719	5696.3506	0.0000	1.0000
F81IG	16	-2832.4556	5697.5998	0.0000	1.0000
JC69	11	-2856.9644	5736.2609	0.0000	1.0000
F81	14	-2854.6600	5737.8504	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2774.3145	5588.6289	0.3292	0.3292
GTRI	20	-2775.1079	5590.2158	0.1489	0.4782
GTRIG	21	-2774.2670	5590.5340	0.1270	0.6052
K2PG	13	-2782.6756	5591.3511	0.0844	0.6896
SYMG	17	-2779.0651	5592.1303	0.0572	0.7468
K2PI	13	-2783.0662	5592.1324	0.0571	0.8039
HKYG	16	-2780.1992	5592.3983	0.0500	0.8539
SYMI	17	-2779.4427	5592.8854	0.0392	0.8931
K2PIG	14	-2782.5644	5593.1288	0.0347	0.9278
HKYI	16	-2780.7681	5593.5361	0.0283	0.9561
SYMIG	18	-2778.9345	5593.8691	0.0240	0.9801
HKYIG	17	-2780.1181	5594.2361	0.0199	1.0000
GTR	19	-2797.7008	5633.4016	0.0000	1.0000
SYM	16	-2802.5767	5637.1534	0.0000	1.0000
K2P	12	-2806.6417	5637.2835	0.0000	1.0000
HKY	15	-2803.7481	5637.4962	0.0000	1.0000
JC69G	12	-2834.6665	5693.3330	0.0000	1.0000
JC69I	12	-2834.8903	5693.7805	0.0000	1.0000
JC69IG	13	-2834.5248	5695.0497	0.0000	1.0000
F81G	15	-2832.5709	5695.1417	0.0000	1.0000
F81I	15	-2832.8719	5695.7438	0.0000	1.0000
F81IG	16	-2832.4556	5696.9112	0.0000	1.0000
JC69	11	-2856.9644	5735.9288	0.0000	1.0000
F81	14	-2854.6600	5737.3201	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2782.6756	5652.3643	0.5824	0.5824
K2PI	13	-2783.0662	5653.1456	0.3941	0.9765
K2PIG	14	-2782.5644	5658.8354	0.0229	0.9995
HKYG	16	-2780.1992	5667.4915	0.0003	0.9998
HKYI	16	-2780.7681	5668.6293	0.0002	0.9999
SYMG	17	-2779.0651	5671.9168	0.0000	1.0000
SYMI	17	-2779.4427	5672.6719	0.0000	1.0000
HKYIG	17	-2780.1181	5674.0226	0.0000	1.0000
SYMIG	18	-2778.9345	5678.3489	0.0000	1.0000
GTRG	20	-2774.3145	5682.4954	0.0000	1.0000
GTRI	20	-2775.1079	5684.0822	0.0000	1.0000
GTRIG	21	-2774.2670	5689.0938	0.0000	1.0000
K2P	12	-2806.6417	5693.6033	0.0000	1.0000
HKY	15	-2803.7481	5707.8960	0.0000	1.0000
SYM	16	-2802.5767	5712.2466	0.0000	1.0000
GTR	19	-2797.7008	5722.5748	0.0000	1.0000
JC69G	12	-2834.6665	5749.6529	0.0000	1.0000
JC69I	12	-2834.8903	5750.1004	0.0000	1.0000
JC69IG	13	-2834.5248	5756.0629	0.0000	1.0000
F81G	15	-2832.5709	5765.5416	0.0000	1.0000
F81I	15	-2832.8719	5766.1437	0.0000	1.0000
F81IG	16	-2832.4556	5772.0043	0.0000	1.0000
JC69	11	-2856.9644	5787.5553	0.0000	1.0000
F81	14	-2854.6600	5803.0266	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/chr4_11159.phylip = 1 - 807;
 Partition Dummy = 1:591-tomodel-phylip/chr4_11159.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/chr4_11159.phylip = 1 - 807;
 Partition Dummy = 1:591-tomodel-phylip/chr4_11159.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 (K2PG)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr4_11159.phylip = 1 - 807;
 Partition Dummy = 1:591-tomodel-phylip/chr4_11159.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

