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

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

Minimum AIC  model: GTRG
Minimum AICc model: K2PG
Minimum BIC  model: K2PG


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
K2PG	57	-1246.3004	2629.8826	0.3826	0.3826
HKYG	60	-1242.5708	2631.1914	0.1989	0.5814
K2PIG	58	-1245.5811	2631.3459	0.1841	0.7655
HKYIG	61	-1241.9155	2632.8453	0.0870	0.8525
K2PI	57	-1248.3469	2633.9754	0.0494	0.9019
GTRG	64	-1238.1908	2634.4177	0.0396	0.9415
HKYI	60	-1244.3771	2634.8040	0.0327	0.9742
GTRIG	65	-1237.3946	2635.8762	0.0191	0.9933
GTRI	64	-1240.0806	2638.1973	0.0060	0.9993
SYMG	61	-1247.2213	2643.4569	0.0004	0.9997
SYMIG	62	-1246.5149	2645.0299	0.0002	0.9999
SYMI	61	-1248.9335	2646.8812	0.0001	1.0000
K2P	56	-1269.7076	2673.8152	0.0000	1.0000
HKY	59	-1266.5561	2676.2185	0.0000	1.0000
GTR	63	-1262.8432	2680.6935	0.0000	1.0000
SYM	60	-1271.3929	2688.8356	0.0000	1.0000
F81G	59	-1286.6560	2716.4183	0.0000	1.0000
F81IG	60	-1286.0035	2718.0568	0.0000	1.0000
JC69G	56	-1292.3454	2719.0907	0.0000	1.0000
F81I	59	-1288.4183	2719.9429	0.0000	1.0000
JC69IG	57	-1291.6497	2720.5812	0.0000	1.0000
JC69I	56	-1294.0958	2722.5916	0.0000	1.0000
F81	58	-1309.7505	2759.6848	0.0000	1.0000
JC69	55	-1315.6703	2762.8790	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRG	64	-1238.1908	2604.3816	0.2582	0.2582
GTRIG	65	-1237.3946	2604.7893	0.2106	0.4687
HKYG	60	-1242.5708	2605.1415	0.1766	0.6453
HKYIG	61	-1241.9155	2605.8310	0.1251	0.7704
K2PG	57	-1246.3004	2606.6009	0.0851	0.8555
K2PIG	58	-1245.5811	2607.1622	0.0643	0.9198
GTRI	64	-1240.0806	2608.1612	0.0390	0.9588
HKYI	60	-1244.3771	2608.7542	0.0290	0.9878
K2PI	57	-1248.3469	2610.6937	0.0110	0.9988
SYMG	61	-1247.2213	2616.4426	0.0006	0.9994
SYMIG	62	-1246.5149	2617.0299	0.0005	0.9999
SYMI	61	-1248.9335	2619.8670	0.0001	1.0000
HKY	59	-1266.5561	2651.1121	0.0000	1.0000
K2P	56	-1269.7076	2651.4152	0.0000	1.0000
GTR	63	-1262.8432	2651.6863	0.0000	1.0000
SYM	60	-1271.3929	2662.7858	0.0000	1.0000
F81G	59	-1286.6560	2691.3119	0.0000	1.0000
F81IG	60	-1286.0035	2692.0070	0.0000	1.0000
F81I	59	-1288.4183	2694.8365	0.0000	1.0000
JC69G	56	-1292.3454	2696.6907	0.0000	1.0000
JC69IG	57	-1291.6497	2697.2995	0.0000	1.0000
JC69I	56	-1294.0958	2700.1916	0.0000	1.0000
F81	58	-1309.7505	2735.5011	0.0000	1.0000
JC69	55	-1315.6703	2741.3406	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1246.3004	2825.1851	0.8009	0.8009
K2PI	57	-1248.3469	2829.2780	0.1035	0.9044
K2PIG	58	-1245.5811	2829.5812	0.0889	0.9933
HKYG	60	-1242.5708	2835.2302	0.0053	0.9986
HKYI	60	-1244.3771	2838.8429	0.0009	0.9994
HKYIG	61	-1241.9155	2839.7545	0.0005	1.0000
GTRG	64	-1238.1908	2849.8095	0.0000	1.0000
SYMG	61	-1247.2213	2850.3661	0.0000	1.0000
GTRI	64	-1240.0806	2853.5891	0.0000	1.0000
SYMI	61	-1248.9335	2853.7904	0.0000	1.0000
GTRIG	65	-1237.3946	2854.0520	0.0000	1.0000
SYMIG	62	-1246.5149	2854.7881	0.0000	1.0000
K2P	56	-1269.7076	2866.1646	0.0000	1.0000
HKY	59	-1266.5561	2877.3660	0.0000	1.0000
SYM	60	-1271.3929	2892.8745	0.0000	1.0000
GTR	63	-1262.8432	2893.2794	0.0000	1.0000
JC69G	56	-1292.3454	2911.4401	0.0000	1.0000
JC69I	56	-1294.0958	2914.9410	0.0000	1.0000
JC69IG	57	-1291.6497	2915.8837	0.0000	1.0000
F81G	59	-1286.6560	2917.5657	0.0000	1.0000
F81I	59	-1288.4183	2921.0904	0.0000	1.0000
F81IG	60	-1286.0035	2922.0957	0.0000	1.0000
JC69	55	-1315.6703	2952.2552	0.0000	1.0000
F81	58	-1309.7505	2957.9201	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/chr8_4947.phylip = 1 - 342;
 Partition Dummy = 1:183-phylip/chr8_4947.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 (K2PG)]

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


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

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

