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

Input data from file "261-aggressive-phylip/13190_hg18.phylip" (29 taxa, 446 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
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
GTRG	64	-1135.6217	2421.0807	0.4573	0.4573
GTRI	64	-1136.2622	2422.3618	0.2410	0.6983
GTRIG	65	-1135.2010	2422.9809	0.1768	0.8752
HKYG	60	-1143.0584	2425.1298	0.0604	0.9356
HKYI	60	-1143.3860	2425.7850	0.0435	0.9791
HKYIG	61	-1142.7878	2427.2735	0.0207	0.9998
F81G	59	-1150.7334	2437.8087	0.0001	0.9999
F81I	59	-1151.0481	2438.4382	0.0001	1.0000
F81IG	60	-1150.4830	2439.9789	0.0000	1.0000
GTR	63	-1147.5799	2442.2697	0.0000	1.0000
HKY	59	-1155.8664	2448.0747	0.0000	1.0000
K2PG	57	-1158.7616	2448.5644	0.0000	1.0000
K2PI	57	-1159.0493	2449.1398	0.0000	1.0000
K2PIG	58	-1158.4833	2450.6514	0.0000	1.0000
SYMG	61	-1155.3124	2452.3227	0.0000	1.0000
SYMI	61	-1155.6501	2452.9982	0.0000	1.0000
SYMIG	62	-1155.0017	2454.4002	0.0000	1.0000
JC69G	56	-1165.4944	2459.4002	0.0000	1.0000
JC69I	56	-1165.7699	2459.9511	0.0000	1.0000
F81	58	-1163.2168	2460.1183	0.0000	1.0000
JC69IG	57	-1165.2460	2461.5332	0.0000	1.0000
K2P	56	-1172.0359	2472.4831	0.0000	1.0000
SYM	60	-1168.3136	2475.6402	0.0000	1.0000
JC69	55	-1178.4879	2482.7707	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1135.6217	2399.2435	0.4633	0.4633
GTRIG	65	-1135.2010	2400.4019	0.2596	0.7229
GTRI	64	-1136.2622	2400.5245	0.2442	0.9671
HKYG	60	-1143.0584	2406.1168	0.0149	0.9820
HKYI	60	-1143.3860	2406.7720	0.0107	0.9928
HKYIG	61	-1142.7878	2407.5756	0.0072	1.0000
F81G	59	-1150.7334	2419.4668	0.0000	1.0000
F81I	59	-1151.0481	2420.0963	0.0000	1.0000
F81IG	60	-1150.4830	2420.9660	0.0000	1.0000
GTR	63	-1147.5799	2421.1598	0.0000	1.0000
HKY	59	-1155.8664	2429.7328	0.0000	1.0000
K2PG	57	-1158.7616	2431.5232	0.0000	1.0000
K2PI	57	-1159.0493	2432.0986	0.0000	1.0000
SYMG	61	-1155.3124	2432.6248	0.0000	1.0000
K2PIG	58	-1158.4833	2432.9666	0.0000	1.0000
SYMI	61	-1155.6501	2433.3003	0.0000	1.0000
SYMIG	62	-1155.0017	2434.0033	0.0000	1.0000
F81	58	-1163.2168	2442.4336	0.0000	1.0000
JC69G	56	-1165.4944	2442.9888	0.0000	1.0000
JC69I	56	-1165.7699	2443.5397	0.0000	1.0000
JC69IG	57	-1165.2460	2444.4920	0.0000	1.0000
K2P	56	-1172.0359	2456.0717	0.0000	1.0000
SYM	60	-1168.3136	2456.6273	0.0000	1.0000
JC69	55	-1178.4879	2466.9758	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1143.0584	2652.1360	0.5505	0.5505
HKYI	60	-1143.3860	2652.7912	0.3967	0.9472
HKYIG	61	-1142.7878	2657.6950	0.0342	0.9813
F81G	59	-1150.7334	2661.3856	0.0054	0.9867
GTRG	64	-1135.6217	2661.6639	0.0047	0.9914
F81I	59	-1151.0481	2662.0151	0.0039	0.9954
GTRI	64	-1136.2622	2662.9449	0.0025	0.9978
K2PG	57	-1158.7616	2665.2414	0.0008	0.9986
K2PI	57	-1159.0493	2665.8167	0.0006	0.9992
GTRIG	65	-1135.2010	2666.9227	0.0003	0.9996
F81IG	60	-1150.4830	2666.9851	0.0003	0.9999
K2PIG	58	-1158.4833	2670.7851	0.0000	0.9999
HKY	59	-1155.8664	2671.6516	0.0000	1.0000
JC69G	56	-1165.4944	2672.6067	0.0000	1.0000
JC69I	56	-1165.7699	2673.1576	0.0000	1.0000
JC69IG	57	-1165.2460	2678.2102	0.0000	1.0000
GTR	63	-1147.5799	2679.4799	0.0000	1.0000
F81	58	-1163.2168	2680.2521	0.0000	1.0000
SYMG	61	-1155.3124	2682.7443	0.0000	1.0000
SYMI	61	-1155.6501	2683.4197	0.0000	1.0000
K2P	56	-1172.0359	2685.6896	0.0000	1.0000
SYMIG	62	-1155.0017	2688.2231	0.0000	1.0000
JC69	55	-1178.4879	2692.4933	0.0000	1.0000
SYM	60	-1168.3136	2702.6464	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13190_hg18.phylip = 1 - 446;
 Partition Dummy = 1:261-aggressive-phylip/13190_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13190_hg18.phylip = 1 - 446;
 Partition Dummy = 1:261-aggressive-phylip/13190_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=gamma;
 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) shapepr=Uniform(0.1,50.0);
END;


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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/13190_hg18.phylip = 1 - 446;
 Partition Dummy = 1:261-aggressive-phylip/13190_hg18.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

