Bayesian models for spatio-temporal assessment of disease

This thesis has contributed to the advancement of knowledge in disease modelling by addressing interesting and crucial issues relevant to modelling health data over space and time. The research has led to the increased understanding of spatial scales, temporal scales, and spatial smoothing for modelling diseases, in terms of their methodology and applications. This research is of particular significance to researchers seeking to employ statistical modelling techniques over space and time in various disciplines. A broad class of statistical models are employed to assess what impact of spatial and temporal scales have on simulated and real data.

Spatial and temporal patterns of locally-acquired dengue transmission in northern Queensland, Australia, 1993-2012

BACKGROUND: Dengue has been a major public health concern in Australia since it re-emerged in Queensland in 1992-1993. We explored spatio-temporal characteristics of locally-acquired dengue cases in northern tropical Queensland, Australia during the period 1993-2012. METHODS: Locally-acquired notified cases of dengue were collected for northern tropical Queensland from 1993 to 2012. Descriptive spatial and temporal analyses were conducted using geographic information system tools and geostatistical techniques. RESULTS: 2,398 locally-acquired dengue cases were recorded in northern tropical Queensland during the study period. The areas affected by the dengue cases exhibited spatial and temporal variation over the study period. Notified cases of dengue occurred more frequently in autumn. Mapping of dengue by statistical local areas (census units) reveals the presence of substantial spatio-temporal variation over time and place. Statistically significant differences in dengue incidence rates among males and females (with more cases in females) (χ(2) = 15.17, d.f.  = 1, p<0.01). Differences were observed among age groups, but these were not statistically significant. There was a significant positive spatial autocorrelation of dengue incidence for the four sub-periods, with the Moran's I statistic ranging from 0.011 to 0.463 (p<0.01). Semi-variogram analysis and smoothed maps created from interpolation techniques indicate that the pattern of spatial autocorrelation was not homogeneous across the northern Queensland. CONCLUSIONS: Tropical areas are potential high-risk areas for mosquito-borne diseases such as dengue. This study demonstrated that the locally-acquired dengue cases have exhibited a spatial and temporal variation over the past twenty years in northern tropical Queensland, Australia. Therefore, this study provides an impetus for further investigation of clusters and risk factors in these high-risk areas.

Temporal characteristics of melanopsin inputs to the human pupil light reflex

Rods, cones and melanopsin containing intrinsically photosensitive retinal ganglion cells (ipRGCs) operate in concert to regulate pupil diameter. The temporal properties of intrinsic ipRGC signalling are distinct to those of rods and cones, including longer latencies and sustained signalling after light offset. We examined whether the melanopsin mediated post-illumination pupil response (PIPR) and pupil constriction were dependent upon the inter-stimulus interval (ISI) between successive light pulses and the temporal frequency of sinusoidal light stimuli. Melanopsin excitation was altered by variation of stimulus wavelength (464 nm and 638 nm lights) and irradiance (11.4 and 15.2 log photons cm(-2) s(-1)). We found that 6s PIPR amplitude was independent of ISI and temporal frequency for all melanopsin excitation levels, indicating complete summation. In contrast to the PIPR, the maximum pupil constriction increased with increasing ISI with high and low melanopsin excitation, but time to minimum diameter was slower with high melanopsin excitation only. This melanopsin response to briefly presented pulses (16 and 100 ms) slows the temporal response of the maximum pupil constriction. We also demonstrate that high melanopsin excitation attenuates the phasic peak-trough pupil amplitude compared to conditions with low melanopsin excitation, indicating an interaction between inner and outer retinal inputs to the pupil light reflex. We infer that outer retina summation is important for rapidly controlling pupil diameter in response to short timescale fluctuations in illumination and may occur at two potential sites, one that is presynaptic to extrinsic photoreceptor input to ipRGCs, or another within the pupil control pathway if ipRGCs have differential temporal tuning to extrinsic and intrinsic signalling.

Inferring lung cancer risk factor patterns through joint Bayesian spatio-temporal analysis

