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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
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
HKYG	60	-1243.0236	2632.4731	0.4592	0.4592
HKYIG	61	-1242.5662	2634.5382	0.1635	0.6227
K2PG	57	-1248.4965	2634.6072	0.1580	0.7806
HKYI	60	-1245.0312	2636.4885	0.0617	0.8423
K2PIG	58	-1248.0403	2636.6110	0.0580	0.9003
GTRG	64	-1239.4840	2637.4442	0.0382	0.9385
K2PI	57	-1250.3673	2638.3489	0.0243	0.9628
SYMG	61	-1245.1101	2639.6260	0.0128	0.9757
GTRIG	65	-1239.1042	2639.7525	0.0121	0.9877
GTRI	64	-1241.3902	2641.2566	0.0057	0.9934
SYMIG	62	-1244.7139	2641.8351	0.0043	0.9977
SYMI	61	-1246.8186	2643.0431	0.0023	1.0000
F81G	59	-1264.0051	2671.4778	0.0000	1.0000
JC69G	56	-1269.3386	2673.3961	0.0000	1.0000
F81IG	60	-1263.5905	2673.6070	0.0000	1.0000
HKY	59	-1265.7910	2675.0497	0.0000	1.0000
JC69IG	57	-1268.8476	2675.3094	0.0000	1.0000
F81I	59	-1266.0194	2675.5063	0.0000	1.0000
K2P	56	-1271.1626	2677.0440	0.0000	1.0000
JC69I	56	-1271.1911	2677.1011	0.0000	1.0000
GTR	63	-1262.6418	2680.7142	0.0000	1.0000
SYM	60	-1267.9783	2682.3826	0.0000	1.0000
F81	58	-1286.2566	2713.0437	0.0000	1.0000
JC69	55	-1291.2542	2714.3525	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1243.0236	2606.0471	0.3328	0.3328
GTRG	64	-1239.4840	2606.9680	0.2100	0.5428
HKYIG	61	-1242.5662	2607.1324	0.1934	0.7362
GTRIG	65	-1239.1042	2608.2084	0.1129	0.8492
HKYI	60	-1245.0312	2610.0625	0.0447	0.8938
GTRI	64	-1241.3902	2610.7804	0.0312	0.9251
K2PG	57	-1248.4965	2610.9930	0.0281	0.9531
K2PIG	58	-1248.0403	2612.0805	0.0163	0.9694
SYMG	61	-1245.1101	2612.2202	0.0152	0.9846
SYMIG	62	-1244.7139	2613.4278	0.0083	0.9929
K2PI	57	-1250.3673	2614.7346	0.0043	0.9972
SYMI	61	-1246.8186	2615.6373	0.0028	1.0000
F81G	59	-1264.0051	2646.0102	0.0000	1.0000
F81IG	60	-1263.5905	2647.1810	0.0000	1.0000
HKY	59	-1265.7910	2649.5820	0.0000	1.0000
F81I	59	-1266.0194	2650.0387	0.0000	1.0000
JC69G	56	-1269.3386	2650.6772	0.0000	1.0000
GTR	63	-1262.6418	2651.2836	0.0000	1.0000
JC69IG	57	-1268.8476	2651.6952	0.0000	1.0000
K2P	56	-1271.1626	2654.3252	0.0000	1.0000
JC69I	56	-1271.1911	2654.3822	0.0000	1.0000
SYM	60	-1267.9783	2655.9566	0.0000	1.0000
F81	58	-1286.2566	2688.5133	0.0000	1.0000
JC69	55	-1291.2542	2692.5085	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1248.4965	2828.9066	0.7770	0.7770
K2PI	57	-1250.3673	2832.6482	0.1197	0.8967
K2PIG	58	-1248.0403	2833.8172	0.0667	0.9634
HKYG	60	-1243.0236	2835.4299	0.0298	0.9932
HKYI	60	-1245.0312	2839.4452	0.0040	0.9972
HKYIG	61	-1242.5662	2840.3382	0.0026	0.9997
SYMG	61	-1245.1101	2845.4260	0.0002	0.9999
SYMI	61	-1246.8186	2848.8431	0.0000	1.0000
SYMIG	62	-1244.7139	2850.4566	0.0000	1.0000
GTRG	64	-1239.4840	2851.6430	0.0000	1.0000
GTRI	64	-1241.3902	2855.4553	0.0000	1.0000
GTRIG	65	-1239.1042	2856.7064	0.0000	1.0000
JC69G	56	-1269.3386	2864.7678	0.0000	1.0000
K2P	56	-1271.1626	2868.4157	0.0000	1.0000
JC69I	56	-1271.1911	2868.4728	0.0000	1.0000
JC69IG	57	-1268.8476	2869.6088	0.0000	1.0000
F81G	59	-1264.0051	2871.5699	0.0000	1.0000
HKY	59	-1265.7910	2875.1417	0.0000	1.0000
F81I	59	-1266.0194	2875.5984	0.0000	1.0000
F81IG	60	-1263.5905	2876.5638	0.0000	1.0000
SYM	60	-1267.9783	2885.3394	0.0000	1.0000
GTR	63	-1262.6418	2892.1355	0.0000	1.0000
JC69	55	-1291.2542	2902.7760	0.0000	1.0000
F81	58	-1286.2566	2910.2499	0.0000	1.0000



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

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

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


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

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


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

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

