A continuous-index hidden Markov jump process for modeling DNA copy number data.
The number of copies of DNA in human cells can be measured using array comparative genomic hybridization (aCGH), which provides intensity ratios of sample to reference DNA at genomic locations corresponding to probes on a microarray. In the present paper, we devise a statistical model, based on a latent continuous-index Markov jump process, that is aimed to capture certain features of aCGH data, including probes that are unevenly long, unevenly spaced, and overlapping. The model has a continuous state space, with 1 state representing a normal copy number of 2, and the rest of the states being either amplifications or deletions. We adopt a Bayesian approach and apply Markov chain Monte Carlo (MCMC) methods for estimating the parameters and the Markov process. The model can be applied to data from both tiling bacterial artificial chromosome arrays and oligonucleotide arrays. We also compare a model with normal distributed noise to a model with t-distributed noise, showing that the latter is more robust to outliers.
Stjernqvist
Susann
S
Centre for Mathematical Sciences, Lund University, Box 118, 22100 Lund, Sweden. susann.stjernqvist@matstat.lu.se
Rydén
Tobias
T
D000465:Algorithms; D001499:Bayes Theorem; D056808:Biostatistics; D001943:Breast Neoplasms / Q000235:genetics; D022202:Chromosomes, Artificial, Bacterial / Q000235:genetics; D055028:Comparative Genomic Hybridization / Q000706:statistics & numerical data*; D004247:DNA / Q000032:analysis* / Q000235:genetics*; D004273:DNA, Neoplasm / Q000032:analysis / Q000235:genetics; D005260:Female; D018628:Gene Dosage; D006801:Humans; D008390:Markov Chains; D015233:Models, Statistical; D009010:Monte Carlo Method; D020411:Oligonucleotide Array Sequence Analysis / Q000706:statistics & numerical data
D004273:DNA, Neoplasm; D004247:DNA