Evaluating sampling strategies and logistic regression methods for modelling complex land cover changes

The role of land cover change as a significant component of global change has become increasingly recognized in recent decades. Large databases measuring land cover change, and the data which can potentially be used to explain the observed changes, are also becoming more commonly available. When developing statistical models to investigate observed changes, it is important to be aware that the chosen sampling strategy and modelling techniques can influence results. We present a comparison of three sampling strategies and two forms of grouped logistic regression models (multinomial and ordinal) in the investigation of patterns of successional change after agricultural land abandonment in Switzerland. Results indicated that both ordinal and nominal transitional change occurs in the landscape and that the use of different sampling regimes and modelling techniques as investigative tools yield different results. Synthesis and applications. Our multimodel inference identified successfully a set of consistently selected indicators of land cover change, which can be used to predict further change, including annual average temperature, the number of already overgrown neighbouring areas of land and distance to historically destructive avalanche sites. This allows for more reliable decision making and planning with respect to landscape management. Although both model approaches gave similar results, ordinal regression yielded more parsimonious models that identified the important predictors of land cover change more efficiently. Thus, this approach is favourable where land cover change pattern can be interpreted as an ordinal process. Otherwise, multinomial logistic regression is a viable alternative.

Bayesian Model Averaging in the Instrumental Variable Regression Model

This paper considers the instrumental variable regression model when there is uncertainty about the set of instruments, exogeneity restrictions, the validity of identifying restrictions and the set of exogenous regressors. This uncertainty can result in a huge number of models. To avoid statistical problems associated with standard model selection procedures, we develop a reversible jump Markov chain Monte Carlo algorithm that allows us to do Bayesian model averaging. The algorithm is very exible and can be easily adapted to analyze any of the di¤erent priors that have been proposed in the Bayesian instrumental variables literature. We show how to calculate the probability of any relevant restriction (e.g. the posterior probability that over-identifying restrictions hold) and discuss diagnostic checking using the posterior distribution of discrepancy vectors. We illustrate our methods in a returns-to-schooling application.

Problèmes d'échantillonnage lors d'une enquête sur la plombémie d'enfants citadin [Sampling problems in a survey of blood lead levels in city children].

A study on lead pollution was carried out on a sample of ca. 300 city children. This paper presents the errors producing bias in the sample. It is emphasized that, in Switzerland, the difference between the Swiss and the migrant population (the latter being mainly Italian and Spanish) must be taken into account.

Patrones de tomate: resistencia variable frente al nematodo Meloidogyne

Los patrones de tomate se han propuesto como una alternativa no química al uso del bromuro de metilo. En este estudio se evaluó la respuesta de resistencia de diez patrones de tomate (comerciales y experimentales) a nematodos del género Meloidogyne mediante un ensayo de campo realizado en un invernadero de plástico cuyo suelo se hallaba infestado por Meloidogyne javanica. Al finalizar la campaña agrícola (marzo a julio), siete de los patrones ensayados mostraron altos niveles de resistencia, uno presentó resistencia intermedia y dos resultaron ser susceptibles al nematodo.

Investigation of Postmortem Gases: A Forensic Imaging and Analytical Chemistry Approach

Background: Distinguishing postmortem gas accumulations in the body due to natural decomposition and other phenomena such as gas embolism can prove a difficult task using purely Multi-Detector Computed Tomography (MDCT). The Radiological Alteration Index (RAI) was created with the intention to be able to identify bodies undergoing the putrefaction process based on the quantity of gas detected within the body. The flaw in this approach is the inability to absolutely determine putrefaction as the origin of gas volumes in cases of moderate alteration. The aim of the current study is to identify percentage compositions of O2, N2, CO2 and the presence of gases such as H2 and H2S within these sampling sites in order to resolve this complication. Materials and methods: All cases investigated in our University Center of Legal Medicine are undergoing a Post-Mortem Computed Tomography (PMCT)-scan before external examination or autopsy as a routine investigation. In the obtained images, areas of gas were characterized as 0, I, II or III based on the amount of gas present according to the RAI (1). The criteria for these characterizations were dependent of the site of gas, for example thoracic and abdominal cavities were graded as I (1 - 3cm gas), II (3 - 5cm gas) and III (>5cm gas). Cases showing gaseous sites with grade II or III were selected for this study. The sampling was performed under CT-guidance to target the regions to be punctured. Luer-lock PTFE syringes equipped with a three-way valve and needles were used to sample the gas directly (2). Gaseous samples were then analysed using gas chromatography coupled to a thermal conductivity detector (GC-TCD). The components present in the samples were expressed as a percentage of the overall gas present. Results: Up to now, we have investigated more than 40 cases using our standardized procedure for sampling and analysis of gas. O2, N2 and CO2 were present in most samples. The following distributions were found to correlate to gas origins of gas embolism/scuba diving accidents, trauma and putrefaction: ? Putrefaction → O2 = 1 - 5%; CO2 > 15%; N2 = 10 - 70%; H2 / H2S / CH4 variable presence ? Gas embolism/Scuba diving accidents → O2 and N2= varying percentages; CO2 > 20% ? Trauma → O2 = small percentage; CO2 < 15%; N2 > 65% H2 and H2S indicated levels of putrefaction along with methane which can also gauge environmental conditions or conditions of body storage/burial. Many cases showing large RAI values (advanced alteration) did reveal a radiological diagnosis which was in concordance with the interpretation of the gas composition. However, in certain cases (gas embolism, scuba divers) radiological interpretation was not possible and only chemical gas analysis was found to lead to the correct diagnosis, meaning that it provided complementary information to the radiological diagnosis. Conclusion: Investigation of postmortem gases is a useful tool to determine origin of gas generation which can aid the diagnosis of the cause of death. Levels of gas can provide information on stage of putrefaction and help to perform essential medico-legal diagnosis such as vital gas embolism.

