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

Input data from file "683-tomodel-phylip/chr4_11166.phylip" (7 taxa, 1024 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	-3149.9143	6339.8287	0.5764	0.5764
GTRI	20	-3150.9153	6341.8305	0.2119	0.7883
GTRIG	21	-3149.9324	6341.8648	0.2083	0.9965
SYMG	17	-3158.5963	6351.1925	0.0020	0.9985
SYMI	17	-3159.5266	6353.0531	0.0008	0.9993
SYMIG	18	-3158.6229	6353.2457	0.0007	1.0000
HKYG	16	-3164.7261	6361.4522	0.0000	1.0000
HKYI	16	-3165.2535	6362.5071	0.0000	1.0000
HKYIG	17	-3164.7639	6363.5278	0.0000	1.0000
K2PG	13	-3170.6082	6367.2164	0.0000	1.0000
K2PI	13	-3171.1968	6368.3935	0.0000	1.0000
K2PIG	14	-3170.6507	6369.3014	0.0000	1.0000
GTR	19	-3176.3333	6390.6665	0.0000	1.0000
SYM	16	-3184.5530	6401.1059	0.0000	1.0000
HKY	15	-3189.6919	6409.3838	0.0000	1.0000
K2P	12	-3195.7619	6415.5239	0.0000	1.0000
F81G	15	-3265.8372	6561.6744	0.0000	1.0000
F81I	15	-3266.2302	6562.4605	0.0000	1.0000
F81IG	16	-3265.8636	6563.7272	0.0000	1.0000
JC69G	12	-3270.1510	6564.3020	0.0000	1.0000
JC69I	12	-3270.6812	6565.3625	0.0000	1.0000
JC69IG	13	-3270.1822	6566.3643	0.0000	1.0000
F81	14	-3287.8580	6603.7160	0.0000	1.0000
JC69	11	-3292.4681	6606.9363	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	13	-3170.6082	6431.3256	0.5359	0.5359
K2PI	13	-3171.1968	6432.5026	0.2975	0.8333
SYMG	17	-3158.5963	6435.0275	0.0842	0.9175
SYMI	17	-3159.5266	6436.8881	0.0332	0.9507
K2PIG	14	-3170.6507	6438.3420	0.0160	0.9668
GTRG	20	-3149.9143	6438.4581	0.0151	0.9819
HKYG	16	-3164.7261	6440.3558	0.0059	0.9878
GTRI	20	-3150.9153	6440.4599	0.0056	0.9933
HKYI	16	-3165.2535	6441.4106	0.0035	0.9968
SYMIG	18	-3158.6229	6442.0122	0.0026	0.9994
GTRIG	21	-3149.9324	6445.4257	0.0005	0.9998
HKYIG	17	-3164.7639	6447.3628	0.0002	1.0000
K2P	12	-3195.7619	6474.7015	0.0000	1.0000
SYM	16	-3184.5530	6480.0095	0.0000	1.0000
HKY	15	-3189.6919	6483.3559	0.0000	1.0000
GTR	19	-3176.3333	6484.3645	0.0000	1.0000
JC69G	12	-3270.1510	6623.4797	0.0000	1.0000
JC69I	12	-3270.6812	6624.5401	0.0000	1.0000
JC69IG	13	-3270.1822	6630.4735	0.0000	1.0000
F81G	15	-3265.8372	6635.6464	0.0000	1.0000
F81I	15	-3266.2302	6636.4326	0.0000	1.0000
F81IG	16	-3265.8636	6642.6308	0.0000	1.0000
JC69	11	-3292.4681	6661.1824	0.0000	1.0000
F81	14	-3287.8580	6672.7566	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/chr4_11166.phylip = 1 - 1024;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11166.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/chr4_11166.phylip = 1 - 1024;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11166.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/chr4_11166.phylip = 1 - 1024;
 Partition Dummy = 1:683-tomodel-phylip/chr4_11166.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

