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

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

Minimum AIC  model: HKYIG
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	-1023.9334	2188.8410	0.3868	0.3868
HKYIG	61	-1022.8239	2189.3835	0.2949	0.6817
HKYI	60	-1024.2294	2189.4331	0.2877	0.9694
GTRG	64	-1021.8305	2195.7770	0.0121	0.9815
GTRIG	65	-1020.5951	2196.1321	0.0101	0.9916
GTRI	64	-1022.1875	2196.4910	0.0084	1.0000
K2PG	57	-1043.5425	2219.8691	0.0000	1.0000
K2PIG	58	-1042.4005	2220.2996	0.0000	1.0000
K2PI	57	-1043.7951	2220.3742	0.0000	1.0000
HKY	59	-1041.7077	2221.6440	0.0000	1.0000
SYMG	61	-1042.9337	2229.6030	0.0000	1.0000
GTR	63	-1040.2594	2229.8251	0.0000	1.0000
SYMIG	62	-1041.7757	2230.0643	0.0000	1.0000
SYMI	61	-1043.1756	2230.0869	0.0000	1.0000
F81G	59	-1054.9055	2248.0397	0.0000	1.0000
F81I	59	-1055.1394	2248.5074	0.0000	1.0000
F81IG	60	-1054.0813	2249.1368	0.0000	1.0000
K2P	56	-1061.0800	2252.2449	0.0000	1.0000
SYM	60	-1060.5513	2262.0767	0.0000	1.0000
JC69G	56	-1070.9686	2272.0223	0.0000	1.0000
JC69I	56	-1071.1950	2272.4750	0.0000	1.0000
JC69IG	57	-1070.0899	2272.9638	0.0000	1.0000
F81	58	-1071.7832	2279.0651	0.0000	1.0000
JC69	55	-1087.9100	2303.2211	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYIG	61	-1022.8239	2167.6478	0.3377	0.3377
HKYG	60	-1023.9334	2167.8667	0.3026	0.6403
HKYI	60	-1024.2294	2168.4589	0.2251	0.8654
GTRIG	65	-1020.5951	2171.1902	0.0574	0.9228
GTRG	64	-1021.8305	2171.6611	0.0454	0.9682
GTRI	64	-1022.1875	2172.3751	0.0318	1.0000
K2PIG	58	-1042.4005	2200.8010	0.0000	1.0000
K2PG	57	-1043.5425	2201.0850	0.0000	1.0000
HKY	59	-1041.7077	2201.4154	0.0000	1.0000
K2PI	57	-1043.7951	2201.5901	0.0000	1.0000
GTR	63	-1040.2594	2206.5187	0.0000	1.0000
SYMIG	62	-1041.7757	2207.5514	0.0000	1.0000
SYMG	61	-1042.9337	2207.8674	0.0000	1.0000
SYMI	61	-1043.1756	2208.3513	0.0000	1.0000
F81G	59	-1054.9055	2227.8111	0.0000	1.0000
F81IG	60	-1054.0813	2228.1626	0.0000	1.0000
F81I	59	-1055.1394	2228.2788	0.0000	1.0000
K2P	56	-1061.0800	2234.1599	0.0000	1.0000
SYM	60	-1060.5513	2241.1025	0.0000	1.0000
JC69G	56	-1070.9686	2253.9373	0.0000	1.0000
JC69IG	57	-1070.0899	2254.1797	0.0000	1.0000
JC69I	56	-1071.1950	2254.3900	0.0000	1.0000
F81	58	-1071.7832	2259.5665	0.0000	1.0000
JC69	55	-1087.9100	2285.8200	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-1023.9334	2408.8362	0.5281	0.5281
HKYI	60	-1024.2294	2409.4283	0.3928	0.9208
HKYIG	61	-1022.8239	2412.6334	0.0791	0.9999
GTRG	64	-1021.8305	2428.6951	0.0000	1.0000
GTRI	64	-1022.1875	2429.4091	0.0000	1.0000
K2PG	57	-1043.5425	2430.0060	0.0000	1.0000
K2PI	57	-1043.7951	2430.5111	0.0000	1.0000
GTRIG	65	-1020.5951	2432.2405	0.0000	1.0000
K2PIG	58	-1042.4005	2433.7381	0.0000	1.0000
HKY	59	-1041.7077	2438.3687	0.0000	1.0000
SYMG	61	-1042.9337	2452.8530	0.0000	1.0000
SYMI	61	-1043.1756	2453.3368	0.0000	1.0000
SYMIG	62	-1041.7757	2456.5531	0.0000	1.0000
K2P	56	-1061.0800	2459.0647	0.0000	1.0000
GTR	63	-1040.2594	2459.5366	0.0000	1.0000
F81G	59	-1054.9055	2464.7644	0.0000	1.0000
F81I	59	-1055.1394	2465.2321	0.0000	1.0000
F81IG	60	-1054.0813	2469.1320	0.0000	1.0000
JC69G	56	-1070.9686	2478.8421	0.0000	1.0000
JC69I	56	-1071.1950	2479.2948	0.0000	1.0000
SYM	60	-1060.5513	2482.0719	0.0000	1.0000
JC69IG	57	-1070.0899	2483.1007	0.0000	1.0000
F81	58	-1071.7832	2492.5036	0.0000	1.0000
JC69	55	-1087.9100	2506.7086	0.0000	1.0000



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

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

BEGIN MRBAYES;
 Charset 261-aggressive-phylip/5848_hg18.phylip = 1 - 410;
 Partition Dummy = 1:261-aggressive-phylip/5848_hg18.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


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

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


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

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

