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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1120.9160	2378.2077	0.5887	0.5887
HKYI	60	-1121.7969	2379.9697	0.2439	0.8326
HKYIG	61	-1121.1523	2381.2643	0.1277	0.9603
GTRG	64	-1118.7307	2384.2424	0.0288	0.9891
GTRI	64	-1120.1374	2387.0558	0.0071	0.9962
GTRIG	65	-1119.4824	2388.3765	0.0036	0.9998
F81I	59	-1131.1043	2396.0122	0.0001	0.9999
F81G	59	-1131.3270	2396.4575	0.0001	1.0000
F81IG	60	-1130.4091	2397.1941	0.0000	1.0000
HKY	59	-1134.0056	2401.8147	0.0000	1.0000
GTR	63	-1132.3434	2408.8490	0.0000	1.0000
F81	58	-1144.2251	2419.6929	0.0000	1.0000
K2PG	57	-1148.6671	2426.0276	0.0000	1.0000
K2PI	57	-1149.4472	2427.5878	0.0000	1.0000
K2PIG	58	-1148.7197	2428.6822	0.0000	1.0000
SYMG	61	-1148.3085	2435.5766	0.0000	1.0000
SYMI	61	-1149.2961	2437.5519	0.0000	1.0000
SYMIG	62	-1148.5557	2438.6665	0.0000	1.0000
JC69I	56	-1156.9524	2440.0599	0.0000	1.0000
JC69G	56	-1157.1617	2440.4787	0.0000	1.0000
JC69IG	57	-1156.2156	2441.1245	0.0000	1.0000
K2P	56	-1161.4126	2448.9804	0.0000	1.0000
SYM	60	-1161.0959	2458.5677	0.0000	1.0000
JC69	55	-1169.9342	2463.4968	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1120.9160	2361.8319	0.5165	0.5165
HKYI	60	-1121.7969	2363.5938	0.2140	0.7305
HKYIG	61	-1121.1523	2364.3047	0.1500	0.8805
GTRG	64	-1118.7307	2365.4613	0.0841	0.9647
GTRI	64	-1120.1374	2368.2748	0.0206	0.9853
GTRIG	65	-1119.4824	2368.9648	0.0146	0.9999
F81I	59	-1131.1043	2380.2086	0.0001	0.9999
F81G	59	-1131.3270	2380.6540	0.0000	1.0000
F81IG	60	-1130.4091	2380.8183	0.0000	1.0000
HKY	59	-1134.0056	2386.0112	0.0000	1.0000
GTR	63	-1132.3434	2390.6868	0.0000	1.0000
F81	58	-1144.2251	2404.4501	0.0000	1.0000
K2PG	57	-1148.6671	2411.3342	0.0000	1.0000
K2PI	57	-1149.4472	2412.8944	0.0000	1.0000
K2PIG	58	-1148.7197	2413.4394	0.0000	1.0000
SYMG	61	-1148.3085	2418.6170	0.0000	1.0000
SYMI	61	-1149.2961	2420.5923	0.0000	1.0000
SYMIG	62	-1148.5557	2421.1114	0.0000	1.0000
JC69I	56	-1156.9524	2425.9047	0.0000	1.0000
JC69G	56	-1157.1617	2426.3234	0.0000	1.0000
JC69IG	57	-1156.2156	2426.4312	0.0000	1.0000
K2P	56	-1161.4126	2434.8252	0.0000	1.0000
SYM	60	-1161.0959	2442.1919	0.0000	1.0000
JC69	55	-1169.9342	2449.8685	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1120.9160	2615.6608	0.6891	0.6891
HKYI	60	-1121.7969	2617.4227	0.2856	0.9747
HKYIG	61	-1121.1523	2622.3640	0.0241	0.9988
F81I	59	-1131.1043	2629.8070	0.0006	0.9994
F81G	59	-1131.3270	2630.2524	0.0005	0.9999
F81IG	60	-1130.4091	2634.6472	0.0001	0.9999
HKY	59	-1134.0056	2635.6096	0.0000	1.0000
GTRG	64	-1118.7307	2636.2122	0.0000	1.0000
GTRI	64	-1120.1374	2639.0256	0.0000	1.0000
GTRIG	65	-1119.4824	2643.9461	0.0000	1.0000
F81	58	-1144.2251	2649.8180	0.0000	1.0000
K2PG	57	-1148.6671	2652.4717	0.0000	1.0000
K2PI	57	-1149.4472	2654.0319	0.0000	1.0000
GTR	63	-1132.3434	2657.2072	0.0000	1.0000
K2PIG	58	-1148.7197	2658.8074	0.0000	1.0000
JC69I	56	-1156.9524	2662.8117	0.0000	1.0000
JC69G	56	-1157.1617	2663.2304	0.0000	1.0000
JC69IG	57	-1156.2156	2667.5686	0.0000	1.0000
K2P	56	-1161.4126	2671.7322	0.0000	1.0000
SYMG	61	-1148.3085	2676.6763	0.0000	1.0000
SYMI	61	-1149.2961	2678.6517	0.0000	1.0000
JC69	55	-1169.9342	2682.5449	0.0000	1.0000
SYMIG	62	-1148.5557	2683.4013	0.0000	1.0000
SYM	60	-1161.0959	2696.0208	0.0000	1.0000



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

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

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


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

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


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

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