Background: Preventing risk factor exposure is vital to reduce the high burden from lung cancer. The leading risk factor for developing lung cancer is tobacco smoking. In Australia, despite apparent success in reducing smoking prevalence, there is limited information on small area patterns and small area temporal trends. We sought to estimate spatio-temporal patterns for lung cancer risk factors using routinely collected population-based cancer data. Methods: The analysis used a Bayesian shared component spatio-temporal model, with male and female lung cancer included separately. The shared component reflected exposure to lung cancer risk factors, and was modelled over 477 statistical local areas (SLAs) and 15 years in Queensland, Australia. Analyses were also run adjusting for area-level socioeconomic disadvantage, Indigenous population composition, or remoteness. Results: Strong spatial patterns were observed in the underlying risk factor exposure for both males (median Relative Risk (RR) across SLAs compared to the Queensland average ranged from 0.48-2.00) and females (median RR range across SLAs 0.53-1.80), with high exposure observed in many remote areas. Strong temporal trends were also observed. Males showed a decrease in the underlying risk across time, while females showed an increase followed by a decrease in the final two years. These patterns were largely consistent across each SLA. The high underlying risk estimates observed among disadvantaged, remote and indigenous areas decreased after adjustment, particularly among females. Conclusion: The modelled underlying exposure appeared to reflect previous smoking prevalence, with a lag period of around 30 years, consistent with the time taken to develop lung cancer. The consistent temporal trends in lung cancer risk factors across small areas support the hypothesis that past interventions have been equally effective across the state. However, this also means that spatial inequalities have remained unaddressed, highlighting the potential for future interventions, particularly among remote areas.

Quantifying the state of populations and effects of disturbances at large spatio-temporal scales: The case of coral populations in the great barrier reef

This project was a step forward in applying statistical methods and models to provide new insights for more informed decision-making at large spatial scales. The model has been designed to address complicated effects of ecological processes that govern the state of populations and uncertainties inherent in large spatio-temporal datasets. Specifically, the thesis contributes to better understanding and management of the Great Barrier Reef.

Hot spot detection and spatio-temporal dynamics of dengue in Queensland, Australia

Dengue has been a major public health concern in Australia since it re-emerged in Queensland in 1992-1993. This study explored spatio-temporal distribution and clustering of locally-acquired dengue cases in Queensland State, Australia and identified target areas for effective interventions. A computerised locally-acquired dengue case dataset was collected from Queensland Health for Queensland from 1993 to 2012. Descriptive spatial and temporal analyses were conducted using geographic information system tools and geostatistical techniques. Dengue hot spots were detected using SatScan method. Descriptive spatial analysis showed that a total of 2,398 locally-acquired dengue cases were recorded in central and northern regions of tropical Queensland. A seasonal pattern was observed with most of the cases occurring in autumn. Spatial and temporal variation of dengue cases was observed in the geographic areas affected by dengue over time. Tropical areas are potential high-risk areas for mosquito-borne diseases such as dengue. This study demonstrated that the locally-acquired dengue cases have exhibited a spatial and temporal variation over the past twenty years in tropical Queensland, Australia. There is a clear evidence for the existence of statistically significant clusters of dengue and these clusters varied over time. These findings enabled us to detect and target dengue clusters suggesting that the use of geospatial information can assist the health authority in planning dengue control activities and it would allow for better design and implementation of dengue management programs.

Spatio-temporal modelling of cancer data in Queensland using Bayesian methods

Cancer is the leading contributor to the disease burden in Australia. This thesis develops and applies Bayesian hierarchical models to facilitate an investigation of the spatial and temporal associations for cancer diagnosis and survival among Queenslanders. The key objectives are to document and quantify the importance of spatial inequalities, explore factors influencing these inequalities, and investigate how spatial inequalities change over time. Existing Bayesian hierarchical models are refined, new models and methods developed, and tangible benefits obtained for cancer patients in Queensland. The versatility of using Bayesian models in cancer control are clearly demonstrated through these detailed and comprehensive analyses.

