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

Input data from file "591-tomodel-phylip/chr6_4409.phylip" (7 taxa, 436 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	-1205.2543	2452.5327	0.7584	0.7584
GTRI	20	-1206.9504	2455.9248	0.1391	0.8975
GTRIG	21	-1206.1752	2456.5822	0.1001	0.9976
HKYG	16	-1215.6670	2464.6323	0.0018	0.9994
HKYI	16	-1217.3312	2467.9608	0.0003	0.9998
HKYIG	17	-1216.6246	2468.7134	0.0002	1.0000
GTR	19	-1220.6611	2481.1492	0.0000	1.0000
SYMG	17	-1224.4979	2484.4599	0.0000	1.0000
SYMI	17	-1225.6372	2486.7385	0.0000	1.0000
SYMIG	18	-1225.0814	2487.8031	0.0000	1.0000
K2PG	13	-1230.6091	2488.0807	0.0000	1.0000
K2PI	13	-1231.6286	2490.1198	0.0000	1.0000
K2PIG	14	-1231.0117	2491.0211	0.0000	1.0000
HKY	15	-1230.3296	2491.8020	0.0000	1.0000
F81G	15	-1237.1780	2505.4988	0.0000	1.0000
F81I	15	-1237.6882	2506.5193	0.0000	1.0000
F81IG	16	-1239.2612	2511.8206	0.0000	1.0000
SYM	16	-1240.2227	2513.7437	0.0000	1.0000
K2P	12	-1245.9565	2516.6506	0.0000	1.0000
JC69G	12	-1251.0606	2526.8588	0.0000	1.0000
JC69I	12	-1251.5420	2527.8216	0.0000	1.0000
F81	14	-1250.9863	2530.9702	0.0000	1.0000
JC69IG	13	-1252.3407	2531.5439	0.0000	1.0000
JC69	11	-1265.3345	2553.2916	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1205.2543	2450.5086	0.7507	0.7507
GTRI	20	-1206.9504	2453.9007	0.1377	0.8884
GTRIG	21	-1206.1752	2454.3503	0.1100	0.9984
HKYG	16	-1215.6670	2463.3340	0.0012	0.9996
HKYI	16	-1217.3312	2466.6624	0.0002	0.9998
HKYIG	17	-1216.6246	2467.2493	0.0002	1.0000
GTR	19	-1220.6611	2479.3223	0.0000	1.0000
SYMG	17	-1224.4979	2482.9958	0.0000	1.0000
SYMI	17	-1225.6372	2485.2743	0.0000	1.0000
SYMIG	18	-1225.0814	2486.1628	0.0000	1.0000
K2PG	13	-1230.6091	2487.2181	0.0000	1.0000
K2PI	13	-1231.6286	2489.2573	0.0000	1.0000
K2PIG	14	-1231.0117	2490.0235	0.0000	1.0000
HKY	15	-1230.3296	2490.6591	0.0000	1.0000
F81G	15	-1237.1780	2504.3559	0.0000	1.0000
F81I	15	-1237.6882	2505.3764	0.0000	1.0000
F81IG	16	-1239.2612	2510.5223	0.0000	1.0000
SYM	16	-1240.2227	2512.4453	0.0000	1.0000
K2P	12	-1245.9565	2515.9130	0.0000	1.0000
JC69G	12	-1251.0606	2526.1212	0.0000	1.0000
JC69I	12	-1251.5420	2527.0840	0.0000	1.0000
F81	14	-1250.9863	2529.9726	0.0000	1.0000
JC69IG	13	-1252.3407	2530.6814	0.0000	1.0000
JC69	11	-1265.3345	2552.6690	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1215.6670	2528.5763	0.7031	0.7031
HKYI	16	-1217.3312	2531.9047	0.1331	0.8362
GTRG	20	-1205.2543	2532.0614	0.1231	0.9593
GTRI	20	-1206.9504	2535.4536	0.0226	0.9818
HKYIG	17	-1216.6246	2536.5692	0.0129	0.9947
GTRIG	21	-1206.1752	2539.9808	0.0023	0.9971
K2PG	13	-1230.6091	2540.2275	0.0021	0.9992
K2PI	13	-1231.6286	2542.2666	0.0007	0.9999
K2PIG	14	-1231.0117	2547.1105	0.0001	1.0000
HKY	15	-1230.3296	2551.8238	0.0000	1.0000
SYMG	17	-1224.4979	2552.3157	0.0000	1.0000
SYMI	17	-1225.6372	2554.5943	0.0000	1.0000
GTR	19	-1220.6611	2556.7975	0.0000	1.0000
SYMIG	18	-1225.0814	2559.5604	0.0000	1.0000
K2P	12	-1245.9565	2564.8447	0.0000	1.0000
F81G	15	-1237.1780	2565.5205	0.0000	1.0000
F81I	15	-1237.6882	2566.5411	0.0000	1.0000
JC69G	12	-1251.0606	2575.0529	0.0000	1.0000
F81IG	16	-1239.2612	2575.7646	0.0000	1.0000
JC69I	12	-1251.5420	2576.0157	0.0000	1.0000
SYM	16	-1240.2227	2577.6876	0.0000	1.0000
JC69IG	13	-1252.3407	2583.6907	0.0000	1.0000
F81	14	-1250.9863	2587.0596	0.0000	1.0000
JC69	11	-1265.3345	2597.5231	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/chr6_4409.phylip = 1 - 436;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4409.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/chr6_4409.phylip = 1 - 436;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4409.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/chr6_4409.phylip = 1 - 436;
 Partition Dummy = 1:591-tomodel-phylip/chr6_4409.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

