Consistent estimation of the fixed effects ordered logit model

The paper considers panel data methods for estimating ordered logit models with individual-specific correlated unobserved heterogeneity. We show that a popular approach is inconsistent, whereas some consistent and efficient estimators are available, including minimum distance and generalized method-of-moment estimators. A Monte Carlo study reveals the good properties of an alternative estimator that has not been considered in econometric applications before, is simple to implement and almost as efficient. An illustrative application based on data from the German Socio-Economic Panel confirms the large negative effect of unemployment on life satisfaction that has been found in the previous literature.

On ANOVA Decompositions of Kernels and Gaussian Random Field Paths

The FANOVA (or “Sobol’-Hoeffding”) decomposition of multivariate functions has been used for high-dimensional model representation and global sensitivity analysis. When the objective function f has no simple analytic form and is costly to evaluate, computing FANOVA terms may be unaffordable due to numerical integration costs. Several approximate approaches relying on Gaussian random field (GRF) models have been proposed to alleviate these costs, where f is substituted by a (kriging) predictor or by conditional simulations. Here we focus on FANOVA decompositions of GRF sample paths, and we notably introduce an associated kernel decomposition into 4 d 4d terms called KANOVA. An interpretation in terms of tensor product projections is obtained, and it is shown that projected kernels control both the sparsity of GRF sample paths and the dependence structure between FANOVA effects. Applications on simulated data show the relevance of the approach for designing new classes of covariance kernels dedicated to high-dimensional kriging.

Mortality after failure of antiretroviral therapy in sub-Saharan Africa

Objective  To assess the outcome of patients who experienced treatment failure with antiretrovirals in sub-Saharan Africa. Methods  Analysis of 11 antiretroviral therapy (ART) programmes in sub-Saharan Africa. World Health Organization (WHO) criteria were used to define treatment failure. All ART-naive patients aged ≥16 who started with a non-nucleoside reverse transcriptase inhibitor (NNRTI)-based regimen and had at least 6 months of follow-up were eligible. For each patient who switched to a second-line regimen, 10 matched patients who remained on a non-failing first-line regimen were selected. Time was measured from the time of switching, from the corresponding time in matched patients, or from the time of treatment failure in patients who remained on a failing regimen. Mortality was analysed using Kaplan–Meier curves and random-effects Cox models. Results  Of 16 591 adult patients starting ART, 382 patients (2.3%) switched to a second-line regimen. Another 323 patients (1.9%) did not switch despite developing immunological or virological failure. Cumulative mortality at 1 year was 4.2% (95% CI 2.2–7.8%) in patients who switched to a second-line regimen and 11.7% (7.3%–18.5%) in patients who remained on a failing first-line regimen, compared to 2.2% (1.6–3.0%) in patients on a non-failing first-line regimen (P < 0.0001). Differences in mortality were not explained by nadir CD4 cell count, age or differential loss to follow up. Conclusions  Many patients who meet criteria for treatment failure do not switch to a second-line regimen and die. There is an urgent need to clarify the reasons why in sub-Saharan Africa many patients remain on failing first-line ART.

Electronic medical record systems, data quality and loss to follow-up: survey of antiretroviral therapy programmes in resource-limited settings

OBJECTIVE: To describe the electronic medical databases used in antiretroviral therapy (ART) programmes in lower-income countries and assess the measures such programmes employ to maintain and improve data quality and reduce the loss of patients to follow-up. METHODS: In 15 countries of Africa, South America and Asia, a survey was conducted from December 2006 to February 2007 on the use of electronic medical record systems in ART programmes. Patients enrolled in the sites at the time of the survey but not seen during the previous 12 months were considered lost to follow-up. The quality of the data was assessed by computing the percentage of missing key variables (age, sex, clinical stage of HIV infection, CD4+ lymphocyte count and year of ART initiation). Associations between site characteristics (such as number of staff members dedicated to data management), measures to reduce loss to follow-up (such as the presence of staff dedicated to tracing patients) and data quality and loss to follow-up were analysed using multivariate logit models. FINDINGS: Twenty-one sites that together provided ART to 50 060 patients were included (median number of patients per site: 1000; interquartile range, IQR: 72-19 320). Eighteen sites (86%) used an electronic database for medical record-keeping; 15 (83%) such sites relied on software intended for personal or small business use. The median percentage of missing data for key variables per site was 10.9% (IQR: 2.0-18.9%) and declined with training in data management (odds ratio, OR: 0.58; 95% confidence interval, CI: 0.37-0.90) and weekly hours spent by a clerk on the database per 100 patients on ART (OR: 0.95; 95% CI: 0.90-0.99). About 10 weekly hours per 100 patients on ART were required to reduce missing data for key variables to below 10%. The median percentage of patients lost to follow-up 1 year after starting ART was 8.5% (IQR: 4.2-19.7%). Strategies to reduce loss to follow-up included outreach teams, community-based organizations and checking death registry data. Implementation of all three strategies substantially reduced losses to follow-up (OR: 0.17; 95% CI: 0.15-0.20). CONCLUSION: The quality of the data collected and the retention of patients in ART treatment programmes are unsatisfactory for many sites involved in the scale-up of ART in resource-limited settings, mainly because of insufficient staff trained to manage data and trace patients lost to follow-up.

