Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

A D-optimal designed population pharmacokinetic study of itraconazole capsules and solution in adults with cystic fibrosis

Background: Oral itraconazole (ITRA) is used for the treatment of allergic bronchopulmonary aspergillosis in patients with cystic fibrosis (CF) because of its antifungal activity against Aspergillus species. ITRA has an active hydroxy-metabolite (OH-ITRA) which has similar antifungal activity. ITRA is a highly lipophilic drug which is available in two different oral formulations, a capsule and an oral solution. It is reported that the oral solution has a 60% higher relative bioavailability. The influence of altered gastric physiology associated with CF on the pharmacokinetics (PK) of ITRA and its metabolite has not been previously evaluated. Objectives: 1) To estimate the population (pop) PK parameters for ITRA and its active metabolite OH-ITRA including relative bioavailability of the parent after administration of the parent by both capsule and solution and 2) to assess the performance of the optimal design. Methods: The study was a cross-over design in which 30 patients received the capsule on the first occasion and 3 days later the solution formulation. The design was constrained to have a maximum of 4 blood samples per occasion for estimation of the popPK of both ITRA and OH-ITRA. The sampling times for the population model were optimized previously using POPT v.2.0.[1] POPT is a series of applications that run under MATLAB and provide an evaluation of the information matrix for a nonlinear mixed effects model given a particular design. In addition it can be used to optimize the design based on evaluation of the determinant of the information matrix. The model details for the design were based on prior information obtained from the literature, which suggested that ITRA may have either linear or non-linear elimination. The optimal sampling times were evaluated to provide information for both competing models for the parent and metabolite and for both capsule and solution simultaneously. Blood samples were assayed by validated HPLC.[2] PopPK modelling was performed using FOCE with interaction under NONMEM, version 5 (level 1.1; GloboMax LLC, Hanover, MD, USA). The PK of ITRA and OH‑ITRA was modelled simultaneously using ADVAN 5. Subsequently three methods were assessed for modelling concentrations less than the LOD (limit of detection). These methods (corresponding to methods 5, 6 & 4 from Beal[3], respectively) were (a) where all values less than LOD were assigned to half of LOD, (b) where the closest missing value that is less than LOD was assigned to half the LOD and all previous (if during absorption) or subsequent (if during elimination) missing samples were deleted, and (c) where the contribution of the expectation of each missing concentration to the likelihood is estimated. The LOD was 0.04 mg/L. The final model evaluation was performed via bootstrap with re-sampling and a visual predictive check. The optimal design and the sampling windows of the study were evaluated for execution errors and for agreement between the observed and predicted standard errors. Dosing regimens were simulated for the capsules and the oral solution to assess their ability to achieve ITRA target trough concentration (Cmin,ss of 0.5-2 mg/L) or a combined Cmin,ss for ITRA and OH-ITRA above 1.5mg/L. Results and Discussion: A total of 241 blood samples were collected and analysed, 94% of them were taken within the defined optimal sampling windows, of which 31% where taken within 5 min of the exact optimal times. Forty six per cent of the ITRA values and 28% of the OH-ITRA values were below LOD. The entire profile after administration of the capsule for five patients was below LOD and therefore the data from this occasion was omitted from estimation. A 2-compartment model with 1st order absorption and elimination best described ITRA PK, with 1st order metabolism of the parent to OH-ITRA. For ITRA the clearance (ClItra/F) was 31.5 L/h; apparent volumes of central and peripheral compartments were 56.7 L and 2090 L, respectively. Absorption rate constants for capsule (kacap) and solution (kasol) were 0.0315 h-1 and 0.125 h-1, respectively. Comparative bioavailability of the capsule was 0.82. There was no evidence of nonlinearity in the popPK of ITRA. No screened covariate significantly improved the fit to the data. The results of the parameter estimates from the final model were comparable between the different methods for accounting for missing data, (M4,5,6)[3] and provided similar parameter estimates. The prospective application of an optimal design was found to be successful. Due to the sampling windows, most of the samples could be collected within the daily hospital routine, but still at times that were near optimal for estimating the popPK parameters. The final model was one of the potential competing models considered in the original design. The asymptotic standard errors provided by NONMEM for the final model and empirical values from bootstrap were similar in magnitude to those predicted from the Fisher Information matrix associated with the D-optimal design. Simulations from the final model showed that the current dosing regimen of 200 mg twice daily (bd) would provide a target Cmin,ss (0.5-2 mg/L) for only 35% of patients when administered as the solution and 31% when administered as capsules. The optimal dosing schedule was 500mg bd for both formulations. The target success for this dosing regimen was 87% for the solution with an NNT=4 compared to capsules. This means, for every 4 patients treated with the solution one additional patient will achieve a target success compared to capsule but at an additional cost of AUD $220 per day. The therapeutic target however is still doubtful and potential risks of these dosing schedules need to be assessed on an individual basis. Conclusion: A model was developed which described the popPK of ITRA and its main active metabolite OH-ITRA in adult CF after administration of both capsule and solution. The relative bioavailability of ITRA from the capsule was 82% that of the solution, but considerably more variable. To incorporate missing data, using the simple Beal method 5 (using half LOD for all samples below LOD) provided comparable results to the more complex but theoretically better Beal method 4 (integration method). The optimal sparse design performed well for estimation of model parameters and provided a good fit to the data.

