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

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

Minimum AIC  model: GTRG
Minimum AICc model: HKYI
Minimum BIC  model: HKYI


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
HKYI	60	-1019.2660	2184.8630	0.3473	0.3473
HKYG	60	-1019.3846	2185.1001	0.3084	0.6557
GTRG	64	-1014.1939	2186.7528	0.1350	0.7907
GTRI	64	-1014.4652	2187.2954	0.1029	0.8936
HKYIG	61	-1019.3844	2188.0757	0.0697	0.9633
GTRIG	65	-1014.1961	2189.8207	0.0291	0.9924
K2PG	57	-1028.4534	2194.4369	0.0029	0.9953
K2PI	57	-1028.7386	2195.0074	0.0022	0.9975
HKY	59	-1026.3648	2196.1060	0.0013	0.9987
K2PIG	58	-1028.4526	2197.3480	0.0007	0.9994
SYMG	61	-1025.2070	2199.7208	0.0002	0.9996
SYMI	61	-1025.4422	2200.1913	0.0002	0.9998
GTR	63	-1022.5044	2200.3324	0.0002	0.9999
SYMIG	62	-1025.2075	2202.7194	0.0000	1.0000
K2P	56	-1034.8329	2204.3040	0.0000	1.0000
SYM	60	-1032.0919	2210.5147	0.0000	1.0000
F81G	59	-1036.2519	2215.8801	0.0000	1.0000
F81I	59	-1036.4355	2216.2473	0.0000	1.0000
F81IG	60	-1036.2365	2218.8040	0.0000	1.0000
JC69G	56	-1044.7224	2224.0830	0.0000	1.0000
JC69I	56	-1044.9190	2224.4763	0.0000	1.0000
F81	58	-1042.9489	2226.3407	0.0000	1.0000
JC69IG	57	-1044.7091	2226.9485	0.0000	1.0000
JC69	55	-1051.2000	2234.1667	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1014.1939	2156.3878	0.3459	0.3459
GTRI	64	-1014.4652	2156.9304	0.2637	0.6096
GTRIG	65	-1014.1961	2158.3921	0.1270	0.7365
HKYI	60	-1019.2660	2158.5320	0.1184	0.8549
HKYG	60	-1019.3846	2158.7692	0.1051	0.9600
HKYIG	61	-1019.3844	2160.7688	0.0387	0.9987
HKY	59	-1026.3648	2170.7296	0.0003	0.9990
K2PG	57	-1028.4534	2170.9067	0.0002	0.9992
GTR	63	-1022.5044	2171.0088	0.0002	0.9995
K2PI	57	-1028.7386	2171.4772	0.0002	0.9997
SYMG	61	-1025.2070	2172.4140	0.0001	0.9998
SYMI	61	-1025.4422	2172.8845	0.0001	0.9999
K2PIG	58	-1028.4526	2172.9052	0.0001	1.0000
SYMIG	62	-1025.2075	2174.4150	0.0000	1.0000
K2P	56	-1034.8329	2181.6657	0.0000	1.0000
SYM	60	-1032.0919	2184.1838	0.0000	1.0000
F81G	59	-1036.2519	2190.5038	0.0000	1.0000
F81I	59	-1036.4355	2190.8710	0.0000	1.0000
F81IG	60	-1036.2365	2192.4730	0.0000	1.0000
JC69G	56	-1044.7224	2201.4447	0.0000	1.0000
JC69I	56	-1044.9190	2201.8380	0.0000	1.0000
F81	58	-1042.9489	2201.8978	0.0000	1.0000
JC69IG	57	-1044.7091	2203.4183	0.0000	1.0000
JC69	55	-1051.2000	2212.4000	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-1019.2660	2388.0920	0.3197	0.3197
HKYG	60	-1019.3846	2388.3292	0.2839	0.6036
K2PG	57	-1028.4534	2388.9887	0.2042	0.8078
K2PI	57	-1028.7386	2389.5592	0.1535	0.9613
HKYIG	61	-1019.3844	2394.1548	0.0154	0.9768
K2PIG	58	-1028.4526	2394.8132	0.0111	0.9879
K2P	56	-1034.8329	2395.9217	0.0064	0.9942
HKY	59	-1026.3648	2396.4636	0.0049	0.9991
GTRG	64	-1014.1939	2401.2518	0.0004	0.9995
GTRI	64	-1014.4652	2401.7944	0.0003	0.9999
SYMG	61	-1025.2070	2405.8000	0.0000	0.9999
SYMI	61	-1025.4422	2406.2705	0.0000	1.0000
GTRIG	65	-1014.1961	2407.0821	0.0000	1.0000
SYMIG	62	-1025.2075	2411.6270	0.0000	1.0000
GTR	63	-1022.5044	2412.0468	0.0000	1.0000
SYM	60	-1032.0919	2413.7438	0.0000	1.0000
JC69G	56	-1044.7224	2415.7007	0.0000	1.0000
JC69I	56	-1044.9190	2416.0940	0.0000	1.0000
F81G	59	-1036.2519	2416.2378	0.0000	1.0000
F81I	59	-1036.4355	2416.6050	0.0000	1.0000
JC69IG	57	-1044.7091	2421.5003	0.0000	1.0000
F81IG	60	-1036.2365	2422.0330	0.0000	1.0000
JC69	55	-1051.2000	2422.8300	0.0000	1.0000
F81	58	-1042.9489	2423.8058	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 183-phylip/chr7_11511.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr7_11511.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 (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11511.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr7_11511.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

BEGIN MRBAYES;
 Charset 183-phylip/chr7_11511.phylip = 1 - 339;
 Partition Dummy = 1:183-phylip/chr7_11511.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=propinv;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0) pinvarpr=Uniform(0.0,1.0);
END;


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

