List taxa names present in the BEAST input file

Description

The LOOCV (Leave-one-out-cross-validation) is a statistical procedure aiming to detect whether some particular sequences, when used as calibration in a tip dating analyse, could lead to systematic bias. Once the XMLs produced with the "TaxaOut" function analyzed with BEAST, the PlotLOOCV function allows a graphical comparison of the estimated vs true sequence ages. Similarly to the PlotDRT function, PlotLOOCV requires uploading, reading and analysing the "Log" files generated by BEAST.

Usage

1
PlotLOOCV(name, burnin = 0.1)

Arguments

name

the name of the original XML file (excluding the .xml extension). In the same folder, should be present the i LogFiles generated from the i XMLs produced with the "TaxaOut" function (i being the number of Taxa in the dataset). The name of the Log files should look like "name.Taxon[i].log".

burnin

the fraction of the first MCMC sampling to exclude from the Log files when computing the parameter estimates distribution (default = 0.1, which means 10

Details

The function works only with a .xml file generated with BEAUti

Value

The function returns the names and order of the taxa present in an input xml file for the BEAST software.

References

Rieux A & Khatchikian, C. Unpublished. Drummond AJ, Suchard MA, Xie D & Rambaut A (2012) Bayesian phylogenetics with BEAUti and the BEAST 1.7. Molecular Biology And Evolution 29: 1969-1973.

Examples

1
2
3
4
5
6
## Not run: 
	# example create with example file
	PlotLOOCV("Flu_BEAST_1.8", burnin = 0.1)
	# produce LOOCV plot 

## End(Not run)