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

Input data from file "591-tomodel-phylip/chr8_4319.phylip" (7 taxa, 747 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	20	-2565.8951	5172.9473	0.4746	0.4746
GTRI	20	-2566.8633	5174.8836	0.1802	0.6548
GTRIG	21	-2565.8956	5175.0656	0.1646	0.8194
HKYG	16	-2571.7247	5176.1946	0.0936	0.9130
HKYI	16	-2572.2677	5177.2806	0.0544	0.9673
HKYIG	17	-2571.7370	5178.3135	0.0324	0.9998
GTR	19	-2574.9760	5188.9975	0.0002	0.9999
HKY	15	-2579.8667	5190.3900	0.0001	1.0000
K2PG	13	-2598.0867	5222.6700	0.0000	1.0000
K2PI	13	-2598.2095	5222.9156	0.0000	1.0000
K2PIG	14	-2598.0201	5224.6140	0.0000	1.0000
SYMG	17	-2595.9552	5226.7499	0.0000	1.0000
SYMI	17	-2596.1516	5227.1428	0.0000	1.0000
SYMIG	18	-2595.9188	5228.7771	0.0000	1.0000
K2P	12	-2605.9107	5236.2465	0.0000	1.0000
SYM	16	-2603.8539	5240.4530	0.0000	1.0000
F81G	15	-2646.1240	5322.9046	0.0000	1.0000
F81I	15	-2646.3767	5323.4100	0.0000	1.0000
F81IG	16	-2646.1102	5324.9656	0.0000	1.0000
F81	14	-2652.7676	5334.1089	0.0000	1.0000
JC69G	12	-2665.7443	5355.9137	0.0000	1.0000
JC69I	12	-2665.8516	5356.1282	0.0000	1.0000
JC69IG	13	-2665.6875	5357.8716	0.0000	1.0000
JC69	11	-2672.6094	5367.5779	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2565.8951	5171.7903	0.4854	0.4854
GTRI	20	-2566.8633	5173.7265	0.1844	0.6698
GTRIG	21	-2565.8956	5173.7911	0.1785	0.8483
HKYG	16	-2571.7247	5175.4494	0.0779	0.9262
HKYI	16	-2572.2677	5176.5354	0.0453	0.9715
HKYIG	17	-2571.7370	5177.4740	0.0283	0.9998
GTR	19	-2574.9760	5187.9521	0.0002	0.9999
HKY	15	-2579.8667	5189.7334	0.0001	1.0000
K2PG	13	-2598.0867	5222.1734	0.0000	1.0000
K2PI	13	-2598.2095	5222.4191	0.0000	1.0000
K2PIG	14	-2598.0201	5224.0403	0.0000	1.0000
SYMG	17	-2595.9552	5225.9104	0.0000	1.0000
SYMI	17	-2596.1516	5226.3033	0.0000	1.0000
SYMIG	18	-2595.9188	5227.8375	0.0000	1.0000
K2P	12	-2605.9107	5235.8215	0.0000	1.0000
SYM	16	-2603.8539	5239.7078	0.0000	1.0000
F81G	15	-2646.1240	5322.2479	0.0000	1.0000
F81I	15	-2646.3767	5322.7534	0.0000	1.0000
F81IG	16	-2646.1102	5324.2204	0.0000	1.0000
F81	14	-2652.7676	5333.5351	0.0000	1.0000
JC69G	12	-2665.7443	5355.4886	0.0000	1.0000
JC69I	12	-2665.8516	5355.7031	0.0000	1.0000
JC69IG	13	-2665.6875	5357.3750	0.0000	1.0000
JC69	11	-2672.6094	5367.2187	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2571.7247	5249.3065	0.6150	0.6150
HKYI	16	-2572.2677	5250.3925	0.3573	0.9723
HKYIG	17	-2571.7370	5255.9471	0.0222	0.9946
HKY	15	-2579.8667	5258.9744	0.0049	0.9995
GTRG	20	-2565.8951	5264.1116	0.0004	0.9998
GTRI	20	-2566.8633	5266.0478	0.0001	1.0000
GTRIG	21	-2565.8956	5270.7285	0.0000	1.0000
GTR	19	-2574.9760	5275.6573	0.0000	1.0000
K2PG	13	-2598.0867	5282.1822	0.0000	1.0000
K2PI	13	-2598.2095	5282.4279	0.0000	1.0000
K2PIG	14	-2598.0201	5288.6652	0.0000	1.0000
K2P	12	-2605.9107	5291.2143	0.0000	1.0000
SYMG	17	-2595.9552	5304.3835	0.0000	1.0000
SYMI	17	-2596.1516	5304.7764	0.0000	1.0000
SYMIG	18	-2595.9188	5310.9267	0.0000	1.0000
SYM	16	-2603.8539	5313.5648	0.0000	1.0000
F81G	15	-2646.1240	5391.4889	0.0000	1.0000
F81I	15	-2646.3767	5391.9943	0.0000	1.0000
F81IG	16	-2646.1102	5398.0774	0.0000	1.0000
F81	14	-2652.7676	5398.1601	0.0000	1.0000
JC69G	12	-2665.7443	5410.8814	0.0000	1.0000
JC69I	12	-2665.8516	5411.0959	0.0000	1.0000
JC69IG	13	-2665.6875	5417.3839	0.0000	1.0000
JC69	11	-2672.6094	5417.9954	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/chr8_4319.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4319.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/chr8_4319.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4319.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 591-tomodel-phylip/chr8_4319.phylip = 1 - 747;
 Partition Dummy = 1:591-tomodel-phylip/chr8_4319.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

