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

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

Minimum AIC  model: SYMIG
Minimum AICc model: SYMIG
Minimum BIC  model: SYMG


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
SYMIG	62	-2344.6647	4843.2605	0.5466	0.5466
SYMG	61	-2346.4767	4843.8235	0.4125	0.9592
GTRIG	65	-2343.1544	4849.5646	0.0234	0.9825
GTRG	64	-2345.0126	4850.1488	0.0175	1.0000
HKYIG	61	-2358.6991	4868.2685	0.0000	1.0000
HKYG	60	-2360.9451	4869.7228	0.0000	1.0000
K2PIG	58	-2366.3522	4874.5308	0.0000	1.0000
K2PG	57	-2369.3197	4877.4966	0.0000	1.0000
GTRI	64	-2371.7389	4903.6013	0.0000	1.0000
SYMI	61	-2376.5741	4904.0184	0.0000	1.0000
HKYI	60	-2398.0776	4943.9879	0.0000	1.0000
K2PI	57	-2407.2179	4953.2930	0.0000	1.0000
F81IG	60	-2434.5774	5016.9875	0.0000	1.0000
F81G	59	-2436.9509	5018.7200	0.0000	1.0000
JC69IG	57	-2445.6636	5030.1843	0.0000	1.0000
JC69G	56	-2447.9643	5031.8388	0.0000	1.0000
F81I	59	-2471.6172	5088.0526	0.0000	1.0000
JC69I	56	-2482.7561	5101.4223	0.0000	1.0000
SYM	60	-2510.9005	5169.6336	0.0000	1.0000
GTR	63	-2507.9287	5172.8727	0.0000	1.0000
HKY	59	-2517.7748	5180.3679	0.0000	1.0000
K2P	56	-2525.7250	5187.3602	0.0000	1.0000
F81	58	-2588.0637	5317.9539	0.0000	1.0000
JC69	55	-2602.7766	5338.5384	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
SYMIG	62	-2344.6647	4813.3295	0.5666	0.5666
SYMG	61	-2346.4767	4814.9533	0.2516	0.8181
GTRIG	65	-2343.1544	4816.3088	0.1277	0.9459
GTRG	64	-2345.0126	4818.0252	0.0541	1.0000
HKYIG	61	-2358.6991	4839.3983	0.0000	1.0000
HKYG	60	-2360.9451	4841.8901	0.0000	1.0000
K2PIG	58	-2366.3522	4848.7044	0.0000	1.0000
K2PG	57	-2369.3197	4852.6394	0.0000	1.0000
GTRI	64	-2371.7389	4871.4777	0.0000	1.0000
SYMI	61	-2376.5741	4875.1481	0.0000	1.0000
HKYI	60	-2398.0776	4916.1552	0.0000	1.0000
K2PI	57	-2407.2179	4928.4358	0.0000	1.0000
F81IG	60	-2434.5774	4989.1548	0.0000	1.0000
F81G	59	-2436.9509	4991.9018	0.0000	1.0000
JC69IG	57	-2445.6636	5005.3271	0.0000	1.0000
JC69G	56	-2447.9643	5007.9287	0.0000	1.0000
F81I	59	-2471.6172	5061.2344	0.0000	1.0000
JC69I	56	-2482.7561	5077.5122	0.0000	1.0000
SYM	60	-2510.9005	5141.8009	0.0000	1.0000
GTR	63	-2507.9287	5141.8573	0.0000	1.0000
HKY	59	-2517.7748	5153.5497	0.0000	1.0000
K2P	56	-2525.7250	5163.4501	0.0000	1.0000
F81	58	-2588.0637	5292.1274	0.0000	1.0000
JC69	55	-2602.7766	5315.5533	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
SYMG	61	-2346.4767	5045.5787	0.7456	0.7456
SYMIG	62	-2344.6647	5047.7356	0.2536	0.9992
GTRG	64	-2345.0126	5059.9928	0.0006	0.9998
GTRIG	65	-2343.1544	5062.0571	0.0002	1.0000
K2PIG	58	-2366.3522	5067.9875	0.0000	1.0000
K2PG	57	-2369.3197	5068.1418	0.0000	1.0000
HKYG	60	-2360.9451	5068.7348	0.0000	1.0000
HKYIG	61	-2358.6991	5070.0237	0.0000	1.0000
SYMI	61	-2376.5741	5105.7735	0.0000	1.0000
GTRI	64	-2371.7389	5113.4453	0.0000	1.0000
HKYI	60	-2398.0776	5142.9998	0.0000	1.0000
K2PI	57	-2407.2179	5143.9382	0.0000	1.0000
F81G	59	-2436.9509	5214.9657	0.0000	1.0000
F81IG	60	-2434.5774	5215.9995	0.0000	1.0000
JC69G	56	-2447.9643	5219.6503	0.0000	1.0000
JC69IG	57	-2445.6636	5220.8295	0.0000	1.0000
F81I	59	-2471.6172	5284.2982	0.0000	1.0000
JC69I	56	-2482.7561	5289.2339	0.0000	1.0000
SYM	60	-2510.9005	5368.6455	0.0000	1.0000
K2P	56	-2525.7250	5375.1717	0.0000	1.0000
HKY	59	-2517.7748	5376.6135	0.0000	1.0000
GTR	63	-2507.9287	5380.0442	0.0000	1.0000
F81	58	-2588.0637	5511.4106	0.0000	1.0000
JC69	55	-2602.7766	5523.4942	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_13064.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr2_13064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_13064.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr2_13064.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 Prset applyto=(1) revmatpr=Dirichlet(1.0,1.0,1.0,1.0,1.0,1.0) statefreqpr=Fixed(Equal) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr2_13064.phylip = 1 - 324;
 Partition Dummy = 1:183-phylip/chr2_13064.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=Fixed(Equal);
END;


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

