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

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

Minimum AIC  model: GTRIG
Minimum AICc model: GTRIG
Minimum BIC  model: GTRG


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	-1757.2424	3678.8049	0.6230	0.6230
GTRG	64	-1759.3388	3679.8251	0.3741	0.9971
HKYIG	61	-1769.4127	3690.6050	0.0017	0.9988
HKYG	60	-1771.3819	3691.4696	0.0011	0.9999
K2PIG	58	-1777.4799	3697.5901	0.0001	0.9999
SYMIG	62	-1771.9300	3698.7375	0.0000	1.0000
K2PG	57	-1779.6970	3699.0219	0.0000	1.0000
SYMG	61	-1774.4625	3700.7046	0.0000	1.0000
GTRI	64	-1772.8588	3706.8651	0.0000	1.0000
HKYI	60	-1788.5963	3725.8984	0.0000	1.0000
SYMI	61	-1788.4090	3728.5976	0.0000	1.0000
K2PI	57	-1797.5737	3734.7753	0.0000	1.0000
F81G	59	-1826.4560	3798.5682	0.0000	1.0000
JC69G	56	-1831.3825	3799.4135	0.0000	1.0000
F81IG	60	-1826.9546	3802.6150	0.0000	1.0000
JC69IG	57	-1831.7086	3803.0451	0.0000	1.0000
F81I	59	-1845.7069	3837.0700	0.0000	1.0000
JC69I	56	-1852.9423	3842.5332	0.0000	1.0000
GTR	63	-1852.3712	3862.7425	0.0000	1.0000
HKY	59	-1859.7685	3865.1933	0.0000	1.0000
SYM	60	-1861.1348	3870.9756	0.0000	1.0000
K2P	56	-1867.9344	3872.5175	0.0000	1.0000
JC69	55	-1919.3287	3972.3497	0.0000	1.0000
F81	58	-1915.9116	3974.4535	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1757.2424	3644.4849	0.7494	0.7494
GTRG	64	-1759.3388	3646.6776	0.2503	0.9997
HKYIG	61	-1769.4127	3660.8254	0.0002	0.9999
HKYG	60	-1771.3819	3662.7637	0.0001	1.0000
SYMIG	62	-1771.9300	3667.8601	0.0000	1.0000
SYMG	61	-1774.4625	3670.9250	0.0000	1.0000
K2PIG	58	-1777.4799	3670.9598	0.0000	1.0000
K2PG	57	-1779.6970	3673.3940	0.0000	1.0000
GTRI	64	-1772.8588	3673.7176	0.0000	1.0000
HKYI	60	-1788.5963	3697.1925	0.0000	1.0000
SYMI	61	-1788.4090	3698.8181	0.0000	1.0000
K2PI	57	-1797.5737	3709.1474	0.0000	1.0000
F81G	59	-1826.4560	3770.9120	0.0000	1.0000
F81IG	60	-1826.9546	3773.9091	0.0000	1.0000
JC69G	56	-1831.3825	3774.7649	0.0000	1.0000
JC69IG	57	-1831.7086	3777.4172	0.0000	1.0000
F81I	59	-1845.7069	3809.4137	0.0000	1.0000
JC69I	56	-1852.9423	3817.8845	0.0000	1.0000
GTR	63	-1852.3712	3830.7425	0.0000	1.0000
HKY	59	-1859.7685	3837.5371	0.0000	1.0000
SYM	60	-1861.1348	3842.2697	0.0000	1.0000
K2P	56	-1867.9344	3847.8689	0.0000	1.0000
F81	58	-1915.9116	3947.8231	0.0000	1.0000
JC69	55	-1919.3287	3948.6574	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRG	64	-1759.3388	3887.0451	0.2849	0.2849
K2PG	57	-1779.6970	3887.4713	0.2302	0.5151
HKYG	60	-1771.3819	3888.1083	0.1674	0.6825
GTRIG	65	-1757.2424	3888.6081	0.1304	0.8129
K2PIG	58	-1777.4799	3888.7928	0.1189	0.9318
HKYIG	61	-1769.4127	3889.9257	0.0675	0.9993
SYMG	61	-1774.4625	3900.0253	0.0004	0.9997
SYMIG	62	-1771.9300	3900.7161	0.0003	1.0000
GTRI	64	-1772.8588	3914.0851	0.0000	1.0000
HKYI	60	-1788.5963	3922.5371	0.0000	1.0000
K2PI	57	-1797.5737	3923.2247	0.0000	1.0000
SYMI	61	-1788.4090	3927.9184	0.0000	1.0000
JC69G	56	-1831.3825	3985.0865	0.0000	1.0000
JC69IG	57	-1831.7086	3991.4945	0.0000	1.0000
F81G	59	-1826.4560	3992.5008	0.0000	1.0000
F81IG	60	-1826.9546	3999.2537	0.0000	1.0000
JC69I	56	-1852.9423	4028.2061	0.0000	1.0000
F81I	59	-1845.7069	4031.0025	0.0000	1.0000
K2P	56	-1867.9344	4058.1904	0.0000	1.0000
HKY	59	-1859.7685	4059.1259	0.0000	1.0000
GTR	63	-1852.3712	4067.3542	0.0000	1.0000
SYM	60	-1861.1348	4067.6142	0.0000	1.0000
JC69	55	-1919.3287	4155.2232	0.0000	1.0000
F81	58	-1915.9116	4165.6562	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_5924.phylip = 1 - 316;
 Partition Dummy = 1:183-phylip/chr3_5924.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_5924.phylip = 1 - 316;
 Partition Dummy = 1:183-phylip/chr3_5924.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 (GTRG)]

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


-------------------------------------------------------------
End of Output

