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

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

Minimum AIC  model: HKYG
Minimum AICc model: HKYG
Minimum BIC  model: HKY


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	-1352.1632	2840.4142	0.3651	0.3651
HKYI	60	-1352.4686	2841.0252	0.2690	0.6341
HKY	59	-1353.8030	2841.1324	0.2550	0.8891
HKYIG	61	-1352.1860	2843.0327	0.0986	0.9877
GTRG	64	-1351.3436	2849.1351	0.0047	0.9923
GTRI	64	-1351.6584	2849.7646	0.0034	0.9958
GTR	63	-1353.0860	2850.0128	0.0030	0.9988
GTRIG	65	-1351.3587	2851.7841	0.0012	1.0000
K2PG	57	-1384.9601	2898.3569	0.0000	1.0000
K2P	56	-1386.4124	2898.7333	0.0000	1.0000
K2PI	57	-1385.2161	2898.8689	0.0000	1.0000
K2PIG	58	-1384.9734	2900.9227	0.0000	1.0000
SYMG	61	-1384.4459	2907.5527	0.0000	1.0000
SYM	60	-1385.8411	2907.7701	0.0000	1.0000
SYMI	61	-1384.7118	2908.0844	0.0000	1.0000
SYMIG	62	-1384.4508	2910.1466	0.0000	1.0000
F81	58	-1395.1015	2921.1788	0.0000	1.0000
F81G	59	-1394.2823	2922.0910	0.0000	1.0000
F81I	59	-1394.5340	2922.5942	0.0000	1.0000
F81IG	60	-1394.2818	2924.6514	0.0000	1.0000
JC69	55	-1422.2724	2967.9360	0.0000	1.0000
JC69G	56	-1421.4185	2968.7455	0.0000	1.0000
JC69I	56	-1421.6695	2969.2474	0.0000	1.0000
JC69IG	57	-1421.4174	2971.2714	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
HKYG	60	-1352.1632	2824.3263	0.3662	0.3662
HKYI	60	-1352.4686	2824.9373	0.2698	0.6361
HKY	59	-1353.8030	2825.6060	0.1931	0.8292
HKYIG	61	-1352.1860	2826.3719	0.1317	0.9609
GTRG	64	-1351.3436	2830.6872	0.0152	0.9761
GTRI	64	-1351.6584	2831.3167	0.0111	0.9872
GTR	63	-1353.0860	2832.1721	0.0072	0.9945
GTRIG	65	-1351.3587	2832.7174	0.0055	1.0000
K2PG	57	-1384.9601	2883.9202	0.0000	1.0000
K2PI	57	-1385.2161	2884.4322	0.0000	1.0000
K2P	56	-1386.4124	2884.8248	0.0000	1.0000
K2PIG	58	-1384.9734	2885.9468	0.0000	1.0000
SYMG	61	-1384.4459	2890.8919	0.0000	1.0000
SYMI	61	-1384.7118	2891.4236	0.0000	1.0000
SYM	60	-1385.8411	2891.6822	0.0000	1.0000
SYMIG	62	-1384.4508	2892.9016	0.0000	1.0000
F81	58	-1395.1015	2906.2029	0.0000	1.0000
F81G	59	-1394.2823	2906.5647	0.0000	1.0000
F81I	59	-1394.5340	2907.0679	0.0000	1.0000
F81IG	60	-1394.2818	2908.5635	0.0000	1.0000
JC69	55	-1422.2724	2954.5447	0.0000	1.0000
JC69G	56	-1421.4185	2954.8370	0.0000	1.0000
JC69I	56	-1421.6695	2955.3389	0.0000	1.0000
JC69IG	57	-1421.4174	2956.8348	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
HKY	59	-1353.8030	3076.1263	0.7123	0.7123
HKYG	60	-1352.1632	3079.0927	0.1616	0.8740
HKYI	60	-1352.4686	3079.7037	0.1191	0.9930
HKYIG	61	-1352.1860	3085.3845	0.0070	1.0000
GTR	63	-1353.0860	3099.6768	0.0000	1.0000
GTRG	64	-1351.3436	3102.4381	0.0000	1.0000
GTRI	64	-1351.6584	3103.0676	0.0000	1.0000
GTRIG	65	-1351.3587	3108.7143	0.0000	1.0000
K2P	56	-1386.4124	3122.6067	0.0000	1.0000
K2PG	57	-1384.9601	3125.9483	0.0000	1.0000
K2PI	57	-1385.2161	3126.4603	0.0000	1.0000
K2PIG	58	-1384.9734	3132.2210	0.0000	1.0000
SYM	60	-1385.8411	3146.4486	0.0000	1.0000
SYMG	61	-1384.4459	3149.9044	0.0000	1.0000
SYMI	61	-1384.7118	3150.4362	0.0000	1.0000
F81	58	-1395.1015	3152.4771	0.0000	1.0000
SYMIG	62	-1384.4508	3156.1602	0.0000	1.0000
F81G	59	-1394.2823	3157.0850	0.0000	1.0000
F81I	59	-1394.5340	3157.5882	0.0000	1.0000
F81IG	60	-1394.2818	3163.3299	0.0000	1.0000
JC69	55	-1422.2724	3188.0806	0.0000	1.0000
JC69G	56	-1421.4185	3192.6189	0.0000	1.0000
JC69I	56	-1421.6695	3193.1209	0.0000	1.0000
JC69IG	57	-1421.4174	3198.8628	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/chr2_8698.phylip = 1 - 516;
 Partition Dummy = 1:183-phylip/chr2_8698.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/chr2_8698.phylip = 1 - 516;
 Partition Dummy = 1:183-phylip/chr2_8698.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 (HKY)]

BEGIN MRBAYES;
 Charset 183-phylip/chr2_8698.phylip = 1 - 516;
 Partition Dummy = 1:183-phylip/chr2_8698.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=equal;
 Prset applyto=(1) statefreqpr=Dirichlet(1.0,1.0,1.0,1.0) tratiopr=Beta(1.0,1.0);
END;


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

