Regression and hypothesis tests for multivariate GNSS state time series

A satellite based observation system can continuously or repeatedly generate a user state vector time series that may contain useful information. One typical example is the collection of International GNSS Services (IGS) station daily and weekly combined solutions. Another example is the epoch-by-epoch kinematic position time series of a receiver derived by a GPS real time kinematic (RTK) technique. Although some multivariate analysis techniques have been adopted to assess the noise characteristics of multivariate state time series, statistic testings are limited to univariate time series. After review of frequently used hypotheses test statistics in univariate analysis of GNSS state time series, the paper presents a number of T-squared multivariate analysis statistics for use in the analysis of multivariate GNSS state time series. These T-squared test statistics have taken the correlation between coordinate components into account, which is neglected in univariate analysis. Numerical analysis was conducted with the multi-year time series of an IGS station to schematically demonstrate the results from the multivariate hypothesis testing in comparison with the univariate hypothesis testing results. The results have demonstrated that, in general, the testing for multivariate mean shifts and outliers tends to reject less data samples than the testing for univariate mean shifts and outliers under the same confidence level. It is noted that neither univariate nor multivariate data analysis methods are intended to replace physical analysis. Instead, these should be treated as complementary statistical methods for a prior or posteriori investigations. Physical analysis is necessary subsequently to refine and interpret the results.

Groundwater visualisation system (GVS) : a software framework for integrated display and interrogation of conceptual hydrogeological models, data and time-series animation

Management of groundwater systems requires realistic conceptual hydrogeological models as a framework for numerical simulation modelling, but also for system understanding and communicating this to stakeholders and the broader community. To help overcome these challenges we developed GVS (Groundwater Visualisation System), a stand-alone desktop software package that uses interactive 3D visualisation and animation techniques. The goal was a user-friendly groundwater management tool that could support a range of existing real-world and pre-processed data, both surface and subsurface, including geology and various types of temporal hydrological information. GVS allows these data to be integrated into a single conceptual hydrogeological model. In addition, 3D geological models produced externally using other software packages, can readily be imported into GVS models, as can outputs of simulations (e.g. piezometric surfaces) produced by software such as MODFLOW or FEFLOW. Boreholes can be integrated, showing any down-hole data and properties, including screen information, intersected geology, water level data and water chemistry. Animation is used to display spatial and temporal changes, with time-series data such as rainfall, standing water levels and electrical conductivity, displaying dynamic processes. Time and space variations can be presented using a range of contouring and colour mapping techniques, in addition to interactive plots of time-series parameters. Other types of data, for example, demographics and cultural information, can also be readily incorporated. The GVS software can execute on a standard Windows or Linux-based PC with a minimum of 2 GB RAM, and the model output is easy and inexpensive to distribute, by download or via USB/DVD/CD. Example models are described here for three groundwater systems in Queensland, northeastern Australia: two unconfined alluvial groundwater systems with intensive irrigation, the Lockyer Valley and the upper Condamine Valley, and the Surat Basin, a large sedimentary basin of confined artesian aquifers. This latter example required more detail in the hydrostratigraphy, correlation of formations with drillholes and visualisation of simulation piezometric surfaces. Both alluvial system GVS models were developed during drought conditions to support government strategies to implement groundwater management. The Surat Basin model was industry sponsored research, for coal seam gas groundwater management and community information and consultation. The “virtual” groundwater systems in these 3D GVS models can be interactively interrogated by standard functions, plus production of 2D cross-sections, data selection from the 3D scene, rear end database and plot displays. A unique feature is that GVS allows investigation of time-series data across different display modes, both 2D and 3D. GVS has been used successfully as a tool to enhance community/stakeholder understanding and knowledge of groundwater systems and is of value for training and educational purposes. Projects completed confirm that GVS provides a powerful support to management and decision making, and as a tool for interpretation of groundwater system hydrological processes. A highly effective visualisation output is the production of short videos (e.g. 2–5 min) based on sequences of camera ‘fly-throughs’ and screen images. Further work involves developing support for multi-screen displays and touch-screen technologies, distributed rendering, gestural interaction systems. To highlight the visualisation and animation capability of the GVS software, links to related multimedia hosted online sites are included in the references.

