Validation of statistical shape model based reconstruction of the proximal femur - A morphology study

Seventeen bones (sixteen cadaveric bones and one plastic bone) were used to validate a method for reconstructing a surface model of the proximal femur from 2D X-ray radiographs and a statistical shape model that was constructed from thirty training surface models. Unlike previously introduced validation studies, where surface-based distance errors were used to evaluate the reconstruction accuracy, here we propose to use errors measured based on clinically relevant morphometric parameters. For this purpose, a program was developed to robustly extract those morphometric parameters from the thirty training surface models (training population), from the seventeen surface models reconstructed from X-ray radiographs, and from the seventeen ground truth surface models obtained either by a CT-scan reconstruction method or by a laser-scan reconstruction method. A statistical analysis was then performed to classify the seventeen test bones into two categories: normal cases and outliers. This classification step depends on the measured parameters of the particular test bone. In case all parameters of a test bone were covered by the training population's parameter ranges, this bone is classified as normal bone, otherwise as outlier bone. Our experimental results showed that statistically there was no significant difference between the morphometric parameters extracted from the reconstructed surface models of the normal cases and those extracted from the reconstructed surface models of the outliers. Therefore, our statistical shape model based reconstruction technique can be used to reconstruct not only the surface model of a normal bone but also that of an outlier bone.

ABCtoolbox: a versatile toolkit for approximate Bayesian computations

Background The estimation of demographic parameters from genetic data often requires the computation of likelihoods. However, the likelihood function is computationally intractable for many realistic evolutionary models, and the use of Bayesian inference has therefore been limited to very simple models. The situation changed recently with the advent of Approximate Bayesian Computation (ABC) algorithms allowing one to obtain parameter posterior distributions based on simulations not requiring likelihood computations. Results Here we present ABCtoolbox, a series of open source programs to perform Approximate Bayesian Computations (ABC). It implements various ABC algorithms including rejection sampling, MCMC without likelihood, a Particle-based sampler and ABC-GLM. ABCtoolbox is bundled with, but not limited to, a program that allows parameter inference in a population genetics context and the simultaneous use of different types of markers with different ploidy levels. In addition, ABCtoolbox can also interact with most simulation and summary statistics computation programs. The usability of the ABCtoolbox is demonstrated by inferring the evolutionary history of two evolutionary lineages of Microtus arvalis. Using nuclear microsatellites and mitochondrial sequence data in the same estimation procedure enabled us to infer sex-specific population sizes and migration rates and to find that males show smaller population sizes but much higher levels of migration than females. Conclusion ABCtoolbox allows a user to perform all the necessary steps of a full ABC analysis, from parameter sampling from prior distributions, data simulations, computation of summary statistics, estimation of posterior distributions, model choice, validation of the estimation procedure, and visualization of the results.

Effect of pulsatility on the mathematical modeling of rotary blood pumps

In this study, the effect of time derivatives of flow rate and rotational speed was investigated on the mathematical modeling of a rotary blood pump (RBP). The basic model estimates the pressure head of the pump as a dependent variable using measured flow and speed as predictive variables. Performance of the model was evaluated by adding time derivative terms for flow and speed. First, to create a realistic working condition, the Levitronix CentriMag RBP was implanted in a sheep. All parameters from the model were physically measured and digitally acquired over a wide range of conditions, including pulsatile speed. Second, a statistical analysis of the different variables (flow, speed, and their time derivatives) based on multiple regression analysis was performed to determine the significant variables for pressure head estimation. Finally, different mathematical models were used to show the effect of time derivative terms on the performance of the models. In order to evaluate how well the estimated pressure head using different models fits the measured pressure head, root mean square error and correlation coefficient were used. The results indicate that inclusion of time derivatives of flow and speed can improve model accuracy, but only minimally.

Automation of shear-wave splitting parameter determination of local earthquakes at Yellowstone : application as indicator of crustal stress and temporal variation

