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

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

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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	-3176.2345	6495.6514	0.7305	0.7305
GTRIG	65	-3176.2243	6498.1342	0.2111	0.9416
SYMG	61	-3182.7562	6501.2401	0.0447	0.9862
K2PG	57	-3189.4140	6504.7415	0.0078	0.9940
HKYG	60	-3186.2868	6505.8344	0.0045	0.9985
HKYIG	61	-3186.2864	6508.3005	0.0013	0.9998
SYMIG	62	-3187.1570	6512.5177	0.0002	1.0000
K2PIG	58	-3193.9800	6516.3138	0.0000	1.0000
GTRI	64	-3187.1401	6517.4627	0.0000	1.0000
HKYI	60	-3194.6562	6522.5732	0.0000	1.0000
SYMI	61	-3195.4981	6526.7239	0.0000	1.0000
K2PI	57	-3202.3043	6530.5220	0.0000	1.0000
SYM	60	-3224.4661	6582.1930	0.0000	1.0000
GTR	63	-3221.1614	6583.0113	0.0000	1.0000
K2P	56	-3231.6714	6586.8247	0.0000	1.0000
HKY	59	-3228.6398	6588.0825	0.0000	1.0000
JC69G	56	-3263.4501	6650.3821	0.0000	1.0000
F81G	59	-3259.8156	6650.4341	0.0000	1.0000
F81IG	60	-3264.5050	6662.2708	0.0000	1.0000
JC69IG	57	-3268.2076	6662.3286	0.0000	1.0000
JC69I	56	-3271.3014	6666.0848	0.0000	1.0000
F81I	59	-3267.9608	6666.7245	0.0000	1.0000
JC69	55	-3304.5342	6730.1276	0.0000	1.0000
F81	58	-3301.5342	6731.4222	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3176.2345	6480.4690	0.7109	0.7109
GTRIG	65	-3176.2243	6482.4487	0.2642	0.9751
SYMG	61	-3182.7562	6487.5123	0.0210	0.9961
HKYG	60	-3186.2868	6492.5736	0.0017	0.9978
K2PG	57	-3189.4140	6492.8280	0.0015	0.9993
HKYIG	61	-3186.2864	6494.5727	0.0006	0.9999
SYMIG	62	-3187.1570	6498.3141	0.0001	1.0000
GTRI	64	-3187.1401	6502.2802	0.0000	1.0000
K2PIG	58	-3193.9800	6503.9600	0.0000	1.0000
HKYI	60	-3194.6562	6509.3123	0.0000	1.0000
SYMI	61	-3195.4981	6512.9961	0.0000	1.0000
K2PI	57	-3202.3043	6518.6085	0.0000	1.0000
GTR	63	-3221.1614	6568.3227	0.0000	1.0000
SYM	60	-3224.4661	6568.9321	0.0000	1.0000
HKY	59	-3228.6398	6575.2796	0.0000	1.0000
K2P	56	-3231.6714	6575.3427	0.0000	1.0000
F81G	59	-3259.8156	6637.6312	0.0000	1.0000
JC69G	56	-3263.4501	6638.9001	0.0000	1.0000
F81IG	60	-3264.5050	6649.0099	0.0000	1.0000
JC69IG	57	-3268.2076	6650.4151	0.0000	1.0000
F81I	59	-3267.9608	6653.9216	0.0000	1.0000
JC69I	56	-3271.3014	6654.6028	0.0000	1.0000
JC69	55	-3304.5342	6719.0683	0.0000	1.0000
F81	58	-3301.5342	6719.0684	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-3189.4140	6744.6748	0.9959	0.9959
SYMG	61	-3182.7562	6757.0326	0.0021	0.9979
HKYG	60	-3186.2868	6757.6755	0.0015	0.9994
K2PIG	58	-3193.9800	6760.2252	0.0004	0.9998
GTRG	64	-3176.2345	6763.2443	0.0001	0.9999
HKYIG	61	-3186.2864	6764.0930	0.0001	1.0000
GTRIG	65	-3176.2243	6769.6424	0.0000	1.0000
K2PI	57	-3202.3043	6770.4553	0.0000	1.0000
SYMIG	62	-3187.1570	6772.2527	0.0000	1.0000
HKYI	60	-3194.6562	6774.4142	0.0000	1.0000
SYMI	61	-3195.4981	6782.5164	0.0000	1.0000
GTRI	64	-3187.1401	6785.0556	0.0000	1.0000
K2P	56	-3231.6714	6822.7711	0.0000	1.0000
SYM	60	-3224.4661	6834.0340	0.0000	1.0000
HKY	59	-3228.6398	6835.9632	0.0000	1.0000
GTR	63	-3221.1614	6846.6797	0.0000	1.0000
JC69G	56	-3263.4501	6886.3286	0.0000	1.0000
F81G	59	-3259.8156	6898.3148	0.0000	1.0000
JC69I	56	-3271.3014	6902.0313	0.0000	1.0000
JC69IG	57	-3268.2076	6902.2619	0.0000	1.0000
F81IG	60	-3264.5050	6914.1118	0.0000	1.0000
F81I	59	-3267.9608	6914.6052	0.0000	1.0000
JC69	55	-3304.5342	6962.0784	0.0000	1.0000
F81	58	-3301.5342	6975.3335	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 183-phylip/chr7_11510.phylip = 1 - 613;
 Partition Dummy = 1:183-phylip/chr7_11510.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 183-phylip/chr7_11510.phylip = 1 - 613;
 Partition Dummy = 1:183-phylip/chr7_11510.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 (K2PG)]

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


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