Condition-based maintenance optimisation without a predetermined strategy structure for a two-component series system

Most existing research on maintenance optimisation for multi-component systems only considers the lifetime distribution of the components. When the condition-based maintenance (CBM) strategy is adopted for multi-component systems, the strategy structure becomes complex due to the large number of component states and their combinations. Consequently, some predetermined maintenance strategy structures are often assumed before the maintenance optimisation of a multi-component system in a CBM context. Developing these predetermined strategy structure needs expert experience and the optimality of these strategies is often not proofed. This paper proposed a maintenance optimisation method that does not require any predetermined strategy structure for a two-component series system. The proposed method is developed based on the semi-Markov decision process (SMDP). A simulation study shows that the proposed method can identify the optimal maintenance strategy adaptively for different maintenance costs and parameters of degradation processes. The optimal maintenance strategy structure is also investigated in the simulation study, which provides reference for further research in maintenance optimisation of multi-component systems.

Exact and approximate Bayesian inference for low count time series models with intractable likelihoods

In this paper we present a new simulation methodology in order to obtain exact or approximate Bayesian inference for models for low-valued count time series data that have computationally demanding likelihood functions. The algorithm fits within the framework of particle Markov chain Monte Carlo (PMCMC) methods. The particle filter requires only model simulations and, in this regard, our approach has connections with approximate Bayesian computation (ABC). However, an advantage of using the PMCMC approach in this setting is that simulated data can be matched with data observed one-at-a-time, rather than attempting to match on the full dataset simultaneously or on a low-dimensional non-sufficient summary statistic, which is common practice in ABC. For low-valued count time series data we find that it is often computationally feasible to match simulated data with observed data exactly. Our particle filter maintains $N$ particles by repeating the simulation until $N+1$ exact matches are obtained. Our algorithm creates an unbiased estimate of the likelihood, resulting in exact posterior inferences when included in an MCMC algorithm. In cases where exact matching is computationally prohibitive, a tolerance is introduced as per ABC. A novel aspect of our approach is that we introduce auxiliary variables into our particle filter so that partially observed and/or non-Markovian models can be accommodated. We demonstrate that Bayesian model choice problems can be easily handled in this framework.

Discussion paper and working paper series : Motivations, expectations and experiences of Australian rural and regional planners

Despite playing an extremely important role in shaping communities, the role and contribution of planners is not widely understood or acknowledged. At the same time, there is a shortage of planners in Australia, especially in non-urban areas. Thus, though an online survey of 185 rural and regional planners, this research explores their motivations, expectations and experiences. Most enjoyed and felt confident in their role, explaining that they valued the relaxed family orientated rural lifestyle and the varied nature of the planning work. Although they sometimes felt isolated, the non-urban location provided quicker progression to senior roles, the ability to engage directly with the community and to see the consequences of their decisions. Only half felt their education had prepared them well for their role, citing gaps in terms of computerised modelling, team leadership and conflict resolution skills. Their feedback centred on providing a more practical course, focussing more on regional planning, and encouraging urban and rural experience placements. As the first study to quantifiably explore rural and regional Australian planners perceptions of their role and challenges, the findings illustrate current experiences, key planning challenges, perceived educational gaps and future priorities.

Vibrational spectroscopic characterization of the phosphate mineral series eosphorite–childrenite–(Mn,Fe)Al(PO4)(OH)2·(H2O)

The phosphate mineral series eosphorite–childrenite–(Mn,Fe)Al(PO4)(OH)2·(H2O) has been studied using a combination of electron probe analysis and vibrational spectroscopy. Eosphorite is the manganese rich mineral with lower iron content in comparison with the childrenite which has higher iron and lower manganese content. The determined formulae of the two studied minerals are: (Mn0.72,Fe0.13,Ca0.01)(Al)1.04(PO4, OHPO3)1.07(OH1.89,F0.02)·0.94(H2O) for SAA-090 and (Fe0.49,Mn0.35,Mg0.06,Ca0.04)(Al)1.03(PO4, OHPO3)1.05(OH)1.90·0.95(H2O) for SAA-072. Raman spectroscopy enabled the observation of bands at 970 cm−1 and 1011 cm−1 assigned to monohydrogen phosphate, phosphate and dihydrogen phosphate units. Differences are observed in the area of the peaks between the two eosphorite minerals. Raman bands at 562 cm−1, 595 cm−1, and 608 cm−1 are assigned to the �4 bending modes of the PO4, HPO4 and H2PO4 units; Raman bands at 405 cm−1, 427 cm−1 and 466 cm−1 are attributed to the �2 modes of these units. Raman bands of the hydroxyl and water stretching modes are observed. Vibrational spectroscopy enabled details of the molecular structure of the eosphorite mineral series to be determined.

