A general method to determine sampling windows for nonlinear mixed effects models with an application to population pharmacokinetic studies

Optimal design methods have been proposed to determine the best sampling times when sparse blood sampling is required in clinical pharmacokinetic studies. However, the optimal blood sampling time points may not be feasible in clinical practice. Sampling windows, a time interval for blood sample collection, have been proposed to provide flexibility in blood sampling times while preserving efficient parameter estimation. Because of the complexity of the population pharmacokinetic models, which are generally nonlinear mixed effects models, there is no analytical solution available to determine sampling windows. We propose a method for determination of sampling windows based on MCMC sampling techniques. The proposed method attains a stationary distribution rapidly and provides time-sensitive windows around the optimal design points. The proposed method is applicable to determine sampling windows for any nonlinear mixed effects model although our work focuses on an application to population pharmacokinetic models.

Simulation-based fully Bayesian experimental design for mixed effects models

In this paper, we present fully Bayesian experimental designs for nonlinear mixed effects models, in which we develop simulation-based optimal design methods to search over both continuous and discrete design spaces. Although Bayesian inference has commonly been performed on nonlinear mixed effects models, there is a lack of research into performing Bayesian optimal design for nonlinear mixed effects models that require searches to be performed over several design variables. This is likely due to the fact that it is much more computationally intensive to perform optimal experimental design for nonlinear mixed effects models than it is to perform inference in the Bayesian framework. In this paper, the design problem is to determine the optimal number of subjects and samples per subject, as well as the (near) optimal urine sampling times for a population pharmacokinetic study in horses, so that the population pharmacokinetic parameters can be precisely estimated, subject to cost constraints. The optimal sampling strategies, in terms of the number of subjects and the number of samples per subject, were found to be substantially different between the examples considered in this work, which highlights the fact that the designs are rather problem-dependent and require optimisation using the methods presented in this paper.

Effect of broadcast baiting on abundance patterns of red imported fire ants (Hymenoptera: Formicidae) and key local ant genera at long-term monitoring sites in Brisbane, Australia

In 2001, the red imported fire ant (Solenopsis invicta Buren) was identified in Brisbane, Australia. An eradication program involving broadcast bait treatment with two insect growth regulators and a metabolic inhibitor began in September of that year and is currently ongoing. To gauge the impacts of these treatments on local ant populations, we examined long-term monitoring data and quantified abundance patterns of S. invicta and common local ant genera using a linear mixed-effects model. For S. invicta, presence in pitfalls reduced over time to zero on every site. Significantly higher numbers of S. invicta workers were collected on high-density polygyne sites, which took longer to disinfest compared with monogyne and low-density polygyne sites. For local ants, nine genus groups of the 10 most common genera analyzed either increased in abundance or showed no significant trend. Five of these genus groups were significantly less abundant at the start of monitoring on high-density polygyne sites compared with monogyne and low-density polygyne sites. The genus Pheidole significantly reduced in abundance over time, suggesting that it was affected by treatment efforts. These results demonstrate that the treatment regime used at the time successfully removed S. invicta from these sites in Brisbane, and that most local ant genera were not seriously impacted by the treatment. These results have important implications for current and future prophylactic treatment efforts, and suggest that native ants remain in treated areas to provide some biological resistance to S. invicta.

Growth defect in Grg5 null mice is associated with reduced Ihh signaling in growth plates

Gene-targeted disruption of Grg5, a mouse homologue of Drosophila groucho (gro), results in postnatal growth retardation in mice. The growth defect, most striking in approximately half of the Grg5 null mice, occurs during the first 4-5 weeks of age, but most mice recover retarded growth later. We used the nonlinear mixed-effects model to fit the growth data of wild-type, heterozygous, and Grg5 null mice. On the basis of preliminary evidence suggesting an interaction between Grg5 and the transcription factor Cbfa1/Runx2, critical for skeletal development, we further investigated the skeleton in the mice. A long bone growth plate defect was identified, which included shorter zones of proliferative and hypertrophic chondrocytes and decreased trabecular bone formation. This decreased trabecular bone formation is likely caused by a reduced recruitment of osteoblasts into the growth plate region of Grg5 null mice. Like the growth defect, the growth plate and trabecular bone abnormality improved as the mice grew older. The growth plate defect was associated with reduced Indian hedgehog expression and signaling. We suggest that Grg5, a transcriptional coregulator, modulates the activities of transcription factors, such as Cbfa1/Runx2 in vivo to affect Ihh expression and the function of long bone growth plates.

Stroke education : promising effects on the health beliefs of those at risk