Efficacy of memantine for agitation in Alzheimer's dementia:a randomised double-blind placebo controlled trial

Background - Agitation in Alzheimer’s disease (AD) is common and associated with poor patient life-quality and carer distress. The best evidence-based pharmacological treatments are antipsychotics which have limited benefits with increased morbidity and mortality. There are no memantine trials in clinically significant agitation but post-hoc analyses in other populations found reduced agitation. We tested the primary hypothesis, memantine is superior to placebo for clinically significant agitation, in patients with moderate-to-severe AD. Methods and Findings - We recruited 153 participants with AD and clinically significant agitation from care-homes or hospitals for a double-blind randomised-controlled trial and 149 people started the trial of memantine versus placebo. The primary outcome was 6 weeks mixed model autoregressive analysis of Cohen-Mansfield Agitation Inventory (CMAI). Secondary outcomes were: 12 weeks CMAI; 6 and 12 weeks Neuropsychiatric symptoms (NPI), Clinical Global Impression Change (CGI-C), Standardised Mini Mental State Examination, Severe Impairment Battery. Using a mixed effects model we found no significant differences in the primary outcome, 6 weeks CMAI, between memantine and placebo (memantine lower -3.0; -8.3 to 2.2, p = 0.26); or 12 weeks CMAI; or CGI-C or adverse events at 6 or 12 weeks. NPI mean difference favoured memantine at weeks 6 (-6.9; -12.2 to -1.6; p = 0.012) and 12 (-9.6; -15.0 to -4.3 p = 0.0005). Memantine was significantly better than placebo for cognition. The main study limitation is that it still remains to be determined whether memantine has a role in milder agitation in AD. Conclusions - Memantine did not improve significant agitation in people with in moderate-to-severe AD. Future studies are urgently needed to test other pharmacological candidates in this group and memantine for neuropsychiatric symptoms.

Mortality and climatic/environmental correlates of Widdringtonia cedargensis trees in the Cederberg, South Africa

Background: Conifer populations appear disproportionately threatened by global change. Most examples are, however, drawn from the northern hemisphere and long-term rates of population decline are not well documented as historical data are often lacking. We use a large and long-term (1931-2013) repeat photography dataset together with environmental data and fire records to account for the decline of the critically endangered Widdringtonia cedarbergensis. Eighty-seven historical and repeat photo-pairs were analysed to establish 20th century changes in W. cedarbergensis demography. A generalized linear mixed-effects model was fitted to determine the relative importance of environmental factors and fire-return interval on mortality for the species. Results: From an initial total of 1313 live trees in historical photographs, 74% had died and only 44 (3.4%) had recruited in the repeat photographs, leaving 387 live individuals. Juveniles (mature adults) had decreased (increased) from 27% (73%) to 8% (92%) over the intervening period. Our model demonstrates that mortality is related to greater fire frequency, higher temperatures, lower elevations, less rocky habitats and aspect (i.e. east-facing slopes had the least mortality). Conclusions: Our results show that W. cedarbergensis populations have declined significantly over the recorded period, with a pronounced decline in the last 30 years. Individuals that established in open habitats at lower, hotter elevations and experienced a greater fire frequency appear to be more vulnerable to mortality than individuals growing within protected, rocky environments at higher, cooler locations with less frequent fires. Climate models predict increasing temperatures for our study area (and likely increases in wildfires). If these predictions are realised, further declines in the species can be expected. Urgent management interventions, including seedling out-planting in fire-protected high elevation sites, reducing fire frequency in higher elevation populations, and assisted migration, should be considered.

Adiponectin and insulin in grey seals during suckling and fasting: relationship with nutritional state and body mass during nursing in mothers and pups

