Data-driven estimates of the number of clusters in multivariate time series

An important problem in unsupervised data clustering is how to determine the number of clusters. Here we investigate how this can be achieved in an automated way by using interrelation matrices of multivariate time series. Two nonparametric and purely data driven algorithms are expounded and compared. The first exploits the eigenvalue spectra of surrogate data, while the second employs the eigenvector components of the interrelation matrix. Compared to the first algorithm, the second approach is computationally faster and not limited to linear interrelation measures.

Hierarchical modeling gave plausible estimates of associations between metabolic syndrome and components of antiretroviral therapy

OBJECTIVE: Hierarchical modeling has been proposed as a solution to the multiple exposure problem. We estimate associations between metabolic syndrome and different components of antiretroviral therapy using both conventional and hierarchical models. STUDY DESIGN AND SETTING: We use discrete time survival analysis to estimate the association between metabolic syndrome and cumulative exposure to 16 antiretrovirals from four drug classes. We fit a hierarchical model where the drug class provides a prior model of the association between metabolic syndrome and exposure to each antiretroviral. RESULTS: One thousand two hundred and eighteen patients were followed for a median of 27 months, with 242 cases of metabolic syndrome (20%) at a rate of 7.5 cases per 100 patient years. Metabolic syndrome was more likely to develop in patients exposed to stavudine, but was less likely to develop in those exposed to atazanavir. The estimate for exposure to atazanavir increased from hazard ratio of 0.06 per 6 months' use in the conventional model to 0.37 in the hierarchical model (or from 0.57 to 0.81 when using spline-based covariate adjustment). CONCLUSION: These results are consistent with trials that show the disadvantage of stavudine and advantage of atazanavir relative to other drugs in their respective classes. The hierarchical model gave more plausible results than the equivalent conventional model.

Non-ignorable loss to follow-up: correcting mortality estimates based on additional outcome ascertainment

Loss to follow-up (LTFU) is a common problem in many epidemiological studies. In antiretroviral treatment (ART) programs for patients with human immunodeficiency virus (HIV), mortality estimates can be biased if the LTFU mechanism is non-ignorable, that is, mortality differs between lost and retained patients. In this setting, routine procedures for handling missing data may lead to biased estimates. To appropriately deal with non-ignorable LTFU, explicit modeling of the missing data mechanism is needed. This can be based on additional outcome ascertainment for a sample of patients LTFU, for example, through linkage to national registries or through survey-based methods. In this paper, we demonstrate how this additional information can be used to construct estimators based on inverse probability weights (IPW) or multiple imputation. We use simulations to contrast the performance of the proposed estimators with methods widely used in HIV cohort research for dealing with missing data. The practical implications of our approach are illustrated using South African ART data, which are partially linkable to South African national vital registration data. Our results demonstrate that while IPWs and proper imputation procedures can be easily constructed from additional outcome ascertainment to obtain valid overall estimates, neglecting non-ignorable LTFU can result in substantial bias. We believe the proposed estimators are readily applicable to a growing number of studies where LTFU is appreciable, but additional outcome data are available through linkage or surveys of patients LTFU. Copyright © 2013 John Wiley & Sons, Ltd.

Use of Ultrasound Backfat Estimates to Form Marketing Groups Prior to Finishing for Feedlot Steers

An experiment was conducted using 95 Continental crossbred steers. The cattle were sorted by ultrasound 160 days before slaughter into a low backfat group (Low BF) and a higher backfat group (High BF). Half of the Low BF and half of the High BF were implanted whereas the other halves were not. Data from the experiment were used in two hypothetical markets. One market was a high yield beef program (HY) that did not allow the use of implants. The second market was a commodity beef program (CM) that allowed the use of implants. The cattle were priced as an unsorted group (ALL) and two sorted groups (Low BF and High BF) within the HY (non-implanted) and CM (implanted) markets. The CM program had a base price of $1.05/lb hot carcass weight (HCW) with a $0.15/lb HCW discount for quality grade (QG) Select and a $0.20/lb HCW discount for yield grade (YG) 4. The HY program used a base price of $1.07/lb HCW with premiums ($/lb HCW) paid for YG £ .9 (.15), 1.0 - 1.4 (.10), and 1.5 - 1.9 (.03). The carcasses were discounted ($/lb HCW) for YG 2.5 - 2.9 (.03), 3.0 - 3.9 (.15), and ³ 4.0 (.35). This data set provides good evidence that the end point at which to sell a group of cattle depends on the particular market. Sorting had an economic advantage over ALL in the HY Low BF and the CM High BF groups. The HY High BF cattle should have been sold sooner due to the discounts recieved for increased YG. The increased YG was directly affected by an increase in BF. Furthermore, the CM Low BF group should have been fed longer to increase the number of carcasses grading Choice.

Spatial heterogeneity and lake morphology affect diffusive greenhouse gas emission estimates of lakes

Most estimates of diffusive flux (F) of methane (CH4) and carbon dioxide (CO2) from lakes are based on single-point flux chamber measurements or on piston velocity (k) modeled from wind speed and single-point measurements of surface water gas concentrations (C-aq). We analyzed spatial variability of F of CH4 and CO2, as well as C-aq and k in 22 European lakes during late summer. F and k were higher in the lake centers, leading to considerable bias when extrapolating single-point chamber measurements to whole-lake estimates. The ratio of our empirical k estimates to wind speed-modeled k was related to lake size and shape, suggesting a lake morphology effect on the relationship between wind speed and k. This indicates that the error inherent to established wind speed models can be reduced by determining k and C-aq at multiple sites on lakes to calibrate wind speed-modeled k to the local system.