Variable viral clearance despite adequate ganciclovir plasma levels during valganciclovir treatment for cytomegalovirus disease in D+/R- transplant recipients.

ABSTRACT: BACKGROUND: Valganciclovir, the oral prodrug of ganciclovir, has been demonstrated equivalent to iv ganciclovir for CMV disease treatment in solid organ transplant recipients. Variability in ganciclovir exposure achieved with valganciclovir could be implicated as a contributing factor for explaining variations in the therapeutic response. This prospective observational study aimed to correlate clinical and cytomegalovirus (CMV) viral load response (DNAemia) with ganciclovir plasma concentrations in patients treated with valganciclovir for CMV infection/disease. METHODS: Seven CMV D+/R- transplant recipients (4 kidney, 2 liver and 1 heart) were treated with valganciclovir (initial dose was 900-1800 mg/day for 3-6.5 weeks, followed by 450-900 mg/day for 2-9 weeks). DNAemia was monitored by real time quantitative PCR and ganciclovir plasma concentration was measured at trough (Ctrough) and 3 h after drug administration (C3h) by HPLC. RESULTS: Four patients presented with CMV syndrome, two had CMV tissue-invasive disease after prophylaxis discontinuation, and one liver recipient was treated pre-emptively for asymptomatic rising CMV viral load 5 weeks post-transplantation in the absence of prophylaxis. CMV DNAemia decreased during the first week of treatment in all recipients except in one patient (median decrease: -1.2 log copies/mL, range: -1.8 to 0) despite satisfactory ganciclovir exposure (AUC0-12 = 48 mg.h/L, range for the 7 patients: 40-118 mg.h/L). Viral clearance was obtained in five patients after a median of time of 34 days (range: 28-82 days). Two patients had recurrent CMV disease despite adequate ganciclovir exposure (65 mg.h/L, range: 44-118 mg.h/L). CONCLUSIONS: Valganciclovir treatment for CMV infection/disease in D+/R- transplant recipients can thus result in variable viral clearance despite adequate ganciclovir plasma concentrations, probably correlating inversely with anti-CMV immune responses after primary infection.

Estimation of dynamic latent variable models using simulated nonparametric moments

Abstract. Given a model that can be simulated, conditional moments at a trial parameter value can be calculated with high accuracy by applying kernel smoothing methods to a long simulation. With such conditional moments in hand, standard method of moments techniques can be used to estimate the parameter. Because conditional moments are calculated using kernel smoothing rather than simple averaging, it is not necessary that the model be simulable subject to the conditioning information that is used to define the moment conditions. For this reason, the proposed estimator is applicable to general dynamic latent variable models. It is shown that as the number of simulations diverges, the estimator is consistent and a higher-order expansion reveals the stochastic difference between the infeasible GMM estimator based on the same moment conditions and the simulated version. In particular, we show how to adjust standard errors to account for the simulations. Monte Carlo results show how the estimator may be applied to a range of dynamic latent variable (DLV) models, and that it performs well in comparison to several other estimators that have been proposed for DLV models.

Decision-theoretic analysis of forensic sampling criteria using Bayesian decision networks

