
Output from mraic.pl version 1.4.4 by Johan Nylander
-------------------------------------------------------------
00:54 02/17/2011

Input data from file "183-phylip/chr3_764.phylip" (29 taxa, 685 characters)

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: HKYI


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
GTRIG	65	-2081.8523	4307.5657	0.5757	0.5757
GTRI	64	-2083.6542	4308.7277	0.3220	0.8976
GTRG	64	-2084.8086	4311.0365	0.1015	0.9992
HKYI	60	-2095.0769	4321.8845	0.0004	0.9996
HKYG	60	-2095.5858	4322.9025	0.0003	0.9999
HKYIG	61	-2095.1242	4324.3897	0.0001	1.0000
F81I	59	-2126.3357	4381.9995	0.0000	1.0000
F81G	59	-2126.6287	4382.5854	0.0000	1.0000
F81IG	60	-2126.1584	4384.0475	0.0000	1.0000
GTR	63	-2123.7733	4386.5321	0.0000	1.0000
SYMIG	62	-2129.1540	4394.8674	0.0000	1.0000
SYMI	61	-2131.1941	4396.5294	0.0000	1.0000
SYMG	61	-2131.6424	4397.4261	0.0000	1.0000
HKY	59	-2138.6291	4406.5861	0.0000	1.0000
K2PG	57	-2142.3860	4409.3175	0.0000	1.0000
K2PIG	58	-2141.8981	4410.7292	0.0000	1.0000
K2PI	57	-2143.4165	4411.3784	0.0000	1.0000
JC69I	56	-2165.7434	4453.6523	0.0000	1.0000
JC69G	56	-2166.3750	4454.9157	0.0000	1.0000
JC69IG	57	-2165.8747	4456.2949	0.0000	1.0000
F81	58	-2166.7692	4460.4713	0.0000	1.0000
SYM	60	-2167.5986	4466.9280	0.0000	1.0000
K2P	56	-2182.6469	4487.4594	0.0000	1.0000
JC69	55	-2206.1166	4532.0266	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-2081.8523	4293.7047	0.6288	0.6288
GTRI	64	-2083.6542	4295.3083	0.2820	0.9108
GTRG	64	-2084.8086	4297.6172	0.0889	0.9997
HKYI	60	-2095.0769	4310.1537	0.0002	0.9998
HKYG	60	-2095.5858	4311.1717	0.0001	0.9999
HKYIG	61	-2095.1242	4312.2484	0.0001	1.0000
F81I	59	-2126.3357	4370.6715	0.0000	1.0000
F81G	59	-2126.6287	4371.2574	0.0000	1.0000
F81IG	60	-2126.1584	4372.3168	0.0000	1.0000
GTR	63	-2123.7733	4373.5466	0.0000	1.0000
SYMIG	62	-2129.1540	4382.3079	0.0000	1.0000
SYMI	61	-2131.1941	4384.3881	0.0000	1.0000
SYMG	61	-2131.6424	4385.2848	0.0000	1.0000
HKY	59	-2138.6291	4395.2581	0.0000	1.0000
K2PG	57	-2142.3860	4398.7721	0.0000	1.0000
K2PIG	58	-2141.8981	4399.7963	0.0000	1.0000
K2PI	57	-2143.4165	4400.8329	0.0000	1.0000
JC69I	56	-2165.7434	4443.4867	0.0000	1.0000
JC69G	56	-2166.3750	4444.7501	0.0000	1.0000
JC69IG	57	-2165.8747	4445.7495	0.0000	1.0000
F81	58	-2166.7692	4449.5384	0.0000	1.0000
SYM	60	-2167.5986	4455.1973	0.0000	1.0000
K2P	56	-2182.6469	4477.2938	0.0000	1.0000
JC69	55	-2206.1166	4522.2333	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYI	60	-2095.0769	4581.9189	0.5158	0.5158
HKYG	60	-2095.5858	4582.9368	0.3100	0.8258
GTRI	64	-2083.6542	4585.1911	0.1004	0.9263
GTRG	64	-2084.8086	4587.5000	0.0317	0.9579
GTRIG	65	-2081.8523	4588.1169	0.0233	0.9812
HKYIG	61	-2095.1242	4588.5430	0.0188	1.0000
F81I	59	-2126.3357	4637.9072	0.0000	1.0000
F81G	59	-2126.6287	4638.4932	0.0000	1.0000
F81IG	60	-2126.1584	4644.0819	0.0000	1.0000
K2PG	57	-2142.3860	4656.9490	0.0000	1.0000
GTR	63	-2123.7733	4658.9000	0.0000	1.0000
K2PI	57	-2143.4165	4659.0098	0.0000	1.0000
SYMI	61	-2131.1941	4660.6827	0.0000	1.0000
SYMG	61	-2131.6424	4661.5793	0.0000	1.0000
HKY	59	-2138.6291	4662.4939	0.0000	1.0000
K2PIG	58	-2141.8981	4662.5026	0.0000	1.0000
SYMIG	62	-2129.1540	4663.1319	0.0000	1.0000
JC69I	56	-2165.7434	4697.1342	0.0000	1.0000
JC69G	56	-2166.3750	4698.3975	0.0000	1.0000
JC69IG	57	-2165.8747	4703.9264	0.0000	1.0000
F81	58	-2166.7692	4712.2447	0.0000	1.0000
SYM	60	-2167.5986	4726.9624	0.0000	1.0000
K2P	56	-2182.6469	4730.9413	0.0000	1.0000
JC69	55	-2206.1166	4771.3513	0.0000	1.0000



To specify default values for best AIC, AICc, or BIC models in MrBayes:

[Mrbayes block for the best AIC model (GTRIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_764.phylip = 1 - 685;
 Partition Dummy = 1:183-phylip/chr3_764.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best AICc model (GTRIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_764.phylip = 1 - 685;
 Partition Dummy = 1:183-phylip/chr3_764.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=6 rates=invgamma;
 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) pinvarpr=Uniform(0.0,1.0);
END;


[Mrbayes block for the best BIC model (HKYI)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_764.phylip = 1 - 685;
 Partition Dummy = 1:183-phylip/chr3_764.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

