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

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

Minimum AIC  model: HKYIG
Minimum AICc model: HKYIG
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
HKYIG	61	-2395.2966	4925.0135	0.5734	0.5734
HKYG	60	-2397.5531	4927.1062	0.2014	0.7748
K2PIG	58	-2400.2756	4927.7341	0.1471	0.9219
K2PG	57	-2402.6473	4930.0809	0.0455	0.9674
SYMIG	62	-2397.6087	4932.0660	0.0169	0.9843
SYMG	61	-2399.8256	4934.0715	0.0062	0.9904
GTRIG	65	-2395.0265	4934.2348	0.0057	0.9961
HKYI	60	-2401.7784	4935.5567	0.0029	0.9991
K2PI	57	-2407.0060	4938.7983	0.0006	0.9997
GTRG	64	-2399.1345	4939.9984	0.0003	1.0000
SYMI	61	-2406.9111	4948.2426	0.0000	1.0000
GTRI	64	-2403.5395	4948.8084	0.0000	1.0000
F81IG	60	-2448.6887	5029.3774	0.0000	1.0000
F81G	59	-2450.4904	5030.5684	0.0000	1.0000
JC69IG	57	-2453.5334	5031.8530	0.0000	1.0000
JC69G	56	-2455.3786	5033.1546	0.0000	1.0000
F81I	59	-2454.7339	5039.0554	0.0000	1.0000
JC69I	56	-2459.7540	5041.9054	0.0000	1.0000
HKY	59	-2462.1611	5053.9097	0.0000	1.0000
K2P	56	-2467.4585	5057.3143	0.0000	1.0000
GTR	63	-2464.2133	5067.7116	0.0000	1.0000
SYM	60	-2468.0412	5068.0825	0.0000	1.0000
F81	58	-2514.9609	5157.1048	0.0000	1.0000
JC69	55	-2520.1350	5160.2862	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-2395.2966	4912.5932	0.6515	0.6515
HKYG	60	-2397.5531	4915.1062	0.1855	0.8370
K2PIG	58	-2400.2756	4916.5511	0.0900	0.9270
SYMIG	62	-2397.6087	4919.2174	0.0237	0.9508
K2PG	57	-2402.6473	4919.2946	0.0228	0.9736
GTRIG	65	-2395.0265	4920.0530	0.0156	0.9893
SYMG	61	-2399.8256	4921.6512	0.0070	0.9963
HKYI	60	-2401.7784	4923.5567	0.0027	0.9990
GTRG	64	-2399.1345	4926.2690	0.0007	0.9997
K2PI	57	-2407.0060	4928.0120	0.0003	1.0000
GTRI	64	-2403.5395	4935.0790	0.0000	1.0000
SYMI	61	-2406.9111	4935.8222	0.0000	1.0000
F81IG	60	-2448.6887	5017.3774	0.0000	1.0000
F81G	59	-2450.4904	5018.9809	0.0000	1.0000
JC69IG	57	-2453.5334	5021.0667	0.0000	1.0000
JC69G	56	-2455.3786	5022.7572	0.0000	1.0000
F81I	59	-2454.7339	5027.4678	0.0000	1.0000
JC69I	56	-2459.7540	5031.5080	0.0000	1.0000
HKY	59	-2462.1611	5042.3222	0.0000	1.0000
K2P	56	-2467.4585	5046.9169	0.0000	1.0000
GTR	63	-2464.2133	5054.4266	0.0000	1.0000
SYM	60	-2468.0412	5056.0825	0.0000	1.0000
F81	58	-2514.9609	5145.9218	0.0000	1.0000
JC69	55	-2520.1350	5150.2699	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-2402.6473	5176.2944	0.6947	0.6947
K2PIG	58	-2400.2756	5178.0597	0.2874	0.9821
K2PI	57	-2407.0060	5185.0119	0.0089	0.9909
HKYG	60	-2397.5531	5185.6323	0.0065	0.9975
HKYIG	61	-2395.2966	5187.6281	0.0024	0.9999
HKYI	60	-2401.7784	5194.0829	0.0001	1.0000
SYMG	61	-2399.8256	5196.6861	0.0000	1.0000
SYMIG	62	-2397.6087	5198.7610	0.0000	1.0000
SYMI	61	-2406.9111	5210.8571	0.0000	1.0000
GTRIG	65	-2395.0265	5213.1230	0.0000	1.0000
GTRG	64	-2399.1345	5214.8302	0.0000	1.0000
GTRI	64	-2403.5395	5223.6402	0.0000	1.0000
JC69G	56	-2455.3786	5275.2483	0.0000	1.0000
JC69IG	57	-2453.5334	5278.0666	0.0000	1.0000
JC69I	56	-2459.7540	5283.9991	0.0000	1.0000
F81G	59	-2450.4904	5284.9982	0.0000	1.0000
F81IG	60	-2448.6887	5287.9036	0.0000	1.0000
F81I	59	-2454.7339	5293.4852	0.0000	1.0000
K2P	56	-2467.4585	5299.4080	0.0000	1.0000
HKY	59	-2462.1611	5308.3395	0.0000	1.0000
SYM	60	-2468.0412	5326.6086	0.0000	1.0000
GTR	63	-2464.2133	5338.4791	0.0000	1.0000
JC69	55	-2520.1350	5398.2522	0.0000	1.0000
F81	58	-2514.9609	5407.4304	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr9_6414.phylip = 1 - 671;
 Partition Dummy = 1:183-phylip/chr9_6414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr9_6414.phylip = 1 - 671;
 Partition Dummy = 1:183-phylip/chr9_6414.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr9_6414.phylip = 1 - 671;
 Partition Dummy = 1:183-phylip/chr9_6414.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

