RANDOM EFFECTS MODELS IN A META-ANALYSIS OF THE ACCURACY OF DIAGNOSTIC TESTS WITHIN A GOLD STANDARD IN THE PRESENCE OF MISSING DATA

In evaluating the accuracy of diagnosis tests, it is common to apply two imperfect tests jointly or sequentially to a study population. In a recent meta-analysis of the accuracy of microsatellite instability testing (MSI) and traditional mutation analysis (MUT) in predicting germline mutations of the mismatch repair (MMR) genes, a Bayesian approach (Chen, Watson, and Parmigiani 2005) was proposed to handle missing data resulting from partial testing and the lack of a gold standard. In this paper, we demonstrate an improved estimation of the sensitivities and specificities of MSI and MUT by using a nonlinear mixed model and a Bayesian hierarchical model, both of which account for the heterogeneity across studies through study-specific random effects. The methods can be used to estimate the accuracy of two imperfect diagnostic tests in other meta-analyses when the prevalence of disease, the sensitivities and/or the specificities of diagnostic tests are heterogeneous among studies. Furthermore, simulation studies have demonstrated the importance of carefully selecting appropriate random effects on the estimation of diagnostic accuracy measurements in this scenario.

A Mechanistic Latent Variable Model for Estimating Drug Concentrations in the Male Genital Tract

The purpose of this study is to develop statistical methodology to facilitate indirect estimation of the concentration of antiretroviral drugs and viral loads in the prostate gland and the seminal vesicle. The differences in antiretroviral drug concentrations in these organs may lead to suboptimal concentrations in one gland compared to the other. Suboptimal levels of the antiretroviral drugs will not be able to fully suppress the virus in that gland, lead to a source of sexually transmissible virus and increase the chance of selecting for drug resistant virus. This information may be useful selecting antiretroviral drug regimen that will achieve optimal concentrations in most of male genital tract glands. Using fractionally collected semen ejaculates, Lundquist (1949) measured levels of surrogate markers in each fraction that are uniquely produced by specific male accessory glands. To determine the original glandular concentrations of the surrogate markers, Lundquist solved a simultaneous series of linear equations. This method has several limitations. In particular, it does not yield a unique solution, it does not address measurement error, and it disregards inter-subject variability in the parameters. To cope with these limitations, we developed a mechanistic latent variable model based on the physiology of the male genital tract and surrogate markers. We employ a Bayesian approach and perform a sensitivity analysis with regard to the distributional assumptions on the random effects and priors. The model and Bayesian approach is validated on experimental data where the concentration of a drug should be (biologically) differentially distributed between the two glands. In this example, the Bayesian model-based conclusions are found to be robust to model specification and this hierarchical approach leads to more scientifically valid conclusions than the original methodology. In particular, unlike existing methods, the proposed model based approach was not affected by a common form of outliers.

HHV-8 seroprevalence: a global view

BACKGROUND Human herpes virus 8 (HHV-8) is the underlying infectious cause of Kaposi sarcoma (KS) and other proliferative diseases; that is, primary effusion lymphoma and multicentric Castleman disease. In regions with high HHV-8 seroprevalence in the general population, KS accounts for a major burden of disease. Outside these endemic regions, HHV-8 prevalence is high in men who have sex with men (MSM) and in migrants from endemic regions. We aim to conduct a systematic literature review and meta-analysis in order 1) to define the global distribution of HHV-8 seroprevalence (primary objective) and 2) to identify risk factors for HHV-8 infection, with a focus on HIV status (secondary objective). METHODS/DESIGN We will include observational studies reporting data on seroprevalence of HHV-8 in children and/or adults from any region in the world. Case reports and case series as well as any studies with fewer than 50 participants will be excluded. We will search MEDLINE, EMBASE, and relevant conference proceedings without language restriction. Two reviewers will independently screen the identified studies and extract data on study characteristics and quality, study population, risk factors, and reported outcomes, using a standardized form. For the primary objective we will pool the data using a fully bayesian approach for meta-analysis, with random effects at the study level. For the secondary objective (association of HIV and HHV-8) we aim to pool odds ratios for the association of HIV and HHV-8 using a fully bayesian approach for meta-analysis, with random effects at the study level. Sub-group analyses and meta-regression analyses will be used to explore sources of heterogeneity, including factors such as geographical region, calendar years of recruitment, age, gender, ethnicity, socioeconomic status, different risk groups for sexually and parenterally transmitted infections (MSM, sex workers, hemophiliacs, intravenous drug users), comorbidities such as organ transplantation and malaria, test(s) used to measure HHV-8 infection, study design, and study quality. DISCUSSION Using the proposed systematic review and meta-analysis, we aim to better define the global seroprevalence of HHV-8 and its associated risk factors. This will improve the current understanding of HHV-8 epidemiology, and could suggest measures to prevent HHV-8 infection and to reduce its associated cancer burden.

