Optimizing Presence–Absence Surveys For Detecting Population Trends

Presence-absence surveys are a commonly used method for monitoring broad-scale changes in wildlife distributions. However, the lack of power of these surveys for detecting population trends is problematic for their application in wildlife management. Options for improving power include increasing the sampling effort or arbitrarily relaxing the type I error rate. We present an alternative, whereby targeted sampling of particular habitats in the landscape using information from a habitat model increases power. The advantage of this approach is that it does not require a trade-off with either cost or the Pr(type I error) to achieve greater power. We use a demographic model of koala (Phascolarctos cinereus) population dynamics and simulations of the monitoring process to estimate the power to detect a trend in occupancy for a range of strategies, thereby demonstrating that targeting particular habitat qualities can improve power substantially. If the objective is to detect a decline in occupancy, the optimal strategy is to sample high-quality habitats. Alternatively, if the objective is to detect an increase in occupancy, the optimal strategy is to sample intermediate-quality habitats. The strategies with the highest power remained the same under a range of parameter assumptions, although observation error had a strong influence on the optimal strategy. Our approach specifically applies to monitoring for detecting long-term trends in occupancy or abundance. This is a common and important monitoring objective for wildlife managers, and we provide guidelines for more effectively achieving it.

Motor unit number estimation - A Bayesian approach

All muscle contractions are dependent on the functioning of motor units. In diseases such as amyotrophic lateral sclerosis (ALS), progressive loss of motor units leads to gradual paralysis. A major difficulty in the search for a treatment for these diseases has been the lack of a reliable measure of disease progression. One possible measure would be an estimate of the number of surviving motor units. Despite over 30 years of motor unit number estimation (MUNE), all proposed methods have been met with practical and theoretical objections. Our aim is to develop a method of MUNE that overcomes these objections. We record the compound muscle action potential (CMAP) from a selected muscle in response to a graded electrical stimulation applied to the nerve. As the stimulus increases, the threshold of each motor unit is exceeded, and the size of the CMAP increases until a maximum response is obtained. However, the threshold potential required to excite an axon is not a precise value but fluctuates over a small range leading to probabilistic activation of motor units in response to a given stimulus. When the threshold ranges of motor units overlap, there may be alternation where the number of motor units that fire in response to the stimulus is variable. This means that increments in the value of the CMAP correspond to the firing of different combinations of motor units. At a fixed stimulus, variability in the CMAP, measured as variance, can be used to conduct MUNE using the "statistical" or the "Poisson" method. However, this method relies on the assumptions that the numbers of motor units that are firing probabilistically have the Poisson distribution and that all single motor unit action potentials (MUAP) have a fixed and identical size. These assumptions are not necessarily correct. We propose to develop a Bayesian statistical methodology to analyze electrophysiological data to provide an estimate of motor unit numbers. Our method of MUNE incorporates the variability of the threshold, the variability between and within single MUAPs, and baseline variability. Our model not only gives the most probable number of motor units but also provides information about both the population of units and individual units. We use Markov chain Monte Carlo to obtain information about the characteristics of individual motor units and about the population of motor units and the Bayesian information criterion for MUNE. We test our method of MUNE on three subjects. Our method provides a reproducible estimate for a patient with stable but severe ALS. In a serial study, we demonstrate a decline in the number of motor unit numbers with a patient with rapidly advancing disease. Finally, with our last patient, we show that our method has the capacity to estimate a larger number of motor units.

Assessing the Conservation Status and Threats to Priority Plants: a threat assessment approach and case study in South-East Queensland, Australia

