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

Input data from file "261-aggressive-phylip/7050_hg18.phylip" (29 taxa, 494 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	-3130.8399	6409.0738	0.7293	0.7293
GTRIG	65	-3130.5123	6411.0714	0.2686	0.9979
SYMG	61	-3141.0700	6421.6493	0.0014	0.9993
SYMIG	62	-3140.5921	6423.3095	0.0006	0.9999
HKYG	60	-3144.9039	6426.7132	0.0001	1.0000
HKYIG	61	-3144.8172	6429.1436	0.0000	1.0000
GTRI	64	-3153.1291	6453.6522	0.0000	1.0000
K2PG	57	-3164.4772	6458.1195	0.0000	1.0000
K2PIG	58	-3164.4772	6460.6877	0.0000	1.0000
HKYI	60	-3164.1785	6465.2623	0.0000	1.0000
SYMI	61	-3164.8816	6469.2724	0.0000	1.0000
K2PI	57	-3184.2882	6497.7416	0.0000	1.0000
F81G	59	-3206.7297	6547.7728	0.0000	1.0000
F81IG	60	-3207.0422	6550.9897	0.0000	1.0000
F81I	59	-3227.9194	6590.1522	0.0000	1.0000
JC69G	56	-3239.9517	6606.5122	0.0000	1.0000
JC69IG	57	-3239.9518	6609.0686	0.0000	1.0000
JC69I	56	-3257.1325	6640.8736	0.0000	1.0000
GTR	63	-3283.5576	6711.8687	0.0000	1.0000
SYM	60	-3291.9884	6720.8821	0.0000	1.0000
HKY	59	-3304.2013	6742.7159	0.0000	1.0000
K2P	56	-3323.0861	6772.7810	0.0000	1.0000
F81	58	-3368.9155	6869.5644	0.0000	1.0000
JC69	55	-3394.9849	6914.0338	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3130.8399	6389.6798	0.6615	0.6615
GTRIG	65	-3130.5123	6391.0247	0.3377	0.9992
SYMG	61	-3141.0700	6404.1401	0.0005	0.9997
SYMIG	62	-3140.5921	6405.1842	0.0003	1.0000
HKYG	60	-3144.9039	6409.8079	0.0000	1.0000
HKYIG	61	-3144.8172	6411.6344	0.0000	1.0000
GTRI	64	-3153.1291	6434.2583	0.0000	1.0000
K2PG	57	-3164.4772	6442.9543	0.0000	1.0000
K2PIG	58	-3164.4772	6444.9544	0.0000	1.0000
HKYI	60	-3164.1785	6448.3570	0.0000	1.0000
SYMI	61	-3164.8816	6451.7631	0.0000	1.0000
K2PI	57	-3184.2882	6482.5765	0.0000	1.0000
F81G	59	-3206.7297	6531.4595	0.0000	1.0000
F81IG	60	-3207.0422	6534.0844	0.0000	1.0000
F81I	59	-3227.9194	6573.8388	0.0000	1.0000
JC69G	56	-3239.9517	6591.9035	0.0000	1.0000
JC69IG	57	-3239.9518	6593.9035	0.0000	1.0000
JC69I	56	-3257.1325	6626.2649	0.0000	1.0000
GTR	63	-3283.5576	6693.1152	0.0000	1.0000
SYM	60	-3291.9884	6703.9768	0.0000	1.0000
HKY	59	-3304.2013	6726.4025	0.0000	1.0000
K2P	56	-3323.0861	6758.1723	0.0000	1.0000
F81	58	-3368.9155	6853.8311	0.0000	1.0000
JC69	55	-3394.9849	6899.9699	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3130.8399	6658.6421	0.5928	0.5928
SYMG	61	-3141.0700	6660.4947	0.2348	0.8276
HKYG	60	-3144.9039	6661.9600	0.1128	0.9404
GTRIG	65	-3130.5123	6664.1895	0.0370	0.9774
SYMIG	62	-3140.5921	6665.7414	0.0170	0.9945
HKYIG	61	-3144.8172	6667.9890	0.0055	1.0000
K2PG	57	-3164.4772	6682.4988	0.0000	1.0000
K2PIG	58	-3164.4772	6688.7015	0.0000	1.0000
HKYI	60	-3164.1785	6700.5091	0.0000	1.0000
GTRI	64	-3153.1291	6703.2205	0.0000	1.0000
SYMI	61	-3164.8816	6708.1178	0.0000	1.0000
K2PI	57	-3184.2882	6722.1210	0.0000	1.0000
F81G	59	-3206.7297	6779.4091	0.0000	1.0000
F81IG	60	-3207.0422	6786.2365	0.0000	1.0000
F81I	59	-3227.9194	6821.7884	0.0000	1.0000
JC69G	56	-3239.9517	6827.2455	0.0000	1.0000
JC69IG	57	-3239.9518	6833.4480	0.0000	1.0000
JC69I	56	-3257.1325	6861.6069	0.0000	1.0000
SYM	60	-3291.9884	6956.1290	0.0000	1.0000
GTR	63	-3283.5576	6957.8749	0.0000	1.0000
HKY	59	-3304.2013	6974.3521	0.0000	1.0000
K2P	56	-3323.0861	6993.5143	0.0000	1.0000
F81	58	-3368.9155	7097.5781	0.0000	1.0000
JC69	55	-3394.9849	7131.1093	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/7050_hg18.phylip = 1 - 494;
 Partition Dummy = 1:261-aggressive-phylip/7050_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/7050_hg18.phylip = 1 - 494;
 Partition Dummy = 1:261-aggressive-phylip/7050_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/7050_hg18.phylip = 1 - 494;
 Partition Dummy = 1:261-aggressive-phylip/7050_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

