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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-1450.6229	2943.0180	0.4394	0.4394
GTRI	20	-1451.1399	2944.0520	0.2620	0.7014
GTRIG	21	-1450.5867	2945.1268	0.1531	0.8545
SYMG	17	-1455.6560	2946.5950	0.0735	0.9280
SYMI	17	-1456.1894	2947.6619	0.0431	0.9711
SYMIG	18	-1455.5869	2948.6109	0.0268	0.9979
K2PG	13	-1464.4558	2955.6683	0.0008	0.9987
HKYG	16	-1461.9904	2957.1188	0.0004	0.9990
K2PI	13	-1465.3273	2957.4114	0.0003	0.9994
K2PIG	14	-1464.3161	2957.5072	0.0003	0.9997
HKYI	16	-1462.8256	2958.7893	0.0002	0.9998
HKYIG	17	-1461.8500	2958.9830	0.0002	1.0000
GTR	19	-1469.6393	2978.8786	0.0000	1.0000
SYM	16	-1474.0901	2981.3182	0.0000	1.0000
K2P	12	-1483.0023	2990.6520	0.0000	1.0000
HKY	15	-1480.6359	2992.2739	0.0000	1.0000
JC69G	12	-1499.9985	3024.6444	0.0000	1.0000
F81G	15	-1497.6225	3026.2471	0.0000	1.0000
JC69IG	13	-1499.8637	3026.4841	0.0000	1.0000
JC69I	12	-1500.9811	3026.6095	0.0000	1.0000
F81IG	16	-1497.4926	3028.1233	0.0000	1.0000
F81I	15	-1498.5776	3028.1573	0.0000	1.0000
JC69	11	-1517.1824	3056.9113	0.0000	1.0000
F81	14	-1514.6771	3058.2292	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1450.6229	2941.2458	0.4464	0.4464
GTRI	20	-1451.1399	2942.2799	0.2662	0.7126
GTRIG	21	-1450.5867	2943.1733	0.1703	0.8828
SYMG	17	-1455.6560	2945.3120	0.0584	0.9413
SYMI	17	-1456.1894	2946.3789	0.0343	0.9756
SYMIG	18	-1455.5869	2947.1739	0.0230	0.9986
K2PG	13	-1464.4558	2954.9116	0.0005	0.9991
HKYG	16	-1461.9904	2955.9807	0.0003	0.9994
K2PIG	14	-1464.3161	2956.6322	0.0002	0.9996
K2PI	13	-1465.3273	2956.6547	0.0002	0.9998
HKYI	16	-1462.8256	2957.6512	0.0001	0.9999
HKYIG	17	-1461.8500	2957.7000	0.0001	1.0000
GTR	19	-1469.6393	2977.2786	0.0000	1.0000
SYM	16	-1474.0901	2980.1802	0.0000	1.0000
K2P	12	-1483.0023	2990.0047	0.0000	1.0000
HKY	15	-1480.6359	2991.2718	0.0000	1.0000
JC69G	12	-1499.9985	3023.9971	0.0000	1.0000
F81G	15	-1497.6225	3025.2450	0.0000	1.0000
JC69IG	13	-1499.8637	3025.7273	0.0000	1.0000
JC69I	12	-1500.9811	3025.9622	0.0000	1.0000
F81IG	16	-1497.4926	3026.9852	0.0000	1.0000
F81I	15	-1498.5776	3027.1552	0.0000	1.0000
JC69	11	-1517.1824	3056.3647	0.0000	1.0000
F81	14	-1514.6771	3057.3542	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-1464.4558	3009.5708	0.6594	0.6594
K2PI	13	-1465.3273	3011.3139	0.2758	0.9353
K2PIG	14	-1464.3161	3015.4960	0.0341	0.9694
SYMG	17	-1455.6560	3016.7894	0.0179	0.9872
SYMI	17	-1456.1894	3017.8563	0.0105	0.9977
SYMIG	18	-1455.5869	3022.8559	0.0009	0.9985
HKYG	16	-1461.9904	3023.2537	0.0007	0.9992
HKYI	16	-1462.8256	3024.9241	0.0003	0.9996
GTRG	20	-1450.6229	3025.3370	0.0002	0.9998
GTRI	20	-1451.1399	3026.3710	0.0001	1.0000
HKYIG	17	-1461.8500	3029.1775	0.0000	1.0000
GTRIG	21	-1450.5867	3031.4691	0.0000	1.0000
K2P	12	-1483.0023	3040.4594	0.0000	1.0000
SYM	16	-1474.0901	3047.4531	0.0000	1.0000
HKY	15	-1480.6359	3054.3402	0.0000	1.0000
GTR	19	-1469.6393	3057.1652	0.0000	1.0000
JC69G	12	-1499.9985	3074.4518	0.0000	1.0000
JC69I	12	-1500.9811	3076.4169	0.0000	1.0000
JC69IG	13	-1499.8637	3080.3866	0.0000	1.0000
F81G	15	-1497.6225	3088.3134	0.0000	1.0000
F81I	15	-1498.5776	3090.2236	0.0000	1.0000
F81IG	16	-1497.4926	3094.2581	0.0000	1.0000
JC69	11	-1517.1824	3102.6149	0.0000	1.0000
F81	14	-1514.6771	3116.2180	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/chr21_1554.phylip = 1 - 495;
 Partition Dummy = 1:591-tomodel-phylip/chr21_1554.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/chr21_1554.phylip = 1 - 495;
 Partition Dummy = 1:591-tomodel-phylip/chr21_1554.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 (K2PG)]

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


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

