
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
11:10 08/01/2011

Input data from file "261-aggressive-phylip/11790_hg18.phylip" (29 taxa, 1172 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYIG


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
GTRIG	65	-6897.4008	13932.5592	0.9936	0.9936
GTRG	64	-6903.7283	13942.9723	0.0054	0.9990
HKYIG	61	-6908.7924	13946.3993	0.0010	1.0000
HKYG	60	-6913.8469	13954.2825	0.0000	1.0000
GTRI	64	-6953.4414	14042.3986	0.0000	1.0000
HKYI	60	-6965.8157	14058.2200	0.0000	1.0000
SYMIG	62	-6966.0565	14063.1572	0.0000	1.0000
SYMG	61	-6973.3543	14075.5231	0.0000	1.0000
K2PIG	58	-6981.6664	14085.4819	0.0000	1.0000
K2PG	57	-6986.9060	14093.7473	0.0000	1.0000
SYMI	61	-7014.2246	14157.2635	0.0000	1.0000
K2PI	57	-7029.2672	14178.4698	0.0000	1.0000
F81IG	60	-7136.7394	14400.0675	0.0000	1.0000
F81G	59	-7141.7777	14407.9222	0.0000	1.0000
JC69IG	57	-7174.3189	14468.5731	0.0000	1.0000
F81I	59	-7180.4355	14485.2379	0.0000	1.0000
JC69G	56	-7185.4995	14488.7246	0.0000	1.0000
JC69I	56	-7223.5715	14564.8685	0.0000	1.0000
GTR	63	-7328.5071	14790.2921	0.0000	1.0000
HKY	59	-7340.7494	14805.8657	0.0000	1.0000
SYM	60	-7391.2539	14909.0964	0.0000	1.0000
K2P	56	-7403.3842	14924.4940	0.0000	1.0000
F81	58	-7542.8017	15207.7525	0.0000	1.0000
JC69	55	-7588.5212	15292.5621	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-6897.4008	13924.8015	0.9945	0.9945
GTRG	64	-6903.7283	13935.4565	0.0048	0.9994
HKYIG	61	-6908.7924	13939.5848	0.0006	1.0000
HKYG	60	-6913.8469	13947.6938	0.0000	1.0000
GTRI	64	-6953.4414	14034.8828	0.0000	1.0000
HKYI	60	-6965.8157	14051.6313	0.0000	1.0000
SYMIG	62	-6966.0565	14056.1130	0.0000	1.0000
SYMG	61	-6973.3543	14068.7087	0.0000	1.0000
K2PIG	58	-6981.6664	14079.3328	0.0000	1.0000
K2PG	57	-6986.9060	14087.8119	0.0000	1.0000
SYMI	61	-7014.2246	14150.4491	0.0000	1.0000
K2PI	57	-7029.2672	14172.5344	0.0000	1.0000
F81IG	60	-7136.7394	14393.4788	0.0000	1.0000
F81G	59	-7141.7777	14401.5553	0.0000	1.0000
JC69IG	57	-7174.3189	14462.6377	0.0000	1.0000
F81I	59	-7180.4355	14478.8710	0.0000	1.0000
JC69G	56	-7185.4995	14482.9990	0.0000	1.0000
JC69I	56	-7223.5715	14559.1429	0.0000	1.0000
GTR	63	-7328.5071	14783.0141	0.0000	1.0000
HKY	59	-7340.7494	14799.4988	0.0000	1.0000
SYM	60	-7391.2539	14902.5077	0.0000	1.0000
K2P	56	-7403.3842	14918.7685	0.0000	1.0000
F81	58	-7542.8017	15201.6034	0.0000	1.0000
JC69	55	-7588.5212	15287.0424	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYIG	61	-6908.7924	14248.6393	0.7771	0.7771
HKYG	60	-6913.8469	14251.6818	0.1697	0.9468
GTRIG	65	-6897.4008	14254.1219	0.0501	0.9969
GTRG	64	-6903.7283	14259.7104	0.0031	1.0000
HKYI	60	-6965.8157	14355.6194	0.0000	1.0000
GTRI	64	-6953.4414	14359.1367	0.0000	1.0000
SYMIG	62	-6966.0565	14370.2340	0.0000	1.0000
K2PIG	58	-6981.6664	14373.1879	0.0000	1.0000
K2PG	57	-6986.9060	14376.6006	0.0000	1.0000
SYMG	61	-6973.3543	14377.7631	0.0000	1.0000
SYMI	61	-7014.2246	14459.5036	0.0000	1.0000
K2PI	57	-7029.2672	14461.3230	0.0000	1.0000
F81IG	60	-7136.7394	14697.4668	0.0000	1.0000
F81G	59	-7141.7777	14700.4769	0.0000	1.0000
JC69IG	57	-7174.3189	14751.4263	0.0000	1.0000
JC69G	56	-7185.4995	14766.7212	0.0000	1.0000
F81I	59	-7180.4355	14777.7926	0.0000	1.0000
JC69I	56	-7223.5715	14842.8651	0.0000	1.0000
HKY	59	-7340.7494	15098.4203	0.0000	1.0000
GTR	63	-7328.5071	15102.2015	0.0000	1.0000
K2P	56	-7403.3842	15202.4906	0.0000	1.0000
SYM	60	-7391.2539	15206.4958	0.0000	1.0000
F81	58	-7542.8017	15495.4584	0.0000	1.0000
JC69	55	-7588.5212	15565.6980	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11790_hg18.phylip = 1 - 1172;
 Partition Dummy = 1:261-aggressive-phylip/11790_hg18.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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11790_hg18.phylip = 1 - 1172;
 Partition Dummy = 1:261-aggressive-phylip/11790_hg18.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=Dirichlet(1.0,1.0,1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/11790_hg18.phylip = 1 - 1172;
 Partition Dummy = 1:261-aggressive-phylip/11790_hg18.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;


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

