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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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	-2964.8216	5970.6972	0.6528	0.6528
GTRIG	21	-2964.8158	5972.7925	0.2290	0.8818
GTRI	20	-2966.5306	5974.1151	0.1182	1.0000
HKYG	16	-2980.1380	5992.9552	0.0000	1.0000
HKYIG	17	-2980.1423	5995.0495	0.0000	1.0000
HKYI	16	-2981.4174	5995.5140	0.0000	1.0000
GTR	19	-2988.4394	6015.8311	0.0000	1.0000
HKY	15	-3001.8629	6034.3242	0.0000	1.0000
SYMG	17	-3003.3089	6041.3828	0.0000	1.0000
K2PG	13	-3007.8976	6042.2479	0.0000	1.0000
SYMIG	18	-3003.3130	6043.4820	0.0000	1.0000
K2PIG	14	-3007.9029	6044.3289	0.0000	1.0000
SYMI	17	-3004.8238	6044.4126	0.0000	1.0000
K2PI	13	-3009.5314	6045.5156	0.0000	1.0000
SYM	16	-3024.9504	6082.5800	0.0000	1.0000
K2P	12	-3030.0805	6084.5485	0.0000	1.0000
F81G	15	-3063.3814	6157.3612	0.0000	1.0000
F81I	15	-3063.8209	6158.2403	0.0000	1.0000
F81IG	16	-3063.3937	6159.4665	0.0000	1.0000
JC69G	12	-3082.9513	6190.2902	0.0000	1.0000
F81	14	-3081.3704	6191.2638	0.0000	1.0000
JC69I	12	-3083.6516	6191.6907	0.0000	1.0000
JC69IG	13	-3082.9641	6192.3809	0.0000	1.0000
JC69	11	-3101.5775	6225.4826	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2964.8216	5969.6432	0.6447	0.6447
GTRIG	21	-2964.8158	5971.6317	0.2385	0.8833
GTRI	20	-2966.5306	5973.0612	0.1167	1.0000
HKYG	16	-2980.1380	5992.2760	0.0000	1.0000
HKYIG	17	-2980.1423	5994.2845	0.0000	1.0000
HKYI	16	-2981.4174	5994.8348	0.0000	1.0000
GTR	19	-2988.4394	6014.8788	0.0000	1.0000
HKY	15	-3001.8629	6033.7257	0.0000	1.0000
SYMG	17	-3003.3089	6040.6178	0.0000	1.0000
K2PG	13	-3007.8976	6041.7952	0.0000	1.0000
SYMIG	18	-3003.3130	6042.6260	0.0000	1.0000
SYMI	17	-3004.8238	6043.6476	0.0000	1.0000
K2PIG	14	-3007.9029	6043.8059	0.0000	1.0000
K2PI	13	-3009.5314	6045.0629	0.0000	1.0000
SYM	16	-3024.9504	6081.9009	0.0000	1.0000
K2P	12	-3030.0805	6084.1609	0.0000	1.0000
F81G	15	-3063.3814	6156.7627	0.0000	1.0000
F81I	15	-3063.8209	6157.6418	0.0000	1.0000
F81IG	16	-3063.3937	6158.7873	0.0000	1.0000
JC69G	12	-3082.9513	6189.9026	0.0000	1.0000
F81	14	-3081.3704	6190.7408	0.0000	1.0000
JC69I	12	-3083.6516	6191.3031	0.0000	1.0000
JC69IG	13	-3082.9641	6191.9281	0.0000	1.0000
JC69	11	-3101.5775	6225.1550	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2964.8216	6063.7804	0.7082	0.7082
GTRI	20	-2966.5306	6067.1984	0.1282	0.8364
HKYG	16	-2980.1380	6067.5858	0.1056	0.9420
HKYI	16	-2981.4174	6070.1446	0.0294	0.9714
GTRIG	21	-2964.8158	6070.4758	0.0249	0.9963
HKYIG	17	-2980.1423	6074.3012	0.0037	1.0000
K2PG	13	-3007.8976	6102.9844	0.0000	1.0000
GTR	19	-2988.4394	6104.3091	0.0000	1.0000
HKY	15	-3001.8629	6104.3287	0.0000	1.0000
K2PI	13	-3009.5314	6106.2521	0.0000	1.0000
K2PIG	14	-3007.9029	6109.7020	0.0000	1.0000
SYMG	17	-3003.3089	6120.6345	0.0000	1.0000
SYMI	17	-3004.8238	6123.6643	0.0000	1.0000
SYMIG	18	-3003.3130	6127.3495	0.0000	1.0000
K2P	12	-3030.0805	6140.6433	0.0000	1.0000
SYM	16	-3024.9504	6157.2107	0.0000	1.0000
F81G	15	-3063.3814	6227.3657	0.0000	1.0000
F81I	15	-3063.8209	6228.2448	0.0000	1.0000
F81IG	16	-3063.3937	6234.0971	0.0000	1.0000
JC69G	12	-3082.9513	6246.3849	0.0000	1.0000
JC69I	12	-3083.6516	6247.7854	0.0000	1.0000
JC69IG	13	-3082.9641	6253.1173	0.0000	1.0000
F81	14	-3081.3704	6256.6369	0.0000	1.0000
JC69	11	-3101.5775	6276.9305	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_5887.phylip = 1 - 818;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5887.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_5887.phylip = 1 - 818;
 Partition Dummy = 1:683-tomodel-phylip/chr3_5887.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 (GTRG)]

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


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

