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

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

Minimum AIC  model: K2PIG
Minimum AICc model: K2PIG
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
K2PIG	58	-1038.5578	2221.2802	0.9294	0.9294
K2PI	57	-1043.6531	2228.4046	0.0264	0.9557
HKYIG	61	-1037.5040	2228.5246	0.0248	0.9806
SYMIG	62	-1036.4308	2229.5479	0.0149	0.9955
K2PG	57	-1046.1369	2233.3721	0.0022	0.9977
HKYI	60	-1042.3172	2235.0079	0.0010	0.9986
GTRIG	65	-1034.6055	2235.5670	0.0007	0.9994
SYMI	61	-1041.4722	2236.4611	0.0005	0.9998
HKYG	60	-1044.6945	2239.7624	0.0001	0.9999
SYMG	61	-1043.7403	2240.9973	0.0000	1.0000
GTRI	64	-1039.6328	2242.3710	0.0000	1.0000
GTRG	64	-1041.7323	2246.5701	0.0000	1.0000
JC69IG	57	-1064.1427	2269.3838	0.0000	1.0000
F81IG	60	-1062.4748	2275.3230	0.0000	1.0000
JC69I	56	-1068.9439	2275.9449	0.0000	1.0000
JC69G	56	-1071.2392	2280.5356	0.0000	1.0000
F81I	59	-1066.8162	2280.8886	0.0000	1.0000
F81G	59	-1069.0136	2285.2835	0.0000	1.0000
K2P	56	-1133.4076	2404.8723	0.0000	1.0000
HKY	59	-1130.2346	2407.7254	0.0000	1.0000
GTR	63	-1125.1723	2410.2270	0.0000	1.0000
SYM	60	-1130.2696	2410.9127	0.0000	1.0000
JC69	55	-1156.4734	2447.9875	0.0000	1.0000
F81	58	-1153.0304	2450.2254	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	58	-1038.5578	2193.1156	0.7303	0.7303
SYMIG	62	-1036.4308	2196.8617	0.1122	0.8425
HKYIG	61	-1037.5040	2197.0079	0.1043	0.9468
GTRIG	65	-1034.6055	2199.2110	0.0347	0.9814
K2PI	57	-1043.6531	2201.3062	0.0122	0.9936
HKYI	60	-1042.3172	2204.6345	0.0023	0.9959
SYMI	61	-1041.4722	2204.9444	0.0020	0.9979
K2PG	57	-1046.1369	2206.2737	0.0010	0.9989
GTRI	64	-1039.6328	2207.2656	0.0006	0.9995
HKYG	60	-1044.6945	2209.3889	0.0002	0.9997
SYMG	61	-1043.7403	2209.4806	0.0002	0.9999
GTRG	64	-1041.7323	2211.4647	0.0001	1.0000
JC69IG	57	-1064.1427	2242.2855	0.0000	1.0000
F81IG	60	-1062.4748	2244.9496	0.0000	1.0000
JC69I	56	-1068.9439	2249.8878	0.0000	1.0000
F81I	59	-1066.8162	2251.6324	0.0000	1.0000
JC69G	56	-1071.2392	2254.4785	0.0000	1.0000
F81G	59	-1069.0136	2256.0273	0.0000	1.0000
GTR	63	-1125.1723	2376.3446	0.0000	1.0000
HKY	59	-1130.2346	2378.4692	0.0000	1.0000
K2P	56	-1133.4076	2378.8152	0.0000	1.0000
SYM	60	-1130.2696	2380.5393	0.0000	1.0000
F81	58	-1153.0304	2422.0608	0.0000	1.0000
JC69	55	-1156.4734	2422.9468	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-1038.5578	2408.3204	0.8960	0.8960
K2PI	57	-1043.6531	2412.8006	0.0954	0.9914
K2PG	57	-1046.1369	2417.7680	0.0080	0.9993
HKYIG	61	-1037.5040	2423.3440	0.0005	0.9998
SYMIG	62	-1036.4308	2426.9082	0.0001	0.9999
HKYI	60	-1042.3172	2427.2601	0.0001	1.0000
SYMI	61	-1041.4722	2431.2805	0.0000	1.0000
HKYG	60	-1044.6945	2432.0146	0.0000	1.0000
SYMG	61	-1043.7403	2435.8167	0.0000	1.0000
GTRIG	65	-1034.6055	2440.3888	0.0000	1.0000
GTRI	64	-1039.6328	2444.7329	0.0000	1.0000
GTRG	64	-1041.7323	2448.9320	0.0000	1.0000
JC69IG	57	-1064.1427	2453.7798	0.0000	1.0000
JC69I	56	-1068.9439	2457.6717	0.0000	1.0000
JC69G	56	-1071.2392	2462.2624	0.0000	1.0000
F81IG	60	-1062.4748	2467.5752	0.0000	1.0000
F81I	59	-1066.8162	2470.5476	0.0000	1.0000
F81G	59	-1069.0136	2474.9425	0.0000	1.0000
K2P	56	-1133.4076	2586.5991	0.0000	1.0000
HKY	59	-1130.2346	2597.3844	0.0000	1.0000
SYM	60	-1130.2696	2603.1649	0.0000	1.0000
GTR	63	-1125.1723	2610.1015	0.0000	1.0000
JC69	55	-1156.4734	2627.0203	0.0000	1.0000
F81	58	-1153.0304	2637.2656	0.0000	1.0000



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

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

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


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

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


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

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

