
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:40 02/17/2011

Input data from file "183-phylip/chr1_5460.phylip" (29 taxa, 690 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
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
GTRIG	65	-3712.4021	7568.5542	0.7254	0.7254
GTRG	64	-3714.6047	7570.5213	0.2713	0.9967
HKYIG	61	-3723.1164	7580.2774	0.0021	0.9988
HKYG	60	-3724.8519	7581.3413	0.0012	1.0000
K2PIG	58	-3747.5217	7621.8897	0.0000	1.0000
K2PG	57	-3749.5249	7623.5117	0.0000	1.0000
SYMIG	62	-3746.5555	7629.5704	0.0000	1.0000
GTRI	64	-3744.2678	7629.8476	0.0000	1.0000
SYMG	61	-3748.2778	7630.6001	0.0000	1.0000
HKYI	60	-3749.6235	7630.8846	0.0000	1.0000
K2PI	57	-3772.0363	7668.5346	0.0000	1.0000
SYMI	61	-3771.2061	7676.4568	0.0000	1.0000
F81IG	60	-3817.5219	7766.6813	0.0000	1.0000
F81G	59	-3819.9858	7769.2097	0.0000	1.0000
JC69IG	57	-3833.3802	7791.2224	0.0000	1.0000
JC69G	56	-3836.0243	7794.1339	0.0000	1.0000
F81I	59	-3841.0027	7811.2436	0.0000	1.0000
JC69I	56	-3857.2245	7836.5342	0.0000	1.0000
GTR	63	-3952.2724	8043.4266	0.0000	1.0000
HKY	59	-3961.0784	8051.3948	0.0000	1.0000
K2P	56	-3985.3348	8092.7549	0.0000	1.0000
SYM	60	-3983.8659	8099.3692	0.0000	1.0000
F81	58	-4045.7084	8218.2631	0.0000	1.0000
JC69	55	-4062.1106	8243.9373	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3712.4021	7554.8042	0.7679	0.7679
GTRG	64	-3714.6047	7557.2093	0.2307	0.9986
HKYIG	61	-3723.1164	7568.2328	0.0009	0.9996
HKYG	60	-3724.8519	7569.7037	0.0004	1.0000
K2PIG	58	-3747.5217	7611.0434	0.0000	1.0000
K2PG	57	-3749.5249	7613.0497	0.0000	1.0000
GTRI	64	-3744.2678	7616.5356	0.0000	1.0000
SYMIG	62	-3746.5555	7617.1111	0.0000	1.0000
SYMG	61	-3748.2778	7618.5555	0.0000	1.0000
HKYI	60	-3749.6235	7619.2471	0.0000	1.0000
K2PI	57	-3772.0363	7658.0725	0.0000	1.0000
SYMI	61	-3771.2061	7664.4122	0.0000	1.0000
F81IG	60	-3817.5219	7755.0437	0.0000	1.0000
F81G	59	-3819.9858	7757.9716	0.0000	1.0000
JC69IG	57	-3833.3802	7780.7604	0.0000	1.0000
JC69G	56	-3836.0243	7784.0486	0.0000	1.0000
F81I	59	-3841.0027	7800.0055	0.0000	1.0000
JC69I	56	-3857.2245	7826.4489	0.0000	1.0000
GTR	63	-3952.2724	8030.5448	0.0000	1.0000
HKY	59	-3961.0784	8040.1567	0.0000	1.0000
K2P	56	-3985.3348	8082.6696	0.0000	1.0000
SYM	60	-3983.8659	8087.7317	0.0000	1.0000
F81	58	-4045.7084	8207.4168	0.0000	1.0000
JC69	55	-4062.1106	8234.2213	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3724.8519	7841.9052	0.7719	0.7719
HKYIG	61	-3723.1164	7844.9710	0.1667	0.9385
GTRG	64	-3714.6047	7847.5576	0.0457	0.9842
GTRIG	65	-3712.4021	7849.6892	0.0158	1.0000
K2PG	57	-3749.5249	7871.6411	0.0000	1.0000
K2PIG	58	-3747.5217	7874.1715	0.0000	1.0000
HKYI	60	-3749.6235	7891.4486	0.0000	1.0000
SYMG	61	-3748.2778	7895.2937	0.0000	1.0000
SYMIG	62	-3746.5555	7898.3859	0.0000	1.0000
GTRI	64	-3744.2678	7906.8838	0.0000	1.0000
K2PI	57	-3772.0363	7916.6640	0.0000	1.0000
SYMI	61	-3771.2061	7941.1504	0.0000	1.0000
F81G	59	-3819.9858	8025.6364	0.0000	1.0000
F81IG	60	-3817.5219	8027.2452	0.0000	1.0000
JC69G	56	-3836.0243	8038.1033	0.0000	1.0000
JC69IG	57	-3833.3802	8039.3518	0.0000	1.0000
F81I	59	-3841.0027	8067.6703	0.0000	1.0000
JC69I	56	-3857.2245	8080.5037	0.0000	1.0000
HKY	59	-3961.0784	8307.8215	0.0000	1.0000
GTR	63	-3952.2724	8316.3564	0.0000	1.0000
K2P	56	-3985.3348	8336.7244	0.0000	1.0000
SYM	60	-3983.8659	8359.9332	0.0000	1.0000
F81	58	-4045.7084	8470.5449	0.0000	1.0000
JC69	55	-4062.1106	8483.7393	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 183-phylip/chr1_5460.phylip = 1 - 690;
 Partition Dummy = 1:183-phylip/chr1_5460.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 183-phylip/chr1_5460.phylip = 1 - 690;
 Partition Dummy = 1:183-phylip/chr1_5460.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 (HKYG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr1_5460.phylip = 1 - 690;
 Partition Dummy = 1:183-phylip/chr1_5460.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