On the order of the fractional Laplacian in determining the spatio-temporal evolution of a space-fractional model of cardiac electrophysiology

Space-fractional operators have been used with success in a variety of practical applications to describe transport processes in media characterised by spatial connectivity properties and high structural heterogeneity altering the classical laws of diffusion. This study provides a systematic investigation of the spatio-temporal effects of a space-fractional model in cardiac electrophysiology. We consider a simplified model of electrical pulse propagation through cardiac tissue, namely the monodomain formulation of the Beeler-Reuter cell model on insulated tissue fibres, and obtain a space-fractional modification of the model by using the spectral definition of the one-dimensional continuous fractional Laplacian. The spectral decomposition of the fractional operator allows us to develop an efficient numerical method for the space-fractional problem. Particular attention is paid to the role played by the fractional operator in determining the solution behaviour and to the identification of crucial differences between the non-fractional and the fractional cases. We find a positive linear dependence of the depolarization peak height and a power law decay of notch and dome peak amplitudes for decreasing orders of the fractional operator. Furthermore, we establish a quadratic relationship in conduction velocity, and quantify the increasingly wider action potential foot and more pronounced dispersion of action potential duration, as the fractional order is decreased. A discussion of the physiological interpretation of the presented findings is made.

Calculation of the speed of walking of individuals with transfemoral amputation fitted with bone-anchored prosthesis

Typically, the walking ability of individuals with a transfemoral amputation (TFA) can be represented by the speed of walking (SofW) obtained in experimental settings. Recent developments in portable kinetic systems allow assessing the level of activity of TFA during actual daily living outside the confined space of a gait lab. Unfortunately, only minimal spatio-temporal characteristics could be extracted from the kinetic data including the cadence and the duration on gait cycles. Therefore, there is a need for a way to use some of these characteristics to assess the instantaneous speed of walking during daily living. The purpose of the study was to compare several methods to determine SofW using minimal spatial gait characteristics.

Functional outcome of transfemoral amputees fitted with an osseointegrated fixation : temporal gait characteristics

The purpose of this study was to characterise the functional outcome of 12 transfemoral amputees fitted with osseointegrated fixation using temporal gait characteristics. The objectives were (A) to present the cadence, duration of gait cycle, support and swing phases with an emphasis on the stride-to-stride and participant-to-participant variability, and (B) to compare these temporal variables with normative data extracted from the literature focusing on transfemoral amputees fitted with a socket and able-bodied participants. The temporal variables were extracted from the load applied on the residuum during straight level walking, which was collected at 200 Hz by a transducer. A total of 613 strides were assessed. The cadence (46±4 strides/min), the duration of the gait cycle (1.29±0.11 s), support (0.73±0.07 s, 57±3% of CG) and swing (0.56±0.07 s, 43±3% of GC) phases of the participants were 2% quicker, 3%, 6% shorter and 1% longer than transfemoral amputees using a socket as well as 11% slower, 9%, 6% and 13% longer than able-bodied, respectively. All combined, the results indicated that the fitting of an osseointegrated fixation has enabled this group of amputees to restore their locomotion with a highly functional level. Further longitudinal and cross-sectional studies would be required to confirm these outcomes. Nonetheless, the data presented can be used as benchmark for future comparisons. It can also be used as input in generic algorithms using templates of patterns of loading to recognise activities of daily living and to detect falls.

Spatial and temporal analysis of Barmah Forest virus disease in Queensland, Australia