bstract: During the Regional Forest Agreement (RFA) process in south-east Queensland, the conservation status of, and threats to, priority vascular plant taxa in the region was assessed. Characteristics of biology, demography and distribution were used to assess the species' intrinsic risk of extinction. In contrast, the threats to the taxa (their extrinsic risk of extinction) were assessed using a decision-support protocol for setting conservation targets for taxa lacking population viability analyses and habitat modelling data. Disturbance processes known or suspected to be adversely affecting the taxa were evaluated for their intensity, extent and time-scale. Expert opinion was used to provide much of the data and to assess the recommended protection areas. Five categories of intrinsic risk of extinction were recognised for the 105 priority taxa: critically endangered (43 taxa); endangered (29); vulnerable (21); rare (10); and presumed extinct (2). Only 6 of the 103 extant taxa were found to be adequately reserved and the majority were considered inadequately protected to survive the current regimes of threatening processes affecting them. Data were insufficient to calculate a protection target for one extant taxon. Over half of the taxa require all populations to be conserved as well as active management to alleviate threatening processes. The most common threats to particular taxa were competition from weeds or native species, inappropriate fire regimes, agricultural clearing, forestry, grazing by native or feral species, drought, urban development, illegal collection of plants, and altered hydrology. Apart from drought and competition from native species, these disturbances are largely influenced or initiated by human actions. Therefore, as well as increased protection of most of the taxa, active management interventions are necessary to reduce the effects of threatening processes and to enable the persistence of the taxa.

Bayesian population pharmacokinetic analysis of sirolimus

Objective: It is usual that data collected from routine clinical care is sparse and unable to support the more complex pharmacokinetic (PK) models that may have been reported in previous rich data studies. Informative priors may be a pre-requisite for model development. The aim of this study was to estimate the population PK parameters of sirolimus using a fully Bayesian approach with informative priors. Methods: Informative priors including prior mean and precision of the prior mean were elicited from previous published studies using a meta-analytic technique. Precision of between-subject variability was determined by simulations from a Wishart distribution using MATLAB (version 6.5). Concentration-time data of sirolimus retrospectively collected from kidney transplant patients were analysed using WinBUGS (version 1.3). The candidate models were either one- or two-compartment with first order absorption and first order elimination. Model discrimination was based on computation of the posterior odds supporting the model. Results: A total of 315 concentration-time points were obtained from 25 patients. Most data were clustered at trough concentrations with range of 1.6 to 77 hours post-dose. Using informative priors, either a one- or two-compartment model could be used to describe the data. When a one-compartment model was applied, information was gained from the data for the value of apparent clearance (CL/F = 18.5 L/h), and apparent volume of distribution (V/F = 1406 L) but no information was gained about the absorption rate constant (ka). When a two-compartment model was fitted to the data, the data were informative about CL/F, apparent inter-compartmental clearance, and apparent volume of distribution of the peripheral compartment (13.2 L/h, 20.8 L/h, and 579 L, respectively). The posterior distribution of the volume distribution of central compartment and ka were the same as priors. The posterior odds for the two-compartment model was 8.1, indicating the data supported the two-compartment model. Conclusion: The use of informative priors supported the choice of a more complex and informative model that would otherwise have not been supported by the sparse data.

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.

Paediatric population pharmacokinetics of itraconazole and its active metabolite hydroxy-itraconazole in cystic fibrosis and bone marrow transplant patients

