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

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

Minimum AIC  model: GTRIG
Minimum AICc model: K2PIG
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
K2PIG	58	-1479.5925	3094.5183	0.3201	0.3201
GTRIG	65	-1470.5869	3095.9001	0.1604	0.4805
HKYIG	61	-1476.2080	3095.9658	0.1552	0.6357
K2PG	57	-1481.7967	3096.2188	0.1368	0.7724
GTRG	64	-1472.3579	3096.6239	0.1117	0.8841
HKYG	60	-1478.2807	3097.3568	0.0774	0.9616
SYMIG	62	-1476.8132	3099.9465	0.0212	0.9828
SYMG	61	-1478.8310	3101.2118	0.0113	0.9940
K2PI	57	-1485.7971	3104.2195	0.0025	0.9965
GTRI	64	-1476.4808	3104.8696	0.0018	0.9984
HKYI	60	-1482.2806	3105.3567	0.0014	0.9998
SYMI	61	-1482.7466	3109.0430	0.0002	1.0000
JC69IG	57	-1503.3624	3139.3502	0.0000	1.0000
JC69G	56	-1505.4197	3140.7719	0.0000	1.0000
F81IG	60	-1500.1822	3141.1599	0.0000	1.0000
F81G	59	-1502.1724	3142.4015	0.0000	1.0000
JC69I	56	-1509.3718	3148.6762	0.0000	1.0000
F81I	59	-1506.1147	3150.2860	0.0000	1.0000
GTR	63	-1510.9620	3171.0300	0.0000	1.0000
HKY	59	-1517.9599	3173.9765	0.0000	1.0000
K2P	56	-1522.4377	3174.8079	0.0000	1.0000
SYM	60	-1518.3618	3177.5190	0.0000	1.0000
JC69	55	-1545.5712	3218.3973	0.0000	1.0000
F81	58	-1541.7839	3218.9011	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-1470.5869	3071.1739	0.5085	0.5085
GTRG	64	-1472.3579	3072.7158	0.2352	0.7437
HKYIG	61	-1476.2080	3074.4160	0.1005	0.8442
K2PIG	58	-1479.5925	3075.1850	0.0684	0.9126
HKYG	60	-1478.2807	3076.5613	0.0344	0.9470
K2PG	57	-1481.7967	3077.5935	0.0205	0.9676
SYMIG	62	-1476.8132	3077.6265	0.0202	0.9877
SYMG	61	-1478.8310	3079.6619	0.0073	0.9950
GTRI	64	-1476.4808	3080.9616	0.0038	0.9988
HKYI	60	-1482.2806	3084.5613	0.0006	0.9995
K2PI	57	-1485.7971	3085.5942	0.0004	0.9999
SYMI	61	-1482.7466	3087.4932	0.0001	1.0000
F81IG	60	-1500.1822	3120.3644	0.0000	1.0000
JC69IG	57	-1503.3624	3120.7249	0.0000	1.0000
F81G	59	-1502.1724	3122.3448	0.0000	1.0000
JC69G	56	-1505.4197	3122.8394	0.0000	1.0000
F81I	59	-1506.1147	3130.2294	0.0000	1.0000
JC69I	56	-1509.3718	3130.7436	0.0000	1.0000
GTR	63	-1510.9620	3147.9240	0.0000	1.0000
HKY	59	-1517.9599	3153.9198	0.0000	1.0000
SYM	60	-1518.3618	3156.7235	0.0000	1.0000
K2P	56	-1522.4377	3156.8753	0.0000	1.0000
F81	58	-1541.7839	3199.5678	0.0000	1.0000
JC69	55	-1545.5712	3201.1424	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
K2PG	57	-1481.7967	3306.9300	0.6802	0.6802
K2PIG	58	-1479.5925	3308.5449	0.3034	0.9836
K2PI	57	-1485.7971	3314.9307	0.0125	0.9960
HKYG	60	-1478.2807	3317.9682	0.0027	0.9988
HKYIG	61	-1476.2080	3319.8463	0.0011	0.9998
SYMG	61	-1478.8310	3325.0922	0.0001	0.9999
HKYI	60	-1482.2806	3325.9681	0.0000	1.0000
SYMIG	62	-1476.8132	3327.0802	0.0000	1.0000
GTRG	64	-1472.3579	3330.2165	0.0000	1.0000
GTRIG	65	-1470.5869	3332.6980	0.0000	1.0000
SYMI	61	-1482.7466	3332.9235	0.0000	1.0000
GTRI	64	-1476.4808	3338.4622	0.0000	1.0000
JC69G	56	-1505.4197	3348.1524	0.0000	1.0000
JC69IG	57	-1503.3624	3350.0614	0.0000	1.0000
JC69I	56	-1509.3718	3356.0567	0.0000	1.0000
F81G	59	-1502.1724	3359.7282	0.0000	1.0000
F81IG	60	-1500.1822	3361.7713	0.0000	1.0000
F81I	59	-1506.1147	3367.6128	0.0000	1.0000
K2P	56	-1522.4377	3382.1884	0.0000	1.0000
HKY	59	-1517.9599	3391.3032	0.0000	1.0000
SYM	60	-1518.3618	3398.1304	0.0000	1.0000
GTR	63	-1510.9620	3401.4012	0.0000	1.0000
JC69	55	-1545.5712	3422.4321	0.0000	1.0000
F81	58	-1541.7839	3432.9277	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_24903.phylip = 1 - 413;
 Partition Dummy = 1:183-phylip/chr3_24903.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 (K2PIG)]

BEGIN MRBAYES;
 Charset 183-phylip/chr3_24903.phylip = 1 - 413;
 Partition Dummy = 1:183-phylip/chr3_24903.phylip;
 Set partition = Dummy;
 Lset applyto=(1) nst=2 rates=invgamma;
 Prset applyto=(1) statefreqpr=Fixed(Equal) tratiopr=Beta(1.0,1.0) shapepr=Uniform(0.1,50.0) pinvarpr=Uniform(0.0,1.0);
END;


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

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

