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

Input data from file "183-phylip/chrZ_3358.phylip" (29 taxa, 738 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	64	-3634.5208	7409.4041	0.6327	0.6327
GTRIG	65	-3633.9428	7410.6535	0.3388	0.9714
HKYG	60	-3642.8375	7416.4873	0.0183	0.9898
HKYIG	61	-3642.2324	7417.6542	0.0102	1.0000
K2PG	57	-3657.7518	7439.2271	0.0000	1.0000
K2PIG	58	-3656.8034	7439.6863	0.0000	1.0000
SYMIG	62	-3653.3399	7442.2532	0.0000	1.0000
SYMG	61	-3655.9137	7445.0167	0.0000	1.0000
GTRI	64	-3662.6865	7465.7355	0.0000	1.0000
HKYI	60	-3668.3870	7467.5865	0.0000	1.0000
K2PI	57	-3678.9183	7481.5602	0.0000	1.0000
SYMI	61	-3677.4465	7488.0824	0.0000	1.0000
GTR	63	-3730.9955	7599.9554	0.0000	1.0000
F81IG	60	-3734.7830	7600.3783	0.0000	1.0000
HKY	59	-3739.0915	7606.6255	0.0000	1.0000
F81G	59	-3739.4362	7607.3149	0.0000	1.0000
JC69IG	57	-3747.9006	7619.5248	0.0000	1.0000
JC69G	56	-3752.0390	7625.4524	0.0000	1.0000
K2P	56	-3758.4773	7638.3290	0.0000	1.0000
SYM	60	-3755.5715	7641.9555	0.0000	1.0000
F81I	59	-3761.4803	7651.4031	0.0000	1.0000
JC69I	56	-3769.9731	7661.3207	0.0000	1.0000
F81	58	-3830.4690	7787.0176	0.0000	1.0000
JC69	55	-3845.7401	7810.5124	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-3634.5208	7397.0416	0.5959	0.5959
GTRIG	65	-3633.9428	7397.8856	0.3908	0.9867
HKYG	60	-3642.8375	7405.6749	0.0080	0.9946
HKYIG	61	-3642.2324	7406.4649	0.0054	1.0000
K2PG	57	-3657.7518	7429.5036	0.0000	1.0000
K2PIG	58	-3656.8034	7429.6068	0.0000	1.0000
SYMIG	62	-3653.3399	7430.6798	0.0000	1.0000
SYMG	61	-3655.9137	7433.8274	0.0000	1.0000
GTRI	64	-3662.6865	7453.3730	0.0000	1.0000
HKYI	60	-3668.3870	7456.7741	0.0000	1.0000
K2PI	57	-3678.9183	7471.8367	0.0000	1.0000
SYMI	61	-3677.4465	7476.8931	0.0000	1.0000
GTR	63	-3730.9955	7587.9910	0.0000	1.0000
F81IG	60	-3734.7830	7589.5659	0.0000	1.0000
HKY	59	-3739.0915	7596.1830	0.0000	1.0000
F81G	59	-3739.4362	7596.8724	0.0000	1.0000
JC69IG	57	-3747.9006	7609.8013	0.0000	1.0000
JC69G	56	-3752.0390	7616.0779	0.0000	1.0000
K2P	56	-3758.4773	7628.9546	0.0000	1.0000
SYM	60	-3755.5715	7631.1431	0.0000	1.0000
F81I	59	-3761.4803	7640.9606	0.0000	1.0000
JC69I	56	-3769.9731	7651.9462	0.0000	1.0000
F81	58	-3830.4690	7776.9380	0.0000	1.0000
JC69	55	-3845.7401	7801.4801	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKYG	60	-3642.8375	7681.9116	0.9236	0.9236
HKYIG	61	-3642.2324	7687.3054	0.0623	0.9858
GTRG	64	-3634.5208	7691.6940	0.0069	0.9928
K2PG	57	-3657.7518	7691.9284	0.0062	0.9990
K2PIG	58	-3656.8034	7696.6355	0.0006	0.9995
GTRIG	65	-3633.9428	7697.1420	0.0005	1.0000
SYMG	61	-3655.9137	7714.6680	0.0000	1.0000
SYMIG	62	-3653.3399	7716.1243	0.0000	1.0000
HKYI	60	-3668.3870	7733.0107	0.0000	1.0000
K2PI	57	-3678.9183	7734.2615	0.0000	1.0000
GTRI	64	-3662.6865	7748.0254	0.0000	1.0000
SYMI	61	-3677.4465	7757.7337	0.0000	1.0000
F81IG	60	-3734.7830	7865.8026	0.0000	1.0000
HKY	59	-3739.0915	7867.8157	0.0000	1.0000
F81G	59	-3739.4362	7868.5051	0.0000	1.0000
JC69IG	57	-3747.9006	7872.2261	0.0000	1.0000
JC69G	56	-3752.0390	7873.8988	0.0000	1.0000
GTR	63	-3730.9955	7878.0394	0.0000	1.0000
K2P	56	-3758.4773	7886.7754	0.0000	1.0000
SYM	60	-3755.5715	7907.3797	0.0000	1.0000
JC69I	56	-3769.9731	7909.7671	0.0000	1.0000
F81I	59	-3761.4803	7912.5933	0.0000	1.0000
F81	58	-3830.4690	8043.9668	0.0000	1.0000
JC69	55	-3845.7401	8054.6971	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 183-phylip/chrZ_3358.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chrZ_3358.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 183-phylip/chrZ_3358.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chrZ_3358.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 183-phylip/chrZ_3358.phylip = 1 - 738;
 Partition Dummy = 1:183-phylip/chrZ_3358.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