Allowable carbon emissions lowered by multiple climate targets

Climate targets are designed to inform policies that would limit the magnitude and impacts of climate change caused by anthropogenic emissions of greenhouse gases and other substances. The target that is currently recognized by most world governments1 places a limit of two degrees Celsius on the global mean warming since preindustrial times. This would require large sustained reductions in carbon dioxide emissions during the twenty-first century and beyond2, 3, 4. Such a global temperature target, however, is not sufficient to control many other quantities, such as transient sea level rise5, ocean acidification6, 7 and net primary production on land8, 9. Here, using an Earth system model of intermediate complexity (EMIC) in an observation-informed Bayesian approach, we show that allowable carbon emissions are substantially reduced when multiple climate targets are set. We take into account uncertainties in physical and carbon cycle model parameters, radiative efficiencies10, climate sensitivity11 and carbon cycle feedbacks12, 13 along with a large set of observational constraints. Within this framework, we explore a broad range of economically feasible greenhouse gas scenarios from the integrated assessment community14, 15, 16, 17 to determine the likelihood of meeting a combination of specific global and regional targets under various assumptions. For any given likelihood of meeting a set of such targets, the allowable cumulative emissions are greatly reduced from those inferred from the temperature target alone. Therefore, temperature targets alone are unable to comprehensively limit the risks from anthropogenic emissions.

The dating of dipterocarp tree rings: establishing a record of carbon cycling and climatic change in the tropics

In a first step to obtain a proxy record of past climatic events (including the El Ni (n) over tildeo-Southern Oscillation) in the normally aseasonal tropical environment of Sabah, a radial segment from a recently fallen dipterocarp (Shorea Superba) was radiocarbon dated and subjected to carbon isotope analysis. The high-precision radiocarbon results fell into the ambiguous modern plateau where several calibrated dates can exist for each sample. Dating was achieved by wiggle matching using a Bayesian approach to calibration. Using the defined growth characteristics of Shorea superba, probability density distributions were calculated and improbable dates rejected. It was found that the tree most likely started growing around AD 1660-1685. A total of 173 apparent growth increments were measured and, therefore, it could be determined that the tree formed one ring approximately every two years. Stable carbon isotope values were obtained from resin-extracted wholewood from each ring. Carbon cycling is evident in the `juvenile effect', resulting from the assimilation of respired carbon dioxide and lower light levels below the canopy, and in the `anthropogenic effect' caused by increased industrial activity in the late-nineteenth and twentieth centuries. This study demonstrates that palaeoenvironmental information can be obtained from trees growing in aseasonal environments, where climatic conditions prevent the formation of well-defined annual rings.

Search for contact interactions and large extra dimensions in dilepton events from pp collisions at sqrts=7 TeV with the ATLAS detector

A search for nonresonant new phenomena, originating from either contact interactions or large extra spatial dimensions, has been carried out using events with two isolated electrons or muons. These events, produced at the LHC in proton-proton collisions at root s = 7 TeV, were recorded by the ATLAS detector. The data sample, collected throughout 2011, corresponds to an integrated luminosity of 4.9 and 5.0 fb(-1) in the e(+)e(-) and mu(+)mu(-) channels, respectively. No significant deviations from the Standard Model expectation are observed. Using a Bayesian approach, 95% confidence level lower limits ranging from 9.0 to 13.9 TeV are placed on the energy scale of llqq contact interactions in the left-left isoscalar model. Lower limits ranging from 2.4 to 3.9 TeV are also set on the string scale in large extra dimension models. After combining these limits with results from a similar search in the diphoton channel, slightly more stringent limits are obtained.

