Statistics and Probability Commons^™

Copy Number Variants In Candidate Genes Are Genetic Modifiers Of Hirschsprung Disease, Qian Jiang, Yen Yi Ho, Li Hao, Courtney Nichols Berrios, Aravinda Chakravarti

Faculty Publications

Hirschsprung disease (HSCR) is a neurocristopathy characterized by absence of intramural ganglion cells along variable lengths of the gastrointestinal tract. The HSCR phenotype is highly variable with respect to gender, length of aganglionosis, familiality and the presence of additional anomalies. By molecular genetic analysis, a minimum of 11 neuro-developmental genes (RET, GDNF, NRTN, SOX10, EDNRB, EDN3, ECE1, ZFHX1B, PHOX2B, KIAA1279, TCF4) are known to harbor rare, high-penetrance mutations that confer a large risk to the bearer. In addition, two other genes (RET, NRG1) harbor common, low-penetrance polymorphisms that contribute only partially to risk and can act as genetic modifiers. To …

Component Extraction Of Complex Biomedical Signal And Performance Analysis Based On Different Algorithm, Hemant Pasusangai Kasturiwale

Johns Hopkins University, Dept. of Biostatistics Working Papers

Biomedical signals can arise from one or many sources including heart ,brains and endocrine systems. Multiple sources poses challenge to researchers which may have contaminated with artifacts and noise. The Biomedical time series signal are like electroencephalogram(EEG),electrocardiogram(ECG),etc The morphology of the cardiac signal is very important in most of diagnostics based on the ECG. The diagnosis of patient is based on visual observation of recorded ECG,EEG,etc, may not be accurate. To achieve better understanding , PCA (Principal Component Analysis) and ICA algorithms helps in analyzing ECG signals . The immense scope in the field of biomedical-signal processing Independent Component Analysis( …

Climate Change And Community Dynamics: A Hierarchical Bayesian Model Of Resource-Driven Changes In A Desert Rodent Community, Glenda M. Yenni

All Graduate Plan B and other Reports, Spring 1920 to Spring 2023

Predicting effects of climate change on species persistence often assumes that those species are responding to abiotic effects alone. However, biotic interactions between community members may affect species’ ability to respond to abiotic changes. Latent Gaussian models of resource availability using precipitation and NDVI and accounting for spatial autocorrelation and rodent group-level uncertainty in the process are developed to detect differences in seasons, groups, and the experimental removal of one group. Precipitation and NDVI have overall positive effects on rodent energy use as expected, but meaningful differences were detected. Differences in the importance of seasonality when the dominant group was …

Analysis Of Morris Water Maze Data With Bayesian Statistical Methods, Maxym V. Myroshnychenko, Anton Westveld, Jefferson Kinney

Festival of Communities: UG Symposium (Posters)

Neuroscientists commonly use a Morris Water Maze to assess learning in rodents. In his kind of a maze, the subjects learn to swim toward a platform hidden in opaque water as they orient themselves according to the cues on the walls. This protocol presents a challenge to statistical analysis, because an artificial cut-off must be set for those experimental subjects that do not reach the platform so as they do not drown from exhaustion. This fact leads to the data being right censored. In our experimental data, which compares learning in rodents that have chemically induced symptoms of schizophrenia to …

A Generalized Approach For Testing The Association Of A Set Of Predictors With An Outcome: A Gene Based Test, Benjamin A. Goldstein, Alan E. Hubbard, Lisa F. Barcellos

U.C. Berkeley Division of Biostatistics Working Paper Series

In many analyses, one has data on one level but desires to draw inference on another level. For example, in genetic association studies, one observes units of DNA referred to as SNPs, but wants to determine whether genes that are comprised of SNPs are associated with disease. While there are some available approaches for addressing this issue, they usually involve making parametric assumptions and are not easily generalizable. A statistical test is proposed for testing the association of a set of variables with an outcome of interest. No assumptions are made about the functional form relating the variables to the …

Assessing The Necessity Of Chimpanzee Experimentation, Andrew Knight

Experimentation Collection

No abstract provided.

A Guide To Defining And Implementing Protocols For The Welfare Assessment Of Laboratory Animals: Eleventh Report Of The Bvaawf/Frame/Rspca/Ufaw Joint Working Group On Refinement, P. Hawkins, D. B. Morton, O. Burman, N. Dennison, P. Honess, M. Jennings, S. Lane, V. Middleton, J. V. Roughan, S. Wells, K. Westwood

Research Methodology and Laboratory Animals Collection

The refinement of husbandry and procedures to reduce animal suffering and improve welfare is an essential component of humane science. Successful refinement depends upon the ability to assess animal welfare effectively, and detect any signs of pain or distress as rapidly as possible, so that any suffering can be alleviated. This document provides practical guidance on setting up and operating effective protocols for the welfare assessment of animals used in research and testing. It sets out general principles for more objective observation of animals, recognizing and assessing indicators of pain or distress and tailoring these to individual projects. Systems for …

Inflated Type I Error Rates When Using Aggregation Methods To Analyze Rare Variants In The 1000 Genomes Project Exon Sequencing Data In Unrelated Individuals: Summary Results From Group 7 At Genetic Analysis Workshop 17, Nathan L. Tintle, Hugues Aschard, Inchi Hu, Nora Nock, Haitian Wang, Elizabeth Pugh

Faculty Work Comprehensive List

As part of Genetic Analysis Workshop 17 (GAW17), our group considered the application of novel and standard approaches to the analysis of genotype-phenotype association in next-generation sequencing data. Our group identified a major issue in the analysis of the GAW17 next-generation sequencing data: type I error and false-positive report probability rates higher than those expected based on empirical type I error levels (as high as 90%). Two main causes emerged: population stratification and long-range correlation (gametic phase disequilibrium) between rare variants. Population stratification was expected because of the diverse sample. Correlation between rare variants was attributable to both random causes …

Identification Of Genetic Association Of Multiple Rare Variants Using Collapsing Methods, Yan V. Sun, Yun Ju Sung, Nathan L. Tintle, Andreas Ziegler

Faculty Work Comprehensive List

Next-generation sequencing technology allows investigation of both common and rare variants in humans. Exomes are sequenced on the population level or in families to further study the genetics of human diseases. Genetic Analysis Workshop 17 (GAW17) provided exomic data from the 1000 Genomes Project and simulated phenotypes. These data enabled evaluations of existing and newly developed statistical methods for rare variant sequence analysis for which standard statistical methods fail because of the rareness of the alleles. Various alternative approaches have been proposed that overcome the rareness problem by combining multiple rare variants within a gene. These approaches are termed collapsing …

Computational Methods Of Hidden Markov Models With Respect To Cpg Island Prediction In Dna Sequences, Roberto Angel Ortega

Open Access Theses & Dissertations

Hidden Markov models (HMM's) are a specific case of Markov models where, contrary to Markov chains, the observer is unaware of what state the model was in when the symbol is observed. Like Markov chains, HMM's assume that the future state of a sequence is dependent only on the current state of the sequence. The parameters associated with HMM's are transition and emission probabilities, where transition probabilities are associated with the probability of transitioning from one state to another, and emission probabilities are the probabilities associated with observing a symbol given it came from a specific state.

The structure of …

Evaluating Methods For The Analysis Of Rare Variants In Sequence Data, Alexander Luedtke, Scott Powers, Ashley Petersen, Alexandra Sitarik, Airat Bekmetjev, Nathan L. Tintle

Faculty Work Comprehensive List

A number of rare variant statistical methods have been proposed for analysis of the impending wave of next-generation sequencing data. To date, there are few direct comparisons of these methods on real sequence data. Furthermore, there is a strong need for practical advice on the proper analytic strategies for rare variant analysis. We compare four recently proposed rare variant methods (combined multivariate and collapsing, weighted sum, proportion regression, and cumulative minor allele test) on simulated phenotype and next-generation sequencing data as part of Genetic Analysis Workshop 17. Overall, we find that all analyzed methods have serious practical limitations on identifying …

Evaluating Methods For Combining Rare Variant Data In Pathway-Based Tests Of Genetic Association, Ashley Petersen, Alexandra Sitarik, Alexander Luedtke, Scott Powers, Airat Bekmetjev, Nathan L. Tintle

Faculty Work Comprehensive List

Analyzing sets of genes in genome-wide association studies is a relatively new approach that aims to capitalize on biological knowledge about the interactions of genes in biological pathways. This approach, called pathway analysis or gene set analysis, has not yet been applied to the analysis of rare variants. Applying pathway analysis to rare variants offers two competing approaches. In the first approach rare variant statistics are used to generate p-values for each gene (e.g., combined multivariate collapsing [CMC] or weighted-sum [WS]) and the gene-level p-values are combined using standard pathway analysis methods (e.g., gene set enrichment analysis or …

Identifying Rare Variants From Exome Scans: The Gaw17 Experience, Saurabh Ghosh, Heike Bickeboller, Julia Bailey, Joan E. Bailey-Wilson, Rita Cantor, Robert Culverhouse, Warwick Daw, Anita L. Destefano, Corinne D. Engelman, Anthony Hinrichs, Jeanine Houwing-Duistermaat, Inke R. Konig, Jack Kent, Nan Laird, Nathan Pankratz, Andrew Paterson, Elizabeth Pugh, Brian Suarez, Yan Sun, Alun Thomas, Nathan L. Tintle, Xiaofeng Zhu, Andreas Ziegler, Jean W. Maccluer, Laura Almasy

Faculty Work Comprehensive List

Genetic Analysis Workshop 17 (GAW17) provided a platform for evaluating existing statistical genetic methods and for developing novel methods to analyze rare variants that modulate complex traits. In this article, we present an overview of the 1000 Genomes Project exome data and simulated phenotype data that were distributed to GAW17 participants for analyses, the different issues addressed by the participants, and the process of preparation of manuscripts resulting from the discussions during the workshop

Weighted Scores Method For Regression Models With Dependent Data, Aristidis K. Nikoloulopoulos, Harry Joe, N. Rao Chaganty

Mathematics & Statistics Faculty Publications

There are copula-based statistical models in the literature for regression with dependent data such as clustered and longitudinal overdispersed counts, for which parameter estimation and inference are straightforward. For situations where the main interest is in the regression and other univariate parameters and not the dependence, we propose a "weighted scores method", which is based on weighting score functions of the univariate margins. The weight matrices are obtained initially fitting a discretized multivariate normal distribution, which admits a wide range of dependence. The general methodology is applied to negative binomial regression models. Asymptotic and small-sample efficiency calculations show that our …

Distributional Properties Of Inversions And Segmentation Algorithms For Rna Sequences, Sameera Dhananjaya Viswakula

Open Access Theses & Dissertations

Ribonucleic acid (RNA) is a long single stranded molecule made up of four types of nucleotide bases: Adenine (A), Cytosine(C), Guanine (G) and Uracil (U). It folds back on itself and forms C-G and A-U complementary base pairs. The set of such hydrogen-bonded pairs in an RNA molecule is called its secondary structure. Knowing the secondary structure of RNA is useful for understanding its biological function. Prediction of RNA secondary structure from the nucleotide sequence has been an important bioinformatics problem for over two decades.

The work in this thesis is motivated by the need to improve the secondary structure …

Class Discovery And Prediction Of Tumor With Microarray Data, Bo Liu

All Graduate Theses, Dissertations, and Other Capstone Projects

Current microarray technology is able take a single tissue sample to construct an Affymetrix oglionucleotide array containing (estimated) expression levels of thousands of different genes for that tissue. The objective is to develop a more systematic approach to cancer classification based on Affymetrix oglionucleotide microarrays. For this purpose, I studied published colon cancer microarray data. Colon cancer, with 655,000 deaths worldwide per year, has become the fourth most common form of cancer in the United States and the third leading cause of cancer - related death in the Western world. This research has been focuses in two areas: class discovery, …

A Novel Device For Cell-Cell Electrofusion, Justin T. Stewart

USF Tampa Graduate Theses and Dissertations

Cell transplantation therapy is a potentially powerful tool and can be used to replace defective cells with healthy cells. This offers the possibility of alleviating the destructive symptoms for many diseases such as Parkinson's disease, Alzheimer's disease, stroke, spinal cord trauma, Type I diabetes and many more. While there are many diseases that could be positively impacted from cell transplantation therapy, the focus of this research is insulin dependent, Type I Diabetes.

The Islets of Langerhans are composed of various types of cells located in the pancreas and are responsible for a variety of biochemical functions. Specifically, the beta Islet …

Analysis Of Differential Gene Expression And Alternative Splicing In The Liver And Gastrointestinal Tract In The Lactating Rat, Antony Thomas Athippozhy

University of Kentucky Doctoral Dissertations

Rat exon microarrays were utilized to detect changes in mRNA expression and alternative splicing in the liver, duodenum, jejunum, and ileum of the lactating rat when compared to age-matched virgin controls. Analysis of data at the level of gene expression revealed differential expression of genes involved in cholesterol biosynthesis in each tissue examined, suggesting increased Sterol Response Element Binding Protein activity. We also detected decreased mRNA from components of the T-cell signaling pathway in the jejunum and ileum. We characterized expression of solute carrier and adenosine triphosphate binding cassette proteins. In addition to characterizing genes by pathway, we have also …

Estimating The Reproductive Numbers For The 2008-2009 Cholera Outbreaks In Zimbabwe, Zindoga Mukandavire, Shu Liao, Jin Wang, Holly Gaff, David L. Smith, J. Glenn Morris Jr.

Biological Sciences Faculty Publications

Cholera remains an important global cause of morbidity and mortality, capable of causing periodic epidemic disease. Beginning in August 2008, a major cholera epidemic occurred in Zimbabwe, with 98,585 reported cases and 4,287 deaths. The dynamics of such outbreaks, particularly in nonestuarine regions, are not well understood. We explored the utility of mathematical models in understanding transmission dynamics of cholera and in assessing the magnitude of interventions necessary to control epidemic disease. Weekly data on reported cholera cases were obtained from the Zimbabwe Ministry of Health and Child Welfare (MoHCW) for the period from November 13, 2008 to July 31, …

Computational Network Analysis Of The Anatomical And Genetic Organizations In The Mouse Brain, Shuiwang Ji

Computer Science Faculty Publications

Motivation: The mammalian central nervous system (CNS) generates high-level behavior and cognitive functions. Elucidating the anatomical and genetic organizations in the CNS is a key step toward understanding the functional brain circuitry. The CNS contains an enormous number of cell types, each with unique gene expression patterns. Therefore, it is of central importance to capture the spatial expression patterns in the brain. Currently, genome-wide atlas of spatial expression patterns in the mouse brain has been made available, and the data are in the form of aligned 3D data arrays. The sheer volume and complexity of these data pose significant challenges …

Minimum Description Length Measures Of Evidence For Enrichment, Zhenyu Yang, David R. Bickel

COBRA Preprint Series

In order to functionally interpret differentially expressed genes or other discovered features, researchers seek to detect enrichment in the form of overrepresentation of discovered features associated with a biological process. Most enrichment methods treat the p-value as the measure of evidence using a statistical test such as the binomial test, Fisher's exact test or the hypergeometric test. However, the p-value is not interpretable as a measure of evidence apart from adjustments in light of the sample size. As a measure of evidence supporting one hypothesis over the other, the Bayes factor (BF) overcomes this drawback of the p-value but lacks …

Reconstructability Analysis Of Epistasis, Martin Zwick

Complex Systems Faculty Publications and Presentations

The literature on epistasis describes various methods to detect epistatic interactions and to classify different types of epistasis. Reconstructability analysis (RA) has recently been used to detect epistasis in genomic data. This paper shows that RA offers a classification of types of epistasis at three levels of resolution (variable-based models without loops, variable-based models with loops, state-based models). These types can be defined by the simplest RA structures that model the data without information loss; a more detailed classification can be defined by the information content of multiple candidate structures. The RA classification can be augmented with structures from related …

Approximate Approach To The Das Model Of Fractional Logistic Population Growth, S. Das, P. K. Gupta, K. Vishal

Applications and Applied Mathematics: An International Journal (AAM)

In this article, the analytical method, Homotopy perturbation method (HPM) has been successfully implemented for solving nonlinear logistic model of fractional order. The fractional derivatives are described in the Caputo sense. Using initial value, the explicit solutions of population size for different particular cases have been derived. Numerical results show that the method is extremely efficient to solve this complicated biological model.

A Bayesian Shared Component Model For Genetic Association Studies, Juan J. Abellan, Carlos Abellan, Juan R. Gonzalez

COBRA Preprint Series

We present a novel approach to address genome association studies between single nucleotide polymorphisms (SNPs) and disease. We propose a Bayesian shared component model to tease out the genotype information that is common to cases and controls from the one that is specific to cases only. This allows to detect the SNPs that show the strongest association with the disease. The model can be applied to case-control studies with more than one disease. In fact, we illustrate the use of this model with a dataset of 23,418 SNPs from a case-control study by The Welcome Trust Case Control Consortium (2007) …

Minimum Description Length And Empirical Bayes Methods Of Identifying Snps Associated With Disease, Ye Yang, David R. Bickel

COBRA Preprint Series

The goal of determining which of hundreds of thousands of SNPs are associated with disease poses one of the most challenging multiple testing problems. Using the empirical Bayes approach, the local false discovery rate (LFDR) estimated using popular semiparametric models has enjoyed success in simultaneous inference. However, the estimated LFDR can be biased because the semiparametric approach tends to overestimate the proportion of the non-associated single nucleotide polymorphisms (SNPs). One of the negative consequences is that, like conventional p-values, such LFDR estimates cannot quantify the amount of information in the data that favors the null hypothesis of no disease-association.

We …

Microct Study Of Critical-Sized Defects Restored With A New Rhbmp-2 Carrier, Joanne S. Kim

Loma Linda University Electronic Theses, Dissertations & Projects

Introduction: Regenerating bone in critical-sized maxillofacial osseous defects remains a challenge. The utilization of recombinant human bone morphogenetic protein-2 (rhBMP-2) to form new bone has been shown to be a promising alternative to autogenous bone grafts, and the development of a carrier matrix to effectively deliver the rhBMP-2 has been identified as the critical factor necessary for the successful clinical application of rhBMP-2. The objective of this study was to perform a microCT analysis to examine the density of reconstructed bone using different rhBMP-2 carriers in critical-sized defects.

Material and Methods: Critical-sized defects (approximately 2.5 cm) were created in the …

A Perturbation Method For Inference On Regularized Regression Estimates, Jessica Minnier, Lu Tian, Tianxi Cai

Harvard University Biostatistics Working Paper Series

No abstract provided.

The Strength Of Statistical Evidence For Composite Hypotheses: Inference To The Best Explanation, David R. Bickel

COBRA Preprint Series

A general function to quantify the weight of evidence in a sample of data for one hypothesis over another is derived from the law of likelihood and from a statistical formalization of inference to the best explanation. For a fixed parameter of interest, the resulting weight of evidence that favors one composite hypothesis over another is the likelihood ratio using the parameter value consistent with each hypothesis that maximizes the likelihood function over the parameter of interest. Since the weight of evidence is generally only known up to a nuisance parameter, it is approximated by replacing the likelihood function with …

Powerful Snp Set Analysis For Case-Control Genome Wide Association Studies, Michael C. Wu, Peter Kraft, Michael P. Epstein, Deanne M. Taylor, Stephen J. Chanock, David J. Hunter, Xihong Lin

Harvard University Biostatistics Working Paper Series

No abstract provided.

Survival Prediction For Brain Tumor Patients Using Gene Expression Data, Vinicius Bonato

Dissertations and Theses (Open Access)

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. …