Background: The effect of patient education on reducing stroke has had mixed effects, raising questions about how to achieve optimal benefit. Because past evaluations have typically lacked an appropriate theoretical base, the design of past research may have missed important effects. --------- Method: This study used a social cognitive framework to identify variables that might change in response to education. A mixed design was used to evaluate two approaches to an intervention, both of which included education. Fifty seniors completed a measure of stroke knowledge and beliefs twice: before and after an intervention that was either standard (educational brochure plus activities that were not about stroke) or enhanced (educational brochure plus activities designed to enhance beliefs about stroke). Outcome measures were health beliefs, intention to exercise to reduce stroke, and stroke knowledge. --------- Results: Selected beliefs changed significantly over time but not differentially across conditions. Beliefs that changed were (a) perceived susceptibility to stroke and (b) perceived benefit of exercise to reduce risk. Benefit beliefs, in particular, were strongly and positively associated with intention to exercise. -------- Conclusion: Findings suggest that basic approaches to patient education may influence health beliefs. More effective stroke prevention programs may result from continued consideration of the role of health beliefs in such programs.

A pseudo-marginal sequential Monte Carlo algorithm for random effects models in Bayesian sequential design

A computationally efficient sequential Monte Carlo algorithm is proposed for the sequential design of experiments for the collection of block data described by mixed effects models. The difficulty in applying a sequential Monte Carlo algorithm in such settings is the need to evaluate the observed data likelihood, which is typically intractable for all but linear Gaussian models. To overcome this difficulty, we propose to unbiasedly estimate the likelihood, and perform inference and make decisions based on an exact-approximate algorithm. Two estimates are proposed: using Quasi Monte Carlo methods and using the Laplace approximation with importance sampling. Both of these approaches can be computationally expensive, so we propose exploiting parallel computational architectures to ensure designs can be derived in a timely manner. We also extend our approach to allow for model uncertainty. This research is motivated by important pharmacological studies related to the treatment of critically ill patients.

Systematic review and meta-analysis of the effects of exercise for those with cancer-related lymphedema

Objective: To evaluate the effects of exercise on cancer-related lymphedema and related symptoms, and to determine the need for those with lymphedema to wear compression during exercise. Data Sources: CINAHL, Cochrane, Ebscohost, MEDLINE, Pubmed, ProQuest Health and Medical Complete, ProQuest Nursing and Allied Health Source, Science Direct and SPORTDiscus databases were searched for trials published prior to 1 January, 2015. Study Selection: Randomised and non-randomised, controlled trials, and single group pre-post studies published in English-language were included. Twenty-one (exercise) and four (compression and exercise) studies met inclusion criteria. Data Extraction: Data was extracted into tabular format using predefined data fields by one reviewer and assessed for accuracy by a second reviewer. Study quality was evaluated using the Effective Public Health Practice Project assessment tool. Data Synthesis: Data was pooled using a random effects model to assess the effects of acute and long-term exercise on lymphedema and lymphedema-associated symptoms, with subgroup analyses for exercise mode and intervention length. There was no effect of exercise (acute or intervention) on lymphedema or associated symptoms with standardised mean differences from all analyses ranging between −0.2 and 0.1 (p-values ≥0.22). Findings from subgroup analyses for exercise mode (aerobic, resistance, mixed, other) and intervention duration (>12 weeks or ≤12 weeks) were consistent with these findings; that is, no effect on lymphedema or associated symptoms. There were too few studies evaluating the effect of compression during regular exercise to conduct a meta-analysis. Conclusions: Individuals with secondary lymphedema can safely participate in progressive, regular exercise without experiencing a worsening of lymphedema or related-symptoms. However, the results also do not suggest any improvements will occur in lymphedema. At present, there is insufficient evidence to support or refute the current clinical recommendation to wear compression garments during regular exercise.

Using Trauma Injury Severity Score (TRISS) variables to predict length of hospital stay following trauma in New Zealand

Aim – To develop and assess the predictive capabilities of a statistical model that relates routinely collected Trauma Injury Severity Score (TRISS) variables to length of hospital stay (LOS) in survivors of traumatic injury. Method – Retrospective cohort study of adults who sustained a serious traumatic injury, and who survived until discharge from Auckland City, Middlemore, Waikato, or North Shore Hospitals between 2002 and 2006. Cubic-root transformed LOS was analysed using two-level mixed-effects regression models. Results – 1498 eligible patients were identified, 1446 (97%) injured from a blunt mechanism and 52 (3%) from a penetrating mechanism. For blunt mechanism trauma, 1096 (76%) were male, average age was 37 years (range: 15-94 years), and LOS and TRISS score information was available for 1362 patients. Spearman’s correlation and the median absolute prediction error between LOS and the original TRISS model was ρ=0.31 and 10.8 days, respectively, and between LOS and the final multivariable two-level mixed-effects regression model was ρ=0.38 and 6.0 days, respectively. Insufficient data were available for the analysis of penetrating mechanism models. Conclusions – Neither the original TRISS model nor the refined model has sufficient ability to accurately or reliably predict LOS. Additional predictor variables for LOS and other indicators for morbidity need to be considered.