Switching to second-line antiretroviral therapy in resource-limited settings: comparison of programmes with and without viral load monitoring

BACKGROUND: In high-income countries, viral load is routinely measured to detect failure of antiretroviral therapy (ART) and guide switching to second-line ART. Viral load monitoring is not generally available in resource-limited settings. We examined switching from nonnucleoside reverse transcriptase inhibitor (NNRTI)-based first-line regimens to protease inhibitor-based regimens in Africa, South America and Asia. DESIGN AND METHODS: Multicohort study of 17 ART programmes. All sites monitored CD4 cell count and had access to second-line ART and 10 sites monitored viral load. We compared times to switching, CD4 cell counts at switching and obtained adjusted hazard ratios for switching (aHRs) with 95% confidence intervals (CIs) from random-effects Weibull models. RESULTS: A total of 20 113 patients, including 6369 (31.7%) patients from 10 programmes with access to viral load monitoring, were analysed; 576 patients (2.9%) switched. Low CD4 cell counts at ART initiation were associated with switching in all programmes. Median time to switching was 16.3 months [interquartile range (IQR) 10.1-26.6] in programmes with viral load monitoring and 21.8 months (IQR 14.0-21.8) in programmes without viral load monitoring (P < 0.001). Median CD4 cell counts at switching were 161 cells/microl (IQR 77-265) in programmes with viral load monitoring and 102 cells/microl (44-181) in programmes without viral load monitoring (P < 0.001). Switching was more common in programmes with viral load monitoring during months 7-18 after starting ART (aHR 1.38; 95% CI 0.97-1.98), similar during months 19-30 (aHR 0.97; 95% CI 0.58-1.60) and less common during months 31-42 (aHR 0.29; 95% CI 0.11-0.79). CONCLUSION: In resource-limited settings, switching to second-line regimens tends to occur earlier and at higher CD4 cell counts in ART programmes with viral load monitoring compared with programmes without viral load monitoring.

Pre-therapy inflammation and coagulation activation and long-term CD4 count responses to the initiation of antiretroviral therapy

OBJECTIVES Pre-antiretroviral therapy (ART) inflammation and coagulation activation predict clinical outcomes in HIV-positive individuals. We assessed whether pre-ART inflammatory marker levels predicted the CD4 count response to ART. METHODS Analyses were based on data from the Strategic Management of Antiretroviral Therapy (SMART) trial, an international trial evaluating continuous vs. interrupted ART, and the Flexible Initial Retrovirus Suppressive Therapies (FIRST) trial, evaluating three first-line ART regimens with at least two drug classes. For this analysis, participants had to be ART-naïve or off ART at randomization and (re)starting ART and have C-reactive protein (CRP), interleukin-6 (IL-6) and D-dimer measured pre-ART. Using random effects linear models, we assessed the association between each of the biomarker levels, categorized as quartiles, and change in CD4 count from ART initiation to 24 months post-ART. Analyses adjusted for CD4 count at ART initiation (baseline), study arm, follow-up time and other known confounders. RESULTS Overall, 1084 individuals [659 from SMART (26% ART naïve) and 425 from FIRST] met the eligibility criteria, providing 8264 CD4 count measurements. Seventy-five per cent of individuals were male with the mean age of 42 years. The median (interquartile range) baseline CD4 counts were 416 (350-530) and 100 (22-300) cells/μL in SMART and FIRST, respectively. All of the biomarkers were inversely associated with baseline CD4 count in FIRST but not in SMART. In adjusted models, there was no clear relationship between changing biomarker levels and mean change in CD4 count post-ART (P for trend: CRP, P = 0.97; IL-6, P = 0.25; and D-dimer, P = 0.29). CONCLUSIONS Pre-ART inflammation and coagulation activation do not predict CD4 count response to ART and appear to influence the risk of clinical outcomes through other mechanisms than blunting long-term CD4 count gain.