Using Bayes' rule to define the value of evidence from syndromic surveillance

In this work we propose the adoption of a statistical framework used in the evaluation of forensic evidence as a tool for evaluating and presenting circumstantial "evidence" of a disease outbreak from syndromic surveillance. The basic idea is to exploit the predicted distributions of reported cases to calculate the ratio of the likelihood of observing n cases given an ongoing outbreak over the likelihood of observing n cases given no outbreak. The likelihood ratio defines the Value of Evidence (V). Using Bayes' rule, the prior odds for an ongoing outbreak are multiplied by V to obtain the posterior odds. This approach was applied to time series on the number of horses showing clinical respiratory symptoms or neurological symptoms. The separation between prior beliefs about the probability of an outbreak and the strength of evidence from syndromic surveillance offers a transparent reasoning process suitable for supporting decision makers. The value of evidence can be translated into a verbal statement, as often done in forensics or used for the production of risk maps. Furthermore, a Bayesian approach offers seamless integration of data from syndromic surveillance with results from predictive modeling and with information from other sources such as disease introduction risk assessments.

Agglomeration effects of informal sector: evidence from Cambodia

The presence of a large informal sector in developing economies poses the question of whether informal activity produces agglomeration externalities. This paper uses data on all the nonfarm establishments and enterprises in Cambodia to estimate the impact of informal agglomeration on the regional economic performance of formal and informal firms. We develop a Bayesian approach for a spatial autoregressive model with an endogenous explanatory variable to address endogeneity and spatial dependence. We find a significantly positive effect of informal agglomeration, where informal firms gain more strongly than formal firms. Calculating the spatial marginal effects of increased agglomeration, we demonstrate that more accessible regions are more likely than less accessible regions to benefit strongly from informal agglomeration.

Nonparametric Message Passing Methods for Cooperative Localization and Tracking

The objective of this thesis is the development of cooperative localization and tracking algorithms using nonparametric message passing techniques. In contrast to the most well-known techniques, the goal is to estimate the posterior probability density function (PDF) of the position of each sensor. This problem can be solved using Bayesian approach, but it is intractable in general case. Nevertheless, the particle-based approximation (via nonparametric representation), and an appropriate factorization of the joint PDFs (using message passing methods), make Bayesian approach acceptable for inference in sensor networks. The well-known method for this problem, nonparametric belief propagation (NBP), can lead to inaccurate beliefs and possible non-convergence in loopy networks. Therefore, we propose four novel algorithms which alleviate these problems: nonparametric generalized belief propagation (NGBP) based on junction tree (NGBP-JT), NGBP based on pseudo-junction tree (NGBP-PJT), NBP based on spanning trees (NBP-ST), and uniformly-reweighted NBP (URW-NBP). We also extend NBP for cooperative localization in mobile networks. In contrast to the previous methods, we use an optional smoothing, provide a novel communication protocol, and increase the efficiency of the sampling techniques. Moreover, we propose novel algorithms for distributed tracking, in which the goal is to track the passive object which cannot locate itself. In particular, we develop distributed particle filtering (DPF) based on three asynchronous belief consensus (BC) algorithms: standard belief consensus (SBC), broadcast gossip (BG), and belief propagation (BP). Finally, the last part of this thesis includes the experimental analysis of some of the proposed algorithms, in which we found that the results based on real measurements are very similar with the results based on theoretical models.

Estrategias innovadoras para la gestión del riesgo de sequía en sistemas ganaderos extensivos

