Hierarchical models for 2D presence/absence data having ambiguous zeroes: With a biogeographical case study on dingo behaviour

This dissertation is primarily an applied statistical modelling investigation, motivated by a case study comprising real data and real questions. Theoretical questions on modelling and computation of normalization constants arose from pursuit of these data analytic questions. The essence of the thesis can be described as follows. Consider binary data observed on a two-dimensional lattice. A common problem with such data is the ambiguity of zeroes recorded. These may represent zero response given some threshold (presence) or that the threshold has not been triggered (absence). Suppose that the researcher wishes to estimate the effects of covariates on the binary responses, whilst taking into account underlying spatial variation, which is itself of some interest. This situation arises in many contexts and the dingo, cypress and toad case studies described in the motivation chapter are examples of this. Two main approaches to modelling and inference are investigated in this thesis. The first is frequentist and based on generalized linear models, with spatial variation modelled by using a block structure or by smoothing the residuals spatially. The EM algorithm can be used to obtain point estimates, coupled with bootstrapping or asymptotic MLE estimates for standard errors. The second approach is Bayesian and based on a three- or four-tier hierarchical model, comprising a logistic regression with covariates for the data layer, a binary Markov Random field (MRF) for the underlying spatial process, and suitable priors for parameters in these main models. The three-parameter autologistic model is a particular MRF of interest. Markov chain Monte Carlo (MCMC) methods comprising hybrid Metropolis/Gibbs samplers is suitable for computation in this situation. Model performance can be gauged by MCMC diagnostics. Model choice can be assessed by incorporating another tier in the modelling hierarchy. This requires evaluation of a normalization constant, a notoriously difficult problem. Difficulty with estimating the normalization constant for the MRF can be overcome by using a path integral approach, although this is a highly computationally intensive method. Different methods of estimating ratios of normalization constants (N Cs) are investigated, including importance sampling Monte Carlo (ISMC), dependent Monte Carlo based on MCMC simulations (MCMC), and reverse logistic regression (RLR). I develop an idea present though not fully developed in the literature, and propose the Integrated mean canonical statistic (IMCS) method for estimating log NC ratios for binary MRFs. The IMCS method falls within the framework of the newly identified path sampling methods of Gelman & Meng (1998) and outperforms ISMC, MCMC and RLR. It also does not rely on simplifying assumptions, such as ignoring spatio-temporal dependence in the process. A thorough investigation is made of the application of IMCS to the three-parameter Autologistic model. This work introduces background computations required for the full implementation of the four-tier model in Chapter 7. Two different extensions of the three-tier model to a four-tier version are investigated. The first extension incorporates temporal dependence in the underlying spatio-temporal process. The second extensions allows the successes and failures in the data layer to depend on time. The MCMC computational method is extended to incorporate the extra layer. A major contribution of the thesis is the development of a fully Bayesian approach to inference for these hierarchical models for the first time. Note: The author of this thesis has agreed to make it open access but invites people downloading the thesis to send her an email via the 'Contact Author' function.

Positive correlation of male and female Bactrocera tryoni (Froggatt) (Diptera: Tephritidae) catches in orange-ammonia traps

Populations of the Queensland fruit fly, Bactrocera tryoni, are routinely monitored using cue-lure, a male-only attractant. Such monitoring provides no information about females and there is little information available to show if male and female B. tryoni numbers are correlated in the field. Using a data set of 1 148 weekly clearances of orange-ammonia baited traps, which catch both males and females, the correlation between male and female numbers was tested for 48 weeks of the year (four weeks each month) and for the combined data set. Weekly male and female trap catches were almost entirely highly correlated, regardless of mean population size or time of year. For the whole year, the correlation between male and female numbers was r = 0.722, significant at p<0.001. Results suggest that changes in the number if male B. tryoni, as detected through cue-lure sampling, will reflect changes in numbers of female B. tryoni.

Comparison of standard image segmentation methods for segmentation of brain tumors from 2D MR images

In the analysis of medical images for computer-aided diagnosis and therapy, segmentation is often required as a preliminary step. Medical image segmentation is a complex and challenging task due to the complex nature of the images. The brain has a particularly complicated structure and its precise segmentation is very important for detecting tumors, edema, and necrotic tissues in order to prescribe appropriate therapy. Magnetic Resonance Imaging is an important diagnostic imaging technique utilized for early detection of abnormal changes in tissues and organs. It possesses good contrast resolution for different tissues and is, thus, preferred over Computerized Tomography for brain study. Therefore, the majority of research in medical image segmentation concerns MR images. As the core juncture of this research a set of MR images have been segmented using standard image segmentation techniques to isolate a brain tumor from the other regions of the brain. Subsequently the resultant images from the different segmentation techniques were compared with each other and analyzed by professional radiologists to find the segmentation technique which is the most accurate. Experimental results show that the Otsu’s thresholding method is the most suitable image segmentation method to segment a brain tumor from a Magnetic Resonance Image.

Correlation of the expression of human kallikrein-related peptidases 4 and 7 with the prognosis in oral squamous cell carcinoma

Background Very few articles have been written about the expression of kallikreins (KLK4 and KLK7) in oral cancers. Therefore, the purpose of this study was to examine and report on their prognostic potential. Methods Eighty archival blocks of primary oral cancers were sectioned and stained for KLK4 and KLK7 by immunohistochemistry. The percentage and the intensity of malignant keratinocyte staining were correlated with patient survival using Cox regression analysis. Results Both kallikreins were expressed strongly in the majority of tumor cells in 68 of 80 cases: these were mostly moderately or poorly differentiated neoplasms. Staining was particularly intense at the infiltrating front. Patients with intense staining had significantly shorter overall survival (p < .05). Conclusion This is the first observation on the patient survival influenced by kallikrein expression in oral carcinoma. The findings are consistent with those for carcinomas at other sites, in particular the prostate and ovary. KLK4 and/or KLK7 immunohistochemistry seems to have diagnostic and prognostic potential in this disease.