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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKYG


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
HKYG	16	-1219.0124	2471.0630	0.2987	0.2987
HKYI	16	-1219.0318	2471.1017	0.2930	0.5917
GTRI	20	-1215.5374	2472.6901	0.1324	0.7241
GTRG	20	-1215.5637	2472.7427	0.1290	0.8531
HKYIG	17	-1219.0198	2473.2097	0.1021	0.9552
GTRIG	21	-1215.5702	2474.9208	0.0434	0.9987
F81I	15	-1226.9680	2484.8503	0.0003	0.9990
F81G	15	-1226.9979	2484.9102	0.0003	0.9992
HKY	15	-1227.1913	2485.2968	0.0002	0.9995
GTR	19	-1223.7501	2486.9590	0.0001	0.9996
K2PI	13	-1230.1556	2487.0018	0.0001	0.9997
F81IG	16	-1226.9961	2487.0305	0.0001	0.9998
K2PG	13	-1230.1710	2487.0327	0.0001	0.9999
K2PIG	14	-1230.1771	2489.1527	0.0000	0.9999
SYMI	17	-1227.2654	2489.7011	0.0000	1.0000
SYMG	17	-1227.3450	2489.8601	0.0000	1.0000
SYMIG	18	-1227.3449	2492.0002	0.0000	1.0000
F81	14	-1234.8873	2498.5731	0.0000	1.0000
JC69I	12	-1237.3527	2499.2963	0.0000	1.0000
JC69G	12	-1237.4065	2499.4039	0.0000	1.0000
K2P	12	-1238.0683	2500.7275	0.0000	1.0000
JC69IG	13	-1237.4148	2501.5203	0.0000	1.0000
SYM	16	-1235.0139	2503.0659	0.0000	1.0000
JC69	11	-1245.1229	2512.7449	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1219.0124	2470.0248	0.2682	0.2682
HKYI	16	-1219.0318	2470.0636	0.2630	0.5312
GTRI	20	-1215.5374	2471.0747	0.1587	0.6899
GTRG	20	-1215.5637	2471.1273	0.1545	0.8444
HKYIG	17	-1219.0198	2472.0395	0.0979	0.9424
GTRIG	21	-1215.5702	2473.1405	0.0565	0.9988
F81I	15	-1226.9680	2483.9360	0.0003	0.9991
F81G	15	-1226.9979	2483.9959	0.0002	0.9993
HKY	15	-1227.1913	2484.3826	0.0002	0.9995
GTR	19	-1223.7501	2485.5002	0.0001	0.9997
F81IG	16	-1226.9961	2485.9923	0.0001	0.9998
K2PI	13	-1230.1556	2486.3111	0.0001	0.9998
K2PG	13	-1230.1710	2486.3420	0.0001	0.9999
K2PIG	14	-1230.1771	2488.3542	0.0000	0.9999
SYMI	17	-1227.2654	2488.5309	0.0000	1.0000
SYMG	17	-1227.3450	2488.6900	0.0000	1.0000
SYMIG	18	-1227.3449	2490.6899	0.0000	1.0000
F81	14	-1234.8873	2497.7746	0.0000	1.0000
JC69I	12	-1237.3527	2498.7054	0.0000	1.0000
JC69G	12	-1237.4065	2498.8130	0.0000	1.0000
K2P	12	-1238.0683	2500.1366	0.0000	1.0000
JC69IG	13	-1237.4148	2500.8296	0.0000	1.0000
SYM	16	-1235.0139	2502.0277	0.0000	1.0000
JC69	11	-1245.1229	2512.2458	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1219.0124	2538.7195	0.4126	0.4126
HKYI	16	-1219.0318	2538.7583	0.4047	0.8172
K2PI	13	-1230.1556	2542.1256	0.0751	0.8924
K2PG	13	-1230.1710	2542.1564	0.0740	0.9664
HKYIG	17	-1219.0198	2545.0277	0.0176	0.9840
F81I	15	-1226.9680	2548.3373	0.0034	0.9874
F81G	15	-1226.9979	2548.3972	0.0033	0.9906
K2PIG	14	-1230.1771	2548.4621	0.0032	0.9938
HKY	15	-1227.1913	2548.7838	0.0027	0.9965
JC69I	12	-1237.3527	2550.2264	0.0013	0.9978
JC69G	12	-1237.4065	2550.3340	0.0012	0.9990
K2P	12	-1238.0683	2551.6576	0.0006	0.9997
F81IG	16	-1226.9961	2554.6870	0.0001	0.9998
JC69IG	13	-1237.4148	2556.6440	0.0001	0.9999
GTRI	20	-1215.5374	2556.9431	0.0000	0.9999
GTRG	20	-1215.5637	2556.9957	0.0000	0.9999
F81	14	-1234.8873	2557.8825	0.0000	1.0000
JC69	11	-1245.1229	2559.4734	0.0000	1.0000
SYMI	17	-1227.2654	2561.5190	0.0000	1.0000
SYMG	17	-1227.3450	2561.6781	0.0000	1.0000
GTRIG	21	-1215.5702	2563.3023	0.0000	1.0000
GTR	19	-1223.7501	2567.0752	0.0000	1.0000
SYMIG	18	-1227.3449	2567.9714	0.0000	1.0000
SYM	16	-1235.0139	2570.7224	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6327.phylip = 1 - 541;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6327.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6327.phylip = 1 - 541;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6327.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_6327.phylip = 1 - 541;
 Partition Dummy = 1:591-tomodel-phylip/chr7_6327.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

