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

Input data from file "261-aggressive-phylip/11284_hg18.phylip" (29 taxa, 1180 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	-3019.6139	6174.6896	0.7564	0.7564
GTRIG	65	-3019.6271	6176.9561	0.2436	1.0000
HKYG	60	-3039.4125	6205.3666	0.0000	1.0000
HKYIG	61	-3039.4207	6207.6071	0.0000	1.0000
GTR	63	-3070.7584	6274.7426	0.0000	1.0000
GTRI	64	-3070.7645	6276.9908	0.0000	1.0000
F81G	59	-3079.6279	6283.5771	0.0000	1.0000
F81IG	60	-3079.6410	6285.8235	0.0000	1.0000
HKY	59	-3089.5630	6303.4474	0.0000	1.0000
HKYI	60	-3089.5691	6305.6797	0.0000	1.0000
F81	58	-3124.0431	6370.1914	0.0000	1.0000
F81I	59	-3127.0966	6378.5147	0.0000	1.0000
SYMG	61	-3133.3509	6395.4675	0.0000	1.0000
SYMIG	62	-3133.3641	6397.7219	0.0000	1.0000
K2PG	57	-3146.9545	6413.8021	0.0000	1.0000
K2PIG	58	-3146.9675	6416.0403	0.0000	1.0000
JC69G	56	-3171.6617	6461.0082	0.0000	1.0000
JC69IG	57	-3173.5410	6466.9750	0.0000	1.0000
SYM	60	-3176.9209	6480.3834	0.0000	1.0000
SYMI	61	-3176.9274	6482.6204	0.0000	1.0000
K2P	56	-3192.9894	6503.6636	0.0000	1.0000
K2PI	57	-3192.9956	6505.8843	0.0000	1.0000
JC69	55	-3218.0600	6551.6004	0.0000	1.0000
JC69I	56	-3219.2978	6556.2804	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3019.6139	6167.2277	0.7336	0.7336
GTRIG	65	-3019.6271	6169.2542	0.2664	1.0000
HKYG	60	-3039.4125	6198.8251	0.0000	1.0000
HKYIG	61	-3039.4207	6200.8414	0.0000	1.0000
GTR	63	-3070.7584	6267.5168	0.0000	1.0000
GTRI	64	-3070.7645	6269.5289	0.0000	1.0000
F81G	59	-3079.6279	6277.2557	0.0000	1.0000
F81IG	60	-3079.6410	6279.2819	0.0000	1.0000
HKY	59	-3089.5630	6297.1260	0.0000	1.0000
HKYI	60	-3089.5691	6299.1381	0.0000	1.0000
F81	58	-3124.0431	6364.0862	0.0000	1.0000
F81I	59	-3127.0966	6372.1933	0.0000	1.0000
SYMG	61	-3133.3509	6388.7019	0.0000	1.0000
SYMIG	62	-3133.3641	6390.7282	0.0000	1.0000
K2PG	57	-3146.9545	6407.9091	0.0000	1.0000
K2PIG	58	-3146.9675	6409.9350	0.0000	1.0000
JC69G	56	-3171.6617	6455.3235	0.0000	1.0000
JC69IG	57	-3173.5410	6461.0820	0.0000	1.0000
SYM	60	-3176.9209	6473.8418	0.0000	1.0000
SYMI	61	-3176.9274	6475.8548	0.0000	1.0000
K2P	56	-3192.9894	6497.9788	0.0000	1.0000
K2PI	57	-3192.9956	6499.9912	0.0000	1.0000
JC69	55	-3218.0600	6546.1200	0.0000	1.0000
JC69I	56	-3219.2978	6550.5956	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-3019.6139	6491.9170	0.9687	0.9687
GTRIG	65	-3019.6271	6499.0167	0.0278	0.9965
HKYG	60	-3039.4125	6503.2213	0.0034	0.9999
HKYIG	61	-3039.4207	6510.3109	0.0001	1.0000
F81G	59	-3079.6279	6576.5786	0.0000	1.0000
F81IG	60	-3079.6410	6583.6781	0.0000	1.0000
GTR	63	-3070.7584	6587.1328	0.0000	1.0000
GTRI	64	-3070.7645	6594.2182	0.0000	1.0000
HKY	59	-3089.5630	6596.4489	0.0000	1.0000
HKYI	60	-3089.5691	6603.5343	0.0000	1.0000
F81	58	-3124.0431	6658.3358	0.0000	1.0000
F81I	59	-3127.0966	6671.5162	0.0000	1.0000
K2PG	57	-3146.9545	6697.0855	0.0000	1.0000
SYMG	61	-3133.3509	6698.1713	0.0000	1.0000
K2PIG	58	-3146.9675	6704.1847	0.0000	1.0000
SYMIG	62	-3133.3641	6705.2709	0.0000	1.0000
JC69G	56	-3171.6617	6739.4266	0.0000	1.0000
JC69IG	57	-3173.5410	6750.2584	0.0000	1.0000
SYM	60	-3176.9209	6778.2380	0.0000	1.0000
K2P	56	-3192.9894	6782.0819	0.0000	1.0000
SYMI	61	-3176.9274	6785.3242	0.0000	1.0000
K2PI	57	-3192.9956	6789.1676	0.0000	1.0000
JC69	55	-3218.0600	6825.1498	0.0000	1.0000
JC69I	56	-3219.2978	6834.6987	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/11284_hg18.phylip = 1 - 1180;
 Partition Dummy = 1:261-aggressive-phylip/11284_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/11284_hg18.phylip = 1 - 1180;
 Partition Dummy = 1:261-aggressive-phylip/11284_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/11284_hg18.phylip = 1 - 1180;
 Partition Dummy = 1:261-aggressive-phylip/11284_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

