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

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

Minimum AIC  model: GTRI
Minimum AICc model: GTRI
Minimum BIC  model: HKYI


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	-1515.6637	3072.4044	0.2831	0.2831
GTRG	20	-1515.6754	3072.4278	0.2798	0.5630
HKYI	16	-1520.4860	3073.6658	0.1507	0.7137
HKYG	16	-1520.5987	3073.8913	0.1346	0.8483
GTRIG	21	-1515.6640	3074.5142	0.0986	0.9469
HKYIG	17	-1520.4858	3075.7531	0.0531	0.9999
GTR	19	-1525.4881	3089.9493	0.0000	1.0000
HKY	15	-1530.8886	3092.3886	0.0000	1.0000
F81I	15	-1534.9125	3100.4365	0.0000	1.0000
F81G	15	-1535.0758	3100.7632	0.0000	1.0000
F81IG	16	-1534.9126	3102.5191	0.0000	1.0000
F81	14	-1545.3104	3119.1551	0.0000	1.0000
K2PI	13	-1567.6350	3161.7326	0.0000	1.0000
K2PG	13	-1568.1264	3162.7152	0.0000	1.0000
K2PIG	14	-1567.6352	3163.8047	0.0000	1.0000
SYMI	17	-1566.0960	3166.9736	0.0000	1.0000
SYMG	17	-1566.4647	3167.7111	0.0000	1.0000
SYMIG	18	-1566.0958	3169.0662	0.0000	1.0000
K2P	12	-1579.7574	3183.9107	0.0000	1.0000
JC69I	12	-1580.2363	3184.8686	0.0000	1.0000
JC69G	12	-1580.6384	3185.6728	0.0000	1.0000
JC69IG	13	-1580.2363	3186.9352	0.0000	1.0000
SYM	16	-1577.7350	3188.1640	0.0000	1.0000
JC69	11	-1591.6740	3205.6827	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-1515.6637	3071.3275	0.2985	0.2985
GTRG	20	-1515.6754	3071.3508	0.2950	0.5935
HKYI	16	-1520.4860	3072.9719	0.1312	0.7247
HKYG	16	-1520.5987	3073.1974	0.1172	0.8419
GTRIG	21	-1515.6640	3073.3280	0.1098	0.9517
HKYIG	17	-1520.4858	3074.9715	0.0483	0.9999
GTR	19	-1525.4881	3088.9762	0.0000	1.0000
HKY	15	-1530.8886	3091.7772	0.0000	1.0000
F81I	15	-1534.9125	3099.8250	0.0000	1.0000
F81G	15	-1535.0758	3100.1517	0.0000	1.0000
F81IG	16	-1534.9126	3101.8253	0.0000	1.0000
F81	14	-1545.3104	3118.6208	0.0000	1.0000
K2PI	13	-1567.6350	3161.2701	0.0000	1.0000
K2PG	13	-1568.1264	3162.2527	0.0000	1.0000
K2PIG	14	-1567.6352	3163.2703	0.0000	1.0000
SYMI	17	-1566.0960	3166.1920	0.0000	1.0000
SYMG	17	-1566.4647	3166.9295	0.0000	1.0000
SYMIG	18	-1566.0958	3168.1915	0.0000	1.0000
K2P	12	-1579.7574	3183.5148	0.0000	1.0000
JC69I	12	-1580.2363	3184.4726	0.0000	1.0000
JC69G	12	-1580.6384	3185.2768	0.0000	1.0000
JC69IG	13	-1580.2363	3186.4727	0.0000	1.0000
SYM	16	-1577.7350	3187.4701	0.0000	1.0000
JC69	11	-1591.6740	3205.3481	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-1520.4860	3147.9457	0.5181	0.5181
HKYG	16	-1520.5987	3148.1712	0.4629	0.9810
HKYIG	17	-1520.4858	3154.6312	0.0183	0.9993
HKY	15	-1530.8886	3162.0651	0.0004	0.9998
GTRI	20	-1515.6637	3165.0447	0.0001	0.9999
GTRG	20	-1515.6754	3165.0681	0.0001	1.0000
F81I	15	-1534.9125	3170.1130	0.0000	1.0000
F81G	15	-1535.0758	3170.4396	0.0000	1.0000
GTRIG	21	-1515.6640	3171.7311	0.0000	1.0000
F81IG	16	-1534.9126	3176.7990	0.0000	1.0000
GTR	19	-1525.4881	3178.0075	0.0000	1.0000
F81	14	-1545.3104	3184.2228	0.0000	1.0000
K2PI	13	-1567.6350	3222.1863	0.0000	1.0000
K2PG	13	-1568.1264	3223.1689	0.0000	1.0000
K2PIG	14	-1567.6352	3228.8724	0.0000	1.0000
K2P	12	-1579.7574	3239.7451	0.0000	1.0000
JC69I	12	-1580.2363	3240.7030	0.0000	1.0000
JC69G	12	-1580.6384	3241.5072	0.0000	1.0000
SYMI	17	-1566.0960	3245.8516	0.0000	1.0000
SYMG	17	-1566.4647	3246.5891	0.0000	1.0000
JC69IG	13	-1580.2363	3247.3889	0.0000	1.0000
SYMIG	18	-1566.0958	3252.5370	0.0000	1.0000
JC69	11	-1591.6740	3256.8926	0.0000	1.0000
SYM	16	-1577.7350	3262.4439	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/chr7_5163.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5163.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/chr7_5163.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5163.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr7_5163.phylip = 1 - 801;
 Partition Dummy = 1:591-tomodel-phylip/chr7_5163.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

