Semiparametric Latent Variable Regression Models for Spatio-temporal Modeling of Mobile Source Particles in the Greater Boston Area

Traffic particle concentrations show considerable spatial variability within a metropolitan area. We consider latent variable semiparametric regression models for modeling the spatial and temporal variability of black carbon and elemental carbon concentrations in the greater Boston area. Measurements of these pollutants, which are markers of traffic particles, were obtained from several individual exposure studies conducted at specific household locations as well as 15 ambient monitoring sites in the city. The models allow for both flexible, nonlinear effects of covariates and for unexplained spatial and temporal variability in exposure. In addition, the different individual exposure studies recorded different surrogates of traffic particles, with some recording only outdoor concentrations of black or elemental carbon, some recording indoor concentrations of black carbon, and others recording both indoor and outdoor concentrations of black carbon. A joint model for outdoor and indoor exposure that specifies a spatially varying latent variable provides greater spatial coverage in the area of interest. We propose a penalised spline formation of the model that relates to generalised kringing of the latent traffic pollution variable and leads to a natural Bayesian Markov Chain Monte Carlo algorithm for model fitting. We propose methods that allow us to control the degress of freedom of the smoother in a Bayesian framework. Finally, we present results from an analysis that applies the model to data from summer and winter separately

Analyzing Panel Count Data with Informative Observation Times

In this paper, we study panel count data with informative observation times. We assume nonparametric and semiparametric proportional rate models for the underlying recurrent event process, where the form of the baseline rate function is left unspecified and a subject-specific frailty variable inflates or deflates the rate function multiplicatively. The proposed models allow the recurrent event processes and observation times to be correlated through their connections with the unobserved frailty; moreover, the distributions of both the frailty variable and observation times are considered as nuisance parameters. The baseline rate function and the regression parameters are estimated by maximizing a conditional likelihood function of observed event counts and solving estimation equations. Large sample properties of the proposed estimators are studied. Numerical studies demonstrate that the proposed estimation procedures perform well for moderate sample sizes. An application to a bladder tumor study is presented to illustrate the use of the proposed methods.

Checking Assumptions in Latent Class Regression Models via a Markov Chain Monte Carlo Estimation Approach: An Application to Depression and Socio-Economic Status

Latent class regression models are useful tools for assessing associations between covariates and latent variables. However, evaluation of key model assumptions cannot be performed using methods from standard regression models due to the unobserved nature of latent outcome variables. This paper presents graphical diagnostic tools to evaluate whether or not latent class regression models adhere to standard assumptions of the model: conditional independence and non-differential measurement. An integral part of these methods is the use of a Markov Chain Monte Carlo estimation procedure. Unlike standard maximum likelihood implementations for latent class regression model estimation, the MCMC approach allows us to calculate posterior distributions and point estimates of any functions of parameters. It is this convenience that allows us to provide the diagnostic methods that we introduce. As a motivating example we present an analysis focusing on the association between depression and socioeconomic status, using data from the Epidemiologic Catchment Area study. We consider a latent class regression analysis investigating the association between depression and socioeconomic status measures, where the latent variable depression is regressed on education and income indicators, in addition to age, gender, and marital status variables. While the fitted latent class regression model yields interesting results, the model parameters are found to be invalid due to the violation of model assumptions. The violation of these assumptions is clearly identified by the presented diagnostic plots. These methods can be applied to standard latent class and latent class regression models, and the general principle can be extended to evaluate model assumptions in other types of models.

MODIFIED TEST STATISTICS BY INTER-VOXEL VARIANCE SHRINKAGE WITH AN APPLICATION TO fMRI