Sampling issues represent a topic of ongoing interest to the forensic science community essentially because of their crucial role in laboratory planning and working protocols. For this purpose, forensic literature described thorough (Bayesian) probabilistic sampling approaches. These are now widely implemented in practice. They allow, for instance, to obtain probability statements that parameters of interest (e.g., the proportion of a seizure of items that present particular features, such as an illegal substance) satisfy particular criteria (e.g., a threshold or an otherwise limiting value). Currently, there are many approaches that allow one to derive probability statements relating to a population proportion, but questions on how a forensic decision maker - typically a client of a forensic examination or a scientist acting on behalf of a client - ought actually to decide about a proportion or a sample size, remained largely unexplored to date. The research presented here intends to address methodology from decision theory that may help to cope usefully with the wide range of sampling issues typically encountered in forensic science applications. The procedures explored in this paper enable scientists to address a variety of concepts such as the (net) value of sample information, the (expected) value of sample information or the (expected) decision loss. All of these aspects directly relate to questions that are regularly encountered in casework. Besides probability theory and Bayesian inference, the proposed approach requires some additional elements from decision theory that may increase the efforts needed for practical implementation. In view of this challenge, the present paper will emphasise the merits of graphical modelling concepts, such as decision trees and Bayesian decision networks. These can support forensic scientists in applying the methodology in practice. How this may be achieved is illustrated with several examples. The graphical devices invoked here also serve the purpose of supporting the discussion of the similarities, differences and complementary aspects of existing Bayesian probabilistic sampling criteria and the decision-theoretic approach proposed throughout this paper.

Large multi-gene phylogenetic trees of the grasses (Poaceae): progress towards complete tribal and generic level sampling.

In this paper we included a very broad representation of grass family diversity (84% of tribes and 42% of genera). Phylogenetic inference was based on three plastid DNA regions rbcL, matK and trnL-F, using maximum parsimony and Bayesian methods. Our results resolved most of the subfamily relationships within the major clades (BEP and PACCMAD), which had previously been unclear, such as, among others the: (i) BEP and PACCMAD sister relationship, (ii) composition of clades and the sister-relationship of Ehrhartoideae and Bambusoideae + Pooideae, (iii) paraphyly of tribe Bambuseae, (iv) position of Gynerium as sister to Panicoideae, (v) phylogenetic position of Micrairoideae. With the presence of a relatively large amount of missing data, we were able to increase taxon sampling substantially in our analyses from 107 to 295 taxa. However, bootstrap support and to a lesser extent Bayesian inference posterior probabilities were generally lower in analyses involving missing data than those not including them. We produced a fully resolved phylogenetic summary tree for the grass family at subfamily level and indicated the most likely relationships of all included tribes in our analysis.

Prevalence and determinants of QT interval prolongation in medical inpatients.

BACKGROUND: QT interval prolongation carries an increased risk of torsade de pointes and death. AIM: We sought to determine the prevalence of QT prolongation in medical inpatients and to identify determinants of this condition. METHODS: We enrolled consecutive patients who were admitted to the internal medicine ward and who had an electrocardiogram performed within 24 h of admission. We collected information on baseline patient characteristics and the use of QT-prolonging drugs. Two blinded readers manually measured the QT intervals. QT intervals were corrected for heart rate using the traditional Bazett formula and the linear regression-based Framingham formula. We used logistic regression to identify patient characteristics and drugs that were independently associated with QTc prolongation. RESULTS: Of 537 inpatients, 22.3% had a prolonged QTc based on the Bazett formula. The adjusted odds for QTc prolongation based on the Bazett correction were significantly higher in patients who had liver disease (OR 2.9, 95% CI: 1.5-5.6), hypokalaemia (OR 3.3, 95% CI: 1.9-5.6) and who were taking ≥1 QT-prolonging drug at admission (OR 1.7, 95% CI: 1.1-2.6). Overall, 50.8% of patients with QTc prolongation received additional QT-prolonging drugs during hospitalisation. CONCLUSIONS: The prevalence of QTc prolongation was high among medical inpatients but depended on the method used to correct for heart rate. The use of QT-prolonging drugs, hypokalaemia and liver disease increased the risk of QTc prolongation. Many patients with QTc prolongation received additional QT-prolonging drugs during hospitalisation, further increasing the risk of torsade de pointes and death.

Variable phenotypic expressivity in a Swiss family with autosomal dominant retinitis pigmentosa due to mutation T494m in the Prpf3 gene

