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

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

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


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	-2137.1089	4315.3911	0.4841	0.4841
GTRG	20	-2137.4545	4316.0823	0.3427	0.8268
GTRIG	21	-2137.0774	4317.4470	0.1732	1.0000
HKYI	16	-2154.1061	4340.9678	0.0000	1.0000
HKYG	16	-2154.5902	4341.9359	0.0000	1.0000
HKYIG	17	-2154.8369	4344.5250	0.0000	1.0000
GTR	19	-2171.3583	4381.7766	0.0000	1.0000
SYMI	17	-2184.1367	4403.1245	0.0000	1.0000
SYMG	17	-2184.4230	4403.6972	0.0000	1.0000
HKY	15	-2187.1421	4404.9499	0.0000	1.0000
SYMIG	18	-2184.4230	4405.7986	0.0000	1.0000
K2PI	13	-2191.4319	4409.3673	0.0000	1.0000
K2PG	13	-2191.9133	4410.3301	0.0000	1.0000
K2PIG	14	-2191.9134	4412.4085	0.0000	1.0000
F81I	15	-2191.4563	4413.5784	0.0000	1.0000
F81G	15	-2192.0201	4414.7060	0.0000	1.0000
F81IG	16	-2192.0202	4416.7960	0.0000	1.0000
SYM	16	-2215.5302	4463.8159	0.0000	1.0000
JC69I	12	-2221.8636	4468.1582	0.0000	1.0000
JC69G	12	-2222.3014	4469.0338	0.0000	1.0000
K2P	12	-2223.2931	4471.0172	0.0000	1.0000
JC69IG	13	-2222.3014	4471.1063	0.0000	1.0000
F81	14	-2222.8680	4474.3177	0.0000	1.0000
JC69	11	-2252.9328	4528.2298	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2137.1089	4314.2179	0.4790	0.4790
GTRG	20	-2137.4545	4314.9091	0.3391	0.8181
GTRIG	21	-2137.0774	4316.1547	0.1819	1.0000
HKYI	16	-2154.1061	4340.2122	0.0000	1.0000
HKYG	16	-2154.5902	4341.1804	0.0000	1.0000
HKYIG	17	-2154.8369	4343.6738	0.0000	1.0000
GTR	19	-2171.3583	4380.7167	0.0000	1.0000
SYMI	17	-2184.1367	4402.2733	0.0000	1.0000
SYMG	17	-2184.4230	4402.8460	0.0000	1.0000
HKY	15	-2187.1421	4404.2842	0.0000	1.0000
SYMIG	18	-2184.4230	4404.8460	0.0000	1.0000
K2PI	13	-2191.4319	4408.8638	0.0000	1.0000
K2PG	13	-2191.9133	4409.8267	0.0000	1.0000
K2PIG	14	-2191.9134	4411.8268	0.0000	1.0000
F81I	15	-2191.4563	4412.9127	0.0000	1.0000
F81G	15	-2192.0201	4414.0402	0.0000	1.0000
F81IG	16	-2192.0202	4416.0404	0.0000	1.0000
SYM	16	-2215.5302	4463.0604	0.0000	1.0000
JC69I	12	-2221.8636	4467.7273	0.0000	1.0000
JC69G	12	-2222.3014	4468.6028	0.0000	1.0000
K2P	12	-2223.2931	4470.5863	0.0000	1.0000
JC69IG	13	-2222.3014	4470.6029	0.0000	1.0000
F81	14	-2222.8680	4473.7360	0.0000	1.0000
JC69	11	-2252.9328	4527.8657	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRI	20	-2137.1089	4406.2696	0.5610	0.5610
GTRG	20	-2137.4545	4406.9608	0.3970	0.9580
GTRIG	21	-2137.0774	4412.8091	0.0213	0.9793
HKYI	16	-2154.1061	4413.8536	0.0127	0.9920
HKYG	16	-2154.5902	4414.8218	0.0078	0.9998
HKYIG	17	-2154.8369	4421.9178	0.0002	1.0000
GTR	19	-2171.3583	4468.1658	0.0000	1.0000
K2PI	13	-2191.4319	4468.6975	0.0000	1.0000
K2PG	13	-2191.9133	4469.6603	0.0000	1.0000
HKY	15	-2187.1421	4473.3230	0.0000	1.0000
K2PIG	14	-2191.9134	4476.2630	0.0000	1.0000
SYMI	17	-2184.1367	4480.5173	0.0000	1.0000
SYMG	17	-2184.4230	4481.0900	0.0000	1.0000
F81I	15	-2191.4563	4481.9515	0.0000	1.0000
F81G	15	-2192.0201	4483.0791	0.0000	1.0000
SYMIG	18	-2184.4230	4487.6926	0.0000	1.0000
F81IG	16	-2192.0202	4489.6818	0.0000	1.0000
JC69I	12	-2221.8636	4522.9584	0.0000	1.0000
JC69G	12	-2222.3014	4523.8339	0.0000	1.0000
K2P	12	-2223.2931	4525.8173	0.0000	1.0000
JC69IG	13	-2222.3014	4530.4365	0.0000	1.0000
SYM	16	-2215.5302	4536.7018	0.0000	1.0000
F81	14	-2222.8680	4538.1723	0.0000	1.0000
JC69	11	-2252.9328	4578.4941	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/chr2_7963.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7963.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/chr2_7963.phylip = 1 - 737;
 Partition Dummy = 1:591-tomodel-phylip/chr2_7963.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 (GTRI)]

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


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

