
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
12:01 08/01/2011

Input data from file "261-aggressive-phylip/7924_hg18.phylip" (29 taxa, 269 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	-974.2559	2097.1023	0.8944	0.8944
HKYIG	61	-972.0335	2102.6082	0.0570	0.9514
K2PG	57	-979.0974	2103.5312	0.0359	0.9873
K2PI	57	-981.0242	2107.3849	0.0052	0.9925
SYMIG	62	-973.0494	2108.0211	0.0038	0.9964
HKYG	60	-976.8880	2108.9683	0.0024	0.9987
HKYI	60	-977.7954	2110.7830	0.0010	0.9997
GTRIG	65	-970.9815	2114.2289	0.0002	0.9999
SYMG	61	-978.1617	2114.8645	0.0001	1.0000
SYMI	61	-980.0604	2118.6619	0.0000	1.0000
GTRG	64	-976.3806	2121.5455	0.0000	1.0000
GTRI	64	-978.3877	2125.5598	0.0000	1.0000
JC69IG	57	-1005.5345	2156.4056	0.0000	1.0000
JC69G	56	-1009.5287	2161.1706	0.0000	1.0000
F81IG	60	-1003.6649	2162.5221	0.0000	1.0000
JC69I	56	-1011.1099	2164.3330	0.0000	1.0000
F81G	59	-1007.6622	2167.1999	0.0000	1.0000
F81I	59	-1009.2471	2170.3697	0.0000	1.0000
K2P	56	-1038.6714	2219.4560	0.0000	1.0000
HKY	59	-1036.6740	2225.2235	0.0000	1.0000
SYM	60	-1035.3553	2225.9030	0.0000	1.0000
GTR	63	-1033.4905	2232.3177	0.0000	1.0000
JC69	55	-1065.6334	2270.1870	0.0000	1.0000
F81	58	-1064.1672	2276.9248	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
K2PIG	58	-974.2559	2064.5118	0.6311	0.6311
HKYIG	61	-972.0335	2066.0671	0.2900	0.9211
SYMIG	62	-973.0494	2070.0988	0.0386	0.9597
GTRIG	65	-970.9815	2071.9629	0.0152	0.9749
K2PG	57	-979.0974	2072.1947	0.0135	0.9885
HKYG	60	-976.8880	2073.7760	0.0061	0.9946
HKYI	60	-977.7954	2075.5907	0.0025	0.9971
K2PI	57	-981.0242	2076.0484	0.0020	0.9991
SYMG	61	-978.1617	2078.3235	0.0006	0.9997
GTRG	64	-976.3806	2080.7612	0.0002	0.9999
SYMI	61	-980.0604	2082.1209	0.0001	1.0000
GTRI	64	-978.3877	2084.7755	0.0000	1.0000
JC69IG	57	-1005.5345	2125.0691	0.0000	1.0000
F81IG	60	-1003.6649	2127.3298	0.0000	1.0000
JC69G	56	-1009.5287	2131.0574	0.0000	1.0000
F81G	59	-1007.6622	2133.3243	0.0000	1.0000
JC69I	56	-1011.1099	2134.2198	0.0000	1.0000
F81I	59	-1009.2471	2136.4941	0.0000	1.0000
K2P	56	-1038.6714	2189.3428	0.0000	1.0000
SYM	60	-1035.3553	2190.7107	0.0000	1.0000
HKY	59	-1036.6740	2191.3479	0.0000	1.0000
GTR	63	-1033.4905	2192.9811	0.0000	1.0000
JC69	55	-1065.6334	2241.2668	0.0000	1.0000
F81	58	-1064.1672	2244.3344	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PIG	58	-974.2559	2273.0051	0.8689	0.8689
K2PG	57	-979.0974	2277.0933	0.1125	0.9814
K2PI	57	-981.0242	2280.9469	0.0164	0.9978
HKYIG	61	-972.0335	2285.3445	0.0018	0.9996
HKYG	60	-976.8880	2289.4586	0.0002	0.9999
HKYI	60	-977.7954	2291.2734	0.0001	1.0000
SYMIG	62	-973.0494	2292.9709	0.0000	1.0000
SYMG	61	-978.1617	2297.6009	0.0000	1.0000
SYMI	61	-980.0604	2301.3983	0.0000	1.0000
GTRIG	65	-970.9815	2305.6192	0.0000	1.0000
GTRG	64	-976.3806	2310.8227	0.0000	1.0000
GTRI	64	-978.3877	2314.8370	0.0000	1.0000
JC69IG	57	-1005.5345	2329.9676	0.0000	1.0000
JC69G	56	-1009.5287	2332.3612	0.0000	1.0000
JC69I	56	-1011.1099	2335.5237	0.0000	1.0000
F81IG	60	-1003.6649	2343.0125	0.0000	1.0000
F81G	59	-1007.6622	2345.4123	0.0000	1.0000
F81I	59	-1009.2471	2348.5821	0.0000	1.0000
K2P	56	-1038.6714	2390.6466	0.0000	1.0000
HKY	59	-1036.6740	2403.4359	0.0000	1.0000
SYM	60	-1035.3553	2406.3934	0.0000	1.0000
GTR	63	-1033.4905	2419.4479	0.0000	1.0000
JC69	55	-1065.6334	2438.9759	0.0000	1.0000
F81	58	-1064.1672	2452.8276	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 261-aggressive-phylip/7924_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/7924_hg18.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 261-aggressive-phylip/7924_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/7924_hg18.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 261-aggressive-phylip/7924_hg18.phylip = 1 - 269;
 Partition Dummy = 1:261-aggressive-phylip/7924_hg18.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

