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

Input data from file "183-phylip/chr5_11292.phylip" (29 taxa, 815 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	-3041.2136	6223.8825	0.5705	0.5705
GTRG	64	-3042.7825	6224.6583	0.3871	0.9576
HKYIG	61	-3048.9647	6229.9745	0.0271	0.9847
HKYG	60	-3050.7184	6231.1451	0.0151	0.9998
GTRI	64	-3050.3949	6239.8832	0.0002	1.0000
HKYI	60	-3058.5527	6246.8137	0.0000	1.0000
K2PIG	58	-3075.2077	6275.4683	0.0000	1.0000
K2PG	57	-3077.1932	6277.1210	0.0000	1.0000
SYMIG	62	-3071.3729	6277.1342	0.0000	1.0000
SYMG	61	-3073.3041	6278.6533	0.0000	1.0000
K2PI	57	-3085.3192	6293.3729	0.0000	1.0000
SYMI	61	-3081.5350	6295.1151	0.0000	1.0000
F81IG	60	-3114.4367	6358.5816	0.0000	1.0000
F81G	59	-3115.8912	6359.1599	0.0000	1.0000
GTR	63	-3117.8981	6372.5338	0.0000	1.0000
F81I	59	-3123.0422	6373.4619	0.0000	1.0000
HKY	59	-3124.8650	6377.1075	0.0000	1.0000
JC69IG	57	-3136.1055	6394.9454	0.0000	1.0000
JC69G	56	-3137.9047	6396.2315	0.0000	1.0000
JC69I	56	-3145.5263	6411.4747	0.0000	1.0000
SYM	60	-3147.5888	6424.8858	0.0000	1.0000
K2P	56	-3152.3584	6425.1390	0.0000	1.0000
F81	58	-3188.0467	6501.1462	0.0000	1.0000
JC69	55	-3210.5074	6539.1308	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3041.2136	6212.4272	0.6246	0.6246
GTRG	64	-3042.7825	6213.5649	0.3536	0.9782
HKYIG	61	-3048.9647	6219.9293	0.0147	0.9929
HKYG	60	-3050.7184	6221.4368	0.0069	0.9998
GTRI	64	-3050.3949	6228.7898	0.0002	1.0000
HKYI	60	-3058.5527	6237.1054	0.0000	1.0000
K2PIG	58	-3075.2077	6266.4154	0.0000	1.0000
SYMIG	62	-3071.3729	6266.7459	0.0000	1.0000
K2PG	57	-3077.1932	6268.3865	0.0000	1.0000
SYMG	61	-3073.3041	6268.6082	0.0000	1.0000
K2PI	57	-3085.3192	6284.6384	0.0000	1.0000
SYMI	61	-3081.5350	6285.0699	0.0000	1.0000
F81IG	60	-3114.4367	6348.8734	0.0000	1.0000
F81G	59	-3115.8912	6349.7824	0.0000	1.0000
GTR	63	-3117.8981	6361.7961	0.0000	1.0000
F81I	59	-3123.0422	6364.0844	0.0000	1.0000
HKY	59	-3124.8650	6367.7301	0.0000	1.0000
JC69IG	57	-3136.1055	6386.2109	0.0000	1.0000
JC69G	56	-3137.9047	6387.8094	0.0000	1.0000
JC69I	56	-3145.5263	6403.0526	0.0000	1.0000
SYM	60	-3147.5888	6415.1776	0.0000	1.0000
K2P	56	-3152.3584	6416.7168	0.0000	1.0000
F81	58	-3188.0467	6492.0933	0.0000	1.0000
JC69	55	-3210.5074	6531.0149	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3050.7184	6503.6281	0.8281	0.8281
HKYIG	61	-3048.9647	6506.8238	0.1675	0.9956
GTRG	64	-3042.7825	6514.5690	0.0035	0.9991
GTRIG	65	-3041.2136	6518.1344	0.0006	0.9997
HKYI	60	-3058.5527	6519.2967	0.0003	1.0000
GTRI	64	-3050.3949	6529.7939	0.0000	1.0000
K2PG	57	-3077.1932	6536.4682	0.0000	1.0000
K2PIG	58	-3075.2077	6539.2003	0.0000	1.0000
K2PI	57	-3085.3192	6552.7201	0.0000	1.0000
SYMG	61	-3073.3041	6555.5026	0.0000	1.0000
SYMIG	62	-3071.3729	6558.3436	0.0000	1.0000
SYMI	61	-3081.5350	6571.9644	0.0000	1.0000
F81G	59	-3115.8912	6627.2705	0.0000	1.0000
F81IG	60	-3114.4367	6631.0647	0.0000	1.0000
F81I	59	-3123.0422	6641.5725	0.0000	1.0000
HKY	59	-3124.8650	6645.2182	0.0000	1.0000
JC69G	56	-3137.9047	6651.1879	0.0000	1.0000
JC69IG	57	-3136.1055	6654.2926	0.0000	1.0000
GTR	63	-3117.8981	6658.0970	0.0000	1.0000
JC69I	56	-3145.5263	6666.4311	0.0000	1.0000
K2P	56	-3152.3584	6680.0954	0.0000	1.0000
SYM	60	-3147.5888	6697.3689	0.0000	1.0000
F81	58	-3188.0467	6764.8782	0.0000	1.0000
JC69	55	-3210.5074	6789.6902	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/chr5_11292.phylip = 1 - 815;
 Partition Dummy = 1:183-phylip/chr5_11292.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/chr5_11292.phylip = 1 - 815;
 Partition Dummy = 1:183-phylip/chr5_11292.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/chr5_11292.phylip = 1 - 815;
 Partition Dummy = 1:183-phylip/chr5_11292.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

