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

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

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


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	-2788.1649	5715.4379	0.5902	0.5902
GTRIG	65	-2787.3619	5716.1945	0.4043	0.9945
GTRI	64	-2792.8530	5724.8141	0.0054	1.0000
SYMG	61	-2802.2588	5736.5762	0.0000	1.0000
SYMIG	62	-2801.3120	5737.0261	0.0000	1.0000
SYMI	61	-2806.7831	5745.6248	0.0000	1.0000
HKYG	60	-2808.5613	5746.8437	0.0000	1.0000
HKYIG	61	-2807.7620	5747.5825	0.0000	1.0000
HKYI	60	-2812.7504	5755.2220	0.0000	1.0000
K2PIG	58	-2820.3360	5765.7370	0.0000	1.0000
K2PG	57	-2822.7006	5768.1473	0.0000	1.0000
K2PI	57	-2825.0239	5772.7939	0.0000	1.0000
F81G	59	-2867.2084	5861.8067	0.0000	1.0000
F81IG	60	-2867.1018	5863.9248	0.0000	1.0000
F81I	59	-2870.4915	5868.3729	0.0000	1.0000
JC69G	56	-2883.3525	5887.1383	0.0000	1.0000
JC69IG	57	-2882.8536	5888.4532	0.0000	1.0000
GTR	63	-2876.0799	5888.9117	0.0000	1.0000
JC69I	56	-2886.1188	5892.6710	0.0000	1.0000
SYM	60	-2889.6447	5909.0106	0.0000	1.0000
HKY	59	-2895.1298	5917.6496	0.0000	1.0000
K2P	56	-2907.9537	5936.3407	0.0000	1.0000
F81	58	-2950.6985	6026.4619	0.0000	1.0000
JC69	55	-2966.5106	6051.1478	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-2788.1649	5704.3297	0.5463	0.5463
GTRIG	65	-2787.3619	5704.7239	0.4486	0.9950
GTRI	64	-2792.8530	5713.7060	0.0050	1.0000
SYMG	61	-2802.2588	5726.5177	0.0000	1.0000
SYMIG	62	-2801.3120	5726.6240	0.0000	1.0000
SYMI	61	-2806.7831	5735.5663	0.0000	1.0000
HKYG	60	-2808.5613	5737.1225	0.0000	1.0000
HKYIG	61	-2807.7620	5737.5240	0.0000	1.0000
HKYI	60	-2812.7504	5745.5009	0.0000	1.0000
K2PIG	58	-2820.3360	5756.6721	0.0000	1.0000
K2PG	57	-2822.7006	5759.4012	0.0000	1.0000
K2PI	57	-2825.0239	5764.0479	0.0000	1.0000
F81G	59	-2867.2084	5852.4168	0.0000	1.0000
F81IG	60	-2867.1018	5854.2036	0.0000	1.0000
F81I	59	-2870.4915	5858.9830	0.0000	1.0000
GTR	63	-2876.0799	5878.1597	0.0000	1.0000
JC69G	56	-2883.3525	5878.7050	0.0000	1.0000
JC69IG	57	-2882.8536	5879.7072	0.0000	1.0000
JC69I	56	-2886.1188	5884.2377	0.0000	1.0000
SYM	60	-2889.6447	5899.2895	0.0000	1.0000
HKY	59	-2895.1298	5908.2597	0.0000	1.0000
K2P	56	-2907.9537	5927.9074	0.0000	1.0000
F81	58	-2950.6985	6017.3970	0.0000	1.0000
JC69	55	-2966.5106	6043.0211	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-2788.1649	6005.2552	0.9027	0.9027
GTRIG	65	-2787.3619	6010.3513	0.0706	0.9733
SYMG	61	-2802.2588	6013.3372	0.0159	0.9892
GTRI	64	-2792.8530	6014.6314	0.0083	0.9975
SYMIG	62	-2801.3120	6018.1455	0.0014	0.9989
HKYG	60	-2808.5613	6019.2402	0.0008	0.9997
SYMI	61	-2806.7831	6022.3858	0.0002	0.9999
HKYIG	61	-2807.7620	6024.3436	0.0001	1.0000
K2PG	57	-2822.7006	6027.4130	0.0000	1.0000
HKYI	60	-2812.7504	6027.6185	0.0000	1.0000
K2PIG	58	-2820.3360	6029.3858	0.0000	1.0000
K2PI	57	-2825.0239	6032.0596	0.0000	1.0000
F81G	59	-2867.2084	6129.8324	0.0000	1.0000
F81IG	60	-2867.1018	6136.3213	0.0000	1.0000
F81I	59	-2870.4915	6136.3987	0.0000	1.0000
JC69G	56	-2883.3525	6142.0148	0.0000	1.0000
JC69I	56	-2886.1188	6147.5475	0.0000	1.0000
JC69IG	57	-2882.8536	6147.7189	0.0000	1.0000
GTR	63	-2876.0799	6174.3832	0.0000	1.0000
SYM	60	-2889.6447	6181.4071	0.0000	1.0000
HKY	59	-2895.1298	6185.6753	0.0000	1.0000
K2P	56	-2907.9537	6191.2172	0.0000	1.0000
F81	58	-2950.6985	6290.1107	0.0000	1.0000
JC69	55	-2966.5106	6301.6290	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/chr12_5837.phylip = 1 - 814;
 Partition Dummy = 1:183-phylip/chr12_5837.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/chr12_5837.phylip = 1 - 814;
 Partition Dummy = 1:183-phylip/chr12_5837.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 (GTRG)]

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


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