La vulnerabilidad de los sistemas ganaderos de pastoreo pone en evidencia la necesidad de herramientas para evaluar y mitigar los efectos de la sequía. El avance en la teledetección ha despertado el interés por explotar potenciales aplicaciones, y está dando lugar a un intenso desarrollo de innovaciones en distintos campos. Una de estas áreas es la gestión del riesgo climático, en donde la utilización de índices de vegetación permite la evaluación de la sequía. En esta investigación, se analiza el impacto de la sequía y se evalúa el potencial de nuevas tecnologías como la teledetección para la gestión del riesgo de sequía en sistemas de ganadería extensiva. Para ello, se desarrollan tres aplicaciones: (i) evaluar el impacto económico de la sequía en una explotación ganadera extensiva de la dehesa de Andalucía, (ii) elaborar mapas de vulnerabilidad a la sequía en pastos de Chile y (iii) diseñar y evaluar el potencial de un seguro indexado para sequía en pastos en la región de Coquimbo en Chile. En la primera aplicación, se diseña un modelo dinámico y estocástico que integra aspectos climáticos, ecológicos, agronómicos y socioeconómicos para evaluar el riesgo de sequía. El modelo simula una explotación ganadera tipo de la dehesa de Andalucía para el período 1999-2010. El método de Análisis Histórico y la simulación de MonteCarlo se utilizan para identificar los principales factores de riesgo de la explotación, entre los que destacan, los periodos de inicios del verano e inicios de invierno. Los resultados muestran la existencia de un desfase temporal entre el riesgo climático y riesgo económico, teniendo este último un periodo de duración más extenso en el tiempo. También, revelan que la intensidad, frecuencia y duración son tres atributos cruciales que determinan el impacto económico de la sequía. La estrategia de reducción de la carga ganadera permite aminorar el riesgo, pero conlleva una disminución en el margen bruto de la explotación. La segunda aplicación está dedicada a la elaboración de mapas de vulnerabilidad a la sequia en pastos de Chile. Para ello, se propone y desarrolla un índice de riesgo económico (IRESP) sencillo de interpretar y replicable, que integra factores de riesgo y estrategias de adaptación para obtener una medida del Valor en Riesgo, es decir, la máxima pérdida esperada en un año con un nivel de significación del 5%.La representación espacial del IRESP pone en evidencia patrones espaciales y diferencias significativas en la vulnerabilidad a la sequía a lo largo de Chile. Además, refleja que la vulnerabilidad no siempre esta correlacionada con el riesgo climático y demuestra la importancia de considerar las estrategias de adaptación. Las medidas de autocorrelación espacial revelan que el riesgo sistémico es considerablemente mayor en el sur que en el resto de zonas. Los resultados demuestran que el IRESP transmite información pertinente y, que los mapas de vulnerabilidad pueden ser una herramienta útil en el diseño de políticas y toma de decisiones para la gestión del riesgo de sequía. La tercera aplicación evalúa el potencial de un seguro indexado para sequía en pastos en la región de Coquimbo en Chile. Para lo cual, se desarrolla un modelo estocástico para estimar la prima actuarialmente justa del seguro y se proponen y evalúan pautas alternativas para mejorar el diseño del contrato. Se aborda el riesgo base, el principal problema de los seguros indexados identificado en la literatura y, que está referido a la correlación imperfecta del índice con las pérdidas de la explotación. Para ello, se sigue un enfoque bayesiano que permite evaluar el impacto en el riesgo base de las pautas de diseño propuestas: i) una zonificación por clúster que considera aspectos espacio-temporales, ii) un período de garantía acotado a los ciclos fenológicos del pasto y iii) umbral de garantía. Los resultados muestran que tanto la zonificación como el periodo de garantía reducen el riesgo base considerablemente. Sin embargo, el umbral de garantía tiene un efecto ambiguo sobre el riesgo base. Por otra parte, la zonificación por clúster contribuye a aminorar el riesgo sistémico que enfrentan las aseguradoras. Estos resultados han puesto de manifiesto que un buen diseño de contrato puede tener un doble dividendo, por un lado aumentar su utilidad y, por otro, reducir el coste del seguro. Un diseño de contrato eficiente junto con los avances en la teledetección y un adecuado marco institucional son los pilares básicos para el buen funcionamiento de un programa de seguro. Las nuevas tecnologías ofrecen un importante potencial para la innovación en la gestión del riesgo climático. Los avances en este campo pueden proporcionar importantes beneficios sociales en los países en desarrollo y regiones vulnerables, donde las herramientas para gestionar eficazmente los riesgos sistémicos como la sequía pueden ser de gran ayuda para el desarrollo. The vulnerability of grazing livestock systems highlights the need for tools to assess and mitigate the adverse impact of drought. The recent and rapid progress in remote sensing has awakened an interest for tapping into potential applications, triggering intensive efforts to develop innovations in a number of spheres. One of these areas is climate risk management, where the use of vegetation indices facilitates assessment of drought. This research analyzes drought impacts and evaluates the potential of new technologies such as remote sensing to manage drought risk in extensive livestock systems. Three essays in drought risk management are developed to: (i) assess the economic impact of drought on a livestock farm in the Andalusian Dehesa, (ii) build drought vulnerability maps in Chilean grazing lands, and (iii) design and evaluate the potential of an index insurance policy to address the risk of drought in grazing lands in Coquimbo, Chile. In the first essay, a dynamic and stochastic farm model is designed combining climate, agronomic, socio-economic and ecological aspects to assess drought risk. The model is developed to simulate a representative livestock farm in the Dehesa of Andalusia for the time period 1999-2010. Burn analysis and MonteCarlo simulation methods are used to identify the significance of various risk sources at the farm. Most notably, early summer and early winter are identified as periods of peak risk. Moreover, there is a significant time lag between climate and economic risk and this later last longer than the former. It is shown that intensity, frequency and duration of the drought are three crucial attributes that shape the economic impact of drought. Sensitivity analysis is conducted to assess the sustainability of farm management strategies and demonstrates that lowering the stocking rate reduces farmer exposure to drought risk but entails a reduction in the expected gross margin. The second essay, mapping drought vulnerability in Chilean grazing lands, proposes and builds an index of economic risk (IRESP) that is replicable and simple to interpret. This methodology integrates risk factors and adaptation strategies to deliver information on Value at Risk, maximum expected losses at 5% significance level. Mapping IRESP provides evidence about spatial patterns and significant differences in drought vulnerability across Chilean grazing lands. Spatial autocorrelation measures reveal that systemic risk is considerably larger in the South as compared to Northern or Central Regions. Furthermore, it is shown that vulnerability is not necessarily correlated with climate risk and that adaptation strategies do matter. These results show that IRESP conveys relevant information and that vulnerability maps may be useful tools to assess policy design and decision-making in drought risk management. The third essay develops a stochastic model to estimate the actuarially fair premium and evaluates the potential of an indexed insurance policy to manage drought risk in Coquimbo, a relevant livestock farming region of Chile. Basis risk refers to the imperfect correlation of the index and farmer loses and is identified in the literature as a main limitation of index insurance. A Bayesian approach is proposed to assess the impact on basis risk of alternative guidelines in contract design: i) A cluster zoning that considers space-time aspects, ii) A guarantee period bounded to fit phenological cycles, and iii) the triggering index threshold. Results show that both the proposed zoning and guarantee period considerably reduces basis risk. However, the triggering index threshold has an ambiguous effect on basis risk. On the other hand, cluster zoning contributes to ameliorate systemic risk faced by the insurer. These results highlighted that adequate contract design is important and may result in double dividend. On the one hand, increasing farmers’ utility and, secondly, reducing the cost of insurance. An efficient contract design coupled with advances in remote sensing and an appropriate institutional framework are the basis for an efficient operation of an insurance program. The new technologies offer significant potential for innovation in climate risk managements. Progress in this field is capturing increasing attention and may provide important social gains in developing countries and vulnerable regions where the tools to efficiently manage systemic risks, such as drought, may be a means to foster development.

