Marginal reversible jump Markov chain Monte Carlo with application to motor unit number estimation

Motor unit number estimation (MUNE) is a method which aims to provide a quantitative indicator of progression of diseases that lead to loss of motor units, such as motor neurone disease. However the development of a reliable, repeatable and fast real-time MUNE method has proved elusive hitherto. Ridall et al. (2007) implement a reversible jump Markov chain Monte Carlo (RJMCMC) algorithm to produce a posterior distribution for the number of motor units using a Bayesian hierarchical model that takes into account biological information about motor unit activation. However we find that the approach can be unreliable for some datasets since it can suffer from poor cross-dimensional mixing. Here we focus on improved inference by marginalising over latent variables to create the likelihood. In particular we explore how this can improve the RJMCMC mixing and investigate alternative approaches that utilise the likelihood (e.g. DIC (Spiegelhalter et al., 2002)). For this model the marginalisation is over latent variables which, for a larger number of motor units, is an intractable summation over all combinations of a set of latent binary variables whose joint sample space increases exponentially with the number of motor units. We provide a tractable and accurate approximation for this quantity and also investigate simulation approaches incorporated into RJMCMC using results of Andrieu and Roberts (2009).

Surgical treatment approaches and reimbursement costs of surgical site infections post hip arthroplasty in Australia : a retrospective analysis

BACKGROUND: The treatment for deep surgical site infection (SSI) following primary total hip arthroplasty (THA) varies internationally and it is at present unclear which treatment approaches are used in Australia. The aim of this study is to identify current treatment approaches in Queensland, Australia, show success rates and quantify the costs of different treatments. METHODS: Data for patients undergoing primary THA and treatment for infection between January 2006 and December 2009 in Queensland hospitals were extracted from routinely used hospital databases. Records were linked with pathology information to confirm positive organisms. Diagnosis and treatment of infection was determined using ICD-10-AM and ACHI codes, respectively. Treatment costs were estimated based on AR-DRG cost accounting codes assigned to each patient hospital episode. RESULTS: A total of n=114 patients with deep surgical site infection were identified. The majority of patients (74%) were first treated with debridement, antibiotics and implant retention (DAIR), which was successful in eradicating the infection in 60.3% of patients with an average cost of $13,187. The remaining first treatments were 1-stage revision, successful in 89.7% with average costs of $27,006, and 2-stage revisions, successful in 92.9% of cases with average costs of $42,772. Multiple treatments following 'failed DAIR' cost on average $29,560, for failed 1-stage revision were $24,357, for failed 2-stage revision were $70,381 and were $23,805 for excision arthroplasty. CONCLUSIONS: As treatment costs in Australia are high primary prevention is important and the economics of competing treatment choices should be carefully considered. These currently vary greatly across international settings.

Pseudo-marginal algorithms with multiple CPUs

The emergence of pseudo-marginal algorithms has led to improved computational efficiency for dealing with complex Bayesian models with latent variables. Here an unbiased estimator of the likelihood replaces the true likelihood in order to produce a Bayesian algorithm that remains on the marginal space of the model parameter (with latent variables integrated out), with a target distribution that is still the correct posterior distribution. Very efficient proposal distributions can be developed on the marginal space relative to the joint space of model parameter and latent variables. Thus psuedo-marginal algorithms tend to have substantially better mixing properties. However, for pseudo-marginal approaches to perform well, the likelihood has to be estimated rather precisely. This can be difficult to achieve in complex applications. In this paper we propose to take advantage of multiple central processing units (CPUs), that are readily available on most standard desktop computers. Here the likelihood is estimated independently on the multiple CPUs, with the ultimate estimate of the likelihood being the average of the estimates obtained from the multiple CPUs. The estimate remains unbiased, but the variability is reduced. We compare and contrast two different technologies that allow the implementation of this idea, both of which require a negligible amount of extra programming effort. The superior performance of this idea over the standard approach is demonstrated on simulated data from a stochastic volatility model.

