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

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

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


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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-3368.0818	6776.1635	0.3357	0.3357
GTRI	20	-3368.5059	6777.0118	0.2197	0.5553
GTRIG	21	-3367.9296	6777.8591	0.1438	0.6991
HKYI	16	-3373.1016	6778.2032	0.1211	0.8202
HKYG	16	-3373.1159	6778.2318	0.1193	0.9395
HKYIG	17	-3372.7961	6779.5922	0.0605	1.0000
HKY	15	-3402.7836	6835.5672	0.0000	1.0000
GTR	19	-3399.0867	6836.1735	0.0000	1.0000
F81I	15	-3436.1423	6902.2845	0.0000	1.0000
F81G	15	-3436.9141	6903.8282	0.0000	1.0000
F81IG	16	-3436.1283	6904.2567	0.0000	1.0000
K2PG	13	-3455.1359	6936.2717	0.0000	1.0000
K2PI	13	-3455.6188	6937.2375	0.0000	1.0000
K2PIG	14	-3455.0473	6938.0947	0.0000	1.0000
SYMG	17	-3452.7877	6939.5754	0.0000	1.0000
SYMI	17	-3453.5468	6941.0936	0.0000	1.0000
SYMIG	18	-3452.7646	6941.5291	0.0000	1.0000
F81	14	-3463.5955	6955.1909	0.0000	1.0000
K2P	12	-3486.3940	6996.7881	0.0000	1.0000
SYM	16	-3484.9483	7001.8966	0.0000	1.0000
JC69I	12	-3505.8029	7035.6058	0.0000	1.0000
JC69G	12	-3506.1319	7036.2638	0.0000	1.0000
JC69IG	13	-3505.6849	7037.3698	0.0000	1.0000
JC69	11	-3534.4105	7090.8209	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	16	-3373.1016	6857.3702	0.4931	0.4931
HKYG	16	-3373.1159	6857.3988	0.4861	0.9791
HKYIG	17	-3372.7961	6863.7071	0.0207	0.9999
GTRG	20	-3368.0818	6875.1223	0.0001	1.0000
GTRI	20	-3368.5059	6875.9705	0.0000	1.0000
GTRIG	21	-3367.9296	6881.7658	0.0000	1.0000
HKY	15	-3402.7836	6909.7863	0.0000	1.0000
GTR	19	-3399.0867	6930.1843	0.0000	1.0000
F81I	15	-3436.1423	6976.5036	0.0000	1.0000
F81G	15	-3436.9141	6978.0473	0.0000	1.0000
F81IG	16	-3436.1283	6983.4237	0.0000	1.0000
K2PG	13	-3455.1359	7000.5949	0.0000	1.0000
K2PI	13	-3455.6188	7001.5607	0.0000	1.0000
K2PIG	14	-3455.0473	7007.3658	0.0000	1.0000
SYMG	17	-3452.7877	7023.6903	0.0000	1.0000
F81	14	-3463.5955	7024.4620	0.0000	1.0000
SYMI	17	-3453.5468	7025.2086	0.0000	1.0000
SYMIG	18	-3452.7646	7030.5920	0.0000	1.0000
K2P	12	-3486.3940	7056.1633	0.0000	1.0000
SYM	16	-3484.9483	7081.0636	0.0000	1.0000
JC69I	12	-3505.8029	7094.9811	0.0000	1.0000
JC69G	12	-3506.1319	7095.6391	0.0000	1.0000
JC69IG	13	-3505.6849	7101.6930	0.0000	1.0000
JC69	11	-3534.4105	7145.2482	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/chr2_26639.phylip = 1 - 1041;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26639.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/chr2_26639.phylip = 1 - 1041;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26639.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 (HKYI)]

BEGIN MRBAYES;
 Charset 591-tomodel-phylip/chr2_26639.phylip = 1 - 1041;
 Partition Dummy = 1:591-tomodel-phylip/chr2_26639.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