Diversificação de espécies e da morfologia em serpentes da família Viperidae: padrões e processos

A diversidade de espécies e fenotípica pode variar consideravelmente entre grupos taxonômicos e ao longo do tempo em uma mesma linhagem. O estudo de tais variações tornou-se um dos principais objetivos da biologia evolutiva fornecendo informações importantes a respeito dos possíveis mecanismos que regulam a biodiversidade. Dessa forma, o objetivo geral da presente tese foi investigar os padrões da diversificação de espécies e da morfologia em um grupo cosmopolita de serpentes, a família Viperidae, e os potenciais processos subjacentes. Primeiramente, (1) reconstruímos as relações filogenéticas e estimamos os tempos de divergência entre as linhagens da família Viperidae utilizando uma abordagem Bayesiana. (2) Aplicando um método recentemente desenvolvido (BAMM), exploramos como as taxas de especiação e extinção variaram ao longo da radiação do grupo inferindo os possíveis processos reguladores. Por fim, (3) analisamos se a evolução do tamanho do corpo e as taxas de especiação variam nos diferentes habitats ocupados pelos viperídeos (terrestres vs arborícola). Nesta tese geramos a filogenia molecular de viperídeos mais completa até o momento utilizando sequências para 11 genes mitocondriais e nucleares abrangendo 79% das espécies viventes (264 terminais) e todos com exceção de um gênero. De maneira geral, foi possível obter relações filogenéticas robustas para o grupo com a maioria dos gêneros sendo monofilética. Os tempos de divergência obtidos indicam que os viperídeos começaram a diversificar em meados do Paleoceno tardio/meio do Eoceno inferindo idades um pouco mais tardias que o encontrado em estudos anteriores. Durante a radiação do grupo, um aumento nas taxas de especiação parece ter ocorrido durante a diversificação dos crotalíneos (pit vipers) em decorrência não só da evolução das fossetas loreais mas também como resultado de mudanças geológicas e climáticas na Ásia e da invasão do novo mundo. Após este rápido aumento inicial, as taxas de especiação desaceleraram em direção ao presente. Por fim, os resultados aqui apresentados indicam que apesar dos habitats arborícolas limitarem a evolução morfológica nos viperídeos, a evolução da arborealidade parece não afetar as taxas de especiação que permanecem similares entre linhagens arborícolas e terrestres. Isto sugere dois cenários: (1) a especiação acontece de forma independente das mudanças morfológicas nos viperídeos; ou (2) o isolamento geográfico seria um mecanismo importante na diversificação de linhagens arborícolas contrabalançando decréscimos nas oportunidades de especiação possivelmente relacionados às pressões seletivas impostas pelo ambiente arborícola. A presente tese contribui para entendermos mais sobre como evoluíram os viperídeos ao longo dos seus ∼50 milhões de anos. Além de propor cenários e hipóteses a serem futuramente explorados com os viperídeos, elaboramos uma discussão ampla e conceitual a respeito dos possíveis mecanismos por trás da diversificação de espécies e da morfologia que poderiam também ser contemplados para outros grupos de organismos. Portanto, a presente tese contribui não só para entendermos os mecanismos que geram e mantém a diversidade de serpentes, mas também para enriquecer a discussão dos mecanismos que geram e mantém a biodiversidade como um todo