Barmah Forest virus (BFV) disease is one of the most widespread mosquito-borne diseases in Australia. The number of outbreaks and the incidence rate of BFV in Australia have attracted growing concerns about the spatio-temporal complexity and underlying risk factors of BFV disease. A large number of notifications has been recorded continuously in Queensland since 1992. Yet, little is known about the spatial and temporal characteristics of the disease. I aim to use notification data to better understand the effects of climatic, demographic, socio-economic and ecological risk factors on the spatial epidemiology of BFV disease transmission, develop predictive risk models and forecast future disease risks under climate change scenarios. Computerised data files of daily notifications of BFV disease and climatic variables in Queensland during 1992-2008 were obtained from Queensland Health and Australian Bureau of Meteorology, respectively. Projections on climate data for years 2025, 2050 and 2100 were obtained from Council of Scientific Industrial Research Organisation. Data on socio-economic, demographic and ecological factors were also obtained from relevant government departments as follows: 1) socio-economic and demographic data from Australian Bureau of Statistics; 2) wetlands data from Department of Environment and Resource Management and 3) tidal readings from Queensland Department of Transport and Main roads. Disease notifications were geocoded and spatial and temporal patterns of disease were investigated using geostatistics. Visualisation of BFV disease incidence rates through mapping reveals the presence of substantial spatio-temporal variation at statistical local areas (SLA) over time. Results reveal high incidence rates of BFV disease along coastal areas compared to the whole area of Queensland. A Mantel-Haenszel Chi-square analysis for trend reveals a statistically significant relationship between BFV disease incidence rates and age groups (ƒÓ2 = 7587, p<0.01). Semi-variogram analysis and smoothed maps created from interpolation techniques indicate that the pattern of spatial autocorrelation was not homogeneous across the state. A cluster analysis was used to detect the hot spots/clusters of BFV disease at a SLA level. Most likely spatial and space-time clusters are detected at the same locations across coastal Queensland (p<0.05). The study demonstrates heterogeneity of disease risk at a SLA level and reveals the spatial and temporal clustering of BFV disease in Queensland. Discriminant analysis was employed to establish a link between wetland classes, climate zones and BFV disease. This is because the importance of wetlands in the transmission of BFV disease remains unclear. The multivariable discriminant modelling analyses demonstrate that wetland types of saline 1, riverine and saline tidal influence were the most significant risk factors for BFV disease in all climate and buffer zones, while lacustrine, palustrine, estuarine and saline 2 and saline 3 wetlands were less important. The model accuracies were 76%, 98% and 100% for BFV risk in subtropical, tropical and temperate climate zones, respectively. This study demonstrates that BFV disease risk varied with wetland class and climate zone. The study suggests that wetlands may act as potential breeding habitats for BFV vectors. Multivariable spatial regression models were applied to assess the impact of spatial climatic, socio-economic and tidal factors on the BFV disease in Queensland. Spatial regression models were developed to account for spatial effects. Spatial regression models generated superior estimates over a traditional regression model. In the spatial regression models, BFV disease incidence shows an inverse relationship with minimum temperature, low tide and distance to coast, and positive relationship with rainfall in coastal areas whereas in whole Queensland the disease shows an inverse relationship with minimum temperature and high tide and positive relationship with rainfall. This study determines the most significant spatial risk factors for BFV disease across Queensland. Empirical models were developed to forecast the future risk of BFV disease outbreaks in coastal Queensland using existing climatic, socio-economic and tidal conditions under climate change scenarios. Logistic regression models were developed using BFV disease outbreak data for the existing period (2000-2008). The most parsimonious model had high sensitivity, specificity and accuracy and this model was used to estimate and forecast BFV disease outbreaks for years 2025, 2050 and 2100 under climate change scenarios for Australia. Important contributions arising from this research are that: (i) it is innovative to identify high-risk coastal areas by creating buffers based on grid-centroid and the use of fine-grained spatial units, i.e., mesh blocks; (ii) a spatial regression method was used to account for spatial dependence and heterogeneity of data in the study area; (iii) it determined a range of potential spatial risk factors for BFV disease; and (iv) it predicted the future risk of BFV disease outbreaks under climate change scenarios in Queensland, Australia. In conclusion, the thesis demonstrates that the distribution of BFV disease exhibits a distinct spatial and temporal variation. Such variation is influenced by a range of spatial risk factors including climatic, demographic, socio-economic, ecological and tidal variables. The thesis demonstrates that spatial regression method can be applied to better understand the transmission dynamics of BFV disease and its risk factors. The research findings show that disease notification data can be integrated with multi-factorial risk factor data to develop build-up models and forecast future potential disease risks under climate change scenarios. This thesis may have implications in BFV disease control and prevention programs in Queensland.