Shear-wave splitting can be a useful technique for determining crustal stress fields in volcanic settings and temporal variations associated with activity. Splitting parameters were determined for a subset of local earthquakes recorded from 2000-2010 at Yellowstone. Analysis was automated using an unsupervised cluster analysis technique to determine optimum splitting parameters from 270 analysis windows for each event. Six stations clearly exhibit preferential fast polarization values sub-orthogonal to the direction of minimum horizontal compression. Yellowstone deformation results in a local crustal stress field differing from the regional field dominated by NE-SW extension, and fast directions reflect this difference rotating around the caldera maintaining perpendicularity to the rim. One station exhibits temporal variations concordant with identified periods of caldera subsidence and uplift. From splitting measurements, we calculated a crustal anisotropy of ~17-23% and crack density ~0.12-0.17 possibly resulting from stress-aligned fluid filled microcracks in the upper crust and an active hydrothermal system.

Heart rate response determines long term exercise capacity after heart transplantation

BACKGROUND: Exercise capacity after heart transplantation (HTx) remains limited despite normal left ventricular systolic function of the allograft. Various clinical and haemodynamic parameters are predictive of exercise capacity following HTx. However, the predictive significance of chronotropic competence has not been demonstrated unequivocally despite its immediate relevance for cardiac output. AIMS: This study assesses the predictive value of various clinical and haemodynamic parameters for exercise capacity in HTx recipients with complete chronotropic competence evolving within the first 6 postoperative months. METHODS: 51 patients were enrolled in this exercise study. Patients were included when at least >6 months after HTx and without negative chronotropic medication or factors limiting exercise capacity such as significant transplant vasculopathy or allograft rejection. Clinical parameters were obtained by chart review, haemodynamic parameters from current cardiac catheterisation, and exercise capacity was assessed by treadmill stress testing. A stepwise multiple regression model analysed the proportion of the variance explained by the predictive parameters. RESULTS: The mean age of these 51 HTx recipients was 55.4 +/- 13.2 yrs on inclusion, 42 pts were male and the mean time interval after cardiac transplantation was 5.1 +/- 2.8 yrs. Five independent predictors explained 47.5% of the variance observed for peak exercise capacity (adjusted R2 = 0.475). In detail, heart rate response explained 31.6%, male gender 5.2%, age 4.1%, pulmonary vascular resistance 3.7%, and body-mass index 2.9%. CONCLUSION: Heart rate response is one of the most important predictors of exercise capacity in HTx recipients with complete chronotropic competence and without relevant transplant vasculopathy or acute allograft rejection.

Predicting the population impact of chlamydia screening programmes: comparative mathematical modelling study

BACKGROUND: Published individual-based, dynamic sexual network modelling studies reach different conclusions about the population impact of screening for Chlamydia trachomatis. The objective of this study was to conduct a direct comparison of the effect of organised chlamydia screening in different models. METHODS: Three models simulating population-level sexual behaviour, chlamydia transmission, screening and partner notification were used. Parameters describing a hypothetical annual opportunistic screening program in 16-24 year olds were standardised, whereas other parameters from the three original studies were retained. Model predictions of the change in chlamydia prevalence were compared under a range of scenarios. RESULTS: Initial overall chlamydia prevalence rates were similar in women but not men and there were age and sex-specific differences between models. The number of screening tests carried out was comparable in all models but there were large differences in the predicted impact of screening. After 10 years of screening, the predicted reduction in chlamydia prevalence in women aged 16-44 years ranged from 4% to 85%. Screening men and women had a greater impact than screening women alone in all models. There were marked differences between models in assumptions about treatment seeking and sexual behaviour before the start of the screening intervention. CONCLUSIONS: Future models of chlamydia transmission should be fitted to both incidence and prevalence data. This meta-modelling study provides essential information for explaining differences between published studies and increasing the utility of individual-based chlamydia transmission models for policy making.

Measurement of the vascular input function in mice for DCE-MRI

