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

Input data from file "591-tomodel-phylip/chr2_23281.phylip" (7 taxa, 754 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	20	-2962.6330	5966.4120	0.5491	0.5491
GTRIG	21	-2962.2045	5967.6713	0.2925	0.8416
GTRI	20	-2963.8766	5968.8991	0.1583	1.0000
HKYI	16	-2977.2873	5987.3127	0.0000	1.0000
HKYG	16	-2977.4759	5987.6899	0.0000	1.0000
HKYIG	17	-2976.6132	5988.0580	0.0000	1.0000
SYMG	17	-2988.4790	6011.7895	0.0000	1.0000
SYMI	17	-2988.4802	6011.7919	0.0000	1.0000
SYMIG	18	-2987.7164	6012.3633	0.0000	1.0000
K2PI	13	-2993.4092	6013.3103	0.0000	1.0000
K2PG	13	-2993.4999	6013.4917	0.0000	1.0000
K2PIG	14	-2992.6756	6013.9195	0.0000	1.0000
GTR	19	-3010.4222	6059.8798	0.0000	1.0000
HKY	15	-3030.0335	6090.7174	0.0000	1.0000
SYM	16	-3043.1271	6118.9923	0.0000	1.0000
K2P	12	-3050.8890	6126.1991	0.0000	1.0000
F81I	15	-3062.5186	6155.6876	0.0000	1.0000
F81IG	16	-3062.5211	6157.7804	0.0000	1.0000
F81G	15	-3064.8353	6160.3209	0.0000	1.0000
JC69I	12	-3074.7583	6173.9376	0.0000	1.0000
JC69IG	13	-3074.7551	6176.0020	0.0000	1.0000
JC69G	12	-3076.8922	6178.2054	0.0000	1.0000
F81	14	-3112.4728	6253.5140	0.0000	1.0000
JC69	11	-3127.0062	6276.3683	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2962.6330	5965.2660	0.5396	0.5396
GTRIG	21	-2962.2045	5966.4090	0.3047	0.8444
GTRI	20	-2963.8766	5967.7532	0.1556	1.0000
HKYI	16	-2977.2873	5986.5746	0.0000	1.0000
HKYG	16	-2977.4759	5986.9518	0.0000	1.0000
HKYIG	17	-2976.6132	5987.2265	0.0000	1.0000
SYMG	17	-2988.4790	6010.9580	0.0000	1.0000
SYMI	17	-2988.4802	6010.9604	0.0000	1.0000
SYMIG	18	-2987.7164	6011.4327	0.0000	1.0000
K2PI	13	-2993.4092	6012.8184	0.0000	1.0000
K2PG	13	-2993.4999	6012.9998	0.0000	1.0000
K2PIG	14	-2992.6756	6013.3512	0.0000	1.0000
GTR	19	-3010.4222	6058.8443	0.0000	1.0000
HKY	15	-3030.0335	6090.0670	0.0000	1.0000
SYM	16	-3043.1271	6118.2541	0.0000	1.0000
K2P	12	-3050.8890	6125.7780	0.0000	1.0000
F81I	15	-3062.5186	6155.0372	0.0000	1.0000
F81IG	16	-3062.5211	6157.0422	0.0000	1.0000
F81G	15	-3064.8353	6159.6705	0.0000	1.0000
JC69I	12	-3074.7583	6173.5166	0.0000	1.0000
JC69IG	13	-3074.7551	6175.5101	0.0000	1.0000
JC69G	12	-3076.8922	6177.7843	0.0000	1.0000
F81	14	-3112.4728	6252.9456	0.0000	1.0000
JC69	11	-3127.0062	6276.0125	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	20	-2962.6330	6057.7738	0.5519	0.5519
GTRI	20	-2963.8766	6060.2610	0.1591	0.7110
HKYI	16	-2977.2873	6060.5809	0.1356	0.8466
HKYG	16	-2977.4759	6060.9581	0.1123	0.9589
GTRIG	21	-2962.2045	6063.5422	0.0308	0.9898
HKYIG	17	-2976.6132	6065.8581	0.0097	0.9994
K2PI	13	-2993.4092	6072.9485	0.0003	0.9997
K2PG	13	-2993.4999	6073.1299	0.0003	1.0000
K2PIG	14	-2992.6756	6078.1067	0.0000	1.0000
SYMG	17	-2988.4790	6089.5896	0.0000	1.0000
SYMI	17	-2988.4802	6089.5921	0.0000	1.0000
SYMIG	18	-2987.7164	6094.6898	0.0000	1.0000
GTR	19	-3010.4222	6146.7268	0.0000	1.0000
HKY	15	-3030.0335	6159.4479	0.0000	1.0000
K2P	12	-3050.8890	6181.2827	0.0000	1.0000
SYM	16	-3043.1271	6192.2604	0.0000	1.0000
F81I	15	-3062.5186	6224.4181	0.0000	1.0000
JC69I	12	-3074.7583	6229.0213	0.0000	1.0000
F81G	15	-3064.8353	6229.0514	0.0000	1.0000
F81IG	16	-3062.5211	6231.0485	0.0000	1.0000
JC69G	12	-3076.8922	6233.2890	0.0000	1.0000
JC69IG	13	-3074.7551	6235.6402	0.0000	1.0000
F81	14	-3112.4728	6317.7011	0.0000	1.0000
JC69	11	-3127.0062	6326.8918	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 591-tomodel-phylip/chr2_23281.phylip = 1 - 754;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23281.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 591-tomodel-phylip/chr2_23281.phylip = 1 - 754;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23281.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 591-tomodel-phylip/chr2_23281.phylip = 1 - 754;
 Partition Dummy = 1:591-tomodel-phylip/chr2_23281.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

