
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
21:32 04/04/2011

Input data from file "591-tomodel-phylip/chr1_29891.phylip" (7 taxa, 441 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2P


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	20	-1122.5159	2287.0318	0.3843	0.3843
GTRI	20	-1122.5865	2287.1730	0.3581	0.7424
GTRIG	21	-1122.5196	2289.2446	0.1271	0.8695
GTR	19	-1124.8890	2289.5831	0.1073	0.9768
HKYG	16	-1130.8776	2295.0382	0.0070	0.9838
HKYI	16	-1131.0899	2295.4629	0.0057	0.9895
HKYIG	17	-1130.8877	2297.2223	0.0024	0.9918
HKY	15	-1133.1799	2297.4893	0.0021	0.9939
SYMG	17	-1131.0997	2297.6463	0.0019	0.9958
SYMI	17	-1131.2239	2297.8945	0.0017	0.9975
SYMIG	18	-1131.1045	2299.8298	0.0006	0.9981
K2PG	13	-1136.6146	2300.0816	0.0006	0.9987
SYM	16	-1133.4156	2300.1142	0.0006	0.9992
K2PI	13	-1136.8436	2300.5396	0.0004	0.9997
K2PIG	14	-1136.6213	2302.2285	0.0002	0.9999
K2P	12	-1139.0997	2302.9284	0.0001	1.0000
F81G	15	-1152.3790	2335.8874	0.0000	1.0000
F81I	15	-1152.4639	2336.0572	0.0000	1.0000
F81	14	-1154.4146	2337.8152	0.0000	1.0000
F81IG	16	-1152.3845	2338.0519	0.0000	1.0000
JC69G	12	-1157.8293	2340.3876	0.0000	1.0000
JC69I	12	-1157.9412	2340.6114	0.0000	1.0000
JC69IG	13	-1157.8356	2342.5237	0.0000	1.0000
JC69	11	-1159.9861	2342.5877	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-1122.5159	2285.0318	0.3855	0.3855
GTRI	20	-1122.5865	2285.1730	0.3592	0.7446
GTRIG	21	-1122.5196	2287.0393	0.1413	0.8859
GTR	19	-1124.8890	2287.7779	0.0976	0.9836
HKYG	16	-1130.8776	2293.7551	0.0049	0.9885
HKYI	16	-1131.0899	2294.1798	0.0040	0.9925
HKYIG	17	-1130.8877	2295.7755	0.0018	0.9943
SYMG	17	-1131.0997	2296.1995	0.0014	0.9957
HKY	15	-1133.1799	2296.3598	0.0013	0.9970
SYMI	17	-1131.2239	2296.4477	0.0013	0.9983
SYMIG	18	-1131.1045	2298.2089	0.0005	0.9989
SYM	16	-1133.4156	2298.8312	0.0004	0.9992
K2PG	13	-1136.6146	2299.2291	0.0003	0.9996
K2PI	13	-1136.8436	2299.6871	0.0003	0.9998
K2PIG	14	-1136.6213	2301.2426	0.0001	0.9999
K2P	12	-1139.0997	2302.1994	0.0001	1.0000
F81G	15	-1152.3790	2334.7579	0.0000	1.0000
F81I	15	-1152.4639	2334.9278	0.0000	1.0000
F81IG	16	-1152.3845	2336.7689	0.0000	1.0000
F81	14	-1154.4146	2336.8293	0.0000	1.0000
JC69G	12	-1157.8293	2339.6586	0.0000	1.0000
JC69I	12	-1157.9412	2339.8824	0.0000	1.0000
JC69IG	13	-1157.8356	2341.6713	0.0000	1.0000
JC69	11	-1159.9861	2341.9723	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2P	12	-1139.0997	2351.2679	0.4687	0.4687
K2PG	13	-1136.6146	2352.3867	0.2679	0.7366
K2PI	13	-1136.8436	2352.8447	0.2131	0.9497
HKY	15	-1133.1799	2357.6955	0.0188	0.9685
K2PIG	14	-1136.6213	2358.4892	0.0127	0.9812
HKYG	16	-1130.8776	2359.1799	0.0090	0.9902
HKYI	16	-1131.0899	2359.6046	0.0073	0.9974
SYM	16	-1133.4156	2364.2559	0.0007	0.9981
HKYIG	17	-1130.8877	2365.2892	0.0004	0.9986
GTR	19	-1124.8890	2365.4698	0.0004	0.9989
SYMG	17	-1131.0997	2365.7132	0.0003	0.9993
SYMI	17	-1131.2239	2365.9615	0.0003	0.9996
GTRG	20	-1122.5159	2366.8127	0.0002	0.9998
GTRI	20	-1122.5865	2366.9539	0.0002	1.0000
SYMIG	18	-1131.1045	2371.8117	0.0000	1.0000
GTRIG	21	-1122.5196	2372.9092	0.0000	1.0000
JC69	11	-1159.9861	2386.9518	0.0000	1.0000
JC69G	12	-1157.8293	2388.7272	0.0000	1.0000
JC69I	12	-1157.9412	2388.9510	0.0000	1.0000
F81	14	-1154.4146	2394.0759	0.0000	1.0000
JC69IG	13	-1157.8356	2394.8289	0.0000	1.0000
F81G	15	-1152.3790	2396.0936	0.0000	1.0000
F81I	15	-1152.4639	2396.2635	0.0000	1.0000
F81IG	16	-1152.3845	2402.1936	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 591-tomodel-phylip/chr1_29891.phylip = 1 - 441;
 Partition Dummy = 1:591-tomodel-phylip/chr1_29891.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 591-tomodel-phylip/chr1_29891.phylip = 1 - 441;
 Partition Dummy = 1:591-tomodel-phylip/chr1_29891.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 (K2P)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr1_29891.phylip = 1 - 441;
 Partition Dummy = 1:591-tomodel-phylip/chr1_29891.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0);
END;


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