Raman, infrared and near-infrared spectroscopic characterization of the herderite-hydroxylherderite mineral series

Natural single-crystal specimens of the herderite-hydroxylherderite series from Brazil, with general formula CaBePO4(F,OH), were investigated by electron microprobe, Raman, infrared and near-infrared spectroscopies. The minerals occur as secondary products in granitic pegmatites. Herderite and hydroxylherderite minerals show extensive solid solution formation. The Raman spectra of hydroxylherderite are characterized by bands at around 985 and 998 cm-1, assigned to ν1 symmetric stretching mode of the HOPO33- and PO43- units. Raman bands at around 1085, 1128 and 1138 cm-1 are attributed to both the HOP and PO antisymmetric stretching vibrations. The set of Raman bands observed at 563, 568, 577, 598, 616 and 633 cm-1 are assigned to the ν4 out of plane bending modes of the PO4 and H2PO4 units. The OH Raman stretching vibrations of hydroxylherderite were observed ranging from 3626 cm-1 to 3609 cm-1. The infrared stretching vibrations of hydroxylherderites were observed between 3606 cm-1 and 3599 cm-1. By using a Libowitzky type function, hydrogen bond distances based upon the OH stretching bands were calculated. Characteristic NIR bands at around 6961 and 7054 cm-1 were assigned to the first overtone of the fundamental, whilst NIR bands at 10194 and 10329 cm-1 are assigned to the second overtone of the fundamental OH stretching vibration. Insight into the structure of the herderite-hydroxylherderite series is assessed by vibrational spectroscopy.

Econometric Modelling with Time Series : Specification, Estimation and Testing

This book provides a general framework for specifying, estimating, and testing time series econometric models. Special emphasis is given to estimation by maximum likelihood, but other methods are also discussed, including quasi-maximum likelihood estimation, generalized method of moments estimation, nonparametric estimation, and estimation by simulation. An important advantage of adopting the principle of maximum likelihood as the unifying framework for the book is that many of the estimators and test statistics proposed in econometrics can be derived within a likelihood framework, thereby providing a coherent vehicle for understanding their properties and interrelationships. In contrast to many existing econometric textbooks, which deal mainly with the theoretical properties of estimators and test statistics through a theorem-proof presentation, this book squarely addresses implementation to provide direct conduits between the theory and applied work.

Spatiotemporal model or time series model for assessing city-wide temperature effects on mortality?

Most studies examining the temperature–mortality association in a city used temperatures from one site or the average from a network of sites. This may cause measurement error as temperature varies across a city due to effects such as urban heat islands. We examined whether spatiotemporal models using spatially resolved temperatures produced different associations between temperature and mortality compared with time series models that used non-spatial temperatures. We obtained daily mortality data in 163 areas across Brisbane city, Australia from 2000 to 2004. We used ordinary kriging to interpolate spatial temperature variation across the city based on 19 monitoring sites. We used a spatiotemporal model to examine the impact of spatially resolved temperatures on mortality. Also, we used a time series model to examine non-spatial temperatures using a single site and the average temperature from three sites. We used squared Pearson scaled residuals to compare model fit. We found that kriged temperatures were consistent with observed temperatures. Spatiotemporal models using kriged temperature data yielded slightly better model fit than time series models using a single site or the average of three sites' data. Despite this better fit, spatiotemporal and time series models produced similar associations between temperature and mortality. In conclusion, time series models using non-spatial temperatures were equally good at estimating the city-wide association between temperature and mortality as spatiotemporal models.