Lifecycle costs of ultracapacitors in electric vehicle applications

The pulse power characteristics of ultracapacitors appear well suited to electric vehicle applications, where they may supply the peak power more efficiently than the battery, and can prevent excessive over sizing of the battery pack due to peak power demands. Operation of ultracapacitors in battery electric vehicles is examined for possible improvements in system efficiency, vehicle driving range, battery pack lifetime, and potential reductions in system lifecycle cost. The lifecycle operation of these ultracapacitors is simulated using custom-built, dynamic simulation code constructed in Matlab. Despite apparent gains in system efficiency and driving range, the results strongly suggest that the inclusion of ultracapacitors in the electric vehicle does not make sense from a lifecycle cost perspective. Furthermore, a comparison with results from earlier work shows that this outcome is highly dependant upon the efficiency and cost of the battery under consideration. However, it is likely that the lifecycle cost benefits of ultracapacitors in these electric vehicles would be, at most, marginal and do not justify the additional capital costs and system complexity that would be incurred in the vehicle

Simulated lifecycle costs of ultracapacitors in battery electric vehicles

The pulse power characteristics of ultracapacitors appear well suited to electric vehicle applications, where they may supply the peak power more efficiently than the battery, and can prevent excessive over sizing of the battery pack due to peak power demands. Operation of ultracapacitors in battery electric vehicles (BEVs) is examined for possible improvements in system efficiency, vehicle driving range, battery pack lifetime, and potential reductions in system lifecycle cost. The lifecycle operation of these ultracapacitors is simulated using a custom-built, dynamic simulation code constructed in Matlab. Despite apparent gains in system efficiency and driving range, the lifecycle cost benefits as simulated appear to be marginal, and are heavily influenced by the incremental cost of power components. However, additional factors are identified which, in reality, will drive ultracapacitors towards viability in electric vehicle applications.

The health outcomes and costs of people attending an interdisciplinary chronic disease service in regional Australia : protocol for a longitudinal cohort investigation

Background Rates of chronic disease are escalating around the world. To date health service evaluations have focused on interventions for single chronic diseases. However, evaluations of the effectiveness of new intervention strategies that target single chronic diseases as well as multimorbidity are required, particularly in areas outside major metropolitan centres where access to services, such as specialist care, is difficult and where the retention and recruitment of health professionals affects service provision. Methods This study is a longitudinal investigation with a baseline and three follow-up assessments comparing the health and health costs of people with chronic disease before and after intervention at a chronic disease clinic, in regional Australia. The clinic is led by students under the supervision of health professionals. The study will provide preliminary evidence regarding the effectiveness of the intervention, and evaluate the influence of a range of factors on the health outcomes and costs of the patients attending the clinic. Patients will be evaluated at baseline (intake to the service), and at 3-, 6-, and 12-months after intake to the service. Health will be measured using the SF-36 and health costs will be measured using government and medical record sources. The intervention involves students and health professionals from multiple professions working together to treat patients with programs that include education and exercise therapy programs for back pain, and Healthy Lifestyle programs; as well as individual consultations involving single professions. Discussion Understanding the effect of a range of factors on the health state and health costs of people attending an interdisciplinary clinic will inform health service provision for this clinical group and will determine which factors need to be controlled for in future observational studies. Preliminary evidence regarding changes in health and health costs associated with the intervention will be a platform for future clinical trials of intervention effectiveness. The results will be of interest to teams investigating new chronic disease programs particularly for people with multimorbidity, and in areas outside major metropolitan centres.

Distribution tap changer adjustment to improve small-scale embedded generator penetration and mitigate voltage rise

