
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:15 04/04/2011

Input data from file "683-tomodel-phylip/chr2_7976.phylip" (7 taxa, 797 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	-2189.4650	4411.6275	0.5216	0.5216
HKYI	16	-2190.1710	4413.0395	0.2575	0.7791
HKYIG	17	-2189.4504	4413.6864	0.1863	0.9654
GTRG	20	-2188.5826	4418.2477	0.0190	0.9845
GTRI	20	-2189.3580	4419.7985	0.0088	0.9932
GTRIG	21	-2188.5648	4420.3219	0.0068	1.0000
HKY	15	-2202.9071	4436.4288	0.0000	1.0000
GTR	19	-2201.9429	4442.8638	0.0000	1.0000
K2PG	13	-2248.1601	4522.7851	0.0000	1.0000
SYMG	17	-2244.0123	4522.8101	0.0000	1.0000
K2PI	13	-2248.5226	4523.5100	0.0000	1.0000
SYMI	17	-2244.4479	4523.6815	0.0000	1.0000
K2PIG	14	-2248.1700	4524.8771	0.0000	1.0000
SYMIG	18	-2244.0172	4524.9135	0.0000	1.0000
SYM	16	-2254.8348	4542.3671	0.0000	1.0000
K2P	12	-2259.2492	4542.8964	0.0000	1.0000
F81I	15	-2258.6658	4547.9462	0.0000	1.0000
F81G	15	-2259.1918	4548.9983	0.0000	1.0000
F81IG	16	-2259.1996	4551.0967	0.0000	1.0000
F81	14	-2267.6971	4563.9314	0.0000	1.0000
JC69I	12	-2306.3983	4637.1946	0.0000	1.0000
JC69G	12	-2306.8534	4638.1048	0.0000	1.0000
JC69IG	13	-2306.8657	4640.1964	0.0000	1.0000
JC69	11	-2314.9066	4652.1494	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	16	-2189.4650	4410.9301	0.5133	0.5133
HKYI	16	-2190.1710	4412.3421	0.2534	0.7667
HKYIG	17	-2189.4504	4412.9008	0.1916	0.9583
GTRG	20	-2188.5826	4417.1652	0.0227	0.9810
GTRI	20	-2189.3580	4418.7160	0.0105	0.9915
GTRIG	21	-2188.5648	4419.1296	0.0085	1.0000
HKY	15	-2202.9071	4435.8142	0.0000	1.0000
GTR	19	-2201.9429	4441.8857	0.0000	1.0000
SYMG	17	-2244.0123	4522.0245	0.0000	1.0000
K2PG	13	-2248.1601	4522.3202	0.0000	1.0000
SYMI	17	-2244.4479	4522.8959	0.0000	1.0000
K2PI	13	-2248.5226	4523.0451	0.0000	1.0000
SYMIG	18	-2244.0172	4524.0343	0.0000	1.0000
K2PIG	14	-2248.1700	4524.3400	0.0000	1.0000
SYM	16	-2254.8348	4541.6696	0.0000	1.0000
K2P	12	-2259.2492	4542.4985	0.0000	1.0000
F81I	15	-2258.6658	4547.3316	0.0000	1.0000
F81G	15	-2259.1918	4548.3837	0.0000	1.0000
F81IG	16	-2259.1996	4550.3993	0.0000	1.0000
F81	14	-2267.6971	4563.3943	0.0000	1.0000
JC69I	12	-2306.3983	4636.7966	0.0000	1.0000
JC69G	12	-2306.8534	4637.7069	0.0000	1.0000
JC69IG	13	-2306.8657	4639.7315	0.0000	1.0000
JC69	11	-2314.9066	4651.8131	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2189.4650	4485.8238	0.6538	0.6538
HKYI	16	-2190.1710	4487.2358	0.3227	0.9765
HKYIG	17	-2189.4504	4492.4753	0.0235	1.0000
HKY	15	-2202.9071	4506.0270	0.0000	1.0000
GTRG	20	-2188.5826	4510.7823	0.0000	1.0000
GTRI	20	-2189.3580	4512.3331	0.0000	1.0000
GTRIG	21	-2188.5648	4517.4276	0.0000	1.0000
GTR	19	-2201.9429	4530.8220	0.0000	1.0000
K2PG	13	-2248.1601	4583.1713	0.0000	1.0000
K2PI	13	-2248.5226	4583.8963	0.0000	1.0000
K2PIG	14	-2248.1700	4589.8720	0.0000	1.0000
K2P	12	-2259.2492	4598.6687	0.0000	1.0000
SYMG	17	-2244.0123	4601.5990	0.0000	1.0000
SYMI	17	-2244.4479	4602.4704	0.0000	1.0000
SYMIG	18	-2244.0172	4608.2897	0.0000	1.0000
SYM	16	-2254.8348	4616.5633	0.0000	1.0000
F81I	15	-2258.6658	4617.5444	0.0000	1.0000
F81G	15	-2259.1918	4618.5965	0.0000	1.0000
F81IG	16	-2259.1996	4625.2929	0.0000	1.0000
F81	14	-2267.6971	4628.9263	0.0000	1.0000
JC69I	12	-2306.3983	4692.9669	0.0000	1.0000
JC69G	12	-2306.8534	4693.8771	0.0000	1.0000
JC69IG	13	-2306.8657	4700.5826	0.0000	1.0000
JC69	11	-2314.9066	4703.3025	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 683-tomodel-phylip/chr2_7976.phylip = 1 - 797;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7976.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 683-tomodel-phylip/chr2_7976.phylip = 1 - 797;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7976.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 683-tomodel-phylip/chr2_7976.phylip = 1 - 797;
 Partition Dummy = 1:683-tomodel-phylip/chr2_7976.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

