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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	16	-1001.1599	2035.9151	0.5505	0.5505
HKYIG	17	-1000.8763	2037.5527	0.2428	0.7933
HKYI	16	-1002.6561	2038.9076	0.1233	0.9165
GTRG	20	-999.2397	2040.9721	0.0439	0.9605
GTRIG	21	-998.9627	2042.6753	0.0187	0.9792
GTRI	20	-1000.7561	2044.0049	0.0096	0.9888
F81IG	16	-1005.5210	2044.6374	0.0070	0.9959
K2PG	13	-1010.2256	2047.5094	0.0017	0.9975
F81G	15	-1008.7161	2048.8358	0.0009	0.9984
K2PIG	14	-1009.9238	2049.0721	0.0008	0.9992
K2PI	13	-1011.6563	2050.3707	0.0004	0.9996
F81I	15	-1009.9629	2051.3293	0.0002	0.9998
SYMG	17	-1009.2520	2054.3039	0.0001	0.9999
HKY	15	-1011.7071	2054.8176	0.0000	0.9999
JC69IG	13	-1014.0871	2055.2323	0.0000	1.0000
SYMIG	18	-1008.9442	2055.9061	0.0000	1.0000
SYMI	17	-1010.6197	2057.0393	0.0000	1.0000
JC69G	12	-1017.3829	2059.6701	0.0000	1.0000
JC69I	12	-1018.6043	2062.1129	0.0000	1.0000
K2P	12	-1021.1391	2067.1825	0.0000	1.0000
GTR	19	-1014.0310	2068.3106	0.0000	1.0000
F81	14	-1020.0690	2069.3625	0.0000	1.0000
JC69	11	-1029.0640	2080.8910	0.0000	1.0000
SYM	16	-1024.3729	2082.3411	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-1001.1599	2034.3198	0.5143	0.5143
HKYIG	17	-1000.8763	2035.7527	0.2512	0.7655
HKYI	16	-1002.6561	2037.3123	0.1152	0.8807
GTRG	20	-999.2397	2038.4795	0.0643	0.9450
GTRIG	21	-998.9627	2039.9253	0.0312	0.9761
GTRI	20	-1000.7561	2041.5123	0.0141	0.9903
F81IG	16	-1005.5210	2043.0421	0.0066	0.9968
K2PG	13	-1010.2256	2046.4513	0.0012	0.9980
F81G	15	-1008.7161	2047.4323	0.0007	0.9987
K2PIG	14	-1009.9238	2047.8476	0.0006	0.9993
K2PI	13	-1011.6563	2049.3126	0.0003	0.9996
F81I	15	-1009.9629	2049.9258	0.0002	0.9998
SYMG	17	-1009.2520	2052.5039	0.0001	0.9999
HKY	15	-1011.7071	2053.4141	0.0000	0.9999
SYMIG	18	-1008.9442	2053.8884	0.0000	1.0000
JC69IG	13	-1014.0871	2054.1742	0.0000	1.0000
SYMI	17	-1010.6197	2055.2393	0.0000	1.0000
JC69G	12	-1017.3829	2058.7657	0.0000	1.0000
JC69I	12	-1018.6043	2061.2086	0.0000	1.0000
GTR	19	-1014.0310	2066.0621	0.0000	1.0000
K2P	12	-1021.1391	2066.2781	0.0000	1.0000
F81	14	-1020.0690	2068.1380	0.0000	1.0000
JC69	11	-1029.0640	2080.1280	0.0000	1.0000
SYM	16	-1024.3729	2080.7458	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-1001.1599	2096.4083	0.4206	0.4206
K2PG	13	-1010.2256	2096.8982	0.3292	0.7497
HKYI	16	-1002.6561	2099.4008	0.0942	0.8439
K2PI	13	-1011.6563	2099.7595	0.0787	0.9227
HKYIG	17	-1000.8763	2101.7217	0.0295	0.9522
K2PIG	14	-1009.9238	2102.1751	0.0235	0.9757
JC69IG	13	-1014.0871	2104.6211	0.0069	0.9826
F81IG	16	-1005.5210	2105.1306	0.0054	0.9880
JC69G	12	-1017.3829	2105.3321	0.0049	0.9929
F81G	15	-1008.7161	2105.6403	0.0042	0.9970
JC69I	12	-1018.6043	2107.7750	0.0014	0.9984
F81I	15	-1009.9629	2108.1338	0.0012	0.9996
HKY	15	-1011.7071	2111.6221	0.0002	0.9998
K2P	12	-1021.1391	2112.8445	0.0001	1.0000
GTRG	20	-999.2397	2116.0902	0.0000	1.0000
SYMG	17	-1009.2520	2118.4730	0.0000	1.0000
GTRI	20	-1000.7561	2119.1229	0.0000	1.0000
SYMI	17	-1010.6197	2121.2084	0.0000	1.0000
GTRIG	21	-998.9627	2121.4165	0.0000	1.0000
F81	14	-1020.0690	2122.4655	0.0000	1.0000
JC69	11	-1029.0640	2122.8138	0.0000	1.0000
SYMIG	18	-1008.9442	2123.7380	0.0000	1.0000
GTR	19	-1014.0310	2139.7922	0.0000	1.0000
SYM	16	-1024.3729	2142.8344	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/chr4_17854.phylip = 1 - 358;
 Partition Dummy = 1:591-tomodel-phylip/chr4_17854.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/chr4_17854.phylip = 1 - 358;
 Partition Dummy = 1:591-tomodel-phylip/chr4_17854.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 (HKYG)]

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