Bluetooth and Wi-Fi MAC address based crowd data collection and monitoring : benefits, challenges and enhancement

This paper firstly presents the benefits and critical challenges on the use of Bluetooth and Wi-Fi for crowd data collection and monitoring. The major challenges include antenna characteristics, environment’s complexity and scanning features. Wi-Fi and Bluetooth are compared in this paper in terms of architecture, discovery time, popularity of use and signal strength. Type of antennas used and the environment’s complexity such as trees for outdoor and partitions for indoor spaces highly affect the scanning range. The aforementioned challenges are empirically evaluated by “real” experiments using Bluetooth and Wi-Fi Scanners. The issues related to the antenna characteristics are also highlighted by experimenting with different antenna types. Novel scanning approaches including Overlapped Zones and Single Point Multi-Range detection methods will be then presented and verified by real-world tests. These novel techniques will be applied for location identification of the MAC IDs captured that can extract more information about people movement dynamics.

Dynamic Bayesian network inferencing for non-homogeneous complex systems

Dynamic Bayesian Networks (DBNs) provide a versatile platform for predicting and analysing the behaviour of complex systems. As such, they are well suited to the prediction of complex ecosystem population trajectories under anthropogenic disturbances such as the dredging of marine seagrass ecosystems. However, DBNs assume a homogeneous Markov chain whereas a key characteristics of complex ecosystems is the presence of feedback loops, path dependencies and regime changes whereby the behaviour of the system can vary based on past states. This paper develops a method based on the small world structure of complex systems networks to modularise a non-homogeneous DBN and enable the computation of posterior marginal probabilities given evidence in forwards inference. It also provides an approach for an approximate solution for backwards inference as convergence is not guaranteed for a path dependent system. When applied to the seagrass dredging problem, the incorporation of path dependency can implement conditional absorption and allows release from the zero state in line with environmental and ecological observations. As dredging has a marked global impact on seagrass and other marine ecosystems of high environmental and economic value, using such a complex systems model to develop practical ways to meet the needs of conservation and industry through enhancing resistance and/or recovery is of paramount importance.

Variability in exposure to ambient ultrafine particles in urban schools: Comparative assessment between Australia and Spain

Ambient ultrafine particle number concentrations (PNC) have inhomogeneous spatio-temporal distributions and depend on a number of different urban factors, including background conditions and distant sources. This paper quantitatively compares exposure to ambient ultrafine particles at urban schools in two cities in developed countries, with high insolation climatic conditions, namely Brisbane (Australia) and Barcelona (Spain). The analysis used comprehensive indoor and outdoor air quality measurements at 25 schools in Brisbane and 39 schools in Barcelona. PNC modes were analysed with respect to ambient temperature, land use and urban characteristics, combined with the measured elemental carbon concentrations, NOx (Brisbane) and NO2 (Barcelona). The trends and modes of the quantified weekday average daily cycles of ambient PNC exhibited significant differences between the two cities. PNC increases were observed during traffic rush hours in both cases. However, the mid-day peak was dominant in Brisbane schools and had the highest contribution to total PNC for both indoors and outdoors. In Barcelona, the contribution from traffic was highest for ambient PNC, while the mid-day peak had a slightly higher contribution for indoor concentrations. Analysis of the relationships between PNC and land use characteristics in Barcelona schools showed a moderate correlation with the percentage of road network area and an anti-correlation with the percentage of green area. No statistically significant correlations were found for Brisbane. Overall, despite many similarities between the two cities, school-based exposure patterns were different. The main source of ambient PNC at schools was shown to be traffic in Barcelona and mid-day new particle formation in Brisbane. The mid-day PNC peak in Brisbane could have been driven by the combined effect of background and meteorological conditions, as well as other local/distant sources. The results have implications for urban development, especially in terms of air quality mitigation and management at schools.