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

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

Minimum AIC  model: K2PI
Minimum AICc model: K2PI
Minimum BIC  model: K2PI


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
K2PI	13	-602.4295	1232.4147	0.2706	0.2706
K2PG	13	-602.4640	1232.4836	0.2615	0.5321
K2PIG	14	-602.5032	1234.8089	0.0817	0.6139
SYMI	17	-599.1486	1234.9581	0.0759	0.6897
SYMG	17	-599.1527	1234.9663	0.0756	0.7653
HKYI	16	-600.5453	1235.4457	0.0595	0.8247
HKYG	16	-600.5715	1235.4979	0.0579	0.8827
SYMIG	18	-599.1562	1237.2992	0.0235	0.9062
GTRG	20	-596.8451	1237.3907	0.0225	0.9287
GTRI	20	-596.8820	1237.4644	0.0217	0.9504
HKYIG	17	-600.6003	1237.8614	0.0178	0.9681
JC69G	12	-607.1669	1239.6615	0.0072	0.9753
GTRIG	21	-596.8342	1239.7568	0.0069	0.9822
JC69I	12	-607.3057	1239.9390	0.0063	0.9885
JC69IG	13	-607.0812	1241.7180	0.0026	0.9911
K2P	12	-608.3319	1241.9915	0.0023	0.9934
SYM	16	-604.0054	1242.3659	0.0019	0.9952
F81G	15	-605.4137	1242.8963	0.0014	0.9967
F81I	15	-605.5701	1243.2091	0.0012	0.9979
GTR	19	-601.2859	1243.9051	0.0009	0.9988
HKY	15	-606.1684	1244.4057	0.0007	0.9994
F81IG	16	-605.3243	1245.0036	0.0005	0.9999
JC69	11	-613.0648	1249.2483	0.0001	1.0000
F81	14	-611.1560	1252.1146	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PI	13	-602.4295	1230.8591	0.2047	0.2047
K2PG	13	-602.4640	1230.9280	0.1978	0.4025
SYMI	17	-599.1486	1232.2972	0.0997	0.5022
SYMG	17	-599.1527	1232.3054	0.0993	0.6016
K2PIG	14	-602.5032	1233.0063	0.0700	0.6715
HKYI	16	-600.5453	1233.0907	0.0671	0.7386
HKYG	16	-600.5715	1233.1429	0.0653	0.8040
GTRG	20	-596.8451	1233.6903	0.0497	0.8537
GTRI	20	-596.8820	1233.7640	0.0479	0.9016
SYMIG	18	-599.1562	1234.3123	0.0364	0.9380
HKYIG	17	-600.6003	1235.2006	0.0234	0.9613
GTRIG	21	-596.8342	1235.6683	0.0185	0.9798
JC69G	12	-607.1669	1238.3338	0.0049	0.9847
JC69I	12	-607.3057	1238.6113	0.0042	0.9890
SYM	16	-604.0054	1240.0109	0.0021	0.9911
JC69IG	13	-607.0812	1240.1624	0.0020	0.9930
GTR	19	-601.2859	1240.5718	0.0016	0.9946
K2P	12	-608.3319	1240.6638	0.0015	0.9961
F81G	15	-605.4137	1240.8274	0.0014	0.9975
F81I	15	-605.5701	1241.1401	0.0012	0.9987
HKY	15	-606.1684	1242.3367	0.0007	0.9994
F81IG	16	-605.3243	1242.6486	0.0006	1.0000
JC69	11	-613.0648	1248.1296	0.0000	1.0000
F81	14	-611.1560	1250.3120	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-602.4295	1276.5337	0.4261	0.4261
K2PG	13	-602.4640	1276.6026	0.4117	0.8379
JC69G	12	-607.1669	1280.4950	0.0588	0.8967
JC69I	12	-607.3057	1280.7725	0.0512	0.9478
K2PIG	14	-602.5032	1282.1943	0.0251	0.9730
K2P	12	-608.3319	1282.8250	0.0183	0.9913
JC69IG	13	-607.0812	1285.8370	0.0041	0.9954
JC69	11	-613.0648	1286.7773	0.0025	0.9979
HKYI	16	-600.5453	1289.3055	0.0007	0.9986
HKYG	16	-600.5715	1289.3578	0.0007	0.9993
SYMI	17	-599.1486	1292.0255	0.0002	0.9995
SYMG	17	-599.1527	1292.0337	0.0002	0.9997
F81G	15	-605.4137	1293.5288	0.0001	0.9998
F81I	15	-605.5701	1293.8416	0.0001	0.9999
HKYIG	17	-600.6003	1294.9288	0.0000	0.9999
HKY	15	-606.1684	1295.0382	0.0000	1.0000
SYM	16	-604.0054	1296.2258	0.0000	1.0000
SYMIG	18	-599.1562	1297.5541	0.0000	1.0000
F81IG	16	-605.3243	1298.8635	0.0000	1.0000
F81	14	-611.1560	1299.5000	0.0000	1.0000
GTRG	20	-596.8451	1303.9589	0.0000	1.0000
GTRI	20	-596.8820	1304.0326	0.0000	1.0000
GTR	19	-601.2859	1307.3269	0.0000	1.0000
GTRIG	21	-596.8342	1309.4503	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_29129.phylip = 1 - 248;
 Partition Dummy = 1:591-tomodel-phylip/chr2_29129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_29129.phylip = 1 - 248;
 Partition Dummy = 1:591-tomodel-phylip/chr2_29129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_29129.phylip = 1 - 248;
 Partition Dummy = 1:591-tomodel-phylip/chr2_29129.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

