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

Input data from file "683-tomodel-phylip/chr4_11159.phylip" (7 taxa, 729 characters)

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
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
K2PG	13	-2260.1934	4546.8958	0.3061	0.3061
HKYG	16	-2257.8677	4548.4993	0.1373	0.4433
K2PIG	14	-2260.0455	4548.6792	0.1255	0.5688
GTRG	20	-2253.8070	4548.8005	0.1181	0.6869
SYMG	17	-2257.0885	4549.0378	0.1049	0.7918
HKYIG	17	-2257.7140	4550.2888	0.0561	0.8479
GTRIG	21	-2253.5789	4550.4647	0.0514	0.8993
SYMIG	18	-2256.9139	4550.7911	0.0436	0.9429
K2PI	13	-2262.5397	4551.5886	0.0293	0.9722
HKYI	16	-2260.2016	4553.1673	0.0133	0.9855
SYMI	17	-2259.6429	4554.1465	0.0082	0.9937
GTRI	20	-2256.7351	4554.6566	0.0063	1.0000
K2P	12	-2290.5119	4605.4596	0.0000	1.0000
HKY	15	-2288.4323	4607.5378	0.0000	1.0000
SYM	16	-2288.0138	4608.7917	0.0000	1.0000
GTR	19	-2285.2987	4609.6694	0.0000	1.0000
JC69G	12	-2316.0243	4656.4844	0.0000	1.0000
JC69IG	13	-2315.9228	4658.3546	0.0000	1.0000
F81G	15	-2314.1178	4658.9088	0.0000	1.0000
F81IG	16	-2314.0100	4660.7839	0.0000	1.0000
JC69I	12	-2318.2917	4661.0191	0.0000	1.0000
F81I	15	-2316.4338	4663.5407	0.0000	1.0000
JC69	11	-2343.8919	4710.1520	0.0000	1.0000
F81	14	-2341.9436	4712.4754	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	13	-2260.1934	4546.3867	0.2674	0.2674
GTRG	20	-2253.8070	4547.6141	0.1448	0.4122
HKYG	16	-2257.8677	4547.7353	0.1363	0.5484
K2PIG	14	-2260.0455	4548.0910	0.1141	0.6625
SYMG	17	-2257.0885	4548.1770	0.1093	0.7717
GTRIG	21	-2253.5789	4549.1578	0.0669	0.8387
HKYIG	17	-2257.7140	4549.4280	0.0584	0.8971
SYMIG	18	-2256.9139	4549.8278	0.0479	0.9450
K2PI	13	-2262.5397	4551.0795	0.0256	0.9706
HKYI	16	-2260.2016	4552.4033	0.0132	0.9838
SYMI	17	-2259.6429	4553.2857	0.0085	0.9923
GTRI	20	-2256.7351	4553.4702	0.0077	1.0000
K2P	12	-2290.5119	4605.0238	0.0000	1.0000
HKY	15	-2288.4323	4606.8646	0.0000	1.0000
SYM	16	-2288.0138	4608.0276	0.0000	1.0000
GTR	19	-2285.2987	4608.5975	0.0000	1.0000
JC69G	12	-2316.0243	4656.0486	0.0000	1.0000
JC69IG	13	-2315.9228	4657.8455	0.0000	1.0000
F81G	15	-2314.1178	4658.2356	0.0000	1.0000
F81IG	16	-2314.0100	4660.0199	0.0000	1.0000
JC69I	12	-2318.2917	4660.5833	0.0000	1.0000
F81I	15	-2316.4338	4662.8675	0.0000	1.0000
JC69	11	-2343.8919	4709.7838	0.0000	1.0000
F81	14	-2341.9436	4711.8872	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-2260.1934	4606.0785	0.8777	0.8777
K2PI	13	-2262.5397	4610.7713	0.0840	0.9617
K2PIG	14	-2260.0455	4612.3744	0.0377	0.9994
HKYG	16	-2257.8677	4621.2021	0.0005	0.9999
HKYI	16	-2260.2016	4625.8701	0.0000	0.9999
SYMG	17	-2257.0885	4626.2355	0.0000	1.0000
HKYIG	17	-2257.7140	4627.4865	0.0000	1.0000
SYMI	17	-2259.6429	4631.3442	0.0000	1.0000
SYMIG	18	-2256.9139	4632.4779	0.0000	1.0000
GTRG	20	-2253.8070	4639.4475	0.0000	1.0000
GTRI	20	-2256.7351	4645.3036	0.0000	1.0000
GTRIG	21	-2253.5789	4645.5829	0.0000	1.0000
K2P	12	-2290.5119	4660.1239	0.0000	1.0000
HKY	15	-2288.4323	4675.7397	0.0000	1.0000
SYM	16	-2288.0138	4681.4944	0.0000	1.0000
GTR	19	-2285.2987	4695.8393	0.0000	1.0000
JC69G	12	-2316.0243	4711.1487	0.0000	1.0000
JC69I	12	-2318.2917	4715.6834	0.0000	1.0000
JC69IG	13	-2315.9228	4717.5373	0.0000	1.0000
F81G	15	-2314.1178	4727.1107	0.0000	1.0000
F81I	15	-2316.4338	4731.7426	0.0000	1.0000
F81IG	16	-2314.0100	4733.4867	0.0000	1.0000
JC69	11	-2343.8919	4760.2922	0.0000	1.0000
F81	14	-2341.9436	4776.1706	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_11159.phylip = 1 - 729;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11159.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_11159.phylip = 1 - 729;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11159.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;


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

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr4_11159.phylip = 1 - 729;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11159.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

