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

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

Minimum AIC  model: K2PG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1002.9374	2145.8042	0.4742	0.4742
K2PI	57	-1003.2315	2146.3924	0.3533	0.8275
K2PIG	58	-1002.8705	2148.6859	0.1122	0.9398
SYMG	61	-999.9594	2152.0544	0.0208	0.9606
SYMI	61	-1000.4532	2153.0418	0.0127	0.9733
HKYG	60	-1002.2344	2153.5163	0.0100	0.9833
HKYI	60	-1002.4402	2153.9280	0.0082	0.9915
SYMIG	62	-999.8103	2154.8686	0.0051	0.9966
HKYIG	61	-1002.1584	2156.4523	0.0023	0.9989
GTRG	64	-998.7970	2159.1423	0.0006	0.9995
GTRI	64	-999.3334	2160.2152	0.0004	0.9999
GTRIG	65	-998.6308	2161.9984	0.0001	1.0000
JC69G	56	-1018.7281	2174.3937	0.0000	1.0000
JC69I	56	-1019.0530	2175.0435	0.0000	1.0000
JC69IG	57	-1018.6524	2177.2342	0.0000	1.0000
K2P	56	-1021.7318	2180.4010	0.0000	1.0000
F81G	59	-1017.7287	2181.4417	0.0000	1.0000
F81I	59	-1017.9278	2181.8398	0.0000	1.0000
SYM	60	-1017.5013	2184.0502	0.0000	1.0000
F81IG	60	-1017.6568	2184.3612	0.0000	1.0000
HKY	59	-1020.8053	2187.5947	0.0000	1.0000
GTR	63	-1016.1099	2190.6053	0.0000	1.0000
JC69	55	-1037.1354	2208.2396	0.0000	1.0000
F81	58	-1036.1761	2215.2970	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PG	57	-1002.9374	2119.8747	0.3116	0.3116
K2PI	57	-1003.2315	2120.4630	0.2322	0.5437
K2PIG	58	-1002.8705	2121.7410	0.1225	0.6663
SYMG	61	-999.9594	2121.9189	0.1121	0.7784
SYMI	61	-1000.4532	2122.9064	0.0684	0.8468
SYMIG	62	-999.8103	2123.6206	0.0479	0.8947
HKYG	60	-1002.2344	2124.4687	0.0313	0.9260
HKYI	60	-1002.4402	2124.8804	0.0255	0.9515
GTRG	64	-998.7970	2125.5939	0.0178	0.9694
HKYIG	61	-1002.1584	2126.3168	0.0124	0.9818
GTRI	64	-999.3334	2126.6669	0.0104	0.9922
GTRIG	65	-998.6308	2127.2616	0.0078	1.0000
JC69G	56	-1018.7281	2149.4562	0.0000	1.0000
JC69I	56	-1019.0530	2150.1060	0.0000	1.0000
JC69IG	57	-1018.6524	2151.3047	0.0000	1.0000
F81G	59	-1017.7287	2153.4575	0.0000	1.0000
F81I	59	-1017.9278	2153.8556	0.0000	1.0000
SYM	60	-1017.5013	2155.0026	0.0000	1.0000
F81IG	60	-1017.6568	2155.3136	0.0000	1.0000
K2P	56	-1021.7318	2155.4635	0.0000	1.0000
GTR	63	-1016.1099	2158.2198	0.0000	1.0000
HKY	59	-1020.8053	2159.6106	0.0000	1.0000
JC69	55	-1037.1354	2184.2707	0.0000	1.0000
F81	58	-1036.1761	2188.3521	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1002.9374	2333.4083	0.5535	0.5535
K2PI	57	-1003.2315	2333.9965	0.4125	0.9660
K2PIG	58	-1002.8705	2339.0208	0.0335	0.9994
HKYG	60	-1002.2344	2349.2409	0.0002	0.9996
HKYI	60	-1002.4402	2349.6526	0.0002	0.9998
SYMG	61	-999.9594	2350.4373	0.0001	0.9999
SYMI	61	-1000.4532	2351.4248	0.0001	1.0000
HKYIG	61	-1002.1584	2354.8352	0.0000	1.0000
SYMIG	62	-999.8103	2355.8852	0.0000	1.0000
JC69G	56	-1018.7281	2359.2436	0.0000	1.0000
JC69I	56	-1019.0530	2359.8934	0.0000	1.0000
JC69IG	57	-1018.6524	2364.8383	0.0000	1.0000
K2P	56	-1021.7318	2365.2509	0.0000	1.0000
GTRG	64	-998.7970	2365.3509	0.0000	1.0000
GTRI	64	-999.3334	2366.4239	0.0000	1.0000
GTRIG	65	-998.6308	2370.7648	0.0000	1.0000
F81G	59	-1017.7287	2374.4835	0.0000	1.0000
F81I	59	-1017.9278	2374.8816	0.0000	1.0000
SYM	60	-1017.5013	2379.7748	0.0000	1.0000
F81IG	60	-1017.6568	2380.0858	0.0000	1.0000
HKY	59	-1020.8053	2380.6365	0.0000	1.0000
JC69	55	-1037.1354	2390.3119	0.0000	1.0000
GTR	63	-1016.1099	2394.2306	0.0000	1.0000
F81	58	-1036.1761	2405.6319	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_7559.phylip = 1 - 313;
 Partition Dummy = 1:183-phylip/chr4_7559.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_7559.phylip = 1 - 313;
 Partition Dummy = 1:183-phylip/chr4_7559.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr4_7559.phylip = 1 - 313;
 Partition Dummy = 1:183-phylip/chr4_7559.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