Functional Magnetic Resonance Imaging (fMRI) is a non-invasive technique which is commonly used to quantify changes in blood oxygenation and flow coupled to neuronal activation. One of the primary goals of fMRI studies is to identify localized brain regions where neuronal activation levels vary between groups. Single voxel t-tests have been commonly used to determine whether activation related to the protocol differs across groups. Due to the generally limited number of subjects within each study, accurate estimation of variance at each voxel is difficult. Thus, combining information across voxels in the statistical analysis of fMRI data is desirable in order to improve efficiency. Here we construct a hierarchical model and apply an Empirical Bayes framework on the analysis of group fMRI data, employing techniques used in high throughput genomic studies. The key idea is to shrink residual variances by combining information across voxels, and subsequently to construct an improved test statistic in lieu of the classical t-statistic. This hierarchical model results in a shrinkage of voxel-wise residual sample variances towards a common value. The shrunken estimator for voxelspecific variance components on the group analyses outperforms the classical residual error estimator in terms of mean squared error. Moreover, the shrunken test-statistic decreases false positive rate when testing differences in brain contrast maps across a wide range of simulation studies. This methodology was also applied to experimental data regarding a cognitive activation task.

AN OVERVIEW OF OBSERVATIONAL SLEEP RESEARCH WITH APPLICATION TO SLEEP STAGE TRANSITIONING

In this manuscript we give an overview of observational sleep research with a particular emphasis on sleep stage transitions. Sleep states represent a categorization of sleep electroencephalogram behavior over the night. We postulate that the rate of transitioning between sleep states is an important predictor of health. This claim is evaluated by comparing subjects with sleep disordered breathing to matched controls.

Kernel Estimation of Rate Function for Recurrent Event Data

Recurrent event data are largely characterized by the rate function but smoothing techniques for estimating the rate function have never been rigorously developed or studied in statistical literature. This paper considers the moment and least squares methods for estimating the rate function from recurrent event data. With an independent censoring assumption on the recurrent event process, we study statistical properties of the proposed estimators and propose bootstrap procedures for the bandwidth selection and for the approximation of confidence intervals in the estimation of the occurrence rate function. It is identified that the moment method without resmoothing via a smaller bandwidth will produce curve with nicks occurring at the censoring times, whereas there is no such problem with the least squares method. Furthermore, the asymptotic variance of the least squares estimator is shown to be smaller under regularity conditions. However, in the implementation of the bootstrap procedures, the moment method is computationally more efficient than the least squares method because the former approach uses condensed bootstrap data. The performance of the proposed procedures is studied through Monte Carlo simulations and an epidemiological example on intravenous drug users.

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.

Anaesthesia for castration of piglets: comparison between intranasal and intramuscular application of ketamine, climazolam and azaperone