Reproducibility and accuracy of linear measurements on dental models derived from cone-beam computed tomography compared with digital dental casts.

INTRODUCTION The aim of this study was to determine the reproducibility and accuracy of linear measurements on 2 types of dental models derived from cone-beam computed tomography (CBCT) scans: CBCT images, and Anatomodels (InVivoDental, San Jose, Calif); these were compared with digital models generated from dental impressions (Digimodels; Orthoproof, Nieuwegein, The Netherlands). The Digimodels were used as the reference standard. METHODS The 3 types of digital models were made from 10 subjects. Four examiners repeated 37 linear tooth and arch measurements 10 times. Paired t tests and the intraclass correlation coefficient were performed to determine the reproducibility and accuracy of the measurements. RESULTS The CBCT images showed significantly smaller intraclass correlation coefficient values and larger duplicate measurement errors compared with the corresponding values for Digimodels and Anatomodels. The average difference between measurements on CBCT images and Digimodels ranged from -0.4 to 1.65 mm, with limits of agreement values up to 1.3 mm for crown-width measurements. The average difference between Anatomodels and Digimodels ranged from -0.42 to 0.84 mm with limits of agreement values up to 1.65 mm. CONCLUSIONS Statistically significant differences between measurements on Digimodels and Anatomodels, and between Digimodels and CBCT images, were found. Although the mean differences might be clinically acceptable, the random errors were relatively large compared with corresponding measurements reported in the literature for both Anatomodels and CBCT images, and might be clinically important. Therefore, with the CBCT settings used in this study, measurements made directly on CBCT images and Anatomodels are not as accurate as measurements on Digimodels.

Transmission of Chlamydia trachomatis through sexual partnerships: a comparison between three individual-based models and empirical data

Chlamydia trachomatis is the most common bacterial sexually transmitted infection (STI) in many developed countries. The highest prevalence rates are found among young adults who have frequent partner change rates. Three published individual-based models have incorporated a detailed description of age-specific sexual behaviour in order to quantify the transmission of C. trachomatis in the population and to assess the impact of screening interventions. Owing to varying assumptions about sexual partnership formation and dissolution and the great uncertainty about critical parameters, such models show conflicting results about the impact of preventive interventions. Here, we perform a detailed evaluation of these models by comparing the partnership formation and dissolution dynamics with data from Natsal 2000, a population-based probability sample survey of sexual attitudes and lifestyles in Britain. The data also allow us to describe the dispersion of C. trachomatis infections as a function of sexual behaviour, using the Gini coefficient. We suggest that the Gini coefficient is a useful measure for calibrating infectious disease models that include risk structure and highlight the need to estimate this measure for other STIs.

Metabolomics reveals attenuation of the SLC6A20 kidney transporter in nonhuman primate and mouse models of type 2 diabetes mellitus

To enhance understanding of the metabolic indicators of type 2 diabetes mellitus (T2DM) disease pathogenesis and progression, the urinary metabolomes of well characterized rhesus macaques (normal or spontaneously and naturally diabetic) were examined. High-resolution ultra-performance liquid chromatography coupled with the accurate mass determination of time-of-flight mass spectrometry was used to analyze spot urine samples from normal (n = 10) and T2DM (n = 11) male monkeys. The machine-learning algorithm random forests classified urine samples as either from normal or T2DM monkeys. The metabolites important for developing the classifier were further examined for their biological significance. Random forests models had a misclassification error of less than 5%. Metabolites were identified based on accurate masses (<10 ppm) and confirmed by tandem mass spectrometry of authentic compounds. Urinary compounds significantly increased (p < 0.05) in the T2DM when compared with the normal group included glycine betaine (9-fold), citric acid (2.8-fold), kynurenic acid (1.8-fold), glucose (68-fold), and pipecolic acid (6.5-fold). When compared with the conventional definition of T2DM, the metabolites were also useful in defining the T2DM condition, and the urinary elevations in glycine betaine and pipecolic acid (as well as proline) indicated defective re-absorption in the kidney proximal tubules by SLC6A20, a Na(+)-dependent transporter. The mRNA levels of SLC6A20 were significantly reduced in the kidneys of monkeys with T2DM. These observations were validated in the db/db mouse model of T2DM. This study provides convincing evidence of the power of metabolomics for identifying functional changes at many levels in the omics pipeline.

