A generalized Markov sampler

A recent development of the Markov chain Monte Carlo (MCMC) technique is the emergence of MCMC samplers that allow transitions between different models. Such samplers make possible a range of computational tasks involving models, including model selection, model evaluation, model averaging and hypothesis testing. An example of this type of sampler is the reversible jump MCMC sampler, which is a generalization of the Metropolis-Hastings algorithm. Here, we present a new MCMC sampler of this type. The new sampler is a generalization of the Gibbs sampler, but somewhat surprisingly, it also turns out to encompass as particular cases all of the well-known MCMC samplers, including those of Metropolis, Barker, and Hastings. Moreover, the new sampler generalizes the reversible jump MCMC. It therefore appears to be a very general framework for MCMC sampling. This paper describes the new sampler and illustrates its use in three applications in Computational Biology, specifically determination of consensus sequences, phylogenetic inference and delineation of isochores via multiple change-point analysis.

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.

Sharpening the edge of market driven intangibles and innovations

A review of the extant literature concludes that market-driven intangibles and innovations are increasingly considered to be the most critical firm-specific resources, but also finds a lack of elaboration of which types of these resources are most important. In this paper, we incorporate these observations into a conceptual model and link it to highly developed institutional settings for the model evaluation. From the point of view of firm revenue management, we can anticipate that performance advantages created through deployment of intellectual and relational capital in marketing and innovation are more likely to be superior. In essence, they constitute the integration of organisational intangibles both in cognitive and behavioural level to create an idiosyncratic combination for each firm. Our research findings show feasible paths for sharpening the edge of market-driven intangibles and innovations. We discuss the key results for research and practice.

COOPER-framework:a unified process for non-parametric projects

Practitioners assess performance of entities in increasingly large and complicated datasets. If non-parametric models, such as Data Envelopment Analysis, were ever considered as simple push-button technologies, this is impossible when many variables are available or when data have to be compiled from several sources. This paper introduces by the 'COOPER-framework' a comprehensive model for carrying out non-parametric projects. The framework consists of six interrelated phases: Concepts and objectives, On structuring data, Operational models, Performance comparison model, Evaluation, and Result and deployment. Each of the phases describes some necessary steps a researcher should examine for a well defined and repeatable analysis. The COOPER-framework provides for the novice analyst guidance, structure and advice for a sound non-parametric analysis. The more experienced analyst benefits from a check list such that important issues are not forgotten. In addition, by the use of a standardized framework non-parametric assessments will be more reliable, more repeatable, more manageable, faster and less costly. © 2010 Elsevier B.V. All rights reserved.

Travelling waves in a complex Ginzburg-Landau equation with time-delay feedback

One of the simplest ways to create nonlinear oscillations is the Hopf bifurcation. The spatiotemporal dynamics observed in an extended medium with diffusion (e.g., a chemical reaction) undergoing this bifurcation is governed by the complex Ginzburg-Landau equation, one of the best-studied generic models for pattern formation, where besides uniform oscillations, spiral waves, coherent structures and turbulence are found. The presence of time delay terms in this equation changes the pattern formation scenario, and different kind of travelling waves have been reported. In particular, we study the complex Ginzburg-Landau equation that contains local and global time-delay feedback terms. We focus our attention on plane wave solutions in this model. The first novel result is the derivation of the plane wave solution in the presence of time-delay feedback with global and local contributions. The second and more important result of this study consists of a linear stability analysis of plane waves in that model. Evaluation of the eigenvalue equation does not show stabilisation of plane waves for the parameters studied. We discuss these results and compare to results of other models.

Klímaváltozás és az európai lepkeszúnyogok. Az areamodellezés módszertani kérdései

The future northward expansion of the arthropod vectors of leishmaniasis caused by climate change seems to be essential veterinary and medical problem. Our aim was to build and evaluate a Climate Envelope Model (CEM) to assess the potential effects of climate change on five European sandfly species. The studied species – Phlebotomus ariasi Tonn., P. neglectus Tonn., P. papatasi Scop., P. perfiliewi Parrot, P. perniciosus Newst., P. sergenti Parrot, P. similis Perfiliev, P. tobbi Adler, Theodor et Lourie – are important vectors of the parasite Leishmania infantum or other Leishmania species. The projections were based on REMO regional climate model with European domain. The climate data were available in a 25 km resolution grid for the reference period (1961-90) and two future periods (2011-40, 2041-70). The regional climate model was based on the IPCC SRES A1B scenario. Three types of climatic parameters were used for every month (averaged in the 30-years periods). The model was supported by VBORNET digital area database (distribution maps), ESRI ArcGIS 10 software’s Spatial Analyst module (modeling environment), PAST (calibration of the model with statistical method). Iterative model evaluation was done by summarizing two types of model errors based on an aggregated distribution. The results show that the best model results can be achieved by leaving 5-5 percentiles from the two extrema of the mean temperature, 2-2 percentiles from the two extrema of the minimum temperature, 0 percentile from the minimum of and 8 percentiles from the maximum of the precipitation.

