
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:47 02/17/2011

Input data from file "183-phylip/chr2_6767.phylip" (29 taxa, 452 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: GTRG


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
GTRG	64	-2154.3747	4458.2482	0.6374	0.6374
GTRIG	65	-2153.5754	4459.3787	0.3622	0.9995
GTRI	64	-2161.5530	4472.6047	0.0005	1.0000
HKYIG	61	-2171.2207	4483.8364	0.0000	1.0000
HKYG	60	-2173.1895	4485.1001	0.0000	1.0000
HKYI	60	-2178.4047	4495.5305	0.0000	1.0000
SYMIG	62	-2185.0504	4514.1830	0.0000	1.0000
SYMG	61	-2186.9153	4515.2255	0.0000	1.0000
K2PIG	58	-2191.4920	4516.3988	0.0000	1.0000
K2PG	57	-2193.2873	4517.3563	0.0000	1.0000
SYMI	61	-2194.2826	4529.9600	0.0000	1.0000
K2PI	57	-2200.8507	4532.4832	0.0000	1.0000
GTR	63	-2203.4263	4553.6361	0.0000	1.0000
HKY	59	-2222.0130	4580.0873	0.0000	1.0000
F81G	59	-2225.6541	4587.3695	0.0000	1.0000
F81IG	60	-2224.9465	4588.6143	0.0000	1.0000
F81I	59	-2233.8815	4603.8243	0.0000	1.0000
JC69G	56	-2239.3160	4606.7941	0.0000	1.0000
JC69IG	57	-2238.5833	4607.9483	0.0000	1.0000
SYM	60	-2237.1478	4613.0168	0.0000	1.0000
K2P	56	-2245.3331	4618.8282	0.0000	1.0000
JC69I	56	-2246.9449	4622.0519	0.0000	1.0000
F81	58	-2274.9492	4683.3131	0.0000	1.0000
JC69	55	-2289.9885	4705.5326	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2154.3747	4436.7495	0.5498	0.5498
GTRIG	65	-2153.5754	4437.1507	0.4498	0.9996
GTRI	64	-2161.5530	4451.1060	0.0004	1.0000
HKYIG	61	-2171.2207	4464.4415	0.0000	1.0000
HKYG	60	-2173.1895	4466.3789	0.0000	1.0000
HKYI	60	-2178.4047	4476.8093	0.0000	1.0000
SYMIG	62	-2185.0504	4494.1007	0.0000	1.0000
SYMG	61	-2186.9153	4495.8306	0.0000	1.0000
K2PIG	58	-2191.4920	4498.9840	0.0000	1.0000
K2PG	57	-2193.2873	4500.5746	0.0000	1.0000
SYMI	61	-2194.2826	4510.5651	0.0000	1.0000
K2PI	57	-2200.8507	4515.7015	0.0000	1.0000
GTR	63	-2203.4263	4532.8526	0.0000	1.0000
HKY	59	-2222.0130	4562.0261	0.0000	1.0000
F81G	59	-2225.6541	4569.3083	0.0000	1.0000
F81IG	60	-2224.9465	4569.8931	0.0000	1.0000
F81I	59	-2233.8815	4585.7631	0.0000	1.0000
JC69G	56	-2239.3160	4590.6321	0.0000	1.0000
JC69IG	57	-2238.5833	4591.1666	0.0000	1.0000
SYM	60	-2237.1478	4594.2956	0.0000	1.0000
K2P	56	-2245.3331	4602.6662	0.0000	1.0000
JC69I	56	-2246.9449	4605.8898	0.0000	1.0000
F81	58	-2274.9492	4665.8984	0.0000	1.0000
JC69	55	-2289.9885	4689.9770	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2154.3747	4700.0251	0.9032	0.9032
GTRIG	65	-2153.5754	4704.5400	0.0945	0.9976
HKYG	60	-2173.1895	4713.1999	0.0012	0.9989
GTRI	64	-2161.5530	4714.3817	0.0007	0.9996
HKYIG	61	-2171.2207	4715.3761	0.0004	1.0000
HKYI	60	-2178.4047	4723.6303	0.0000	1.0000
K2PG	57	-2193.2873	4735.0544	0.0000	1.0000
K2PIG	58	-2191.4920	4737.5776	0.0000	1.0000
SYMG	61	-2186.9153	4746.7652	0.0000	1.0000
SYMIG	62	-2185.0504	4749.1490	0.0000	1.0000
K2PI	57	-2200.8507	4750.1813	0.0000	1.0000
SYMI	61	-2194.2826	4761.4998	0.0000	1.0000
GTR	63	-2203.4263	4792.0146	0.0000	1.0000
HKY	59	-2222.0130	4804.7333	0.0000	1.0000
F81G	59	-2225.6541	4812.0155	0.0000	1.0000
F81IG	60	-2224.9465	4816.7140	0.0000	1.0000
JC69G	56	-2239.3160	4820.9983	0.0000	1.0000
JC69IG	57	-2238.5833	4825.6464	0.0000	1.0000
F81I	59	-2233.8815	4828.4703	0.0000	1.0000
K2P	56	-2245.3331	4833.0324	0.0000	1.0000
JC69I	56	-2246.9449	4836.2560	0.0000	1.0000
SYM	60	-2237.1478	4841.1165	0.0000	1.0000
F81	58	-2274.9492	4904.4919	0.0000	1.0000
JC69	55	-2289.9885	4916.2295	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_6767.phylip = 1 - 452;
 Partition Dummy = 1:183-phylip/chr2_6767.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_6767.phylip = 1 - 452;
 Partition Dummy = 1:183-phylip/chr2_6767.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_6767.phylip = 1 - 452;
 Partition Dummy = 1:183-phylip/chr2_6767.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