Recognition memory across the lifespan: the impact of word frequency and study-test interval on estimates of familiarity and recollection

The goal of this study was to investigate recognition memory performance across the lifespan and to determine how estimates of recollection and familiarity contribute to performance. In each of three experiments, participants from five groups from 14 up to 85 years of age (children, young adults, middle-aged adults, young-old adults, and old-old adults) were presented with high- and low-frequency words in a study phase and were tested immediately afterwards and/or after a one day retention interval. The results showed that word frequency and retention interval affected recognition memory performance as well as estimates of recollection and familiarity. Across the lifespan, the trajectory of recognition memory followed an inverse u-shape function that was neither affected by word frequency nor by retention interval. The trajectory of estimates of recollection also followed an inverse u-shape function, and was especially pronounced for low-frequency words. In contrast, estimates of familiarity did not differ across the lifespan. The results indicate that age differences in recognition memory are mainly due to differences in processes related to recollection while the contribution of familiarity-based processes seems to be age-invariant.

High resolution quantitative computed tomography-based assessment of trabecular microstructure and strength estimates by finite-element analysis of the spine, but not DXA, reflects vertebral fracture status in men with glucocorticoid-induced osteoporosis

High-resolution quantitative computed tomography (HRQCT)-based analysis of spinal bone density and microstructure, finite element analysis (FEA), and DXA were used to investigate the vertebral bone status of men with glucocorticoid-induced osteoporosis (GIO). DXA of L1–L3 and total hip, QCT of L1–L3, and HRQCT of T12 were available for 73 men (54.6±14.0years) with GIO. Prevalent vertebral fracture status was evaluated on radiographs using a semi-quantitative (SQ) score (normal=0 to severe fracture=3), and the spinal deformity index (SDI) score (sum of SQ scores of T4 to L4 vertebrae). Thirty-one (42.4%) subjects had prevalent vertebral fractures. Cortical BMD (Ct.BMD) and thickness (Ct.Th), trabecular BMD (Tb.BMD), apparent trabecular bone volume fraction (app.BV/TV), and apparent trabecular separation (app.Tb.Sp) were analyzed by HRQCT. Stiffness and strength of T12 were computed by HRQCT-based nonlinear FEA for axial compression, anterior bending and axial torsion. In logistic regressions adjusted for age, glucocorticoid dose and osteoporosis treatment, Tb.BMD was most closely associated with vertebral fracture status (standardized odds ratio [sOR]: Tb.BMD T12: 4.05 [95% CI: 1.8–9.0], Tb.BMD L1–L3: 3.95 [1.8–8.9]). Strength divided by cross-sectional area for axial compression showed the most significant association with spine fracture status among FEA variables (2.56 [1.29–5.07]). SDI was best predicted by a microstructural model using Ct.Th and app.Tb.Sp (r2=0.57, p<0.001). Spinal or hip DXA measurements did not show significant associations with fracture status or severity. In this cross-sectional study of males with GIO, QCT, HRQCT-based measurements and FEA variables were superior to DXA in discriminating between patients of differing prevalent vertebral fracture status. A microstructural model combining aspects of cortical and trabecular bone reflected fracture severity most accurately.

Pressure-temperature estimates of the lizardite/antigorite transition in high pressure serpentinites

Serpentine minerals in natural samples are dominated by lizardite and antigorite. In spite of numerous laboratory experiments, the stability fields of these species remain poorly constrained. This paper presents petrological observations and the Raman spectroscopy and XRD analyses of natural serpentinites from the Alpine paleo-accretionary wedge. Serpentine varieties were identified from a range of metamorphic pressure and temperature conditions from sub-greenschist (P < 4 kbar, T ~ 200–300 °C) to eclogite facies conditions (P > 20 kbar, T > 460 °C) along a subduction geothermal gradient. We use the observed mineral assemblage in natural serpentinite along with the Tmax estimated by Raman spectroscopy of the carbonaceous matter in associated metasediments to constrain the temperature of the lizardite to antigorite transition at high pressures. We show that below 300 °C, lizardite and locally chrysotile are the dominant species in the mesh texture. Between 320 and 390 °C, lizardite is progressively replaced by antigorite at the grain boundaries through dissolution–precipitation processes in the presence of SiO2 enriched fluids and in the cores of the lizardite mesh. Above 390 °C, under high-grade blueschist to eclogite facies conditions, antigorite is the sole stable serpentine mineral until the onset of secondary olivine crystallization at 460 °C.

Plotting regression coefficients and other estimates in Stata

Graphical presentation of regression results has become increasingly popular in the scientific literature, as graphs are much easier to read than tables in many cases. In Stata such plots can be produced by the -marginsplot- command. However, while -marginsplot- is very versatile and flexible, it has two major limitations: it can only process results left behind by -margins- and it can only handle one set of results at the time. In this article I introduce a new command called -coefplot- that overcomes these limitations. It plots results from any estimation command and combines results from several models into a single graph. The default behavior of -coefplot- is to plot markers for coefficients and horizontal spikes for confidence intervals. However, -coefplot- can also produce various other types of graphs. The capabilities of -coefplot- are illustrated in this article using a series of examples.