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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: K2PIG


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	-3574.7139	7288.6536	0.8655	0.8655
HKYIG	61	-3581.2456	7292.5897	0.1209	0.9865
SYMIG	62	-3582.5326	7297.4383	0.0107	0.9972
K2PIG	58	-3589.2392	7301.7826	0.0012	0.9984
GTRG	64	-3582.4255	7301.7875	0.0012	0.9996
HKYI	60	-3589.0508	7305.9305	0.0002	0.9998
HKYG	60	-3589.0712	7305.9713	0.0002	0.9999
GTRI	64	-3585.4445	7307.8257	0.0001	1.0000
SYMG	61	-3589.9897	7310.0779	0.0000	1.0000
SYMI	61	-3591.2516	7312.6017	0.0000	1.0000
K2PG	57	-3596.8164	7314.6819	0.0000	1.0000
K2PI	57	-3597.3497	7315.7485	0.0000	1.0000
F81IG	60	-3677.5658	7482.9604	0.0000	1.0000
JC69IG	57	-3683.4929	7488.0348	0.0000	1.0000
F81G	59	-3682.2326	7490.0292	0.0000	1.0000
F81I	59	-3686.9862	7499.5366	0.0000	1.0000
JC69G	56	-3691.1517	7501.1021	0.0000	1.0000
JC69I	56	-3695.5121	7509.8229	0.0000	1.0000
GTR	63	-3772.8785	7680.4094	0.0000	1.0000
SYM	60	-3776.4220	7680.6729	0.0000	1.0000
HKY	59	-3782.5762	7690.7165	0.0000	1.0000
K2P	56	-3791.9724	7702.7434	0.0000	1.0000
F81	58	-3877.5369	7878.3779	0.0000	1.0000
JC69	55	-3884.1610	7884.8753	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3574.7139	7279.4278	0.9178	0.9178
HKYIG	61	-3581.2456	7284.4912	0.0730	0.9907
SYMIG	62	-3582.5326	7289.0653	0.0074	0.9982
GTRG	64	-3582.4255	7292.8509	0.0011	0.9993
K2PIG	58	-3589.2392	7294.4784	0.0005	0.9998
HKYI	60	-3589.0508	7298.1016	0.0001	0.9999
HKYG	60	-3589.0712	7298.1424	0.0001	0.9999
GTRI	64	-3585.4445	7298.8891	0.0001	1.0000
SYMG	61	-3589.9897	7301.9794	0.0000	1.0000
SYMI	61	-3591.2516	7304.5032	0.0000	1.0000
K2PG	57	-3596.8164	7307.6329	0.0000	1.0000
K2PI	57	-3597.3497	7308.6995	0.0000	1.0000
F81IG	60	-3677.5658	7475.1316	0.0000	1.0000
JC69IG	57	-3683.4929	7480.9857	0.0000	1.0000
F81G	59	-3682.2326	7482.4651	0.0000	1.0000
F81I	59	-3686.9862	7491.9725	0.0000	1.0000
JC69G	56	-3691.1517	7494.3033	0.0000	1.0000
JC69I	56	-3695.5121	7503.0242	0.0000	1.0000
GTR	63	-3772.8785	7671.7571	0.0000	1.0000
SYM	60	-3776.4220	7672.8440	0.0000	1.0000
HKY	59	-3782.5762	7683.1524	0.0000	1.0000
K2P	56	-3791.9724	7695.9447	0.0000	1.0000
F81	58	-3877.5369	7871.0737	0.0000	1.0000
JC69	55	-3884.1610	7878.3221	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-3589.2392	7578.8958	0.8903	0.8903
HKYIG	61	-3581.2456	7583.6198	0.0839	0.9742
K2PG	57	-3596.8164	7587.1465	0.0144	0.9886
K2PI	57	-3597.3497	7588.2131	0.0084	0.9971
HKYI	60	-3589.0508	7592.3264	0.0011	0.9981
HKYG	60	-3589.0712	7592.3673	0.0011	0.9992
SYMIG	62	-3582.5326	7593.0976	0.0007	0.9999
GTRIG	65	-3574.7139	7598.1714	0.0001	1.0000
SYMG	61	-3589.9897	7601.1080	0.0000	1.0000
SYMI	61	-3591.2516	7603.6318	0.0000	1.0000
GTRG	64	-3582.4255	7606.6907	0.0000	1.0000
GTRI	64	-3585.4445	7612.7289	0.0000	1.0000
JC69IG	57	-3683.4929	7760.4993	0.0000	1.0000
JC69G	56	-3691.1517	7768.9132	0.0000	1.0000
F81IG	60	-3677.5658	7769.3564	0.0000	1.0000
F81G	59	-3682.2326	7771.7862	0.0000	1.0000
JC69I	56	-3695.5121	7777.6340	0.0000	1.0000
F81I	59	-3686.9862	7781.2935	0.0000	1.0000
SYM	60	-3776.4220	7967.0688	0.0000	1.0000
K2P	56	-3791.9724	7970.5546	0.0000	1.0000
HKY	59	-3782.5762	7972.4735	0.0000	1.0000
GTR	63	-3772.8785	7980.6932	0.0000	1.0000
JC69	55	-3884.1610	8148.0282	0.0000	1.0000
F81	58	-3877.5369	8155.4911	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/chr2_24815.phylip = 1 - 996;
 Partition Dummy = 1:183-phylip/chr2_24815.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/chr2_24815.phylip = 1 - 996;
 Partition Dummy = 1:183-phylip/chr2_24815.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 (K2PIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_24815.phylip = 1 - 996;
 Partition Dummy = 1:183-phylip/chr2_24815.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