Análise de portfólio: uma perspectiva bayesiana

Este trabalho tem com objetivo abordar o problema de alocação de ativos (análise de portfólio) sob uma ótica Bayesiana. Para isto foi necessário revisar toda a análise teórica do modelo clássico de média-variância e na sequencia identificar suas deficiências que comprometem sua eficácia em casos reais. Curiosamente, sua maior deficiência não esta relacionado com o próprio modelo e sim pelos seus dados de entrada em especial ao retorno esperado calculado com dados históricos. Para superar esta deficiência a abordagem Bayesiana (modelo de Black-Litterman) trata o retorno esperado como uma variável aleatória e na sequência constrói uma distribuição a priori (baseado no modelo de CAPM) e uma distribuição de verossimilhança (baseado na visão de mercado sob a ótica do investidor) para finalmente aplicar o teorema de Bayes tendo como resultado a distribuição a posteriori. O novo valor esperado do retorno, que emerge da distribuição a posteriori, é que substituirá a estimativa anterior do retorno esperado calculado com dados históricos. Os resultados obtidos mostraram que o modelo Bayesiano apresenta resultados conservadores e intuitivos em relação ao modelo clássico de média-variância.

Genetic analysis of complex demographic scenarios: Spatially expanding populations of the cane toad, Bufo marinus