Animals that fast during breeding and/or development, such as phocids, must regulate energy balance carefully to maximize reproductive fitness and survival probability. Adiponectin, produced by adipose tissue, contributes to metabolic regulation by modulating sensitivity to insulin, increasing fatty acid oxidation by liver and muscle, and promoting adipogenesis and lipid storage in fat tissue. We tested the hypotheses that (1) circulating adiponectin, insulin, or relative adiponectin gene expression is related to nutritional state, body mass, and mass gain in wild gray seal pups; (2) plasma adiponectin or insulin is related to maternal lactation duration, body mass, percentage milk fat, or free fatty acid (FFA) concentration; and (3) plasma adiponectin and insulin are correlated with circulating FFA in females and pups. In pups, plasma adiponectin decreased during suckling (linear mixed-effects model [LME]: T = 4.49; P < 0.001) and the early postweaning fast (LME: T = 3.39; P = 0.004). In contrast, their blubber adiponectin gene expression was higher during the early postweaning fast than early in suckling (LME: T = 2.11; P = 0.046). Insulin levels were significantly higher in early (LME: T = 3.52; P = 0.004) and late (LME: T = 6.99; P < 0.001) suckling than in fasting and, given the effect of nutritional state, were also positively related to body mass (LME: T = 3.58; P = 0.004). Adiponectin and insulin levels did not change during lactation and were unrelated to milk FFA or percentage milk fat in adult females. Our data suggest that adiponectin, in conjunction with insulin, may facilitate fat storage in seals and is likely to be particularly important in the development of blubber reserves in pups.

Wind Farm Facilities in Germany Kill Noctule Bats from Near and Far

Over recent years, it became widely accepted that alternative, renewable energy may come at some risk for wildlife, for example, when wind turbines cause large numbers of bat fatalities. To better assess likely populations effects of wind turbine related wildlife fatalities, we studied the geographical origin of the most common bat species found dead below German wind turbines, the noctule bat (Nyctalus noctula). We measured stable isotope ratios of non-exchangeable hydrogen in fur keratin to separate migrants from local individuals, used a linear mixed-effects model to identify temporal, spatial and biological factors explaining the variance in measured stable isotope ratios and determined the geographical breeding provenance of killed migrants using isoscape origin models. We found that 72% of noctule bat casualties (n = 136) were of local origin, while 28% were long-distance migrants. These findings highlight that bat fatalities at German wind turbines may affect both local and distant populations. Our results indicated a sex and age-specific vulnerability of bats towards lethal accidents at turbines, i.e. a relatively high proportion of killed females were recorded among migratory individuals, whereas more juveniles than adults were recorded among killed bats of local origin. Migratory noctule bats were found to originate from distant populations in the Northeastern parts of Europe. The large catchment areas of German wind turbines and high vulnerability of female and juvenile noctule bats call for immediate action to reduce the negative cross-boundary effects of bat fatalities at wind turbines on local and distant populations. Further, our study highlights the importance of implementing effective mitigation measures and developing species and scale-specific conservation approaches on both national and international levels to protect source populations of bats. The efficacy of local compensatory measures appears doubtful, at least for migrant noctule bats, considering the large geographical catchment areas of German wind turbines for this species.

A Triple-Isotope Approach to Predict the Breeding Origins of European Bats

Despite a commitment by the European Union to protect its migratory bat populations, conservation efforts are hindered by a poor understanding of bat migratory strategies and connectivity between breeding and wintering grounds. Traditional methods like mark-recapture are ineffective to study broad-scale bat migratory patterns. Stable hydrogen isotopes (delta D) have been proven useful in establishing spatial migratory connectivity of animal populations. Before applying this tool, the method was calibrated using bat samples of known origin. Here we established the potential of delta D as a robust geographical tracer of breeding origins of European bats by measuring delta D in hair of five sedentary bat species from 45 locations throughout Europe. The delta D of bat hair strongly correlated with well-established spatial isotopic patterns in mean annual precipitation in Europe, and therefore was highly correlated with latitude. We calculated a linear mixed-effects model, with species as random effect, linking delta D of bat hair to precipitation delta D of the areas of hair growth. This model can be used to predict breeding origins of European migrating bats. We used delta C-13 and delta N-15 to discriminate among potential origins of bats, and found that these isotopes can be used as variables to further refine origin predictions. A triple-isotope approach could thereby pinpoint populations or subpopulations that have distinct origins. Our results further corroborated stable isotope analysis as a powerful method to delineate animal migrations in Europe.

Statistical modeling of protein lysate array data

The protein lysate array is an emerging technology for quantifying the protein concentration ratios in multiple biological samples. It is gaining popularity, and has the potential to answer questions about post-translational modifications and protein pathway relationships. Statistical inference for a parametric quantification procedure has been inadequately addressed in the literature, mainly due to two challenges: the increasing dimension of the parameter space and the need to account for dependence in the data. Each chapter of this thesis addresses one of these issues. In Chapter 1, an introduction to the protein lysate array quantification is presented, followed by the motivations and goals for this thesis work. In Chapter 2, we develop a multi-step procedure for the Sigmoidal models, ensuring consistent estimation of the concentration level with full asymptotic efficiency. The results obtained in this chapter justify inferential procedures based on large-sample approximations. Simulation studies and real data analysis are used to illustrate the performance of the proposed method in finite-samples. The multi-step procedure is simpler in both theory and computation than the single-step least squares method that has been used in current practice. In Chapter 3, we introduce a new model to account for the dependence structure of the errors by a nonlinear mixed effects model. We consider a method to approximate the maximum likelihood estimator of all the parameters. Using the simulation studies on various error structures, we show that for data with non-i.i.d. errors the proposed method leads to more accurate estimates and better confidence intervals than the existing single-step least squares method.