Contemporary New Zealand coefficients for the Trauma Injury Severity Score: TRISS(NZ)

Aims – To develop local contemporary coefficients for the Trauma Injury Severity Score in New Zealand, TRISS(NZ), and to evaluate their performance at predicting survival against the original TRISS coefficients. Methods – Retrospective cohort study of adults who sustained a serious traumatic injury, and who survived until presentation at Auckland City, Middlemore, Waikato, or North Shore Hospitals between 2002 and 2006. Coefficients were estimated using ordinary and multilevel mixed-effects logistic regression models. Results – 1735 eligible patients were identified, 1672 (96%) injured from a blunt mechanism and 63 (4%) from a penetrating mechanism. For blunt mechanism trauma, 1250 (75%) were male and average age was 38 years (range: 15-94 years). TRISS information was available for 1565 patients of whom 204 (13%) died. Area under the Receiver Operating Characteristic (ROC) curves was 0.901 (95%CI: 0.879-0.923) for the TRISS(NZ) model and 0.890 (95% CI: 0.866-0.913) for TRISS (P<0.001). Insufficient data were available to determine coefficients for penetrating mechanism TRISS(NZ) models. Conclusions – Both TRISS models accurately predicted survival for blunt mechanism trauma. However, TRISS(NZ) coefficients were statistically superior to TRISS coefficients. A strong case exists for replacing TRISS coefficients in the New Zealand benchmarking software with these updated TRISS(NZ) estimates.

Difference-based meta-analytic procedures for between-participant and/or within-participant designs : a tutorial review for sports and exercise scientists

The aim of this paper is to provide a contemporary summary of statistical and non-statistical meta-analytic procedures that have relevance to the type of experimental designs often used by sport scientists when examining differences/change in dependent measure(s) as a result of one or more independent manipulation(s). Using worked examples from studies on observational learning in the motor behaviour literature, we adopt a random effects model and give a detailed explanation of the statistical procedures for the three types of raw score difference-based analyses applicable to between-participant, within-participant, and mixed-participant designs. Major merits and concerns associated with these quantitative procedures are identified and agreed methods are reported for minimizing biased outcomes, such as those for dealing with multiple dependent measures from single studies, design variation across studies, different metrics (i.e. raw scores and difference scores), and variations in sample size. To complement the worked examples, we summarize the general considerations required when conducting and reporting a meta-analysis, including how to deal with publication bias, what information to present regarding the primary studies, and approaches for dealing with outliers. By bringing together these statistical and non-statistical meta-analytic procedures, we provide the tools required to clarify understanding of key concepts and principles.

Empirical Bayes method in the study of traffic safety via heterogeneous negative multinomial model

In the study of traffic safety, expected crash frequencies across sites are generally estimated via the negative binomial model, assuming time invariant safety. Since the time invariant safety assumption may be invalid, Hauer (1997) proposed a modified empirical Bayes (EB) method. Despite the modification, no attempts have been made to examine the generalisable form of the marginal distribution resulting from the modified EB framework. Because the hyper-parameters needed to apply the modified EB method are not readily available, an assessment is lacking on how accurately the modified EB method estimates safety in the presence of the time variant safety and regression-to-the-mean (RTM) effects. This study derives the closed form marginal distribution, and reveals that the marginal distribution in the modified EB method is equivalent to the negative multinomial (NM) distribution, which is essentially the same as the likelihood function used in the random effects Poisson model. As a result, this study shows that the gamma posterior distribution from the multivariate Poisson-gamma mixture can be estimated using the NM model or the random effects Poisson model. This study also shows that the estimation errors from the modified EB method are systematically smaller than those from the comparison group method by simultaneously accounting for the RTM and time variant safety effects. Hence, the modified EB method via the NM model is a generalisable method for estimating safety in the presence of the time variant safety and the RTM effects.

Adapting, renewing, subverting : what place the performing arts in the creative industries?