Inferring the spatial expansion dynamics of invading species from molecular data is notoriously difficult due to the complexity of the processes involved. For these demographic scenarios, genetic data obtained from highly variable markers may be profitably combined with specific sampling schemes and information from other sources using a Bayesian approach. The geographic range of the introduced toad Bufo marinus is still expanding in eastern and northern Australia, in each case from isolates established around 1960. A large amount of demographic and historical information is available on both expansion areas. In each area, samples were collected along a transect representing populations of different ages and genotyped at 10 microsatellite loci. Five demographic models of expansion, differing in the dispersal pattern for migrants and founders and in the number of founders, were considered. Because the demographic history is complex, we used an approximate Bayesian method, based on a rejection-regression algorithm. to formally test the relative likelihoods of the five models of expansion and to infer demographic parameters. A stepwise migration-foundation model with founder events was statistically better supported than other four models in both expansion areas. Posterior distributions supported different dynamics of expansion in the studied areas. Populations in the eastern expansion area have a lower stable effective population size and have been founded by a smaller number of individuals than those in the northern expansion area. Once demographically stabilized, populations exchange a substantial number of effective migrants per generation in both expansion areas, and such exchanges are larger in northern than in eastern Australia. The effective number of migrants appears to be considerably lower than that of founders in both expansion areas. We found our inferences to be relatively robust to various assumptions on marker. demographic, and historical features. The method presented here is the only robust, model-based method available so far, which allows inferring complex population dynamics over a short time scale. It also provides the basis for investigating the interplay between population dynamics, drift, and selection in invasive species.

Pharmacokinetic-pharmacodynamic modelling of QT interval prolongation following citalopram overdoses

Aims To develop a pharmacokinetic-pharmacodynamic model describing the time-course of QT interval prolongation after citalopram overdose and to evaluate the effect of charcoal on the relative risk of developing abnormal QT and heart-rate combinations. Methods Plasma concentrations and electrocardiograph (ECG) data from 52 patients after 62 citalopram overdose events were analysed in WinBUGS using a Bayesian approach. The reported doses ranged from 20 to 1700 mg and on 17 of the events a single dose of activated charcoal was administered. The developed pharmacokinetic-pharmacodynamic model was used for predicting the probability of having abnormal combinations of QT-RR, which was assumed to be related to an increased risk for torsade de pointes (TdP). Results The absolute QT interval was related to the observed heart rate with an estimated individual heart-rate correction factor [alpha = 0.36, between-subject coefficient of variation (CV) = 29%]. The heart-rate corrected QT interval was linearly dependent on the predicted citalopram concentration (slope = 40 ms l mg(-1), between-subject CV = 70%) in a hypothetical effect-compartment (half-life of effect-delay = 1.4 h). The heart-rate corrected QT was predicted to be higher in women than in men and to increase with age. Administration of activated charcoal resulted in a pronounced reduction of the QT prolongation and was shown to reduce the risk of having abnormal combinations of QT-RR by approximately 60% for citalopram doses above 600 mg. Conclusion Citalopram caused a delayed lengthening of the QT interval. Administration of activated charcoal was shown to reduce the risk that the QT interval exceeds a previously defined threshold and therefore is expected to reduce the risk of TdP.

Efficiency effects of agricultural economics research in the United States

Allocations of research funds across programs are often made for efficiency reasons. Social science research is shown to have small, lagged but significant effects on U.S. agricultural efficiency when public agricultural R&D and extension are simultaneously taken into account. Farm management and marketing research variables are used to explain variations in estimates of allocative and technical efficiency using a Bayesian approach that incorporates stylized facts concerning lagged research impacts in a way that is less restrictive than popular polynomial distributed lags. Results are reported in terms of means and standard deviations of estimated probability distributions of parameters and long-run total multipliers. Extension is estimated to have a greater impact on both allocative and technical efficiency than either R&D or social science research.