Type
Text
Type
Dissertation
Advisor
Wu, Song | Zhang, Michael Q., Xing, Haipeng | Fang, Yixin.
Date
2012-12-01
Keywords
Statistics
Department
Department of Applied Mathematics and Statistics
Language
en_US
Source
This work is sponsored by the Stony Brook University Graduate School in compliance with the requirements for completion of degree.
Identifier
http://hdl.handle.net/11401/71350
Publisher
The Graduate School, Stony Brook University: Stony Brook, NY.
Format
application/pdf
Abstract
Nowadays, Hidden Markov Model (HMM) has been widely used in analysis of various biological data for both smoothing and clustering. However, characterizing each hidden state by a single distribution, the classical HMM might have some limitations on the data whose hidden state is composed by a mixture of distributions (Heng Lian et al. | 2006). To address this issue, we proposed a new stochastic segmentation model and an associated estimation procedure that has attractive analytical and computational properties. We combined the forward and backward filter together based on Bayes' theorem to calculate the posterior mean and variance. Besides, we developed an expectation-maximization (EM) algorithm to estimate the hyper-parameters. Furthermore, we utilized a bounded complexity mixture (BCMIX) approximation whose computational complexity is linear in sequence length. Another important feature of this segmentation model is that it yields explicit formulas for posterior means and probability of categorical states, which can be used to make inference on both categorical and continuous aspects of the data. Other quantities relating to the posterior distribution that are useful for making confidence assessments of any given segmentation can also be estimated by using our method. We perform intensive simulation studies (1) to compare the Bayes and BCMIX estimates (2) to evaluate the BCMIX estimates in terms of sum square error, Kullback-Leibler divergence and the identification ratio of true segments. We also applied our model on two biological data sets: (1) reduced representation bisulfite sequencing (RRBS) data (A.Molaro et al. | 2011) (2) ENCODE Nimblegen tilled arrays (Sabo et al. | 2006). Our model shows good performance on segmentation of these two sequential data. In RRBS data it can further help identify differential methylation region (DMR) while in microarray data it can discover the DNAsel Hypersensitive Sites (DHSs). | 107 pages
Recommended Citation
Mo, Yifan, "A Stochastic Segmentation Model for Categorical and Continuous Features of various biological sequential" (2012). Stony Brook Theses and Dissertations Collection, 2006-2020 (closed to submissions). 556.
https://commons.library.stonybrook.edu/stony-brook-theses-and-dissertations-collection/556