This paper addresses the voltage rise constraints that are initiated from increased renewable generation resources in low voltage distribution networks. In this paper, an approach which is able to mitigate these voltage rise constraints and allow for increased distributed generator penetration is presented. The proposed approach involves utilizing the distribution transformers static tap changer to reduce the distribution feeder voltage setpoint. The proposed approach is modeled on a generic low voltage distribution network using the PSS SINCAL© simulation software package and is also implemented in a real low voltage distribution network to verify its practicality. Results indicate that this approach can be implemented to mitigate the voltage rise constraint and increase small-scale embedded generator penetration in a high proportion of low voltage feeders while avoiding any substantial network costs.

Loss reduction experiences in electric power distribution companies of Iran

The experiences of the loss reduction projects in electric power distribution companies (EPDCs) of Iran are presented. The loss reduction methods, which are proposed individually by 14 EPDCs, corresponding energy saving (ES), Investment costs (IC), and loss rate reductions are provided. In order to illustrate the effectiveness and performance of the loss reduction methods, three parameters are proposed as energy saving per investment costs (ESIC), energy saving per quantity (ESPQ), and investment costs per quantity (ICPQ). The overall ESIC of 14 EPDC as well as individual average and standard deviation of the EISC for each method is presented and compared. In addition, the average and standard deviation of the ESPQs and ICPQs for the loss reduction methods, individually, are provided and investigated. These parameters are useful for EPDCs that intend to reduce the electric losses in distribution networks as a benchmark and as a background in the planning purposes.

Load profile based determination of distribution feeder configuration by dynamic programming

This paper presents a new method to determine feeder reconfiguration scheme considering variable load profile. The objective function consists of system losses, reliability costs and also switching costs. In order to achieve an optimal solution the proposed method compares these costs dynamically and determines when and how it is reasonable to have a switching operation. The proposed method divides a year into several equal time periods, then using particle swarm optimization (PSO), optimal candidate configurations for each period are obtained. System losses and customer interruption cost of each configuration during each period is also calculated. Then, considering switching cost from a configuration to another one, dynamic programming algorithm (DPA) is used to determine the annual reconfiguration scheme. Several test systems were used to validate the proposed method. The obtained results denote that to have an optimum solution it is necessary to compare operation costs dynamically.

Benefit evaluation for information providing system by WebGIS in view of users' travel costs

This paper estimates the benefit of a plan for information providing system on road administration by WebGIS. The system will reduce travel costs of visitors from their business establishments to a road administration section of a city office. The authors had individual interviews with the visitors at the section of the Ichikawa City Office. Annual total sum of travel costs was estimated at 37 million yen at most. This paper also proposes formulas which expect the frequency of visits or the total sum of travel costs from the spatial distribution of the business establishments without questionnaires.

Accounting for management costs in sensitivity analyses of matrix population models

Traditional sensitivity and elasticity analyses of matrix population models have been used to inform management decisions, but they ignore the economic costs of manipulating vital rates. For example, the growth rate of a population is often most sensitive to changes in adult survival rate, but this does not mean that increasing that rate is the best option for managing the population because it may be much more expensive than other options. To explore how managers should optimize their manipulation of vital rates, we incorporated the cost of changing those rates into matrix population models. We derived analytic expressions for locations in parameter space where managers should shift between management of fecundity and survival, for the balance between fecundity and survival management at those boundaries, and for the allocation of management resources to sustain that optimal balance. For simple matrices, the optimal budget allocation can often be expressed as simple functions of vital rates and the relative costs of changing them. We applied our method to management of the Helmeted Honeyeater (Lichenostomus melanops cassidix; an endangered Australian bird) and the koala (Phascolarctos cinereus) as examples. Our method showed that cost-efficient management of the Helmeted Honeyeater should focus on increasing fecundity via nest protection, whereas optimal koala management should focus on manipulating both fecundity and survival simultaneously. These findings are contrary to the cost-negligent recommendations of elasticity analysis, which would suggest focusing on managing survival in both cases. A further investigation of Helmeted Honeyeater management options, based on an individual-based model incorporating density dependence, spatial structure, and environmental stochasticity, confirmed that fecundity management was the most cost-effective strategy. Our results demonstrate that decisions that ignore economic factors will reduce management efficiency. ©2006 Society for Conservation Biology.