Risk Prediction for Epithelial Ovarian Cancer in 11 United States-Based Case-Control Studies: Incorporation of Epidemiologic Risk Factors and 17 Confirmed Genetic Loci.

Previously developed models for predicting absolute risk of invasive epithelial ovarian cancer have included a limited number of risk factors and have had low discriminatory power (area under the receiver operating characteristic curve (AUC) < 0.60). Because of this, we developed and internally validated a relative risk prediction model that incorporates 17 established epidemiologic risk factors and 17 genome-wide significant single nucleotide polymorphisms (SNPs) using data from 11 case-control studies in the United States (5,793 cases; 9,512 controls) from the Ovarian Cancer Association Consortium (data accrued from 1992 to 2010). We developed a hierarchical logistic regression model for predicting case-control status that included imputation of missing data. We randomly divided the data into an 80% training sample and used the remaining 20% for model evaluation. The AUC for the full model was 0.664. A reduced model without SNPs performed similarly (AUC = 0.649). Both models performed better than a baseline model that included age and study site only (AUC = 0.563). The best predictive power was obtained in the full model among women younger than 50 years of age (AUC = 0.714); however, the addition of SNPs increased the AUC the most for women older than 50 years of age (AUC = 0.638 vs. 0.616). Adapting this improved model to estimate absolute risk and evaluating it in prospective data sets is warranted.

Application Of A Semi-empirical Model For The Evaluation Of Transmission Properties Of Barite Mortar.

The aim of this study was to estimate barite mortar attenuation curves using X-ray spectra weighted by a workload distribution. A semi-empirical model was used for the evaluation of transmission properties of this material. Since ambient dose equivalent, H(⁎)(10), is the radiation quantity adopted by IAEA for dose assessment, the variation of the H(⁎)(10) as a function of barite mortar thickness was calculated using primary experimental spectra. A CdTe detector was used for the measurement of these spectra. The resulting spectra were adopted for estimating the optimized thickness of protective barrier needed for shielding an area in an X-ray imaging facility.

Evaluation Of Characteristic-to-total Spectrum Ratio: Comparison Between Experimental And A Semi-empirical Model.

Primary X-ray spectra were measured in the range of 80-150kV in order to validate a computer program based on a semiempirical model. The ratio between the characteristic and total air Kerma was considered to compare computed results and experimental data. Results show that the experimental spectra have higher first HVL and mean energy than the calculated ones. The ratios between the characteristic and total air Kerma for calculated spectra are in good agreement with experimental results for all filtrations used.

Evaluation of an innovative critical power model in intermittent vertical jump

The aim of this study was to test if the critical power model can be used to determine the critical rest interval (CRI) between vertical jumps. Ten males performed intermittent countermovement jumps on a force platform with different resting periods (4.1 +/- 0.3 s, 5.0 +/- 0.4 s, 5.9 +/- 0.6 s). Jump trials were interrupted when participants could no longer maintain 95% of their maximal jump height. After interruption, number of jumps, total exercise duration and total external work were computed. Time to exhaustion (s) and total external work (J) were used to solve the equation Work = a + b . time. The CRI (corresponding to the shortest resting interval that allowed jump height to be maintained for a long time without fatigue) was determined dividing the average external work needed to jump at a fixed height (J) by b parameter (J/s). in the final session, participants jumped at their calculated CRI. A high coefficient of determination (0.995 +/- 0.007) and the CRI (7.5 +/- 1.6 s) were obtained. In addition, the longer the resting period, the greater the number of jumps (44 13, 71 28, 105 30, 169 53 jumps; p<0.0001), time to exhaustion (179 +/- 50, 351 +/- 120, 610 +/- 141, 1,282 +/- 417 s; p<0.0001) and total external work (28.0 +/- 8.3, 45.0 +/- 16.6, 67.6 +/- 17.8, 111.9 +/- 34.6 kJ; p<0.0001). Therefore, the critical power model may be an alternative approach to determine the CRI during intermittent vertical jumps.

Evaluation of the role of environmental factors in the human gastrointestinal tract on the behaviour of probiotic cultures of Lactobacillus casei Shirota and Lactobacillus casei LC01 by the use of a semi-dynamic in vitro model

This study evaluated the influence of gastrointestinal environmental factors (pH, digestive enzymes, food components, medicaments) on the survival of Lactobacillus casei Shirota and Lactobacillus casei LC01, using a semi-dynamic in vitro model that simulates the transit of microorganisms through the human GIT. The strains were first exposed to different simulated gastric juices for different periods of time (0, 30, 60 and 120 min), and then to simulated intestinal fluids for zero, 120, 180 and 240 min, in a step-wise format. The number of viable cells was determined after each step. The influence of food residues (skim milk) in the fluids and resistance to medicaments commonly used for varied therapeutic purposes (analgesics, antiarrhythmics, antibiotics, antihistaminics, proton pump inhibitors, etc.) were also evaluated. Results indicated that survival of both cultures was pH and time dependent, and digestive enzymes had little influence. Milk components presented a protective effect, and medicaments, especially anti-inflammatory drugs, influenced markedly the viability of the probiotic cultures, indicating that the beneficial effects of the two probiotic cultures to health are dependent of environmental factors encountered in the human gastrointestinal tract.

