title: Robust Analysis of High Throughput Screening (HTS) Assay Data.
abstract: Quantitative high throughput screening (qHTS) assays use cells or tissues to screen thousands of compounds in a short period of time. Data generated from qHTS assays are then evaluated using nonlinear regression models, such as the Hill model, and decisions regarding toxicity are made using the estimates of the parameters of the model. For any given compound, the variability in the observed response may either be constant across dose groups (homoscedasticity) or vary with dose (heteroscedasticity). Since thousands of compounds are simultaneously evaluated in a qHTS assay, it is not practically feasible for an investigator to perform residual analysis to determine the variance structure before performing statistical inferences on each compound. Since it is well-known that the variance structure plays an important role in the analysis of linear and nonlinear regression models it is therefore important to have practically useful and easy to interpret methodology which is robust to the variance structure. Furthermore, given the number of chemicals that are investigated in the qHTS assay, outliers and influential observations are not uncommon. In this article we describe preliminary test estimation (PTE) based methodology which is robust to the variance structure as well as any potential outliers and influential observations. Performance of the proposed methodology is evaluated in terms of false discovery rate (FDR) and power using a simulation study mimicking a real qHTS data. Of the two methods currently in use, our simulations studies suggest that one is extremely conservative with very small power in comparison to the proposed PTE based method whereas the other method is very liberal. In contrast, the proposed PTE based methodology achieves a better control of FDR while maintaining good power. The proposed methodology is illustrated using a data set obtained from the National Toxicology Program (NTP). Additional information, simulation results, data and computer code are available online as supplementary materials.
authors.affiliation: Department of Mathematics and Statistics, Loyola University Chicago, 1032 W Sheridan Rd, Chicago, IL 60660.
authors.forename: Changwon
authors.initials: C
authors.lastname: Lim
authors.forename: Pranab K
authors.initials: PK
authors.lastname: Sen
authors.forename: Shyamal D
authors.initials: SD
authors.lastname: Peddada
keywords: Dose-response study; False discovery rate (FDR); Heteroscedasticity; Hill model; M-estimation procedure; Nonlinear regression model; Power; Toxicology
references.citation: J Biomol Screen. 2009 Dec;14(10):1216-27
references.pmid: 19828774
references.citation: Proc Natl Acad Sci U S A. 2006 Aug 1;103(31):11473-8
references.pmid: 16864780
references.citation: Assay Drug Dev Technol. 2008 Oct;6(5):637-57
references.pmid: 19035846
references.citation: Environ Health Perspect. 2008 Mar;116(3):284-91
references.pmid: 18335092
references.citation: J Stat Plan Inference. 2012 May 1;142(5):1047-1062
references.pmid: 22345900
references.citation: J Biomol Screen. 2007 Aug;12(5):645-55
references.pmid: 17517904