
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
16:12 04/04/2011

Input data from file "683-tomodel-phylip/chr2_4311.phylip" (7 taxa, 780 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.6922	5172.4911	0.6313	0.6313
GTRIG	21	-2565.7015	5174.6220	0.2175	0.8488
GTRI	20	-2567.3690	5175.8447	0.1180	0.9669
HKYG	16	-2573.2977	5179.3083	0.0209	0.9878
HKYIG	17	-2573.3054	5181.4140	0.0073	0.9951
HKYI	16	-2574.7407	5182.1944	0.0049	1.0000
GTR	19	-2589.0259	5217.0519	0.0000	1.0000
HKY	15	-2595.7635	5222.1552	0.0000	1.0000
SYMG	17	-2610.1727	5255.1485	0.0000	1.0000
SYMI	17	-2610.5943	5255.9918	0.0000	1.0000
K2PG	13	-2615.1425	5256.7602	0.0000	1.0000
SYMIG	18	-2610.2132	5257.3253	0.0000	1.0000
K2PI	13	-2615.8117	5258.0986	0.0000	1.0000
K2PIG	14	-2615.1611	5258.8713	0.0000	1.0000
SYM	16	-2631.2727	5295.2584	0.0000	1.0000
K2P	12	-2636.0871	5296.5811	0.0000	1.0000
F81G	15	-2653.5107	5337.6497	0.0000	1.0000
F81I	15	-2653.8938	5338.4159	0.0000	1.0000
F81IG	16	-2653.5388	5339.7906	0.0000	1.0000
F81	14	-2671.9609	5372.4708	0.0000	1.0000
JC69G	12	-2684.3796	5393.1659	0.0000	1.0000
JC69I	12	-2684.5149	5393.4365	0.0000	1.0000
JC69IG	13	-2684.3837	5395.2425	0.0000	1.0000
JC69	11	-2702.6728	5427.6894	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2565.6922	5171.3844	0.6270	0.6270
GTRIG	21	-2565.7015	5173.4030	0.2285	0.8555
GTRI	20	-2567.3690	5174.7380	0.1172	0.9727
HKYG	16	-2573.2977	5178.5953	0.0170	0.9898
HKYIG	17	-2573.3054	5180.6108	0.0062	0.9960
HKYI	16	-2574.7407	5181.4814	0.0040	1.0000
GTR	19	-2589.0259	5216.0519	0.0000	1.0000
HKY	15	-2595.7635	5221.5269	0.0000	1.0000
SYMG	17	-2610.1727	5254.3454	0.0000	1.0000
SYMI	17	-2610.5943	5255.1887	0.0000	1.0000
K2PG	13	-2615.1425	5256.2850	0.0000	1.0000
SYMIG	18	-2610.2132	5256.4265	0.0000	1.0000
K2PI	13	-2615.8117	5257.6234	0.0000	1.0000
K2PIG	14	-2615.1611	5258.3223	0.0000	1.0000
SYM	16	-2631.2727	5294.5454	0.0000	1.0000
K2P	12	-2636.0871	5296.1743	0.0000	1.0000
F81G	15	-2653.5107	5337.0214	0.0000	1.0000
F81I	15	-2653.8938	5337.7876	0.0000	1.0000
F81IG	16	-2653.5388	5339.0776	0.0000	1.0000
F81	14	-2671.9609	5371.9218	0.0000	1.0000
JC69G	12	-2684.3796	5392.7591	0.0000	1.0000
JC69I	12	-2684.5149	5393.0297	0.0000	1.0000
JC69IG	13	-2684.3837	5394.7673	0.0000	1.0000
JC69	11	-2702.6728	5427.3457	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2573.2977	5253.1440	0.7838	0.7838
HKYI	16	-2574.7407	5256.0301	0.1851	0.9690
HKYIG	17	-2573.3054	5259.8188	0.0278	0.9968
GTRG	20	-2565.6922	5264.5703	0.0026	0.9994
GTRI	20	-2567.3690	5267.9238	0.0005	0.9999
GTRIG	21	-2565.7015	5271.2482	0.0001	1.0000
HKY	15	-2595.7635	5291.4163	0.0000	1.0000
GTR	19	-2589.0259	5304.5785	0.0000	1.0000
K2PG	13	-2615.1425	5316.8558	0.0000	1.0000
K2PI	13	-2615.8117	5318.1942	0.0000	1.0000
K2PIG	14	-2615.1611	5323.5524	0.0000	1.0000
SYMG	17	-2610.1727	5333.5534	0.0000	1.0000
SYMI	17	-2610.5943	5334.3967	0.0000	1.0000
SYMIG	18	-2610.2132	5340.2938	0.0000	1.0000
K2P	12	-2636.0871	5352.0858	0.0000	1.0000
SYM	16	-2631.2727	5369.0941	0.0000	1.0000
F81G	15	-2653.5107	5406.9108	0.0000	1.0000
F81I	15	-2653.8938	5407.6770	0.0000	1.0000
F81IG	16	-2653.5388	5413.6263	0.0000	1.0000
F81	14	-2671.9609	5437.1519	0.0000	1.0000
JC69G	12	-2684.3796	5448.6706	0.0000	1.0000
JC69I	12	-2684.5149	5448.9412	0.0000	1.0000
JC69IG	13	-2684.3837	5455.3381	0.0000	1.0000
JC69	11	-2702.6728	5478.5979	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 683-tomodel-phylip/chr2_4311.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4311.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 683-tomodel-phylip/chr2_4311.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4311.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 683-tomodel-phylip/chr2_4311.phylip = 1 - 780;
 Partition Dummy = 1:683-tomodel-phylip/chr2_4311.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