Purpose:to describe the clinical features in a five generations family segregating autosomal dominant retinitis pigmentosa and to identify the causative gene Patient and Methods:Twenty five individuals of a large five-generation family originating from Western Switzerland were ascertained for phenotypic and genotypic characterization. Ophthalmologic evaluations included color vision testing, Goldman perimetry and digital fundus photography. Some patients had autofluorescence (AF) imaging, ocular coherence tomography (OCT) and ISCEV-standard full-field electroretinography (ERG). Blood samples were collected from 10 affected (4 to 70 years of age) and 15 unaffected members after informed consent. DNA was isolated and exons and intron-exons junctions of known adRP genes were sequenced using a Big Dye sequencing kit 1.1. Results:Age of onset of nightblindness and severity of progression of the disease was variable between members of the family. Some patients had early onset of nightblindess aged 3, others at mid-twenties. Most patients had visual acuity above 0.6 for the first 4 decades. Two older patients still had good vision (0.4) in their seventies. Myopia (range: -2 to -5) was noticed in most affected subjects. Fundus findings showed areas of atrophy along the arcades. The AF imaging showed a large high density ring bilaterally. A T494M change was found in exon 11 of PRPF3 gene. The change segregates with the disease in the family. Conclusion: A mutation in the PRPF3 gene is rare compared with other genes causing ADRP. Although a T494M change has been reported, our family is the first one with a variable expressivity. Mutations in PRPF3 gene can cause a variable phenotype of ADRP unlike the previously described Danish and English families. Our report gives a better understanding as to the phenotype/genotype description of ADRP due to PRPF3 mutation.

Birth records from Swiss married couples analyzed over of the past 35 years reveal an aging of first-time mothers by 5.1 years while the interpregnancy interval has shortened

RésuméIntroduction : Le travail de thèse est un article publié dans le journal « Fertility & Sterility » s'intitulant « Birth records from Swiss married couples analyzed over the past 35 years reveal an aging of first-time mothers by 5.1 years while the interpregnancy interval has shortened ».Méthodes : Les données concernant l'âge auquel les femmes mariées donnaient naissance ont été obtenues de l'Office fédéral de la statistique. Nous avons examiné la période allant de 1969 à 2006. Cet intervalle de temps choisi nous a permis de prendre en considération un total de 2'716'370 naissances. L'âge moyen des parturientes à la naissance de leur 1er, 2ème et 3ème enfant a été calculé pour chaque année à l'aide du logiciel Excel. Grâce à ces données on a pu obtenir les intervalles inter-gestationnels théoriques entre le 1er et 2ème enfant, ainsi qu'entre le 2ème et 3ème enfant.Résultats : Nous pouvons constater que l'intervalle inter-gestationnel théorique entre la première et la deuxième grossesse était de 23.2 mois en 1969 et passait à 13.0 mois en 2006. L'intervalle compris entre la deuxième et la troisième grossesse passait de 22.4 mois en 1969 à 7.9 mois en 2006. Notre analyse suggère donc que c'est le facteur social qui exerce un effet plus important que le facteur biologique sur les intervalles inter-gestationnels, car ces intervalles diminuaient malgré l'augmentation de l'âge des mères à la naissance.

Comparative Histopathology of Biomphalaria glabrata, B. tenagophila and B. straminea with Variable Degrees of Resistance to Schistosoma mansoni Miracidia

A comparative histopathological study of three snails species - Biomphalaria glabrata, B. tenagophila and B. straminea - which had been infected with Schistosoma mansoni miracidia revealed similar qualitative features, consisting of areas of sporocyst proliferation and differentiation associated with reactive host reaction, at the time they were actively eliminating great number of cercariae. However, in specimens that were exposed to miracidia but failed to eliminate cercariae later on, different histopathological pictures were observed in different snail species. While B. glabrata exhibited frequent focal (granulomatous) proliferation of amebocytes in several organs, B. tenagophila and B. straminea only rarely showed such reactive changes, suggesting that the mechanism of resistance to miracidial infection probably follows different pathways in the snail species studied

Two monoclonal rat antibodies with specificity for the beta-chain variable region V beta 6 of the murine T-cell receptor.

Two rat monoclonal antibodies (mAbs), 44-22-1 and 46-6B5, which recognize an alloreactive cytotoxic clone, 3F9, have been further tested on a panel of T hybridomas and cytotoxic T-cell clones for binding and functional activities. The mAbs recognized only those cells sharing the expression of the T-cell receptor beta-chain variable region gene V beta 6 with 3F9. All V beta 6+ cells were activated by these mAbs under cross-linking conditions and their antigen-specific activation was blocked by soluble mAb. Furthermore, depletion of 46-6B5+ normal lymph node T cells eliminated all cells expressing the epitope recognized by 44-22-1 and V beta 6 mRNA.