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

Input data from file "683-tomodel-phylip/chr1_31083.phylip" (7 taxa, 1122 characters)

Minimum AIC  model: GTRG
Minimum AICc model: GTRG
Minimum BIC  model: K2PG


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3048.4000	6136.8000	0.6223	0.6223
GTRIG	21	-3048.3050	6138.6099	0.2518	0.8741
GTRI	20	-3050.0176	6140.0353	0.1234	0.9975
SYMG	17	-3057.6109	6149.2218	0.0012	0.9988
SYMIG	18	-3057.5355	6151.0710	0.0005	0.9993
SYMI	17	-3058.9567	6151.9134	0.0003	0.9996
HKYG	16	-3060.2993	6152.5985	0.0002	0.9998
HKYIG	17	-3060.2434	6154.4868	0.0001	0.9999
HKYI	16	-3061.7052	6155.4104	0.0001	1.0000
K2PG	13	-3067.1810	6160.3620	0.0000	1.0000
K2PIG	14	-3067.1341	6162.2681	0.0000	1.0000
K2PI	13	-3068.5293	6163.0587	0.0000	1.0000
GTR	19	-3068.8912	6175.7823	0.0000	1.0000
SYM	16	-3076.6523	6185.3047	0.0000	1.0000
HKY	15	-3079.1853	6188.3706	0.0000	1.0000
K2P	12	-3085.6307	6195.2613	0.0000	1.0000
F81G	15	-3108.2099	6246.4199	0.0000	1.0000
F81IG	16	-3108.1715	6248.3430	0.0000	1.0000
F81I	15	-3109.5685	6249.1369	0.0000	1.0000
JC69G	12	-3114.3468	6252.6936	0.0000	1.0000
JC69IG	13	-3114.3180	6254.6359	0.0000	1.0000
JC69I	12	-3115.6580	6255.3159	0.0000	1.0000
F81	14	-3125.7330	6279.4660	0.0000	1.0000
JC69	11	-3131.5332	6285.0663	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3067.1810	6225.6593	0.7449	0.7449
K2PI	13	-3068.5293	6228.3560	0.1934	0.9383
K2PIG	14	-3067.1341	6232.5883	0.0233	0.9616
HKYG	16	-3060.2993	6232.9644	0.0193	0.9809
SYMG	17	-3057.6109	6234.6105	0.0085	0.9894
HKYI	16	-3061.7052	6235.7763	0.0047	0.9941
GTRG	20	-3048.4000	6237.2573	0.0023	0.9964
SYMI	17	-3058.9567	6237.3022	0.0022	0.9986
HKYIG	17	-3060.2434	6239.8755	0.0006	0.9992
GTRI	20	-3050.0176	6240.4927	0.0004	0.9997
SYMIG	18	-3057.5355	6241.4826	0.0003	0.9999
GTRIG	21	-3048.3050	6244.0901	0.0001	1.0000
K2P	12	-3085.6307	6255.5357	0.0000	1.0000
HKY	15	-3079.1853	6263.7136	0.0000	1.0000
SYM	16	-3076.6523	6265.6705	0.0000	1.0000
GTR	19	-3068.8912	6271.2168	0.0000	1.0000
JC69G	12	-3114.3468	6312.9680	0.0000	1.0000
JC69I	12	-3115.6580	6315.5903	0.0000	1.0000
JC69IG	13	-3114.3180	6319.9332	0.0000	1.0000
F81G	15	-3108.2099	6321.7629	0.0000	1.0000
F81I	15	-3109.5685	6324.4799	0.0000	1.0000
F81IG	16	-3108.1715	6328.7089	0.0000	1.0000
JC69	11	-3131.5332	6340.3179	0.0000	1.0000
F81	14	-3125.7330	6349.7861	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/chr1_31083.phylip = 1 - 1122;
 Partition Dummy = 1:683-tomodel-phylip/chr1_31083.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/chr1_31083.phylip = 1 - 1122;
 Partition Dummy = 1:683-tomodel-phylip/chr1_31083.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 (K2PG)]

BEGIN MRBAYES;
 Charset 683-tomodel-phylip/chr1_31083.phylip = 1 - 1122;
 Partition Dummy = 1:683-tomodel-phylip/chr1_31083.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=gamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0);
END;


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

