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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: HKYG


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	-1165.7141	2482.5394	0.3742	0.3742
GTRI	64	-1166.1733	2483.4578	0.2364	0.6106
HKYG	60	-1172.1350	2484.3799	0.1491	0.7597
GTRIG	65	-1165.4914	2484.8825	0.1160	0.8757
HKYI	60	-1172.8492	2485.8083	0.0730	0.9487
HKYIG	61	-1171.8374	2486.5122	0.0513	1.0000
GTR	63	-1180.3506	2509.0391	0.0000	1.0000
HKY	59	-1187.0968	2511.5908	0.0000	1.0000
K2PG	57	-1204.0471	2540.1106	0.0000	1.0000
K2PI	57	-1204.7803	2541.5770	0.0000	1.0000
SYMG	61	-1199.3929	2541.6232	0.0000	1.0000
K2PIG	58	-1203.5685	2541.8365	0.0000	1.0000
SYMI	61	-1199.8959	2542.6292	0.0000	1.0000
F81G	59	-1202.9326	2543.2625	0.0000	1.0000
SYMIG	62	-1199.0681	2543.7163	0.0000	1.0000
F81I	59	-1203.4658	2544.3288	0.0000	1.0000
F81IG	60	-1202.8613	2545.8326	0.0000	1.0000
F81	58	-1217.6603	2570.0201	0.0000	1.0000
K2P	56	-1220.3678	2570.0834	0.0000	1.0000
SYM	60	-1215.4662	2571.0424	0.0000	1.0000
JC69G	56	-1231.9719	2593.2916	0.0000	1.0000
JC69IG	57	-1231.4368	2594.8900	0.0000	1.0000
JC69I	56	-1233.0366	2595.4210	0.0000	1.0000
JC69	55	-1246.1658	2619.0254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1165.7141	2459.4282	0.4410	0.4410
GTRI	64	-1166.1733	2460.3467	0.2786	0.7195
GTRIG	65	-1165.4914	2460.9828	0.2027	0.9222
HKYG	60	-1172.1350	2464.2700	0.0392	0.9614
HKYIG	61	-1171.8374	2465.6748	0.0194	0.9808
HKYI	60	-1172.8492	2465.6984	0.0192	1.0000
GTR	63	-1180.3506	2486.7011	0.0000	1.0000
HKY	59	-1187.0968	2492.1936	0.0000	1.0000
SYMG	61	-1199.3929	2520.7857	0.0000	1.0000
SYMI	61	-1199.8959	2521.7917	0.0000	1.0000
K2PG	57	-1204.0471	2522.0942	0.0000	1.0000
SYMIG	62	-1199.0681	2522.1362	0.0000	1.0000
K2PIG	58	-1203.5685	2523.1370	0.0000	1.0000
K2PI	57	-1204.7803	2523.5607	0.0000	1.0000
F81G	59	-1202.9326	2523.8653	0.0000	1.0000
F81I	59	-1203.4658	2524.9316	0.0000	1.0000
F81IG	60	-1202.8613	2525.7227	0.0000	1.0000
SYM	60	-1215.4662	2550.9325	0.0000	1.0000
F81	58	-1217.6603	2551.3206	0.0000	1.0000
K2P	56	-1220.3678	2552.7356	0.0000	1.0000
JC69G	56	-1231.9719	2575.9438	0.0000	1.0000
JC69IG	57	-1231.4368	2576.8736	0.0000	1.0000
JC69I	56	-1233.0366	2578.0732	0.0000	1.0000
JC69	55	-1246.1658	2602.3316	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1172.1350	2707.3954	0.6407	0.6407
HKYI	60	-1172.8492	2708.8238	0.3137	0.9545
HKYIG	61	-1171.8374	2712.8522	0.0419	0.9963
GTRG	64	-1165.7141	2718.7619	0.0022	0.9985
GTRI	64	-1166.1733	2719.6804	0.0014	0.9999
GTRIG	65	-1165.4914	2724.3686	0.0001	1.0000
HKY	59	-1187.0968	2731.2668	0.0000	1.0000
GTR	63	-1180.3506	2741.9828	0.0000	1.0000
K2PG	57	-1204.0471	2753.0633	0.0000	1.0000
K2PI	57	-1204.7803	2754.5298	0.0000	1.0000
K2PIG	58	-1203.5685	2758.1582	0.0000	1.0000
F81G	59	-1202.9326	2762.9385	0.0000	1.0000
F81I	59	-1203.4658	2764.0048	0.0000	1.0000
SYMG	61	-1199.3929	2767.9632	0.0000	1.0000
F81IG	60	-1202.8613	2768.8481	0.0000	1.0000
SYMI	61	-1199.8959	2768.9692	0.0000	1.0000
SYMIG	62	-1199.0681	2773.3657	0.0000	1.0000
K2P	56	-1220.3678	2779.6526	0.0000	1.0000
F81	58	-1217.6603	2786.3418	0.0000	1.0000
SYM	60	-1215.4662	2794.0578	0.0000	1.0000
JC69G	56	-1231.9719	2802.8608	0.0000	1.0000
JC69I	56	-1233.0366	2804.9902	0.0000	1.0000
JC69IG	57	-1231.4368	2807.8427	0.0000	1.0000
JC69	55	-1246.1658	2825.1965	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/1632_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/1632_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/1632_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/1632_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 (HKYG)]

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/1632_hg18.phylip = 1 - 425;
 Partition Dummy = 1:261-aggressive-phylip/1632_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 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);
END;


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