The aim of this study was to compare the effect of an anaesthetic combination given either intramuscularly (IM) or intranasally (IN) for castration of piglets. Forty piglets aged 4 to 7 days were randomly assigned to receive a mixture of ketamine 15 mg kg-1, climazolam 1.5 mg kg-1 and azaperone 1.0 mg kg-1, IN or IM, 10 minutes prior to castration. Physiological parameters were measured. Castration was videotaped for evaluation by 3 independent observers using a scoring system. Reaction and vocalization to the skin incision and cutting of spermatic cord was evaluated and scored (0 = no reaction, 16 = strong reaction). The IN group had a significantly higher (P < 0.01) castration score, compared to the IM group. There was an association between castration score and room temperature in the IN group (with temperatures below 18 "C associated with a higher castration scores (P < 0.001). Heart rate was significantly higher 10 minutes after castration in the IN group (P < 0.05). Respiratory rate was significantly higher in the IM group at time points -5, -1, 10, 20 and 30 (P < 0.05).The IN group was walking significantly (P < 0.0001) faster than the IM group. In conclusion, this combination provides effective anaesthesia for routine castration of newborn piglets when administered IM. IN administration provided shorter recovery times but had significantly higher castration scores.

Effects of non-ionic iodinated contrast media on patient heart rate and pressures during intra-cardiac or intra-arterial injection

PURPOSE: To compare the effects on heart rate (HR), on left ventricular (LV) or arterial pressures, and the general safety of a non-ionic low-osmolar contrast medium (CM) and a non-ionic iso-osmolar CM in patients undergoing cardiac angiography (CA) or peripheral intra-arterial digital subtraction angiography (IA-DSA). MATERIALS AND METHODS: Two double-blind, randomized studies were conducted in 216 patients who underwent CA (n=120) or peripheral IA-DSA (n=96). Patients referred for CA received a low-osmolar monomeric CM (iomeprol-350, n=60) or an iso-osmolar dimeric CM (iodixanol-320; n=60). HR and LV peak systolic and end-diastolic pressures were determined before and after the first injection during left and right coronary arteriography and left ventriculography. Monitoring for all types of adverse event (AE) was performed for 24 h following the procedure. t-tests were performed to compare CM for effects on HR. Patients referred for IA-DSA received iomeprol-300 (n=49) or iodixanol-320 (n=47). HR and arterial blood pressure (BP) were evaluated before and after the first 4 injections. Monitoring for AE was performed for 4 h following the procedure. Repeated-measures ANOVA was used to compare mean HR changes across the first 4 injections, whereas changes after the first injection were compared using t-tests. RESULTS: No significant differences were noted between iomeprol and iodixanol in terms of mean changes in HR during left coronary arteriography (p=0.8), right coronary arteriography (p=0.9), and left ventriculography (p=0.8). In patients undergoing IA-DSA, no differences between CM were noted for effects on mean HR after the first injection (p=0.6) or across the first 4 injections (p=0.2). No significant differences (p>0.05) were noted in terms of effects on arterial BP in either study or on LV pressures in patients undergoing CA. Non-serious AE considered possibly CM-related (primarily headache and events affecting the cardiovascular and digestive systems) were reported more frequently by patients undergoing CA and more frequently after iodixanol (14/60 [23.3%] and 2/47 [4.3%]; CA and IA-DSA, respectively) than iomeprol (10/60 [16.7%] and 1/49 [2%], respectively). CONCLUSIONS: Iomeprol and iodixanol are safe and have equally negligible effects on HR and LV pressures or arterial BP during and after selective intra-cardiac injection and peripheral IA-DSA. CLINICAL APPLICATION: Iomeprol and iodixanol are safe and equally well tolerated with regard to cardiac rhythm and clinical preference should be based on diagnostic image quality alone.

Precise pose recovery of distal locking holes from single calibrated fluoroscopic image via a novel variable decomposition approach

This paper presents a novel variable decomposition approach for pose recovery of the distal locking holes using single calibrated fluoroscopic image. The problem is formulated as a model-based optimal fitting process, where the control variables are decomposed into two sets: (a) the angle between the nail axis and its projection on the imaging plane, and (b) the translation and rotation of the geometrical model of the distal locking hole around the nail axis. By using an iterative algorithm to find the optimal values of the latter set of variables for any given value of the former variable, we reduce the multiple-dimensional model-based optimal fitting problem to a one-dimensional search along a finite interval. We report the results of our in vitro experiments, which demonstrate that the accuracy of our approach is adequate for successful distal locking of intramedullary nails.

Indocyanine green plasma disappearance rate during the anhepatic phase of orthotopic liver transplantation

Non-invasive pulse spectrophotometry to measure indocyanine green (ICG) elimination correlates well with the conventional invasive ICG clearance test. Nevertheless, the precision of this method remains unclear for any application, including small-for-size liver remnants. We therefore measured ICG plasma disappearance rate (PDR) during the anhepatic phase of orthotopic liver transplantation using pulse spectrophotometry. Measurements were done in 24 patients. The median PDR after exclusion of two outliers and two patients with inconstant signal was 1.55%/min (95% confidence interval [CI]=0.8-2.2). No correlation with patient age, gender, body mass, blood loss, administration of fresh frozen plasma, norepinephrine dose, postoperative albumin (serum), or difference in pre and post transplant body weight was detected. In conclusion, we found an ICG-PDR different from zero in the anhepatic phase, an overestimation that may arise in particular from a redistribution into the interstitial space. If ICG pulse spectrophotometry is used to measure functional hepatic reserve, the verified average difference from zero (1.55%/min) determined in our study needs to be taken into account.

Occlusion length is a crucial determinant of efficiency and complication rate in thrombectomy for acute ischemic stroke

BACKGROUND AND PURPOSE: Although mechanical thrombectomy (MT) has an encouragingly high recanalization rate in treating stroke, it is associated with severe complications of which the underlying factors have yet to be identified. Because MT is a mechanical approach, the mechanical properties of the thrombus might be crucial for its success. The present study assesses the effect of thrombus length on the in vivo effectiveness and complication rate of MT. MATERIALS AND METHODS: Angiography and embolization of 21 cranial vessels with radiopaque whole-blood thrombi 10, 20, and 40 mm in length (7 occlusions each) were performed in 7 swine. MT was carried out using a distal snarelike device (BCR Roadsaver) with proximal balloon occlusion. A total of 61 retrievals were attempted. RESULTS: In the group of 10-mm occlusions, 77.8% of the attempts achieved complete recanalisation. For longer occlusions, the success rates decreased significantly to 20% of attempts for 20-mm occlusions (odds ratio [OR], 14; 95% confidence interval [CI], 2.2-89.2) and 11.1% for 40-mm occlusions (OR, 28; 95% CI, 3.9-202.2; P < .005). The low success rates were largely due to complications associated with thrombus compaction during retrieval. Similarly, the rate of thromboembolic events increased from 0% in 10-mm occlusions to 14.8% in 40-mm occlusions. CONCLUSIONS: MT using a distal device proved to be a fast, effective, and safe procedure for recanalizing short (10-mm) occlusions in the animal model. However, occlusion length emerged as a crucial determinant for MT with a significant decrease in recanalization success per attempt and increased complication rates. These findings suggest limitations of MT in the clinical application.

Effects of pediatric adiposity on heart rate variability

The relationship between obesity and heart rate variability (HRV) has been studied in adults and adolescents, but is not determined in young pediatrics. The purpose of this study was to assess autonomic activity using HRV in a pediatric population. We hypothesized that obese children would have reduced parasympathetic and increased sympathetic activity compared to age-matched subjects. 42 pediatric subjects (ages 3-5) were classified into 3 groups based on body mass index-for-age; normal, overweight and obese. HRV and respiratory rate were recorded during 3 minute baseline, 2 minute isometric handgrip and 3 minute recovery. HRV was analyzed in the time domain [heart rate (HR), RR interval (RRI) and RRI standard deviation (RRISD)] and frequency domain [low frequency (LF), high frequency (HF) and LF/HF ratio] using repeated measures ANOVA. Spearman’s correlations were used to examine the relations between BMI and HRV at rest. Significant condition effects were found between baseline, exercise and recovery, but these responses were not significantly different between the normal, overweight and obese children. BMI was negatively correlated with LF/HF, while BMI was positively correlated with RRISD, LF, HF and nHF. Our data demonstrate that higher BMI in the pediatric population is correlated with higher parasympathetic and lower sympathetic activity. These findings are contrary to HRV responses observed in adults and adolescents, suggesting complex relationships between age, obesity and autonomic control of the heart. The data supports the concept of an age reliance of HRV and a novel relationship between adiposity and body mass index in 3-5 year olds.

Effect of constant and variable domain glycosylation on pharmacokinetics of therapeutic antibodies in mice

Previous studies on the effect of glycosylation on the elimination rate of antibodies have produced conflicting results. Here, we performed pharmacokinetic studies in mice with two preparations of a monoclonal IgG1 antibody enriched for complex type or high mannose type oligosaccharides at the Fc glycosylation site. No significant difference in the serum half-life was found between the two antibody glycoforms, nor was any difference observed in the serum half-lives of different complex type glycoforms. To evaluate the influence of glycosylation within the variable domain, a second monoclonal antibody, glycosylated in both the Fc and Fv domains, was separated into fractions containing different amounts of Fv-associated sialic acid and administered to mice. Again, no significant difference was found in the clearance rates of variants carrying different amounts of Fv-associated sialic acid or lacking Fv-glycosylation. These results suggest that glycosylation has little or no impact on the pharmacokinetic behavior of these two monoclonal antibodies in mice.