DCE-MRI is an important technique in the study of small animal cancer models because its sensitivity to vascular changes opens the possibility of quantitative assessment of early therapeutic response. However, extraction of physiologically descriptive parameters from DCE-MRI data relies upon measurement of the vascular input function (VIF), which represents the contrast agent concentration time course in the blood plasma. This is difficult in small animal models due to artifacts associated with partial volume, inflow enhancement, and the limited temporal resolution achievable with MR imaging. In this work, the development of a suite of techniques for high temporal resolution, artifact resistant measurement of the VIF in mice is described. One obstacle in VIF measurement is inflow enhancement, which decreases the sensitivity of the MR signal to the presence of contrast agent. Because the traditional techniques used to suppress inflow enhancement degrade the achievable spatiotemporal resolution of the pulse sequence, improvements can be achieved by reducing the time required for the suppression. Thus, a novel RF pulse which provides spatial presaturation contemporaneously with the RF excitation was implemented and evaluated. This maximizes the achievable temporal resolution by removing the additional RF and gradient pulses typically required for suppression of inflow enhancement. A second challenge is achieving the temporal resolution required for accurate characterization of the VIF, which exceeds what can be achieved with conventional imaging techniques while maintaining adequate spatial resolution and tumor coverage. Thus, an anatomically constrained reconstruction strategy was developed that allows for sampling of the VIF at extremely high acceleration factors, permitting capture of the initial pass of the contrast agent in mice. Simulation, phantom, and in vivo validation of all components were performed. Finally, the two components were used to perform VIF measurement in the murine heart. An in vivo study of the VIF reproducibility was performed, and an improvement in the measured injection-to-injection variation was observed. This will lead to improvements in the reliability of quantitative DCE-MRI measurements and increase their sensitivity.

Hepatocellular carcinoma survival in uninsured and underinsured patients.

BACKGROUND: The incidence of hepatitis C virus (HCV) and hepatocellular carcinoma (HCC) is increasing. The purpose of this study is to establish baseline survival in a medically-underserved population and to evaluate the effect of HCV seropositivity on our patient population. MATERIALS AND METHODS: We reviewed clinicopathologic parameters from a prospective tumor registry and medical records from the Harris County Hospital District (HCHD). Outcomes were compared using Kaplan-Meier survival analysis and log-rank tests. RESULTS: A total of 298 HCC patients were identified. The median survival for the entire cohort was 3.4 mo. There was no difference in survival between the HCV seropositive and the HCV seronegative groups (3.6 mo versus 2.6 mo, P = 0.7). Patients with a survival <1 mo had a significant increase in>αfetoprotein (AFP), international normalized ratio (INR), model for end-stage liver disease (MELD) score, and total bilirubin and decrease in albumin compared with patients with a survival ≥ 1 mo. CONCLUSIONS: Survival for HCC patients in the HCHD is extremely poor compared with an anticipated median survival of 7 mo reported in other studies. HCV seropositive patients have no survival advantage over HCV seronegative patients. Poorer liver function at diagnosis appears to be related to shorter survival. Further analysis into variables contributing to decreased survival is needed.

Patient classification system: An integrated method for measuring nursing intensity and optimizing resource allocation