Sunshine, rainfall, humidity and child pneumonia in the tropics : time-series analyses

Few studies have formally examined the relationship between meteorological factors and the incidence of child pneumonia in the tropics, despite the fact that most child pneumonia deaths occur there. We examined the association between four meteorological exposures (rainy days, sunshine, relative humidity, temperature) and the incidence of clinical pneumonia in young children in the Philippines using three time-series methods: correlation of seasonal patterns, distributed lag regression, and case-crossover. Lack of sunshine was most strongly associated with pneumonia in both lagged regression [overall relative risk over the following 60 days for a 1-h increase in sunshine per day was 0·67 (95% confidence interval (CI) 0·51–0·87)] and case-crossover analysis [odds ratio for a 1-h increase in mean daily sunshine 8–14 days earlier was 0·95 (95% CI 0·91–1·00)]. This association is well known in temperate settings but has not been noted previously in the tropics. Further research to assess causality is needed.

Assessing the short-term effects of heatwaves on mortality and morbidity in Brisbane, Australia : comparison of case-crossover and time series analyses

Background Heat-related impacts may have greater public health implications as climate change continues. It is important to appropriately characterize the relationship between heatwave and health outcomes. However, it is unclear whether a case-crossover design can be effectively used to assess the event- or episode-related health effects. This study examined the association between exposure to heatwaves and mortality and emergency hospital admissions (EHAs) from non-external causes in Brisbane, Australia, using both case-crossover and time series analyses approaches. Methods Poisson generalised additive model (GAM) and time-stratified case-crossover analyses were used to assess the short-term impact of heatwaves on mortality and EHAs. Heatwaves exhibited a significant impact on mortality and EHAs after adjusting for air pollution, day of the week, and season. Results For time-stratified case-crossover analysis, odds ratios of mortality and EHAs during heatwaves were 1.62 (95% confidence interval (CI): 1.36–1.94) and 1.22 (95% CI: 1.14–1.30) at lag 1, respectively. Time series GAM models gave similar results. Relative risks of mortality and EHAs ranged from 1.72 (95% CI: 1.40–2.11) to 1.81 (95% CI: 1.56–2.10) and from 1.14 (95% CI: 1.06–1.23) to 1.28 (95% CI: 1.21–1.36) at lag 1, respectively. The risk estimates gradually attenuated after the lag of one day for both case-crossover and time series analyses. Conclusions The risk estimates from both case-crossover and time series models were consistent and comparable. This finding may have implications for future research on the assessment of event- or episode-related (e.g., heatwave) health effects.

Effects of temperature on mortality in Chiang Mai city, Thailand : a time series study

Background The association between temperature and mortality has been examined mainly in North America and Europe. However, less evidence is available in developing countries, especially in Thailand. In this study, we examined the relationship between temperature and mortality in Chiang Mai city, Thailand, during 1999–2008. Method A time series model was used to examine the effects of temperature on cause-specific mortality (non-external, cardiopulmonary, cardiovascular, and respiratory) and age-specific non-external mortality (<=64, 65–74, 75–84, and > =85 years), while controlling for relative humidity, air pollution, day of the week, season and long-term trend. We used a distributed lag non-linear model to examine the delayed effects of temperature on mortality up to 21 days. Results We found non-linear effects of temperature on all mortality types and age groups. Both hot and cold temperatures resulted in immediate increase in all mortality types and age groups. Generally, the hot effects on all mortality types and age groups were short-term, while the cold effects lasted longer. The relative risk of non-external mortality associated with cold temperature (19.35°C, 1st percentile of temperature) relative to 24.7°C (25th percentile of temperature) was 1.29 (95% confidence interval (CI): 1.16, 1.44) for lags 0–21. The relative risk of non-external mortality associated with high temperature (31.7°C, 99th percentile of temperature) relative to 28°C (75th percentile of temperature) was 1.11 (95% CI: 1.00, 1.24) for lags 0–21. Conclusion This study indicates that exposure to both hot and cold temperatures were related to increased mortality. Both cold and hot effects occurred immediately but cold effects lasted longer than hot effects. This study provides useful data for policy makers to better prepare local responses to manage the impact of hot and cold temperatures on population health.

