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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
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
GTRI	20	-1246.5580	2534.6890	0.1938	0.1938
HKYI	16	-1250.8906	2534.7923	0.1840	0.3778
GTRG	20	-1246.6397	2534.8525	0.1786	0.5564
HKYG	16	-1250.9927	2534.9966	0.1662	0.7226
HKYIG	17	-1250.1691	2535.4779	0.1306	0.8532
GTRIG	21	-1245.9518	2535.6371	0.1206	0.9738
SYMI	17	-1253.1476	2541.4348	0.0066	0.9805
SYMG	17	-1253.2688	2541.6772	0.0059	0.9863
SYMIG	18	-1252.5989	2542.4738	0.0040	0.9903
K2PIG	14	-1256.8982	2542.5742	0.0038	0.9941
K2PI	13	-1258.1326	2542.9379	0.0031	0.9972
K2PG	13	-1258.4047	2543.4822	0.0024	0.9996
F81I	15	-1258.9857	2548.8619	0.0002	0.9997
F81G	15	-1259.0649	2549.0203	0.0001	0.9999
F81IG	16	-1258.4564	2549.9240	0.0001	1.0000
JC69I	12	-1265.2745	2555.1246	0.0000	1.0000
JC69G	12	-1265.4184	2555.4124	0.0000	1.0000
JC69IG	13	-1264.7353	2556.1435	0.0000	1.0000
GTR	19	-1264.1483	2567.7172	0.0000	1.0000
HKY	15	-1268.5310	2567.9526	0.0000	1.0000
SYM	16	-1270.8986	2574.8083	0.0000	1.0000
K2P	12	-1276.1004	2576.7765	0.0000	1.0000
F81	14	-1276.2540	2581.2857	0.0000	1.0000
JC69	11	-1282.9043	2588.2948	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1246.5580	2533.1159	0.2171	0.2171
GTRG	20	-1246.6397	2533.2795	0.2001	0.4172
HKYI	16	-1250.8906	2533.7811	0.1557	0.5729
GTRIG	21	-1245.9518	2533.9036	0.1464	0.7194
HKYG	16	-1250.9927	2533.9854	0.1406	0.8599
HKYIG	17	-1250.1691	2534.3383	0.1178	0.9778
SYMI	17	-1253.1476	2540.2951	0.0060	0.9838
SYMG	17	-1253.2688	2540.5376	0.0053	0.9891
SYMIG	18	-1252.5989	2541.1977	0.0038	0.9929
K2PIG	14	-1256.8982	2541.7964	0.0028	0.9957
K2PI	13	-1258.1326	2542.2651	0.0022	0.9980
K2PG	13	-1258.4047	2542.8094	0.0017	0.9997
F81I	15	-1258.9857	2547.9713	0.0001	0.9998
F81G	15	-1259.0649	2548.1298	0.0001	0.9999
F81IG	16	-1258.4564	2548.9128	0.0001	1.0000
JC69I	12	-1265.2745	2554.5490	0.0000	1.0000
JC69G	12	-1265.4184	2554.8368	0.0000	1.0000
JC69IG	13	-1264.7353	2555.4707	0.0000	1.0000
GTR	19	-1264.1483	2566.2966	0.0000	1.0000
HKY	15	-1268.5310	2567.0621	0.0000	1.0000
SYM	16	-1270.8986	2573.7972	0.0000	1.0000
K2P	12	-1276.1004	2576.2009	0.0000	1.0000
F81	14	-1276.2540	2580.5079	0.0000	1.0000
JC69	11	-1282.9043	2587.8086	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PI	13	-1258.1326	2598.4117	0.4632	0.4632
K2PG	13	-1258.4047	2598.9559	0.3528	0.8160
K2PIG	14	-1256.8982	2602.2620	0.0676	0.8836
HKYI	16	-1250.8906	2602.8846	0.0495	0.9331
HKYG	16	-1250.9927	2603.0889	0.0447	0.9778
JC69I	12	-1265.2745	2606.3766	0.0086	0.9864
JC69G	12	-1265.4184	2606.6644	0.0075	0.9939
HKYIG	17	-1250.1691	2607.7607	0.0043	0.9982
JC69IG	13	-1264.7353	2611.6173	0.0006	0.9988
F81I	15	-1258.9857	2612.7559	0.0004	0.9992
F81G	15	-1259.0649	2612.9143	0.0003	0.9995
SYMI	17	-1253.1476	2613.7176	0.0002	0.9997
SYMG	17	-1253.2688	2613.9600	0.0002	0.9999
F81IG	16	-1258.4564	2618.0163	0.0000	1.0000
SYMIG	18	-1252.5989	2618.9391	0.0000	1.0000
GTRI	20	-1246.5580	2619.4953	0.0000	1.0000
GTRG	20	-1246.6397	2619.6588	0.0000	1.0000
GTRIG	21	-1245.9518	2624.6019	0.0000	1.0000
K2P	12	-1276.1004	2628.0285	0.0000	1.0000
HKY	15	-1268.5310	2631.8466	0.0000	1.0000
JC69	11	-1282.9043	2635.3173	0.0000	1.0000
F81	14	-1276.2540	2640.9735	0.0000	1.0000
SYM	16	-1270.8986	2642.9006	0.0000	1.0000
GTR	19	-1264.1483	2648.3570	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_27825.phylip = 1 - 555;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27825.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_27825.phylip = 1 - 555;
 Partition Dummy = 1:591-tomodel-phylip/chr2_27825.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=propinv;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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