Excised Porcine Cornea Integrity Evaluation in an in vitro Model of Iontophoretic Ocular Research

Background/Aims: It is a challenge to adapt traditional in vitro diffusion experiments to ocular tissue. Thus, the aim of this work was to present experimental evidence on the integrity of the porcine cornea, barrier function and maintenance of electrical properties for 6 h of experiment when the tissue is mounted on an inexpensive and easy-to-use in vitro model for ocular iontophoresis. Methods: A modified Franz diffusion cell containing two ports for the insertion of the electrodes and a receiving compartment that does not need gassing with carbogen was used in the studies. Corneal electron transmission microscopy images were obtained, and diffusion experiments with fluorescent markers were performed to examine the integrity of the barrier function. The preservation of the negatively charged corneal epithelium was verified by the determination of the electro-osmotic flow of a hydrophilic and non-ionized molecule. Results: The diffusion cell was able to maintain the temperature, homogenization, porcine epithelial corneal structure integrity, barrier function and electrical characteristics throughout the 6 h of permeation experiment, without requiring CO(2) gassing when the receiving chamber was filled with 25 m M of HEPES buffer solution. Conclusion: The system described here is inexpensive, easy to handle and reliable as an in vitro model for iontophoretic ocular delivery studies. Copyright (C) 2010 S. Karger AG, Basel

Environmental impact evaluation using a cooperative model for implementing EMS (ISO 14001) in small and medium-sized enterprises

Certification of an ISO 14001 Environmental Management System (EMS) is currently an important requirement for those enterprises wishing to sell their products in the context of a global market. The system`s structure is based on environmental impact evaluation (EIE). However, if an erroneous or inadequate methodology is applied, the entire process may be jeopardized. Many methodologies have been developed for making of EIEs, some of them are fairly complex and unsuitable for EMS implementation in an organizational context, principally when small and medium size enterprises (SMEs) are involved. The proposed methodology for EIE is part of a model for implementing EMS. The methodological approach used was a qualitative exploratory research method based upon sources of evidence such as document analyses, semi-structured interviews and participant observations. By adopting a cooperative implementation model based on the theory of system engineering, difficulties relating to implementation of the sub-system were overcome thus encouraging SMEs to implement EMS. (C) 2007 Elsevier Ltd. All rights reserved.

Evaluation of chronic omega-3 fatty acids supplementation on behavioral and neurochemical alterations in 6-hydroxydopamine-lesion model of Parkinson`s disease

Omega-3 polyunsaturated fatty acids (omega-3 PUFAs) have been widely associated to beneficial effects over different neuropathologies, but only a few studies associate them to Parkinson`s disease (PD). Rats were submitted to chronic supplementation (21-90 days of life) with fish oil, rich in omega-3 PUFAs, and were uni- or bilaterally lesioned with 4 mu g of the neurotoxin 6-hydroxydopamine (6-OHDA) in the medial forebrain bundle Although lipid incorporation was evidenced in neuronal membranes, it was not sufficient to compensate motor deficits induced by 6-OHDA. In contrast, omega-3 PUFAs were capable of reducing rotational behavior induced by apomorphine, suggesting neuroprotection over dyskinesia The beneficial effects of omega-3 PUFAs were also evident in the maintenance of thiobarbituric acid reactive substances index from animals lesioned with 6-OHDA similar to levels from SHAM and intact animals. Although omega-3 PUFAs did not modify the tyrosine hydroxylase immunoreactivity in the substantia nigra pars compacta and in the ventral tegmental area, nor the depletion of dopamine (DA) and its metabolites in the striatum, DA turnover was increased after omega-3 PUFAs chronic supplementation Therefore, it is proposed that omega-3 PUFAs action characterizes the adaptation of remaining neurons activity. altering striatal DA turnover without modifying the estimated neuronal population. (C) 2009 Elsevier Ireland Ltd and the Japan Neuroscience Society. All rights reserved

On the evaluation of some differential formulations for the pom-pom constitutive model

In recent years, progress has been made in modelling long chain branched polymers by the introduction of the so-called pompom model. Initially developed by McLeish and Larson (1998), the model has undergone several improvements or alterations, leading to the development of new formulations. Some of these formulations however suffer from certain mathematical defects. The purpose of the present paper is to review some of the formulations of the pom-pom constitutive model, and to investigate their possible mathematical defects. Next, an alternative formulation is proposed, which does not appear to exhibit mathematical defects, and we explore its modelling performance by comparing the predictions with experiments in non-trivial rheometric flows of an LDPE melt. The selected rheometric flows are the double step strain, as well as the large amplitude oscillatory shear experiments. For LAOS experiments, the comparison involves the use of Fourier-transform analysis.