Surgical resection for non-small cell lung cancer : clinical features and outcomes for a consecutive series at an Irish tertiary referral centre

Background: Few patients diagnosed with lung cancer are still alive 5 years after diagnosis. The aim of the current study was to conduct a 10-year review of a consecutive series of patients undergoing curative-intent surgical resection at the largest tertiary referral centre to identify prognostic factors. Methods: Case records of all patients operated on for lung cancer between 1998 and 2008 were reviewed. The clinical features and outcomes of all patients with non-small cell lung cancer (NSCLC) stage I-IV were recorded. Results: A total of 654 patients underwent surgical resection with curative intent during the study period. Median overall survival for the entire cohort was 37 months. The median age at operation was 66 years, with males accounting for 62.7 %. Squamous cell type was the most common histological subtype, and lobectomies were performed in 76.5 % of surgical resections. Pneumonectomy rates decreased significantly in the latter half of the study (25 vs. 16.3 %), while sub-anatomical resection more than doubled (2 vs. 5 %) (p < 0.005). Clinico-pathological characteristics associated with improved survival by univariate analysis include younger age, female sex, smaller tumour size, smoking status, lobectomy, lower T and N status and less advanced pathological stage. Age, gender, smoking status and tumour size, as well as T and N descriptors have emerged as independent prognostic factors by multivariate analysis. Conclusion: We identified several factors that predicted outcome for NSCLC patients undergoing curative-intent surgical resection. Survival rates in our series are comparable to those reported from other thoracic surgery centres. © 2012 Royal Academy of Medicine in Ireland.

Made by motion, Series one

“Made by Motion” is a collaboration between digital artist Paul Van Opdenbosch and performer and choreographer Elise May; a series of studies on captured motion data used to generating experimental visual forms that reverberate in space and time. The project investigates the invisible forces generated by and influencing the movement of a dancer. Along with how the forces can be captured and applied to generating visual outcomes that surpass simple data visualisation, projecting the intent of the performer’s movements. The source or ‘seed’ comes from using an Xsens MVN - Inertial Motion Capture system to capture spontaneous dance movements, with the visual generation conducted through a customised dynamics simulation. In this first series the visual investigation focused on manipulating the movement date at the instance of capture, capture been the recording of three-dimensional movement as ‘seen’ by the hardware and ‘understood’ through the calibration of software. By repositioning the capture hardware on the body we can effectively change how the same sequence of movements is ‘seen’ by the motion capture system thus generating a different visual result from effetely identical movement. The outcomes from the experiments clearly demonstrates the effectiveness of using motion capture hardware as a creative tool to manipulate the perception of the capture subject, in this case been a sequence of dance movements. The creative work exhibited is a cross-section of the experiments conducted in practice with the first animated work (Movement A - Control) using the motion capture hardware in its default ‘normal’ configuration. Following this is the lower body moved to the upper body (Lb-Ub), right arm moved onto the left arm (Ra-La), right leg moved onto the left leg (Rl-Ll) and finally the left leg moved onto a object that is then held in the left hand (Ll-Pf (Lh)).

Targeting histone deacetylases for the treatment of disease : Epigenetics Review Series

The 'histone code' is a well-established hypothesis describing the idea that specific patterns of post-translational modifications to histones act like a molecular 'code' recognized and used by non-histone proteins to regulate specific chromatin functions. One modification, which has received significant attention, is that of histone acetylation. The enzymes that regulate this modification are described as lysine acetyltransferases or KATs, and histone deacetylases or HDACs. Due to their conserved catalytic domain HDACs have been actively targeted as a therapeutic target. The pro-inflammatory environment is increasingly being recognized as a critical element for both degenerative diseases and cancer. The present review will discuss the current knowledge surrounding the clinical potential and current development of histone deacetylases for the treatment of diseases for which a pro-inflammatory environment plays important roles, and the molecular mechanisms by which such inhibitors may play important functions in modulating the pro-inflammatory environment. © 2009 Foundation for Cellular and Molecular Medicine/Blackwell Publishing Ltd.