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

Input data from file "683-tomodel-phylip/chr3_17583.phylip" (7 taxa, 831 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	-2517.4570	5075.9511	0.5897	0.5897
GTRIG	21	-2517.4328	5078.0077	0.2109	0.8006
GTRI	20	-2518.6379	5078.3129	0.1811	0.9817
HKYG	16	-2525.5599	5083.7881	0.0117	0.9934
HKYIG	17	-2525.5553	5085.8634	0.0042	0.9976
HKYI	16	-2527.1326	5086.9335	0.0024	1.0000
GTR	19	-2533.3339	5105.6049	0.0000	1.0000
HKY	15	-2541.5078	5113.6046	0.0000	1.0000
SYMG	17	-2547.5667	5129.8861	0.0000	1.0000
SYMIG	18	-2547.5574	5131.9573	0.0000	1.0000
SYMI	17	-2549.3214	5133.3955	0.0000	1.0000
K2PG	13	-2553.4853	5133.4160	0.0000	1.0000
K2PIG	14	-2553.4827	5135.4801	0.0000	1.0000
K2PI	13	-2555.4165	5137.2785	0.0000	1.0000
SYM	16	-2562.7020	5158.0724	0.0000	1.0000
K2P	12	-2568.8796	5162.1406	0.0000	1.0000
F81G	15	-2597.1112	5224.8113	0.0000	1.0000
F81IG	16	-2597.1110	5226.8904	0.0000	1.0000
F81I	15	-2598.2773	5227.1436	0.0000	1.0000
F81	14	-2611.2313	5250.9773	0.0000	1.0000
JC69G	12	-2618.7341	5261.8496	0.0000	1.0000
JC69IG	13	-2618.7342	5263.9140	0.0000	1.0000
JC69I	12	-2620.1188	5264.6189	0.0000	1.0000
JC69	11	-2633.2930	5288.9084	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	20	-2517.4570	5074.9141	0.5848	0.5848
GTRIG	21	-2517.4328	5076.8655	0.2204	0.8052
GTRI	20	-2518.6379	5077.2759	0.1795	0.9848
HKYG	16	-2525.5599	5083.1198	0.0097	0.9944
HKYIG	17	-2525.5553	5085.1106	0.0036	0.9980
HKYI	16	-2527.1326	5086.2652	0.0020	1.0000
GTR	19	-2533.3339	5104.6678	0.0000	1.0000
HKY	15	-2541.5078	5113.0156	0.0000	1.0000
SYMG	17	-2547.5667	5129.1334	0.0000	1.0000
SYMIG	18	-2547.5574	5131.1149	0.0000	1.0000
SYMI	17	-2549.3214	5132.6427	0.0000	1.0000
K2PG	13	-2553.4853	5132.9705	0.0000	1.0000
K2PIG	14	-2553.4827	5134.9654	0.0000	1.0000
K2PI	13	-2555.4165	5136.8330	0.0000	1.0000
SYM	16	-2562.7020	5157.4041	0.0000	1.0000
K2P	12	-2568.8796	5161.7592	0.0000	1.0000
F81G	15	-2597.1112	5224.2224	0.0000	1.0000
F81IG	16	-2597.1110	5226.2221	0.0000	1.0000
F81I	15	-2598.2773	5226.5546	0.0000	1.0000
F81	14	-2611.2313	5250.4626	0.0000	1.0000
JC69G	12	-2618.7341	5261.4681	0.0000	1.0000
JC69IG	13	-2618.7342	5263.4685	0.0000	1.0000
JC69I	12	-2620.1188	5264.2375	0.0000	1.0000
JC69	11	-2633.2930	5288.5861	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	16	-2525.5599	5158.6819	0.8008	0.8008
HKYI	16	-2527.1326	5161.8272	0.1662	0.9669
HKYIG	17	-2525.5553	5165.3953	0.0279	0.9949
GTRG	20	-2517.4570	5169.3667	0.0038	0.9987
GTRI	20	-2518.6379	5171.7285	0.0012	0.9999
GTRIG	21	-2517.4328	5176.0407	0.0001	1.0000
HKY	15	-2541.5078	5183.8550	0.0000	1.0000
K2PG	13	-2553.4853	5194.3647	0.0000	1.0000
GTR	19	-2533.3339	5194.3978	0.0000	1.0000
K2PI	13	-2555.4165	5198.2272	0.0000	1.0000
K2PIG	14	-2553.4827	5201.0822	0.0000	1.0000
SYMG	17	-2547.5667	5209.4181	0.0000	1.0000
SYMI	17	-2549.3214	5212.9274	0.0000	1.0000
SYMIG	18	-2547.5574	5216.1222	0.0000	1.0000
K2P	12	-2568.8796	5218.4308	0.0000	1.0000
SYM	16	-2562.7020	5232.9661	0.0000	1.0000
F81G	15	-2597.1112	5295.0618	0.0000	1.0000
F81I	15	-2598.2773	5297.3941	0.0000	1.0000
F81IG	16	-2597.1110	5301.7841	0.0000	1.0000
F81	14	-2611.2313	5316.5794	0.0000	1.0000
JC69G	12	-2618.7341	5318.1397	0.0000	1.0000
JC69I	12	-2620.1188	5320.9091	0.0000	1.0000
JC69IG	13	-2618.7342	5324.8627	0.0000	1.0000
JC69	11	-2633.2930	5340.5350	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/chr3_17583.phylip = 1 - 831;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17583.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/chr3_17583.phylip = 1 - 831;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17583.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/chr3_17583.phylip = 1 - 831;
 Partition Dummy = 1:683-tomodel-phylip/chr3_17583.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

