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

Input data from file "591-tomodel-phylip/chr4_13410.phylip" (7 taxa, 784 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	-2235.9405	4512.9818	0.4260	0.4260
GTRG	20	-2236.1317	4513.3642	0.3518	0.7778
GTRIG	21	-2235.7974	4514.8073	0.1710	0.9488
HKYI	16	-2242.9797	4518.6687	0.0248	0.9736
HKYG	16	-2243.3520	4519.4132	0.0171	0.9907
HKYIG	17	-2242.9131	4520.6253	0.0093	1.0000
GTR	19	-2252.4723	4543.9394	0.0000	1.0000
HKY	15	-2259.8912	4550.4073	0.0000	1.0000
SYMI	17	-2280.5194	4595.8378	0.0000	1.0000
SYMG	17	-2281.1378	4597.0745	0.0000	1.0000
SYMIG	18	-2280.5082	4597.9105	0.0000	1.0000
K2PI	13	-2287.2393	4600.9514	0.0000	1.0000
K2PG	13	-2287.8553	4602.1833	0.0000	1.0000
K2PIG	14	-2287.2321	4603.0103	0.0000	1.0000
F81I	15	-2300.5321	4631.6893	0.0000	1.0000
SYM	16	-2299.6488	4632.0068	0.0000	1.0000
F81G	15	-2300.9309	4632.4868	0.0000	1.0000
F81IG	16	-2300.5087	4633.7266	0.0000	1.0000
K2P	12	-2305.4064	4635.2174	0.0000	1.0000
F81	14	-2316.6747	4661.8956	0.0000	1.0000
JC69I	12	-2340.0296	4704.4638	0.0000	1.0000
JC69G	12	-2340.5131	4705.4308	0.0000	1.0000
JC69IG	13	-2340.0220	4706.5166	0.0000	1.0000
JC69	11	-2356.7476	4735.8371	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRI	20	-2235.9405	4511.8809	0.4255	0.4255
GTRG	20	-2236.1317	4512.2633	0.3515	0.7770
GTRIG	21	-2235.7974	4513.5947	0.1806	0.9576
HKYI	16	-2242.9797	4517.9595	0.0204	0.9780
HKYG	16	-2243.3520	4518.7039	0.0140	0.9920
HKYIG	17	-2242.9131	4519.8263	0.0080	1.0000
GTR	19	-2252.4723	4542.9447	0.0000	1.0000
HKY	15	-2259.8912	4549.7823	0.0000	1.0000
SYMI	17	-2280.5194	4595.0389	0.0000	1.0000
SYMG	17	-2281.1378	4596.2755	0.0000	1.0000
SYMIG	18	-2280.5082	4597.0163	0.0000	1.0000
K2PI	13	-2287.2393	4600.4786	0.0000	1.0000
K2PG	13	-2287.8553	4601.7106	0.0000	1.0000
K2PIG	14	-2287.2321	4602.4642	0.0000	1.0000
F81I	15	-2300.5321	4631.0643	0.0000	1.0000
SYM	16	-2299.6488	4631.2975	0.0000	1.0000
F81G	15	-2300.9309	4631.8618	0.0000	1.0000
F81IG	16	-2300.5087	4633.0173	0.0000	1.0000
K2P	12	-2305.4064	4634.8127	0.0000	1.0000
F81	14	-2316.6747	4661.3494	0.0000	1.0000
JC69I	12	-2340.0296	4704.0591	0.0000	1.0000
JC69G	12	-2340.5131	4705.0261	0.0000	1.0000
JC69IG	13	-2340.0220	4706.0439	0.0000	1.0000
JC69	11	-2356.7476	4735.4951	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-2242.9797	4592.5900	0.5778	0.5778
HKYG	16	-2243.3520	4593.3345	0.3982	0.9759
HKYIG	17	-2242.9131	4599.1213	0.0221	0.9980
GTRI	20	-2235.9405	4605.1691	0.0011	0.9991
GTRG	20	-2236.1317	4605.5515	0.0009	1.0000
GTRIG	21	-2235.7974	4611.5473	0.0000	1.0000
HKY	15	-2259.8912	4619.7485	0.0000	1.0000
GTR	19	-2252.4723	4631.5685	0.0000	1.0000
K2PI	13	-2287.2393	4661.1160	0.0000	1.0000
K2PG	13	-2287.8553	4662.3479	0.0000	1.0000
K2PIG	14	-2287.2321	4667.7659	0.0000	1.0000
SYMI	17	-2280.5194	4674.3338	0.0000	1.0000
SYMG	17	-2281.1378	4675.5705	0.0000	1.0000
SYMIG	18	-2280.5082	4680.9757	0.0000	1.0000
K2P	12	-2305.4064	4690.7856	0.0000	1.0000
F81I	15	-2300.5321	4701.0304	0.0000	1.0000
F81G	15	-2300.9309	4701.8280	0.0000	1.0000
SYM	16	-2299.6488	4705.9281	0.0000	1.0000
F81IG	16	-2300.5087	4707.6479	0.0000	1.0000
F81	14	-2316.6747	4726.6511	0.0000	1.0000
JC69I	12	-2340.0296	4760.0320	0.0000	1.0000
JC69G	12	-2340.5131	4760.9990	0.0000	1.0000
JC69IG	13	-2340.0220	4766.6812	0.0000	1.0000
JC69	11	-2356.7476	4786.8036	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/chr4_13410.phylip = 1 - 784;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13410.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/chr4_13410.phylip = 1 - 784;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13410.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/chr4_13410.phylip = 1 - 784;
 Partition Dummy = 1:591-tomodel-phylip/chr4_13410.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

