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

Input data from file "261-aggressive-phylip/813_hg18.phylip" (29 taxa, 1258 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	64	-3897.1512	7929.2763	0.9998	0.9998
GTRIG	65	-3904.3953	7945.9886	0.0002	1.0000
HKYG	60	-3919.2983	7964.7118	0.0000	1.0000
HKYIG	61	-3925.3986	7979.1215	0.0000	1.0000
SYMG	61	-3990.7119	8109.7482	0.0000	1.0000
K2PG	57	-3995.9107	8111.3314	0.0000	1.0000
K2PIG	58	-3999.6080	8120.9241	0.0000	1.0000
SYMIG	62	-3999.3005	8129.1382	0.0000	1.0000
GTRI	64	-4005.1766	8145.3273	0.0000	1.0000
HKYI	60	-4026.8383	8179.7919	0.0000	1.0000
GTR	63	-4025.0132	8182.7801	0.0000	1.0000
F81G	59	-4031.5637	8187.0373	0.0000	1.0000
F81IG	60	-4031.5717	8189.2587	0.0000	1.0000
HKY	59	-4067.9913	8259.8924	0.0000	1.0000
SYMI	61	-4095.0198	8318.3640	0.0000	1.0000
K2PI	57	-4104.6820	8328.8741	0.0000	1.0000
SYM	60	-4108.0898	8342.2949	0.0000	1.0000
JC69G	56	-4115.5334	8348.3823	0.0000	1.0000
JC69IG	57	-4115.1599	8349.8299	0.0000	1.0000
F81I	59	-4127.9572	8379.8242	0.0000	1.0000
K2P	56	-4134.7011	8386.7177	0.0000	1.0000
F81	58	-4146.7970	8415.3021	0.0000	1.0000
JC69I	56	-4210.3282	8537.9721	0.0000	1.0000
JC69	55	-4236.0316	8587.1880	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3897.1512	7922.3023	0.9997	0.9997
GTRIG	65	-3904.3953	7938.7906	0.0003	1.0000
HKYG	60	-3919.2983	7958.5966	0.0000	1.0000
HKYIG	61	-3925.3986	7972.7971	0.0000	1.0000
SYMG	61	-3990.7119	8103.4238	0.0000	1.0000
K2PG	57	-3995.9107	8105.8214	0.0000	1.0000
K2PIG	58	-3999.6080	8115.2160	0.0000	1.0000
SYMIG	62	-3999.3005	8122.6010	0.0000	1.0000
GTRI	64	-4005.1766	8138.3533	0.0000	1.0000
HKYI	60	-4026.8383	8173.6766	0.0000	1.0000
GTR	63	-4025.0132	8176.0263	0.0000	1.0000
F81G	59	-4031.5637	8181.1275	0.0000	1.0000
F81IG	60	-4031.5717	8183.1434	0.0000	1.0000
HKY	59	-4067.9913	8253.9825	0.0000	1.0000
SYMI	61	-4095.0198	8312.0396	0.0000	1.0000
K2PI	57	-4104.6820	8323.3641	0.0000	1.0000
SYM	60	-4108.0898	8336.1796	0.0000	1.0000
JC69G	56	-4115.5334	8343.0668	0.0000	1.0000
JC69IG	57	-4115.1599	8344.3199	0.0000	1.0000
F81I	59	-4127.9572	8373.9144	0.0000	1.0000
K2P	56	-4134.7011	8381.4022	0.0000	1.0000
F81	58	-4146.7970	8409.5940	0.0000	1.0000
JC69I	56	-4210.3282	8532.6565	0.0000	1.0000
JC69	55	-4236.0316	8582.0632	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3897.1512	8251.0881	0.9996	0.9996
HKYG	60	-3919.2983	8266.8333	0.0004	1.0000
GTRIG	65	-3904.3953	8272.7137	0.0000	1.0000
HKYIG	61	-3925.3986	8286.1711	0.0000	1.0000
K2PG	57	-3995.9107	8398.6463	0.0000	1.0000
K2PIG	58	-3999.6080	8413.1782	0.0000	1.0000
SYMG	61	-3990.7119	8416.7978	0.0000	1.0000
SYMIG	62	-3999.3005	8441.1123	0.0000	1.0000
GTRI	64	-4005.1766	8467.1391	0.0000	1.0000
HKYI	60	-4026.8383	8481.9133	0.0000	1.0000
F81G	59	-4031.5637	8484.2269	0.0000	1.0000
F81IG	60	-4031.5717	8491.3801	0.0000	1.0000
GTR	63	-4025.0132	8499.6748	0.0000	1.0000
HKY	59	-4067.9913	8557.0819	0.0000	1.0000
K2PI	57	-4104.6820	8616.1889	0.0000	1.0000
SYMI	61	-4095.0198	8625.4135	0.0000	1.0000
JC69G	56	-4115.5334	8630.7544	0.0000	1.0000
JC69IG	57	-4115.1599	8637.1447	0.0000	1.0000
SYM	60	-4108.0898	8644.4163	0.0000	1.0000
K2P	56	-4134.7011	8669.0898	0.0000	1.0000
F81I	59	-4127.9572	8677.0138	0.0000	1.0000
F81	58	-4146.7970	8707.5561	0.0000	1.0000
JC69I	56	-4210.3282	8820.3441	0.0000	1.0000
JC69	55	-4236.0316	8864.6135	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 261-aggressive-phylip/813_hg18.phylip = 1 - 1258;
 Partition Dummy = 1:261-aggressive-phylip/813_hg18.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 261-aggressive-phylip/813_hg18.phylip = 1 - 1258;
 Partition Dummy = 1:261-aggressive-phylip/813_hg18.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 261-aggressive-phylip/813_hg18.phylip = 1 - 1258;
 Partition Dummy = 1:261-aggressive-phylip/813_hg18.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