Objectives: The aim of the study was to characterise the population pharmacokinetics (popPK) properties of itraconazole (ITRA) and its active metabolite hydroxy-ITRA in a representative paediatric population of cystic fibrosis (CF) and bone marrow transplant (BMT) patients. The goals were to determine the relative bioavailability between the two oral formulations, and to explore improved dosage regimens in these patients. Methods: All paediatric patients with CF taking oral ITRA for the treatment of allergic bronchopulmonary aspergillosis and patients undergoing BMT who were taking ITRA for prophylaxis of any fungal infection were eligible for the study. A minimum of two blood samples were drawn after the capsules and also after switching to oral solution, or vice versa. ITRA and hydroxy-ITRA plasma concentrations were measured by HPLC[1]. A nonlinear mixed-effect modelling approach (NONMEM 5.1.1) was used to describe the PK of ITRA and hydroxy-ITRA simultaneously. Simulations were used to assess dosing strategies in these patients. Results: Forty-nine patients (29CF, 20 BMT) were recruited to the study who provided 227 blood samples for the population analysis. A 1-compartment model with 1st order absorption and elimination best described ITRA kinetics, with 1st order conversion to hydroxy-ITRA. For ITRA, the apparent clearance (ClItra/F) and volume of distribution (Vitra/F) was 35.5L/h and 672L, respectively; the absorption rate constant for the capsule formulation was 0.0901 h-1 and for the oral solution formulation it was 0.959 h-1. The capsule comparative bioavailability (vs. solution) was 0.55. For hydroxy-ITRA, the apparent volume of distribution and clearance were 10.6 L and 5.28 L/h, respectively. Of several screened covariates only allometrically scaled total body weight significantly improved the fit to the data. No difference between the two populations was found. Conclusion: The developed popPK model adequately described the pharmacokinetics of ITRA and hydroxy-ITRA in paediatric patients with CF and patients undergoing BMT. High inter-patient variability confirmed previous data in CF[2], leukaemia and BMT[3] patients. From the population model, simulations showed the standard dose (5 mg/kg/day) needs to be doubled for the solution formulation and even 4 times more given of the capsules to achieve an adequate target therapeutic trough plasma concentration of 0.5 mg/L[4] in these patients.

Use of routine data for determination of the population pharmacokinetics and enteral biovailability of phenytoin in neonates and infants with seizures

Objective: To investigate the population pharmacokinetics and the enteral bioavailability of phenytoin in neonates and infants with seizures. Methods: Data (5 mg kg-1 day-1) from 83 patients were obtained retrospectively from the medical records following written ethical approval. A one-compartment model was fitted to the data using NONMEM with FOCE-interaction. Between-subject variability (BSV) and interoccasion variability (IOV) were modelled exponentially together with a log transform-both-sides exponential residual unexplained variance (RUV) model. Covariates in nested models were screened for significance (X2, 1, 0.01). Model validity was determined by bootstrapping with replacement (N=500 samples) from the dataset. Results: The parameters of final pharmacokinetic were: Clearance (L h-1) = 0.826.(current Weight [kg]/70)0.75.(1+0.0692.(Postnatal age [days]-11)); Volume of distribution (L) = 74.2.(current Weight [kg]/70); Enteral bioavailability = 0.76; Absorption rate constant (h-1) = 0.167. BSV for clearance and volume of distribution were 74.2% and 65.6%, respectively. The IOV in clearance was 54.4%. The RUV was 51.1%. Final model parameters deviated from mean bootstrap estimates by

Learning to play Pac-Man: An evolutionary, rule-based approach

Pac-Man is a well-known, real-time computer game that provides an interesting platform for research. We describe an initial approach to developing an artificial agent that replaces the human to play a simplified version of Pac-Man. The agent is specified as a simple finite state machine and ruleset. with parameters that control the probability of movement by the agent given the constraints of the maze at some instant of time. In contrast to previous approaches, the agent represents a dynamic strategy for playing Pac-Man, rather than a pre-programmed maze-solving method. The agent adaptively "learns" through the application of population-based incremental learning (PBIL) to adjust the agents' parameters. Experimental results are presented that give insight into some of the complexities of the game, as well as highlighting the limitations and difficulties of the representation of the agent.

Towards a public health approach to parenting

There are various parenting, school and personal factors at play in determining a child’s risk of developing serious conduct problems. The temptation is therefore to conclude that “more is better than less”, but we think that has not been convincingly demonstrated. Some large-scale multi-risk-factor reduction approaches that include parenting, school and child-specific interventions with older school-aged children have shown promise but are complex to administer, costly to implement and have yet to show strong long-term outcomes. But in young children (toddler and preschool-aged children) there is strong evidence that social-learning-based parenting programmes are effective with a wide range of families from quite diverse socio-economic and ethnic backgrounds. We choose to focus on such programmes.

A heuristic approach for component scheduling on a high-speed PCB assembly machine

