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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	20	-1217.0552	2475.9446	0.5299	0.5299
GTRIG	21	-1216.7671	2477.5560	0.2367	0.7666
GTRI	20	-1217.9205	2477.6750	0.2231	0.9897
HKYG	16	-1226.1080	2485.3935	0.0047	0.9944
HKYI	16	-1226.7164	2486.6104	0.0026	0.9970
HKYIG	17	-1225.8468	2487.0212	0.0021	0.9991
SYMG	17	-1227.3579	2490.0434	0.0005	0.9995
SYMI	17	-1228.1982	2491.7239	0.0002	0.9997
SYMIG	18	-1227.1328	2491.7526	0.0002	0.9999
GTR	19	-1226.8492	2493.3542	0.0001	1.0000
K2PG	13	-1236.7157	2500.2143	0.0000	1.0000
K2PI	13	-1237.2721	2501.3269	0.0000	1.0000
K2PIG	14	-1236.5303	2501.9658	0.0000	1.0000
HKY	15	-1236.6544	2504.3455	0.0000	1.0000
SYM	16	-1236.4200	2506.0175	0.0000	1.0000
F81G	15	-1241.8986	2514.8339	0.0000	1.0000
F81I	15	-1242.4554	2515.9476	0.0000	1.0000
F81IG	16	-1241.7206	2516.6186	0.0000	1.0000
K2P	12	-1246.5120	2517.6935	0.0000	1.0000
JC69G	12	-1251.1907	2527.0510	0.0000	1.0000
JC69I	12	-1251.7125	2528.0946	0.0000	1.0000
JC69IG	13	-1251.0520	2528.8867	0.0000	1.0000
F81	14	-1251.7370	2532.3792	0.0000	1.0000
JC69	11	-1260.5919	2543.7492	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1217.0552	2474.1105	0.5192	0.5192
GTRIG	21	-1216.7671	2475.5342	0.2548	0.7740
GTRI	20	-1217.9205	2475.8410	0.2186	0.9925
HKYG	16	-1226.1080	2484.2160	0.0033	0.9959
HKYI	16	-1226.7164	2485.4329	0.0018	0.9977
HKYIG	17	-1225.8468	2485.6937	0.0016	0.9993
SYMG	17	-1227.3579	2488.7159	0.0003	0.9996
SYMIG	18	-1227.1328	2490.2656	0.0002	0.9998
SYMI	17	-1228.1982	2490.3964	0.0002	0.9999
GTR	19	-1226.8492	2491.6985	0.0001	1.0000
K2PG	13	-1236.7157	2499.4315	0.0000	1.0000
K2PI	13	-1237.2721	2500.5441	0.0000	1.0000
K2PIG	14	-1236.5303	2501.0607	0.0000	1.0000
HKY	15	-1236.6544	2503.3087	0.0000	1.0000
SYM	16	-1236.4200	2504.8400	0.0000	1.0000
F81G	15	-1241.8986	2513.7972	0.0000	1.0000
F81I	15	-1242.4554	2514.9109	0.0000	1.0000
F81IG	16	-1241.7206	2515.4411	0.0000	1.0000
K2P	12	-1246.5120	2517.0240	0.0000	1.0000
JC69G	12	-1251.1907	2526.3814	0.0000	1.0000
JC69I	12	-1251.7125	2527.4251	0.0000	1.0000
JC69IG	13	-1251.0520	2528.1039	0.0000	1.0000
F81	14	-1251.7370	2531.4740	0.0000	1.0000
JC69	11	-1260.5919	2543.1839	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1226.1080	2550.9632	0.4759	0.4759
HKYI	16	-1226.7164	2552.1801	0.2590	0.7349
K2PG	13	-1236.7157	2553.6636	0.1234	0.8583
K2PI	13	-1237.2721	2554.7762	0.0707	0.9290
HKYIG	17	-1225.8468	2556.6126	0.0282	0.9572
GTRG	20	-1217.0552	2557.5445	0.0177	0.9749
GTRI	20	-1217.9205	2559.2750	0.0075	0.9824
K2PIG	14	-1236.5303	2559.4645	0.0068	0.9892
SYMG	17	-1227.3579	2559.6348	0.0062	0.9954
SYMI	17	-1228.1982	2561.3153	0.0027	0.9981
GTRIG	21	-1216.7671	2563.1399	0.0011	0.9992
SYMIG	18	-1227.1328	2565.3562	0.0004	0.9995
HKY	15	-1236.6544	2565.8842	0.0003	0.9998
K2P	12	-1246.5120	2567.0844	0.0002	1.0000
GTR	19	-1226.8492	2570.9608	0.0000	1.0000
SYM	16	-1236.4200	2571.5872	0.0000	1.0000
F81G	15	-1241.8986	2576.3727	0.0000	1.0000
JC69G	12	-1251.1907	2576.4418	0.0000	1.0000
JC69I	12	-1251.7125	2577.4855	0.0000	1.0000
F81I	15	-1242.4554	2577.4864	0.0000	1.0000
F81IG	16	-1241.7206	2582.1883	0.0000	1.0000
JC69IG	13	-1251.0520	2582.3360	0.0000	1.0000
JC69	11	-1260.5919	2589.0726	0.0000	1.0000
F81	14	-1251.7370	2589.8778	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/chr8_4057.phylip = 1 - 479;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4057.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/chr8_4057.phylip = 1 - 479;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4057.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 (HKYG)]

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