Willingness to Pay for Rural Landscape Improvements: Combining Mixed Logit and Random-Effects Models

This paper reports the findings from a discrete-choice experiment designed to estimate the economic benefits associated with rural landscape improvements in Ireland. Using a mixed logit model, the panel nature of the dataset is exploited to retrieve willingness-to-pay values for every individual in the sample. This departs from customary approaches in which the willingness-to-pay estimates are normally expressed as measures of central tendency of an a priori distribution. Random-effects models for panel data are subsequently used to identify the determinants of the individual-specific willingness-to-pay estimates. In comparison with the standard methods used to incorporate individual-specific variables into the analysis of discrete-choice experiments, the analytical approach outlined in this paper is shown to add considerable explanatory power to the welfare estimates.

Examining the Differential Effects of a Universal SEL Curriculum on Student Functioning on a Dual Continua Model of Mental Health

Thesis (Ph.D.)--University of Washington, 2016-06

Team diversity and leadership: Tests of model of process and direct effects

A model was tested to examine relationships among leadership behaviors, team diversity, and team process measures with team performance and satisfaction at both the team and leader-member levels of analysis. Relationships between leadership behavior and team demographic and cognitive diversity were hypothesized to have both direct effects on organizational outcomes as well as indirect effects through team processes. Leader member differences were investigated to determine the effects of leader-member diversity leader-member exchange quality, individual effectiveness and satisfaction.^ Leadership had little direct effect on team performance, but several strong positive indirect effects through team processes. Demographic Diversity had no impact on team processes, directly impacted only one performance measure, and moderated the leadership to team process relationship.^ Cognitive Diversity had a number of direct and indirect effects on team performance, the net effects uniformly positive, and did not moderate the leadership to team process relationship.^ In sum, the team model suggests a complex combination of leadership behaviors positively impacting team processes, demographic diversity having little impact on team process or performance, cognitive diversity having a positive net impact impact, and team processes having mixed effects on team outcomes.^ At the leader-member level, leadership behaviors were a strong predictor of Leader-Member Exchange (LMX) quality. Leader-member demographic and cognitive dissimilarity were each predictors of LMX quality, but failed to moderate the leader behavior to LMX quality relationship. LMX quality was strongly and positively related to self reported effectiveness and satisfaction.^ The study makes several contributions to the literature. First, it explicitly links leadership and team diversity. Second, demographic and cognitive diversity are conceptualized as distinct and multi-faceted constructs. Third, a methodology for creating an index of categorical demographic and interval cognitive measures is provided so that diversity can be measured in a holistic conjoint fashion. Fourth, the study simultaneously investigates the impact of diversity at the team and leader-member levels of analyses. Fifth, insights into the moderating impact of different forms of team diversity on the leadership to team process relationship are provided. Sixth, this study incorporates a wide range of objective and independent measures to provide a 360$\sp\circ$ assessment of team performance. ^

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.

Genetic effects and genetic relationships among shrunken (sh2) sweet corn lines and F1 hybrids

The objectives of this study were to quantify the components of genetic variance and the genetic effects, and to examine the genetic relationship of inbred lines extracted from various shrunken2(sh2) breeding populations. Ten diverse inbred lines developed from genetic background, were crossed in half diallel. Parents and their F1 hybrids were evaluated at three environments. The parents were genotyped using 20 polymorphic simple sequence repeats (SSR). Agronomic and quality traits were analysed by a mixed linear model according to additive-dominance genetic model. Genetic effects were estimated using an adjusted unbiased prediction method. Additive variance was more important than dominance variance in the expression of traits related to ear aspects (husk ratio and percentage of ear filled) and eating quality (flavour and total soluble solids). For agronomic traits, however, dominance variance was more important than additive variance. The additive genetic correlation between flavour and tenderness was strong (r = 0.84, P <0.01). Flavour, tenderness and kernel colour additive genetic effects were not correlated with yield related traits. Genetic distance (GD), estimated from SSR profiles on the basis of Jaccard's similarity coefficient varied from 0.10 to 0.77 with an average of 0.56. Cluster analysis classified parents according to their pedigree relationships. In most studied traits, F1 performance was not associated with GD.