Depression in Working Adults: Comparing the Costs and Health Outcomes of Working When Ill

Objective Working through a depressive illness can improve mental health but also carries risks and costs from reduced concentration, fatigue, and poor on-the-job performance. However, evidence-based recommendations for managing work attendance decisions, which benefit individuals and employers, are lacking. Therefore, this study has compared the costs and health outcomes of short-term absenteeism versus working while ill (“presenteeism”) amongst employed Australians reporting lifetime major depression. Methods Cohort simulation using state-transition Markov models simulated movement of a hypothetical cohort of workers, reporting lifetime major depression, between health states over one- and five-years according to probabilities derived from a quality epidemiological data source and existing clinical literature. Model outcomes were health service and employment-related costs, and quality-adjusted-life-years (QALYs), captured for absenteeism relative to presenteeism, and stratified by occupation (blue versus white-collar). Results Per employee with depression, absenteeism produced higher mean costs than presenteeism over one- and five-years ($42,573/5-years for absenteeism, $37,791/5-years for presenteeism). However, overlapping confidence intervals rendered differences non-significant. Employment-related costs (lost productive time, job turnover), and antidepressant medication and service use costs of absenteeism and presenteeism were significantly higher for white-collar workers. Health outcomes differed for absenteeism versus presenteeism amongst white-collar workers only. Conclusions Costs and health outcomes for absenteeism and presenteeism were not significantly different; service use costs excepted. Significant variation by occupation type was identified. These findings provide the first occupation-specific cost evidence which can be used by clinicians, employees, and employers to review their management of depression-related work attendance, and may suggest encouraging employees to continue working is warranted.

Integral Formulation of the Problem of Current Distribution in Compulsator Wires of Electromagnetic Launchers and Railguns

Compulsators are power sources of choice for use in electromagnetic launchers and railguns. These devices hold the promise of reducing unit costs of payload to orbit. In an earlier work, the author had calculated the current distribution in compulsator wires by considering the wire to be split into a finite number of separate wires. The present work develops an integral formulation of the problem of current distribution in compulsator wires which leads to an integrodifferential equation. Analytical solutions, including those for the integration constants, are obtained in closed form. The analytical solutions present a much clearer picture of the effect of various input parameters on the cross-sectional current distribution and point to ways in which the desired current density distribution can be achieved. Results are graphically presented and discussed, with particular reference to a 50-kJ compulsator in Bangalore. Finite-element analysis supports the results.

Accelerating Pseudo-Marginal MCMC using Gaussian Processes

Pseudo-marginal methods such as the grouped independence Metropolis-Hastings (GIMH) and Markov chain within Metropolis (MCWM) algorithms have been introduced in the literature as an approach to perform Bayesian inference in latent variable models. These methods replace intractable likelihood calculations with unbiased estimates within Markov chain Monte Carlo algorithms. The GIMH method has the posterior of interest as its limiting distribution, but suffers from poor mixing if it is too computationally intensive to obtain high-precision likelihood estimates. The MCWM algorithm has better mixing properties, but less theoretical support. In this paper we propose to use Gaussian processes (GP) to accelerate the GIMH method, whilst using a short pilot run of MCWM to train the GP. Our new method, GP-GIMH, is illustrated on simulated data from a stochastic volatility and a gene network model.

The impact of home loan key facts sheets on borrowers' comparisons of loan costs

This thesis evaluates the effectiveness of the prescribed design and distribution requirements of the Australian Government's home loan key facts sheets (KFS) aimed at helping borrowers compare loan costs. The findings show that despite effectively improving borrower decision-making, few borrowers were aware of their existence and function. It was also demonstrated that KFS have had limited market impact over the four year window since introduction, likely due to the requirement that KFS provision is not required unless formally requested by a borrower. Recommendations include transferring the burden of disclosure to lenders in the first instance to address this information gap.