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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-522.6068	1080.1384	0.3469	0.3469
HKYI	16	-522.6819	1080.2885	0.3219	0.6688
HKYIG	17	-522.4039	1082.1158	0.1291	0.7979
K2PG	13	-527.6803	1083.2864	0.0719	0.8698
K2PI	13	-527.7604	1083.4468	0.0664	0.9361
K2PIG	14	-527.4247	1085.0835	0.0293	0.9654
GTRG	20	-521.1587	1086.9327	0.0116	0.9770
GTRI	20	-521.2208	1087.0570	0.0109	0.9879
GTRIG	21	-521.1268	1089.3585	0.0035	0.9914
SYMG	17	-526.6747	1090.6575	0.0018	0.9932
F81G	15	-529.0960	1090.7589	0.0017	0.9949
SYMI	17	-526.7907	1090.8894	0.0016	0.9965
F81I	15	-529.2670	1091.1008	0.0014	0.9979
F81IG	16	-528.8114	1092.5475	0.0007	0.9986
SYMIG	18	-526.4174	1092.5522	0.0007	0.9993
JC69G	12	-534.3479	1094.3379	0.0003	0.9996
JC69I	12	-534.5093	1094.6607	0.0002	0.9999
JC69IG	13	-534.0068	1095.9394	0.0001	1.0000
HKY	15	-534.1915	1100.9498	0.0000	1.0000
K2P	12	-539.4665	1104.5751	0.0000	1.0000
GTR	19	-531.9579	1106.0689	0.0000	1.0000
F81	14	-539.3106	1108.8553	0.0000	1.0000
SYM	16	-538.1798	1111.2842	0.0000	1.0000
JC69	11	-544.9713	1113.3248	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-522.6068	1077.2136	0.3465	0.3465
HKYI	16	-522.6819	1077.3638	0.3214	0.6679
HKYIG	17	-522.4039	1078.8077	0.1561	0.8240
K2PG	13	-527.6803	1081.3605	0.0436	0.8676
K2PI	13	-527.7604	1081.5208	0.0402	0.9078
GTRG	20	-521.1587	1082.3174	0.0270	0.9348
GTRI	20	-521.2208	1082.4416	0.0254	0.9602
K2PIG	14	-527.4247	1082.8495	0.0207	0.9809
GTRIG	21	-521.1268	1084.2535	0.0103	0.9911
SYMG	17	-526.6747	1087.3494	0.0022	0.9933
SYMI	17	-526.7907	1087.5813	0.0019	0.9953
F81G	15	-529.0960	1088.1920	0.0014	0.9967
F81I	15	-529.2670	1088.5339	0.0012	0.9979
SYMIG	18	-526.4174	1088.8348	0.0010	0.9989
F81IG	16	-528.8114	1089.6228	0.0007	0.9996
JC69G	12	-534.3479	1092.6958	0.0002	0.9998
JC69I	12	-534.5093	1093.0186	0.0001	0.9999
JC69IG	13	-534.0068	1094.0135	0.0001	1.0000
HKY	15	-534.1915	1098.3830	0.0000	1.0000
GTR	19	-531.9579	1101.9159	0.0000	1.0000
K2P	12	-539.4665	1102.9330	0.0000	1.0000
F81	14	-539.3106	1106.6212	0.0000	1.0000
SYM	16	-538.1798	1108.3595	0.0000	1.0000
JC69	11	-544.9713	1111.9426	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-527.6803	1124.4322	0.4624	0.4624
K2PI	13	-527.7604	1124.5925	0.4268	0.8892
K2PIG	14	-527.4247	1129.2344	0.0419	0.9311
HKYG	16	-522.6068	1130.2249	0.0255	0.9566
HKYI	16	-522.6819	1130.3751	0.0237	0.9803
JC69G	12	-534.3479	1132.4543	0.0084	0.9887
JC69I	12	-534.5093	1132.7771	0.0071	0.9958
HKYIG	17	-522.4039	1135.1322	0.0022	0.9980
JC69IG	13	-534.0068	1137.0852	0.0008	0.9988
F81G	15	-529.0960	1137.8901	0.0006	0.9994
F81I	15	-529.2670	1138.2320	0.0005	0.9998
F81IG	16	-528.8114	1142.6341	0.0001	0.9999
K2P	12	-539.4665	1142.6914	0.0001	0.9999
SYMG	17	-526.6747	1143.6739	0.0000	1.0000
SYMI	17	-526.7907	1143.9058	0.0000	1.0000
HKY	15	-534.1915	1148.0810	0.0000	1.0000
JC69	11	-544.9713	1148.3879	0.0000	1.0000
SYMIG	18	-526.4174	1148.4725	0.0000	1.0000
GTRG	20	-521.1587	1148.5815	0.0000	1.0000
GTRI	20	-521.2208	1148.7057	0.0000	1.0000
F81	14	-539.3106	1153.0061	0.0000	1.0000
GTRIG	21	-521.1268	1153.8308	0.0000	1.0000
SYM	16	-538.1798	1161.3708	0.0000	1.0000
GTR	19	-531.9579	1164.8668	0.0000	1.0000



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

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

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


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

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


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

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

