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

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

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


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	-4079.2689	8297.2045	1.0000	1.0000
GTRG	64	-4091.8347	8320.0650	0.0000	1.0000
HKYIG	61	-4107.7354	8345.0804	0.0000	1.0000
SYMIG	62	-4111.9776	8355.8222	0.0000	1.0000
HKYG	60	-4119.2816	8365.9201	0.0000	1.0000
K2PIG	58	-4126.8827	8376.6301	0.0000	1.0000
SYMG	61	-4124.7995	8379.2087	0.0000	1.0000
K2PG	57	-4138.7824	8398.1900	0.0000	1.0000
GTRI	64	-4153.7825	8443.9606	0.0000	1.0000
HKYI	60	-4181.9327	8491.2222	0.0000	1.0000
SYMI	61	-4186.6599	8502.9294	0.0000	1.0000
K2PI	57	-4202.1195	8524.8642	0.0000	1.0000
GTR	63	-4196.6629	8527.4548	0.0000	1.0000
HKY	59	-4223.3292	8571.7669	0.0000	1.0000
SYM	60	-4229.6859	8586.7287	0.0000	1.0000
F81IG	60	-4239.9196	8607.1960	0.0000	1.0000
K2P	56	-4245.6380	8609.6664	0.0000	1.0000
F81G	59	-4253.8127	8632.7338	0.0000	1.0000
JC69IG	57	-4256.3830	8633.3912	0.0000	1.0000
JC69G	56	-4270.3604	8659.1112	0.0000	1.0000
F81I	59	-4313.5060	8752.1205	0.0000	1.0000
JC69I	56	-4330.2074	8778.8052	0.0000	1.0000
F81	58	-4357.3371	8837.5387	0.0000	1.0000
JC69	55	-4374.3197	8864.7994	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-4079.2689	8288.5378	1.0000	1.0000
GTRG	64	-4091.8347	8311.6694	0.0000	1.0000
HKYIG	61	-4107.7354	8337.4707	0.0000	1.0000
SYMIG	62	-4111.9776	8347.9552	0.0000	1.0000
HKYG	60	-4119.2816	8358.5633	0.0000	1.0000
K2PIG	58	-4126.8827	8369.7655	0.0000	1.0000
SYMG	61	-4124.7995	8371.5990	0.0000	1.0000
K2PG	57	-4138.7824	8391.5647	0.0000	1.0000
GTRI	64	-4153.7825	8435.5650	0.0000	1.0000
HKYI	60	-4181.9327	8483.8654	0.0000	1.0000
SYMI	61	-4186.6599	8495.3198	0.0000	1.0000
K2PI	57	-4202.1195	8518.2390	0.0000	1.0000
GTR	63	-4196.6629	8519.3257	0.0000	1.0000
HKY	59	-4223.3292	8564.6584	0.0000	1.0000
SYM	60	-4229.6859	8579.3719	0.0000	1.0000
F81IG	60	-4239.9196	8599.8392	0.0000	1.0000
K2P	56	-4245.6380	8603.2760	0.0000	1.0000
F81G	59	-4253.8127	8625.6253	0.0000	1.0000
JC69IG	57	-4256.3830	8626.7659	0.0000	1.0000
JC69G	56	-4270.3604	8652.7208	0.0000	1.0000
F81I	59	-4313.5060	8745.0121	0.0000	1.0000
JC69I	56	-4330.2074	8772.4148	0.0000	1.0000
F81	58	-4357.3371	8830.6741	0.0000	1.0000
JC69	55	-4374.3197	8858.6394	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-4079.2689	8611.0836	0.9999	0.9999
GTRG	64	-4091.8347	8629.2530	0.0001	1.0000
HKYIG	61	-4107.7354	8640.1676	0.0000	1.0000
SYMIG	62	-4111.9776	8655.6143	0.0000	1.0000
HKYG	60	-4119.2816	8656.2979	0.0000	1.0000
K2PIG	58	-4126.8827	8657.5756	0.0000	1.0000
SYMG	61	-4124.7995	8674.2959	0.0000	1.0000
K2PG	57	-4138.7824	8674.4126	0.0000	1.0000
GTRI	64	-4153.7825	8753.1486	0.0000	1.0000
HKYI	60	-4181.9327	8781.6000	0.0000	1.0000
SYMI	61	-4186.6599	8798.0166	0.0000	1.0000
K2PI	57	-4202.1195	8801.0869	0.0000	1.0000
GTR	63	-4196.6629	8831.9471	0.0000	1.0000
HKY	59	-4223.3292	8857.4308	0.0000	1.0000
SYM	60	-4229.6859	8877.1065	0.0000	1.0000
K2P	56	-4245.6380	8881.1617	0.0000	1.0000
F81IG	60	-4239.9196	8897.5738	0.0000	1.0000
JC69IG	57	-4256.3830	8909.6138	0.0000	1.0000
F81G	59	-4253.8127	8918.3977	0.0000	1.0000
JC69G	56	-4270.3604	8930.6065	0.0000	1.0000
F81I	59	-4313.5060	9037.7844	0.0000	1.0000
JC69I	56	-4330.2074	9050.3005	0.0000	1.0000
F81	58	-4357.3371	9118.4843	0.0000	1.0000
JC69	55	-4374.3197	9131.5628	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/5049_hg18.phylip = 1 - 1056;
 Partition Dummy = 1:261-aggressive-phylip/5049_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/5049_hg18.phylip = 1 - 1056;
 Partition Dummy = 1:261-aggressive-phylip/5049_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 (GTRIG)]

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


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