Photographs of dental casts or digital models: rating dental arch relationships in bilateral cleft lip and palate

Yardsticks have been developed to measure dental arch relations in cleft lip and palate (CLP) patients as diagnostic proxies for the underlying skeletal relationship. Travelling with plaster casts to compare results between CLP centres is inefficient so the aim of this study was to investigate the reliability of using digital models or photographs of dental casts instead of plaster casts for rating dental arch relationships in children with complete bilateral cleft lip and palate (CBCLP). Dental casts of children with CBCLP (n=20) were included. Plaster casts, digital models and photographs of the plaster casts were available for all the children at 6, 9, and 12 years of age. All three record formats were scored using the bilateral cleft lip and palate (BCLP) yardstick by four observers in random order. No significant differences were found for the BCLP yardstick scores among the three formats. The interobserver weighted kappa scores were between 0.672 and 0.934. Comparison between the formats per observer resulted in weighted kappa scores between 0.692 and 0.885. It is concluded that digital models and photographs of dental casts can be used for rating dental arch relationships in patients with CBCLP. These formats are a reliable alternative for BCLP yardstick assessments on conventional plaster casts.

Predicting spike timing of neocortical pyramidal neurons by simple threshold models

Neurons generate spikes reliably with millisecond precision if driven by a fluctuating current--is it then possible to predict the spike timing knowing the input? We determined parameters of an adapting threshold model using data recorded in vitro from 24 layer 5 pyramidal neurons from rat somatosensory cortex, stimulated intracellularly by a fluctuating current simulating synaptic bombardment in vivo. The model generates output spikes whenever the membrane voltage (a filtered version of the input current) reaches a dynamic threshold. We find that for input currents with large fluctuation amplitude, up to 75% of the spike times can be predicted with a precision of +/-2 ms. Some of the intrinsic neuronal unreliability can be accounted for by a noisy threshold mechanism. Our results suggest that, under random current injection into the soma, (i) neuronal behavior in the subthreshold regime can be well approximated by a simple linear filter; and (ii) most of the nonlinearities are captured by a simple threshold process.

Quantification of atelectatic lung volumes in two different porcine models of ARDS

BACKGROUND: Cyclic recruitment during mechanical ventilation contributes to ventilator associated lung injury. Two different pathomechanisms in acute respiratory distress syndrome (ARDS) are currently discussed: alveolar collapse vs persistent flooding of small airways and alveoli. We compare two different ARDS animal models by computed tomography (CT) to describe different recruitment and derecruitment mechanisms at different airway pressures: (i) lavage-ARDS, favouring alveolar collapse by surfactant depletion; and (ii) oleic acid ARDS, favouring alveolar flooding by capillary leakage. METHODS: In 12 pigs [25 (1) kg], ARDS was randomly induced, either by saline lung lavage or oleic acid (OA) injection, and 3 animals served as controls. A respiratory breathhold manoeuvre without spontaneous breathing at different continuous positive airway pressure (CPAP) was applied in random order (CPAP levels of 5, 10, 15, 30, 35 and 50 cm H(2)O) and spiral-CT scans of the total lung were acquired at each CPAP level (slice thickness=1 mm). In each spiral-CT the volume of total lung parenchyma, tissue, gas, non-aerated, well-aerated, poorly aerated, and over-aerated lung was calculated. RESULTS: In both ARDS models non-aerated lung volume decreased significantly from CPAP 5 to CPAP 50 [oleic acid lung injury (OAI): 346.9 (80.1) to 96.4 (48.8) ml, P<0.001; lavage-ARDS: 245 17.6) to 42.7 (4.8) ml, P<0.001]. In lavage-ARDS poorly aerated lung volume decreased at higher CPAP levels [232 (45.2) at CPAP 10 to 84 (19.4) ml at CPAP 50, P<0.001] whereas in OAI poorly aerated lung volume did not vary at different airway pressures. CONCLUSIONS: In both ARDS models well-aerated and non-aerated lung volume respond to different CPAP levels in a comparable fashion: Thus, a cyclical alveolar collapse seems to be part of the derecruitment process also in the OA-ARDS. In OA-ARDS, the increase in poorly aerated lung volume reflects the specific initial lesion, that is capillary leakage with interstitial and alveolar oedema.