The collect-and-place machine is one of the most widely used placement machines for assembling electronic components on the printed circuit boards (PCBs). Nevertheless, the number of researches concerning the optimisation of the machine performance is very few. This motivates us to study the component scheduling problem for this type of machine with the objective of minimising the total assembly time. The component scheduling problem is an integration of the component sequencing problem, that is, the sequencing of component placements; and the feeder arrangement problem, that is, the assignment of component types to feeders. To solve the component scheduling problem efficiently, a hybrid genetic algorithm is developed in this paper. A numerical example is used to compare the performance of the algorithm with different component grouping approaches and different population sizes.

Typical kernel size and number of sparse random matrices over Galois fields: a statistical physics approach

Using methods of statistical physics, we study the average number and kernel size of general sparse random matrices over GF(q), with a given connectivity profile, in the thermodynamical limit of large matrices. We introduce a mapping of GF(q) matrices onto spin systems using the representation of the cyclic group of order q as the q-th complex roots of unity. This representation facilitates the derivation of the average kernel size of random matrices using the replica approach, under the replica symmetric ansatz, resulting in saddle point equations for general connectivity distributions. Numerical solutions are then obtained for particular cases by population dynamics. Similar techniques also allow us to obtain an expression for the exact and average number of random matrices for any general connectivity profile. We present numerical results for particular distributions.

Meta-analysis of the population and phenotypic expression of CYP2C9/2C19 polymorphism on drug metabolism in different ethnicities

The term "pharmacogenetics" has been defined as the scientific study of inherited factors that affect the human drug response. Many pharmacogenetie studies have been published since 1995 and have focussed on the principal enzyme family involved in drug metabolism, the cytochrome P450 family, particularly cytochrome P4502C9 and 2C19. In order to investigate the pharmacogenetic aspect of pharmacotherapy, the relevant studies describing the association of pharmacogenetic factor(s) in drug responses must be retrieved from existing literature using a systematic review approach. In addition, the estimation of variant allele prevalence for the gene under study between different ethnic populations is important for pharmacogenetic studies. In this thesis, the prevalence of CYP2C9/2C19 alleles between different ethnicities has been estimated through meta-analysis and the population genetic principle. The clinical outcome of CYP2C9/2C19 allelic variation on the pharmacotherapy of epilepsy has been investigated; although many new antiepileptic drugs have been launched into the market, carbamazepine, phenobarbital and phenytoin are still the major agents in the pharmacotherapy of epilepsy. Therefore, phenytoin was chosen as a model AED and the effect of CYP2C9/2C19 genetic polymorphism on phenytoin metabolism was further examined.An estimation of the allele prevalence was undertaken for three CYP2C9/2C19 alleles respectively using a meta-analysis of studies that fit the Hardy-Weinberg equilibrium. The prevalence of CYP2C9*1 is approximately 81%, 96%, 97% and 94% in Caucasian, Chinese, Japanese, African populations respectively; the pooled prevalence of CYP2C19*1 is about 86%, 57%, 58% and 85% in these ethnic populations respectively. However, the studies of association between CYP2C9/2C19 polymorphism and phenytoin metabolism failed to achieve any qualitative or quantitative conclusion. Therefore, mephenytoin metabolism was examined as a probe drug for association between CYP2C19 polymorphism and mephenytoin metabolic ratio. Similarly, analysis of association between CYP2C9 polymorphism and warfarin dose requirement was undertaken.It was confirmed that subjects carrying two mutated CYP2C19 alleles have higher S/R mephenytoin ratio due to deficient CYP2C19 enzyme activity. The studies of warfarin and CYP2C9 polymorphism did not provide a conclusive result due to poor comparability between studies.The genetic polymorphism of drug metabolism enzymes has been studied extensively, however other genetic factors, such as multiple drug resistance genes (MDR) and genes encoding ion channels, which may contribute to variability in function of drug transporters and targets, require more attention in future pharmacogenetic studies of antiepileptic drugs.