A patient classification system was developed integrating a patient acuity instrument with a computerized nursing distribution method based on a linear programming model. The system was designed for real-time measurement of patient acuity (workload) and allocation of nursing personnel to optimize the utilization of resources.^ The acuity instrument was a prototype tool with eight categories of patients defined by patient severity and nursing intensity parameters. From this tool, the demand for nursing care was defined in patient points with one point equal to one hour of RN time. Validity and reliability of the instrument was determined as follows: (1) Content validity by a panel of expert nurses; (2) predictive validity through a paired t-test analysis of preshift and postshift categorization of patients; (3) initial reliability by a one month pilot of the instrument in a practice setting; and (4) interrater reliability by the Kappa statistic.^ The nursing distribution system was a linear programming model using a branch and bound technique for obtaining integer solutions. The objective function was to minimize the total number of nursing personnel used by optimally assigning the staff to meet the acuity needs of the units. A penalty weight was used as a coefficient of the objective function variables to define priorities for allocation of staff.^ The demand constraints were requirements to meet the total acuity points needed for each unit and to have a minimum number of RNs on each unit. Supply constraints were: (1) total availability of each type of staff and the value of that staff member (value was determined relative to that type of staff's ability to perform the job function of an RN (i.e., value for eight hours RN = 8 points, LVN = 6 points); (2) number of personnel available for floating between units.^ The capability of the model to assign staff quantitatively and qualitatively equal to the manual method was established by a thirty day comparison. Sensitivity testing demonstrated appropriate adjustment of the optimal solution to changes in penalty coefficients in the objective function and to acuity totals in the demand constraints.^ Further investigation of the model documented: correct adjustment of assignments in response to staff value changes; and cost minimization by an addition of a dollar coefficient to the objective function. ^

Estimating the mean and standard deviation of concentration distributions using censored samples

Environmental data sets of pollutant concentrations in air, water, and soil frequently include unquantified sample values reported only as being below the analytical method detection limit. These values, referred to as censored values, should be considered in the estimation of distribution parameters as each represents some value of pollutant concentration between zero and the detection limit. Most of the currently accepted methods for estimating the population parameters of environmental data sets containing censored values rely upon the assumption of an underlying normal (or transformed normal) distribution. This assumption can result in unacceptable levels of error in parameter estimation due to the unbounded left tail of the normal distribution. With the beta distribution, which is bounded by the same range of a distribution of concentrations, $\rm\lbrack0\le x\le1\rbrack,$ parameter estimation errors resulting from improper distribution bounds are avoided. This work developed a method that uses the beta distribution to estimate population parameters from censored environmental data sets and evaluated its performance in comparison to currently accepted methods that rely upon an underlying normal (or transformed normal) distribution. Data sets were generated assuming typical values encountered in environmental pollutant evaluation for mean, standard deviation, and number of variates. For each set of model values, data sets were generated assuming that the data was distributed either normally, lognormally, or according to a beta distribution. For varying levels of censoring, two established methods of parameter estimation, regression on normal ordered statistics, and regression on lognormal ordered statistics, were used to estimate the known mean and standard deviation of each data set. The method developed for this study, employing a beta distribution assumption, was also used to estimate parameters and the relative accuracy of all three methods were compared. For data sets of all three distribution types, and for censoring levels up to 50%, the performance of the new method equaled, if not exceeded, the performance of the two established methods. Because of its robustness in parameter estimation regardless of distribution type or censoring level, the method employing the beta distribution should be considered for full development in estimating parameters for censored environmental data sets. ^

Diffusion tensor imaging of rat spinal cord in vivo

Magnetic resonance imaging, with its exquisite soft tissue contrast, is an ideal modality for investigating spinal cord pathology. While conventional MRI techniques are very sensitive for spinal cord pathology, their specificity is somewhat limited. Diffusion MRI is an advanced technique which is a very sensitive and specific indicator of the integrity of white matter tracts. Diffusion imaging has been shown to detect early ischemic changes in white matter, while conventional imaging demonstrates no change. By acquiring the complete apparent diffusion tensor (ADT), tissue diffusion properties can be expressed in terms of quantitative and rotationally invariant parameters. ^ Systematic study of SCI in vivo requires controlled animal models such as the popular rat model. To date, studies of spinal cord using ADT imaging have been performed exclusively in fixed, excised spinal cords, introducing inevitable artifacts and losing the benefits of MRI's noninvasive nature. In vivo imaging reflects the actual in vivo tissue properties, and allows each animal to be imaged at multiple time points, greatly reducing the number of animals required to achieve statistical significance. Because the spinal cord is very small, the available signal-to-noise ratio (SNR) is very low. Prior spin-echo based ADT studies of rat spinal cord have relied on high magnetic field strengths and long imaging times—on the order of 10 hours—for adequate SNR. Such long imaging times are incompatible with in vivo imaging, and are not relevant for imaging the early phases following SCI. Echo planar imaging (EPI) is one of the fastest imaging methods, and is popular for diffusion imaging. However, EPI further lowers the image SNR, and is very sensitive to small imperfections in the magnetic field, such as those introduced by the bony spine. Additionally, The small field-of-view (FOV) needed for spinal cord imaging requires large imaging gradients which generate EPI artifacts. The addition of diffusion gradients introduces yet further artifacts. ^ This work develops a method for rapid EPI-based in vivo diffusion imaging of rat spinal cord. The method involves improving the SNR using an implantable coil; reducing magnetic field inhomogeneities by means of an autoshim, and correcting EPI artifacts by post-processing. New EPI artifacts due to diffusion gradients described, and post-processing correction techniques are developed. ^ These techniques were used to obtain rotationally invariant diffusion parameters from 9 animals in vivo, and were validated using the gold-standard, but slow, spinecho based diffusion sequence. These are the first reported measurements of the ADT in spinal cord in vivo . ^ Many of the techniques described are equally applicable toward imaging of human spinal cord. We anticipate that these techniques will aid in evaluating and optimizing potential therapies, and will lead to improved patient care. ^

Homogeneous reprocessing of GPS, GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from −35 and −38 mm to −12 and −13 mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.

Finite element analysis in asymptomatic, symptomatic, and ruptured abdominal aortic aneurysms: in search of new rupture risk predictors.

OBJECTIVES To compare biomechanical rupture risk parameters of asymptomatic, symptomatic and ruptured abdominal aortic aneurysms (AAA) using finite element analysis (FEA). STUDY DESIGN Retrospective biomechanical single center analysis of asymptomatic, symptomatic, and ruptured AAAs. Comparison of biomechanical parameters from FEA. MATERIALS AND METHODS From 2011 to 2013 computed tomography angiography (CTA) data from 30 asymptomatic, 15 symptomatic, and 15 ruptured AAAs were collected consecutively. FEA was performed according to the successive steps of AAA vessel reconstruction, segmentation and finite element computation. Biomechanical parameters Peak Wall Rupture Risk Index (PWRI), Peak Wall Stress (PWS), and Rupture Risk Equivalent Diameter (RRED) were compared among the three subgroups. RESULTS PWRI differentiated between asymptomatic and symptomatic AAAs (p < .0004) better than PWS (p < .1453). PWRI-dependent RRED was higher in the symptomatic subgroup compared with the asymptomatic subgroup (p < .0004). Maximum AAA external diameters were comparable between the two groups (p < .1355). Ruptured AAAs showed the highest values for external diameter, total intraluminal thrombus volume, PWS, RRED, and PWRI compared with asymptomatic and symptomatic AAAs. In contrast with symptomatic and ruptured AAAs, none of the asymptomatic patients had a PWRI value >1.0. This threshold value might identify patients at imminent risk of rupture. CONCLUSIONS From different FEA derived parameters, PWRI distinguishes most precisely between asymptomatic and symptomatic AAAs. If elevated, this value may represent a negative prognostic factor for asymptomatic AAAs.

Real-Time In Vivo Characterization of Primary Liver Tumors With Diffuse Optical Spectroscopy During Percutaneous Needle Interventions: Feasibility Study in Woodchucks

OBJECTIVE This study presents the first in vivo real-time optical tissue characterization during image-guided percutaneous intervention using near-infrared diffuse optical spectroscopy sensing at the tip of a needle. The goal of this study was to indicate transition boundaries from healthy tissue to tumors, namely, hepatic carcinoma, based on the real-time feedback derived from the optical measurements. MATERIALS AND METHODS Five woodchucks with hepatic carcinoma were used for this study. The woodchucks were imaged with contrast-enhanced cone beam computed tomography with a flat panel detector C-arm system to visualize the carcinoma in the liver. In each animal, 3 insertions were performed, starting from the skin surface toward the hepatic carcinoma under image guidance. In 2 woodchucks, each end point of the insertion was confirmed with pathologic examination of a biopsy sample. While advancing the needle in the animals under image guidance such as fluoroscopy overlaid with cone beam computed tomography slice and ultrasound, optical spectra were acquired at the distal end of the needles. Optical tissue characterization was determined by translating the acquired optical spectra into clinical parameters such as blood, water, lipid, and bile fractions; tissue oxygenation levels; and scattering amplitude related to tissue density. The Kruskal-Wallis test was used to study the difference in the derived clinical parameters from the measurements performed within the healthy tissue and the hepatic carcinoma. Kurtoses were calculated to assess the dispersion of these parameters within the healthy and carcinoma tissues. RESULTS Blood and lipid volume fractions as well as tissue oxygenation and reduced scattering amplitude showed to be significantly different between the healthy part of the liver and the hepatic carcinoma (P < 0.05) being higher in normal liver tissue. A decrease in blood and lipid volume fractions and tissue oxygenation as well as an increase in scattering amplitude were observed when the tip of the needle crossed the margin from the healthy liver tissue to the carcinoma. The kurtosis for each derived clinical parameter was high in the hepatic tumor as compared with that in the healthy liver indicating intracarcinoma variability. CONCLUSIONS Tissue blood content, oxygenation level, lipid content, and tissue density all showed significant differences when the needle tip was guided from the healthy tissue to the carcinoma and can therefore be used to identify tissue boundaries during percutaneous image-guided interventions.

Submaximal cardiopulmonary thresholds on a robotics-assisted tilt table, a cycle and a treadmill: a comparative analysis

BACKGROUND: The robotics-assisted tilt table (RATT), including actuators for tilting and cyclical leg movement, is used for rehabilitation of severely disabled neurological patients. Following further engineering development of the system, i.e. the addition of force sensors and visual bio-feedback, patients can actively participate in exercise testing and training on the device. Peak cardiopulmonary performance parameters were previously investigated, but it also important to compare submaximal parameters with standard devices. The aim of this study was to evaluate the feasibility of the RATT for estimation of submaximal exercise thresholds by comparison with a cycle ergometer and a treadmill. METHODS: 17 healthy subjects randomly performed six maximal individualized incremental exercise tests, with two tests on each of the three exercise modalities. The ventilatory anaerobic threshold (VAT) and respiratory compensation point (RCP) were determined from breath-by-breath data. RESULTS: VAT and RCP on the RATT were lower than the cycle ergometer and the treadmill: oxygen uptake (V'O2) at VAT was [mean (SD)] 1.2 (0.3), 1.5 (0.4) and 1.6 (0.5) L/min, respectively (p < 0.001); V'O2 at RCP was 1.7 (0.4), 2.3 (0.8) and 2.6 (0.9) L/min, respectively (p = 0.001). High correlations for VAT and RCP were found between the RATT vs the cycle ergometer and RATT vs the treadmill (R on the range 0.69-0.80). VAT and RCP demonstrated excellent test-retest reliability for all three devices (ICC from 0.81 to 0.98). Mean differences between the test and retest values on each device were close to zero. The ventilatory equivalent for O2 at VAT for the RATT and cycle ergometer were similar and both were higher than the treadmill. The ventilatory equivalent for CO2 at RCP was similar for all devices. Ventilatory equivalent parameters demonstrated fair-to-excellent reliability and repeatability. CONCLUSIONS: It is feasible to use the RATT for estimation of submaximal exercise thresholds: VAT and RCP on the RATT were lower than the cycle ergometer and the treadmill, but there were high correlations between the RATT vs the cycle ergometer and vs the treadmill. Repeatability and test-retest reliability of all submaximal threshold parameters from the RATT were comparable to those of standard devices.