Creative Industries was adopted as a platform in the 90s by the Blair government in the UK to describe the convergence of the arts, media, communication and information technologies as a newly formed cluster, providing economic and cultural capital for the knowledge economy. The philosophy and rhetoric which has grown around this concept (Leadbeater 2000, Castells 2000, Florida 2000, Caves 2000, Hartley 2000) has been influential in re-contextualising culture and the arts in the 21st century. Where governments and educational institutions have embraced the context of the creative industries, it is having a profound effect on the way arts are being positioned, originally as ‘creative content’ for the new economy. Countries and regions which have actively targeted the Creative Industries as an important economic growth factor in a post-industrial environment are numerous, but it is interesting to note that North and South East Asia and Australia have been at the forefront of developing the Creative Industries in its various guises. It could be argued that the initial phase of Creative Industries concentrated on media and communication technologies to provide commercial outcomes in small incubator business models; developing, for example, products for the games industry. Creative Industries is now entering a second phase of development; one in which the broader palette of the arts, though still not at the forefront of debate, is being re-examined. Both phases of Creative Industries have emphasised creativity and innovation as key drivers in the success and effectiveness of this sector, and although the arts by no means has a monopoly on these drivers, it is where they have an important part to play in the creative industries context. Arguably, the second wave of the creative industries acknowledges to a greater extent that commercialisation works in tandem with government and other support in a complex mixed economic model. In relation to the performing arts, the global market has seen an increase in large-scale cultural events such as festivals which are providing employment for the arts industry and multiplier effects in other parts of the economy. Differentiated product is important in this competitive arena and the use of mediated and digitised environments has been able to increase the amount of arts product available to an international market. This changed environment requires the development of new skills for our artists and producers and has given rise to a reappraisal of approaches to arts training and research in the Higher Degree Education sector (Brown 2007, Cunningham 2006). This paper examines pedagogical changes which took place in the first Creative Industries Faculty in the world at Queensland University of Technology as well as the increased opportunities for leading research initiatives. It concludes with the example of an interdisciplinary artwork produced in a creative industries precinct, exemplifying the convergence of arts and communication technologies and that of artistic practice and research.

Right-angle crash vulnerability of motorcycles at signalized intersections : mixed logit analysis

Motorcycles are particularly vulnerable in right-angle crashes at signalized intersections. The objective of this study is to explore how variations in roadway characteristics, environmental factors, traffic factors, maneuver types, human factors as well as driver demographics influence the right-angle crash vulnerability of motorcycles at intersections. The problem is modeled using a mixed logit model with a binary choice category formulation to differentiate how an at-fault vehicle collides with a not-at-fault motorcycle in comparison to other collision types. The mixed logit formulation allows randomness in the parameters and hence takes into account the underlying heterogeneities potentially inherent in driver behavior, and other unobserved variables. A likelihood ratio test reveals that the mixed logit model is indeed better than the standard logit model. Night time riding shows a positive association with the vulnerability of motorcyclists. Moreover, motorcyclists are particularly vulnerable on single lane roads, on the curb and median lanes of multi-lane roads, and on one-way and two-way road type relative to divided-highway. Drivers who deliberately run red light as well as those who are careless towards motorcyclists especially when making turns at intersections increase the vulnerability of motorcyclists. Drivers appear more restrained when there is a passenger onboard and this has decreased the crash potential with motorcyclists. The presence of red light cameras also significantly decreases right-angle crash vulnerabilities of motorcyclists. The findings of this study would be helpful in developing more targeted countermeasures for traffic enforcement, driver/rider training and/or education, safety awareness programs to reduce the vulnerability of motorcyclists.

Coping resources, coping strategies and adaptation to organisational change : direct or buffering effects?

This study examined the effects of personal and social resources, coping strategies and appraised stress on employees' levels of anxiety and depression. In relation to the effects of resources and coping strategies, two different models were tested. The main effects model proposes that, irrespective of the level of stress, coping resources and coping strategies have direct effects on well-being. In contrast, the buffering model predicts that the buffering effects of coping resources and strategies are only evident at high levels of stress. One hundred lawyers completed a structured self-administered questionnaire that measured their personal and social resources, use of problem-focused and emotion-focused coping strategies, and appraisals of the stressfulness of the situation. Results revealed generally strong support for the main effects model in the prediction of employee levels of anxiety and depression. Lower levels of anxiety were linked to judgements of lower levels of organizational change, greater self-confidence, greater internality of control beliefs and less use of emotion-focused coping strategies. Lower levels of depression in employees were also linked to judgements of lower levels of organizational change, greater use of resources and less appraised stress. There was only limited support for the buffering effects model. Due to the small size of the sample, the findings need to be explored further in other contexts.