Ordering of Omics Features Using Beta Distributions on Montecarlo p -Values

Angela L. Riffo-Campos, Guillermo Ayala and Juan Domingo
Additional contact information
Angela L. Riffo-Campos: Centro de Excelencia de Modelación y Computación Científica, Universidad de La Frontera, Temuco 01145, Chile
Guillermo Ayala: Department of Statistics and Operation Research, Faculty of Mathematics, Universitat de Valencia, 46100 Burjasot, Spain
Juan Domingo: Department of Computer Science, ETSE, Universitat de Valencia, Avda. de la Universidad, s/n, 46100 Burjasot, Spain

Mathematics, 2021, vol. 9, issue 11, 1-18

Abstract: The current trend in genetic research is the study of omics data as a whole, either combining studies or omics techniques. This raises the need for new robust statistical methods that can integrate and order the relevant biological information. A good way to approach the problem is to order the features studied according to the different kinds of data so a key point is to associate good values to the features that permit us a good sorting of them. These values are usually the p -values corresponding to a hypothesis which has been tested for each feature studied. The Montecarlo method is certainly one of the most robust methods for hypothesis testing. However, a large number of simulations is needed to obtain a reliable p -value, so the method becomes computationally infeasible in many situations. We propose a new way to order genes according to their differential features by using a score defined from a beta distribution fitted to the generated p -values. Our approach has been tested using simulated data and colorectal cancer datasets from Infinium methylationEPIC array, Affymetrix gene expression array and Illumina RNA-seq platforms. The results show that this approach allows a proper ordering of genes using a number of simulations much lower than with the Montecarlo method. Furthermore, the score can be interpreted as an estimated p -value and compared with Montecarlo and other approaches like the p -value of the moderated t -tests. We have also identified a new expression pattern of eighteen genes common to all colorectal cancer microarrays, i.e., 21 datasets. Thus, the proposed method is effective for obtaining biological results using different datasets. Our score shows a slightly smaller type I error for small sizes than the Montecarlo p -value. The type II error of Montecarlo p -value is lower than the one obtained with the proposed score and with a moderated p -value, but these differences are highly reduced for larger sample sizes and higher false discovery rates. Similar performances from type I and II errors and the score enable a clear ordering of the features being evaluated.

Keywords: beta distribution; gene ordering; colorectal cancer; Montecarlo method; moderated t -test (search for similar items in EconPapers)
JEL-codes: C (search for similar items in EconPapers)
Date: 2021
References: View references in EconPapers View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/2227-7390/9/11/1307/pdf (application/pdf)
https://www.mdpi.com/2227-7390/9/11/1307/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jmathe:v:9:y:2021:i:11:p:1307-:d:570315

Access Statistics for this article

Mathematics is currently edited by Ms. Emma He

More articles in Mathematics from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().