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

Input data from file "591-tomodel-phylip/chr20_5193.phylip" (7 taxa, 421 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	-1457.5590	2957.2181	0.5001	0.5001
GTRIG	21	-1456.9929	2958.3016	0.2909	0.7909
GTRI	20	-1459.0759	2960.2517	0.1097	0.9007
SYMG	17	-1463.3785	2962.2757	0.0399	0.9405
SYMI	17	-1463.7007	2962.9201	0.0289	0.9694
SYMIG	18	-1462.9632	2963.6280	0.0203	0.9897
HKYG	16	-1466.9934	2967.3333	0.0032	0.9929
HKYI	16	-1467.4461	2968.2388	0.0020	0.9949
K2PG	13	-1470.8981	2968.6905	0.0016	0.9965
HKYIG	17	-1466.7337	2968.9861	0.0014	0.9979
K2PI	13	-1471.1009	2969.0961	0.0013	0.9992
K2PIG	14	-1470.5915	2970.2175	0.0008	1.0000
F81I	15	-1505.7948	3042.7748	0.0000	1.0000
F81G	15	-1505.8287	3042.8425	0.0000	1.0000
F81IG	16	-1505.8288	3045.0041	0.0000	1.0000
JC69I	12	-1510.6917	3046.1481	0.0000	1.0000
JC69G	12	-1510.9814	3046.7274	0.0000	1.0000
GTR	19	-1503.7245	3047.3444	0.0000	1.0000
JC69IG	13	-1510.9814	3048.8572	0.0000	1.0000
SYM	16	-1510.3680	3054.0824	0.0000	1.0000
HKY	15	-1516.3947	3063.9745	0.0000	1.0000
K2P	12	-1520.5246	3065.8138	0.0000	1.0000
F81	14	-1551.8886	3132.8116	0.0000	1.0000
JC69	11	-1557.1805	3137.0066	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1457.5590	2955.1181	0.4958	0.4958
GTRIG	21	-1456.9929	2955.9858	0.3213	0.8171
GTRI	20	-1459.0759	2958.1517	0.1088	0.9259
SYMG	17	-1463.3785	2960.7571	0.0296	0.9555
SYMI	17	-1463.7007	2961.4015	0.0214	0.9769
SYMIG	18	-1462.9632	2961.9265	0.0165	0.9934
HKYG	16	-1466.9934	2965.9867	0.0022	0.9956
HKYI	16	-1467.4461	2966.8923	0.0014	0.9969
HKYIG	17	-1466.7337	2967.4675	0.0010	0.9980
K2PG	13	-1470.8981	2967.7962	0.0009	0.9988
K2PI	13	-1471.1009	2968.2018	0.0007	0.9996
K2PIG	14	-1470.5915	2969.1830	0.0004	1.0000
F81I	15	-1505.7948	3041.5896	0.0000	1.0000
F81G	15	-1505.8287	3041.6573	0.0000	1.0000
F81IG	16	-1505.8288	3043.6575	0.0000	1.0000
JC69I	12	-1510.6917	3045.3834	0.0000	1.0000
GTR	19	-1503.7245	3045.4491	0.0000	1.0000
JC69G	12	-1510.9814	3045.9627	0.0000	1.0000
JC69IG	13	-1510.9814	3047.9628	0.0000	1.0000
SYM	16	-1510.3680	3052.7359	0.0000	1.0000
HKY	15	-1516.3947	3062.7893	0.0000	1.0000
K2P	12	-1520.5246	3065.0491	0.0000	1.0000
F81	14	-1551.8886	3131.7771	0.0000	1.0000
JC69	11	-1557.1805	3136.3611	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1470.8981	3020.3504	0.5231	0.5231
K2PI	13	-1471.1009	3020.7560	0.4271	0.9502
K2PIG	14	-1470.5915	3025.7798	0.0346	0.9848
SYMG	17	-1463.3785	3029.4818	0.0054	0.9903
SYMI	17	-1463.7007	3030.1263	0.0039	0.9942
HKYG	16	-1466.9934	3030.6688	0.0030	0.9972
HKYI	16	-1467.4461	3031.5744	0.0019	0.9991
SYMIG	18	-1462.9632	3034.6939	0.0004	0.9995
GTRG	20	-1457.5590	3035.9707	0.0002	0.9997
HKYIG	17	-1466.7337	3036.1922	0.0002	0.9999
GTRI	20	-1459.0759	3039.0044	0.0000	1.0000
GTRIG	21	-1456.9929	3040.8811	0.0000	1.0000
JC69I	12	-1510.6917	3093.8950	0.0000	1.0000
JC69G	12	-1510.9814	3094.4743	0.0000	1.0000
JC69IG	13	-1510.9814	3100.5170	0.0000	1.0000
F81I	15	-1505.7948	3102.2291	0.0000	1.0000
F81G	15	-1505.8287	3102.2968	0.0000	1.0000
F81IG	16	-1505.8288	3108.3396	0.0000	1.0000
K2P	12	-1520.5246	3113.5607	0.0000	1.0000
SYM	16	-1510.3680	3117.4180	0.0000	1.0000
GTR	19	-1503.7245	3122.2591	0.0000	1.0000
HKY	15	-1516.3947	3123.4288	0.0000	1.0000
JC69	11	-1557.1805	3180.8300	0.0000	1.0000
F81	14	-1551.8886	3188.3740	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/chr20_5193.phylip = 1 - 421;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5193.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/chr20_5193.phylip = 1 - 421;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5193.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/chr20_5193.phylip = 1 - 421;
 Partition Dummy = 1:591-tomodel-phylip/chr20_5193.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

