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

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

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


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	-3264.8658	6672.5567	0.9756	0.9756
GTRG	64	-3269.7567	6679.9314	0.0244	1.0000
HKYG	60	-3282.4777	6695.8160	0.0000	1.0000
HKYIG	61	-3284.6961	6702.6315	0.0000	1.0000
SYMIG	62	-3297.2406	6730.1061	0.0000	1.0000
K2PIG	58	-3302.3873	6730.8988	0.0000	1.0000
SYMG	61	-3300.1770	6733.5931	0.0000	1.0000
GTRI	64	-3297.8825	6736.1829	0.0000	1.0000
K2PG	57	-3306.9610	6737.6886	0.0000	1.0000
HKYI	60	-3313.7707	6758.4020	0.0000	1.0000
SYMI	61	-3323.8845	6781.0082	0.0000	1.0000
K2PI	57	-3331.5904	6786.9474	0.0000	1.0000
GTR	63	-3400.0655	6938.1490	0.0000	1.0000
HKY	59	-3417.9048	6964.2986	0.0000	1.0000
F81G	59	-3428.6622	6985.8134	0.0000	1.0000
F81IG	60	-3429.1137	6989.0879	0.0000	1.0000
JC69IG	57	-3437.0912	6997.9489	0.0000	1.0000
SYM	60	-3437.2160	7005.2925	0.0000	1.0000
JC69G	56	-3444.9371	7011.2902	0.0000	1.0000
K2P	56	-3446.2319	7013.8797	0.0000	1.0000
F81I	59	-3454.8621	7038.2131	0.0000	1.0000
JC69I	56	-3466.2372	7053.8903	0.0000	1.0000
F81	58	-3554.1936	7234.5115	0.0000	1.0000
JC69	55	-3573.5622	7266.1965	0.0000	1.0000



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

Model	df	lnL		AIC		wAIC	cwAIC
GTRIG	65	-3264.8658	6659.7316	0.9800	0.9800
GTRG	64	-3269.7567	6667.5135	0.0200	1.0000
HKYG	60	-3282.4777	6684.9555	0.0000	1.0000
HKYIG	61	-3284.6961	6691.3923	0.0000	1.0000
SYMIG	62	-3297.2406	6718.4811	0.0000	1.0000
K2PIG	58	-3302.3873	6720.7745	0.0000	1.0000
SYMG	61	-3300.1770	6722.3539	0.0000	1.0000
GTRI	64	-3297.8825	6723.7650	0.0000	1.0000
K2PG	57	-3306.9610	6727.9220	0.0000	1.0000
HKYI	60	-3313.7707	6747.5415	0.0000	1.0000
SYMI	61	-3323.8845	6769.7690	0.0000	1.0000
K2PI	57	-3331.5904	6777.1807	0.0000	1.0000
GTR	63	-3400.0655	6926.1311	0.0000	1.0000
HKY	59	-3417.9048	6953.8097	0.0000	1.0000
F81G	59	-3428.6622	6975.3245	0.0000	1.0000
F81IG	60	-3429.1137	6978.2274	0.0000	1.0000
JC69IG	57	-3437.0912	6988.1823	0.0000	1.0000
SYM	60	-3437.2160	6994.4320	0.0000	1.0000
JC69G	56	-3444.9371	7001.8743	0.0000	1.0000
K2P	56	-3446.2319	7004.4638	0.0000	1.0000
F81I	59	-3454.8621	7027.7242	0.0000	1.0000
JC69I	56	-3466.2372	7044.4744	0.0000	1.0000
F81	58	-3554.1936	7224.3872	0.0000	1.0000
JC69	55	-3573.5622	7257.1243	0.0000	1.0000



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

Model	df	lnL		BIC		wBIC	cwBIC
GTRIG	65	-3264.8658	6958.7232	0.6519	0.6519
HKYG	60	-3282.4777	6960.9477	0.2144	0.8663
GTRG	64	-3269.7567	6961.9052	0.1328	0.9991
HKYIG	61	-3284.6961	6971.9844	0.0009	1.0000
K2PIG	58	-3302.3873	6987.5670	0.0000	1.0000
K2PG	57	-3306.9610	6990.1146	0.0000	1.0000
SYMG	61	-3300.1770	7002.9460	0.0000	1.0000
SYMIG	62	-3297.2406	7003.6731	0.0000	1.0000
GTRI	64	-3297.8825	7018.1567	0.0000	1.0000
HKYI	60	-3313.7707	7023.5337	0.0000	1.0000
K2PI	57	-3331.5904	7039.3734	0.0000	1.0000
SYMI	61	-3323.8845	7050.3611	0.0000	1.0000
GTR	63	-3400.0655	7215.9229	0.0000	1.0000
HKY	59	-3417.9048	7225.2020	0.0000	1.0000
F81G	59	-3428.6622	7246.7169	0.0000	1.0000
JC69IG	57	-3437.0912	7250.3749	0.0000	1.0000
F81IG	60	-3429.1137	7254.2196	0.0000	1.0000
JC69G	56	-3444.9371	7259.4670	0.0000	1.0000
K2P	56	-3446.2319	7262.0565	0.0000	1.0000
SYM	60	-3437.2160	7270.4242	0.0000	1.0000
F81I	59	-3454.8621	7299.1165	0.0000	1.0000
JC69I	56	-3466.2372	7302.0671	0.0000	1.0000
F81	58	-3554.1936	7491.1797	0.0000	1.0000
JC69	55	-3573.5622	7510.1172	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/chr13_2821.phylip = 1 - 735;
 Partition Dummy = 1:183-phylip/chr13_2821.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/chr13_2821.phylip = 1 - 735;
 Partition Dummy = 1:183-phylip/chr13_2821.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 (GTRIG)]

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


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

