Statistical analysis of genomic data : a new model for class prediction and inference
View/Open
Authors
Date
2011Supervisor
Type
Access Status
Education Level
Metadata
Show full item recordAbstract
Genomics is a major scientific revolution in this century. Highthroughput genomic data provides an opportunity for identifying genes and SNPs (singlenucleotide polymorphism) that are related to various clinical phenotypes. To deal with the sheer volume of genetic data being produced, it requires advanced methodological development in biostatistics that is lagging behind the technical capability to generate genomic data. SNPs have great importance in biomedical research for comparing regions of the genome between cohorts (such as casecontrol studies). Within a population, SNPs can be assigned a minor allele frequency, the lowest allele frequency at a locus that is observed in a particular population, and be recoded to binary datasets. Therefore, it is important to develop suitable statistical methods for SNPs analysis of genome alteration with the goal of contributing to the understanding of complex human diseases or traits such as mental health.In this thesis, we develop new statistical methodologies for the analysis of schizophrenia genomic data from the WA Genetic Epidemiology Resource (WAGER). The motivation is driven by the schizophrenia class prediction, (i.e. the prediction of individuals’ disease status through their genotype and quantitative traits). In general, individual’s disease status is a nominal variable, while genotypes can be converted into ordinal variables but are of high dimension. Note that the usual nonparametric regression that is developed for continuous variables cannot be applied here. There are some methodologies, such as the treebased logistic Nonparametric Pathwaybased Regression model (NPR) proposed by Wei and Li (2007)available in the literature. However, it is found that this model does not well adapt to the data set that we are analyzing. It is even worse than the (generalized) linear logistic regression model. Using logistic discrimination rule, together with adding quantitative traits, some important results have been obtained. However, some shortcomings remain. Firstly, the generalized linear logistic model has a high type I error rate for schizophrenia classification. Secondly, quantitative traits required for schizophrenia class prediction are performance assessments which demand several hours onsite participation by both assessor and assessee. These traits are generally quite difficult to reach even for a medium size sample. Meanwhile, though the laboratory analyzing cost is high, a person’s genotype can be obtained by merely collecting a drop of blood.Thus, two kinds of nonlinear models are proposed to capture the nonlinear effects in SNP datasets, which are categorical. The main contributions of this thesis are summarized as follows: • Two kinds of nonlinear threshold index logistic regression models are proposed to capture the nonlinear effects by applying the idea of threshold models (Tong (1983, 1990)) which are parametric and therefore applicable to the categorical data. One of the proposed models, which is called the partially linear threshold index logistic regression (PLTILoR) model, is given by log ( P(Yi = 1Xi) 1 − P(Yi = 1Xi) ) = ®TXi + g(¯TXi), (0.1) where Yi is the disease status of the ith person under casecontrol study, taking on values of 1 (case) or 0 (control), Xi is the vector of genotype variables, which is pdimensional, and the superscript T stands for transpose of a vector or matrix. Here, ® and ¯ are pdimensional unknown parameters with ¯ being an index vector used for the reduction of dimension, satisfying k¯k = 1 and ®T¯ = 0 for model identifiability, and g is, therefore, a onedimensional nonlinear function, which is modelled as stepwise linear function through threshold effect (Tong, 1990), given below. g(z) = (b1z + b2)I{z•c} + (b3z + b4)I{z>c}, (0.2) where bi’s and c are unknown parameters to be estimated and IA is an indicator function of the set A. In practice, the first component in model (0.1) could also be nonlinear. In this case, model (0.1) becomes log ( P(Yi = 1Xi) 1 − P(Yi = 1Xi) ) = g1(®TXi) + g2(¯TXi), (0.3) where k®k = 1, k¯k = 1 and ®T¯ = 0 for model identifiability, and g1 and g2 are two onedimensional nonlinear functions which are modelled by stepwise linear functions through threshold effects as follows: gk(z) = (bk1z + bk2)I{z•ck} + (bk3z + bk4)I{z>ck}, k = 1, 2, (0.4) where bki’s and ck’s are unknown parameters to be estimated. Thus, (0.3) and (0.4) form an additive threshold index logistic regression (ATILoR) model. • A maximum likelihood methodology is developed to estimate the unknown parameters in the PLTILoR and ATILoR models. Simulation studies have found that the proposed methodology works well for finite size samples. • Empirical studies of the proposed models applied to the analysis of schizophrenia genomic data from the WA Genetic Epidemiology Resource (WAGER) have shown that ATILoR model is very successful in reducing the type I error rate in schizophrenia classification without even using quantitative traits. It outperforms the generalized linear logistic model that is widely used in the literature.
Department
Collections
Related items
Showing items related by title, author, creator and subject.

Jiang, Z.; Du, C.; Jablensky, A.; Liang, H.; Lu, Z.; Ma, Y.; Teo, Kok Lay (2014)Genetic information, such as single nucleotide polymorphism (SNP) data, has been widely recognized as useful in prediction of disease risk. However, how to model the genetic data that is often categorical in disease class ...

Siva Subramaniam, Nitthiya (2012)The major histocompatibility complex (MHC) is a chromosomal region associated with immune responsiveness in vertebrates. Over four decades many studies have demonstrated important associations between MHC loci and disease ...

Lee, Chee Yang (2009)The main aim of this project has been to provide insights into the genomic organization of the class II region of the ovine major histocompatibility complex (MHC), a chromosomal region containing genes that control adaptive ...