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

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

Minimum AIC  model: GTR
Minimum AICc model: K2P
Minimum BIC  model: K2P


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
K2P	56	-375.9778	925.9362	0.6256	0.6256
K2PI	57	-375.4416	929.7066	0.0950	0.7205
K2PG	57	-375.5112	929.8459	0.0886	0.8091
JC69	55	-380.6635	930.5578	0.0620	0.8711
HKY	59	-370.9550	930.7100	0.0575	0.9286
SYM	60	-369.5377	933.0148	0.0182	0.9468
JC69I	56	-380.1372	934.2551	0.0098	0.9566
JC69G	56	-380.2032	934.3870	0.0091	0.9657
K2PIG	58	-375.4434	934.6491	0.0080	0.9737
HKYI	60	-370.4244	934.7882	0.0075	0.9812
HKYG	60	-370.4881	934.9157	0.0070	0.9883
F81	58	-376.0580	935.8785	0.0043	0.9926
SYMI	61	-369.2004	937.5845	0.0018	0.9944
SYMG	61	-369.2914	937.7665	0.0017	0.9961
GTR	63	-364.4173	938.8345	0.0010	0.9971
JC69IG	57	-380.2023	939.2282	0.0008	0.9979
F81I	59	-375.4962	939.7924	0.0006	0.9985
F81G	59	-375.5683	939.9367	0.0006	0.9991
HKYIG	61	-370.4260	940.0356	0.0005	0.9997
SYMIG	62	-369.2033	942.9426	0.0001	0.9998
GTRI	64	-364.0659	943.7106	0.0001	0.9999
GTRG	64	-364.1630	943.9050	0.0001	0.9999
F81IG	60	-375.4977	944.9349	0.0000	1.0000
GTRIG	65	-364.0665	949.4096	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTR	63	-364.4173	854.8345	0.3741	0.3741
GTRI	64	-364.0659	856.1317	0.1956	0.5697
GTRG	64	-364.1630	856.3261	0.1775	0.7472
GTRIG	65	-364.0665	858.1330	0.0719	0.8191
SYM	60	-369.5377	859.0754	0.0449	0.8640
HKY	59	-370.9550	859.9100	0.0296	0.8936
SYMI	61	-369.2004	860.4008	0.0231	0.9167
SYMG	61	-369.2914	860.5828	0.0211	0.9378
HKYI	60	-370.4244	860.8488	0.0185	0.9563
HKYG	60	-370.4881	860.9763	0.0174	0.9737
SYMIG	62	-369.2033	862.4065	0.0085	0.9822
HKYIG	61	-370.4260	862.8519	0.0068	0.9890
K2P	56	-375.9778	863.9556	0.0039	0.9929
K2PI	57	-375.4416	864.8831	0.0025	0.9954
K2PG	57	-375.5112	865.0224	0.0023	0.9976
K2PIG	58	-375.4434	866.8867	0.0009	0.9985
F81	58	-376.0580	868.1161	0.0005	0.9990
F81I	59	-375.4962	868.9924	0.0003	0.9994
F81G	59	-375.5683	869.1367	0.0003	0.9996
F81IG	60	-375.4977	870.9955	0.0001	0.9998
JC69	55	-380.6635	871.3271	0.0001	0.9999
JC69I	56	-380.1372	872.2745	0.0001	0.9999
JC69G	56	-380.2032	872.4065	0.0001	1.0000
JC69IG	57	-380.2023	874.4047	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	56	-375.9778	1036.1653	0.6431	0.6431
K2PI	57	-375.4416	1040.1680	0.0869	0.7300
K2PG	57	-375.5112	1040.3073	0.0811	0.8110
JC69	55	-380.6635	1040.4616	0.0750	0.8861
HKY	59	-370.9550	1041.3452	0.0482	0.9343
SYM	60	-369.5377	1043.5858	0.0157	0.9501
JC69I	56	-380.1372	1044.4842	0.0100	0.9601
JC69G	56	-380.2032	1044.6162	0.0094	0.9695
K2PIG	58	-375.4434	1045.2468	0.0069	0.9764
HKYI	60	-370.4244	1045.3592	0.0065	0.9829
HKYG	60	-370.4881	1045.4867	0.0061	0.9889
F81	58	-376.0580	1046.4762	0.0037	0.9926
SYMI	61	-369.2004	1047.9864	0.0017	0.9944
SYMG	61	-369.2914	1048.1684	0.0016	0.9960
GTR	63	-364.4173	1048.5705	0.0013	0.9973
JC69IG	57	-380.2023	1049.6896	0.0007	0.9980
F81I	59	-375.4962	1050.4276	0.0005	0.9985
HKYIG	61	-370.4260	1050.4375	0.0005	0.9991
F81G	59	-375.5683	1050.5719	0.0005	0.9995
GTRI	64	-364.0659	1052.9428	0.0001	0.9997
SYMIG	62	-369.2033	1053.0673	0.0001	0.9998
GTRG	64	-364.1630	1053.1372	0.0001	0.9999
F81IG	60	-375.4977	1055.5059	0.0000	1.0000
GTRIG	65	-364.0665	1058.0193	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTR)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3284.phylip = 1 - 160;
 Partition Dummy = 1:183-phylip/chr6_3284.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=equal;
 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);
END;


[Mrbayes block for the best AICc model (K2P)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3284.phylip = 1 - 160;
 Partition Dummy = 1:183-phylip/chr6_3284.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


[Mrbayes block for the best BIC model (K2P)]

BEGIN MRBAYES;
 Charset 183-phylip/chr6_3284.phylip = 1 - 160;
 Partition Dummy = 1:183-phylip/chr6_3284.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

