Search

Background Most cost-effectiveness analyses in the context of cervical cancer prevention involve the use of mathematical models to simulate HPV infection, cervical disease and prevention strategies. However, it is common for professionals who would need to perform these analyses to not be familiar with the models. This work introduces the Online Cost-Effectiveness ANalysis tool, featuring an easy-to-use web interface providing health professionals, researchers and decision makers involved in cervical cancer prevention programmes with a useful instrument to conduct complex cost-effectiveness analyses, which are becoming an essential tool as an approach for supporting decision-making that involves important trade-offs. Results The users can run cost-effectiveness evaluations of cervical cancer prevention strategies without deep knowledge of the underlying mathematical model or any programming language, obtaining the most relevant costs and health outcomes in a user-friendly format. The results provided by the tool are consistent with the existing literature. Conclusions Having such a tool will be an asset to the cervical cancer prevention community, providing researchers with an easy-to-use instrument to conduct cost-effectiveness analyses.

Background High-throughput RNA sequencing (RNA-seq) offers unprecedented power to capture the real dynamics of gene expression. Experimental designs with extensive biological replication present a unique opportunity to exploit this feature and distinguish expression profiles with higher resolution. RNA-seq data analysis methods so far have been mostly applied to data sets with few replicates and their default settings try to provide the best performance under this constraint. These methods are based on two well-known count data distributions: the Poisson and the negative binomial. The way to properly calibrate them with large RNA-seq data sets is not trivial for the non-expert bioinformatics user. Results Here we show that expression profiles produced by extensively-replicated RNA-seq experiments lead to a rich diversity of count data distributions beyond the Poisson and the negative binomial, such as Poisson-Inverse Gaussian or Pólya-Aeppli, which can be captured by a more general family of count data distributions called the Poisson-Tweedie. The flexibility of the Poisson-Tweedie family enables a direct fitting of emerging features of large expression profiles, such as heavy-tails or zero-inflation, without the need to alter a single configuration parameter. We provide a software package for R called implementing a new test for differential expression based on the Poisson-Tweedie family. Using simulations on synthetic and real RNA-seq data we show that yields P-values that are equally or more accurate than competing methods under different configuration parameters. By surveying the tiny fraction of sex-specific gene expression changes in human lymphoblastoid cell lines, we also show that accurately detects differentially expressed genes in a real large RNA-seq data set with improved performance and reproducibility over the previously compared methodologies. Finally, we compared the results with those obtained from microarrays in order to check for reproducibility. Conclusions RNA-seq data with many replicates leads to a handful of count data distributions which can be accurately estimated with the statistical model illustrated in this paper. This method provides a better fit to the underlying biological variability; this may be critical when comparing groups of RNA-seq samples with markedly different count data distributions. The package forms part of the Bioconductor project and it is available for download at http://www.bioconductor.org.