A real time simulation model of glucose-insulin metabolism for type 1 diabetes patients

In this paper, a simulation model of glucose-insulin metabolism for Type 1 diabetes patients is presented. The proposed system is based on the combination of Compartmental Models (CMs) and artificial Neural Networks (NNs). This model aims at the development of an accurate system, in order to assist Type 1 diabetes patients to handle their blood glucose profile and recognize dangerous metabolic states. Data from a Type 1 diabetes patient, stored in a database, have been used as input to the hybrid system. The data contain information about measured blood glucose levels, insulin intake, and description of food intake, along with the corresponding time. The data are passed to three separate CMs, which produce estimations about (i) the effect of Short Acting (SA) insulin intake on blood insulin concentration, (ii) the effect of Intermediate Acting (IA) insulin intake on blood insulin concentration, and (iii) the effect of carbohydrate intake on blood glucose absorption from the gut. The outputs of the three CMs are passed to a Recurrent NN (RNN) in order to predict subsequent blood glucose levels. The RNN is trained with the Real Time Recurrent Learning (RTRL) algorithm. The resulted blood glucose predictions are promising for the use of the proposed model for blood glucose level estimation for Type 1 diabetes patients.

Non-invasive quantitative estimation of bone density in rats throughout the life cycle and in arthritic osteopenia: preliminary results

A longitudinal bone survey was conducted in 86 female Wistar rats in order to assess mineral density kinetics from young age (5 weeks: 115 g) till late adulthood (64 weeks: 586 g). In vivo quantitative radiographic scanning was performed on the caudal vertebrae, taking trabecular mass as the parameter. Measurements were expressed as Relative Optical Density (ROD) units by means of a high resolution densitometric device. Results showed a progressive increase in mineral density throughout the life cycle, with a tendency to level in the higher weight range, indicating that progressive mineral aposition occurs in rats in dependency of age. This phenomenon, however, should be always considered within the context of continuous skeletal growth and related changes typical of this species. Twelve different animals were also examined following induction of articular inflammation with Freund's adjuvant in six of them. Bone survey conducted 12 to 18 days after inoculation revealed a significant (P less than 0.01) reduction in trabecular bone mass of scanned vertebrae in comparison with the weight-matched untreated controls. It is concluded that the in vivo quantitative assessment of bone density illustrated in this report represents a sensitive and useful tool for the long-term survey of naturally occurring or experimentally induced bone changes. Scanning of the same part of the skeleton can be repeated, thereby avoiding sacrifice of the animal and time-consuming preparation of post-mortem material.

Incorporating Latent Variables into Discrete Choice Models - A Simultaneous Estimation Approach Using SEM Software

Integrated choice and latent variable (ICLV) models represent a promising new class of models which merge classic choice models with the structural equation approach (SEM) for latent variables. Despite their conceptual appeal, applications of ICLV models in marketing remain rare. We extend previous ICLV applications by first estimating a multinomial choice model and, second, by estimating hierarchical relations between latent variables. An empirical study on travel mode choice clearly demonstrates the value of ICLV models to enhance the understanding of choice processes. In addition to the usually studied directly observable variables such as travel time, we show how abstract motivations such as power and hedonism as well as attitudes such as a desire for flexibility impact on travel mode choice. Furthermore, we show that it is possible to estimate such a complex ICLV model with the widely available structural equation modeling package Mplus. This finding is likely to encourage more widespread application of this appealing model class in the marketing field.

Prediction of PA12 melt viscosity in Laser Sintering by a Time and Temperature dependent rheological model

Ein auf Basis von Prozessdaten kalibriertes Viskositätsmodell wird vorgeschlagen und zur Vorhersage der Viskosität einer Polyamid 12 (PA12) Kunststoffschmelze als Funktion von Zeit, Temperatur und Schergeschwindigkeit angewandt. Im ersten Schritt wurde das Viskositätsmodell aus experimentellen Daten abgeleitet. Es beruht hauptsächlich auf dem drei-parametrigen Ansatz von Carreau, wobei zwei zusätzliche Verschiebungsfaktoren eingesetzt werden. Die Temperaturabhängigkeit der Viskosität wird mithilfe des Verschiebungsfaktors aT von Arrhenius berücksichtigt. Ein weiterer Verschiebungsfaktor aSC (Structural Change) wird eingeführt, der die Strukturänderung von PA12 als Folge der Prozessbedingungen beim Lasersintern beschreibt. Beobachtet wurde die Strukturänderung in Form einer signifikanten Viskositätserhöhung. Es wurde geschlussfolgert, dass diese Viskositätserhöhung auf einen Molmassenaufbau zurückzuführen ist und als Nachkondensation verstanden werden kann. Abhängig von den Zeit- und Temperaturbedingungen wurde festgestellt, dass die Viskosität als Folge des Molmassenaufbaus exponentiell gegen eine irreversible Grenze strebt. Die Geschwindigkeit dieser Nachkondensation ist zeit- und temperaturabhängig. Es wird angenommen, dass die Pulverbetttemperatur einen Molmassenaufbau verursacht und es damit zur Kettenverlängerung kommt. Dieser fortschreitende Prozess der zunehmenden Kettenlängen setzt molekulare Beweglichkeit herab und unterbindet die weitere Nachkondensation. Der Verschiebungsfaktor aSC drückt diese physikalisch-chemische Modellvorstellung aus und beinhaltet zwei zusätzliche Parameter. Der Parameter aSC,UL entspricht der oberen Viskositätsgrenze, wohingegen k0 die Strukturänderungsrate angibt. Es wurde weiterhin festgestellt, dass es folglich nützlich ist zwischen einer Fließaktivierungsenergie und einer Strukturänderungsaktivierungsenergie für die Berechnung von aT und aSC zu unterscheiden. Die Optimierung der Modellparameter erfolgte mithilfe eines genetischen Algorithmus. Zwischen berechneten und gemessenen Viskositäten wurde eine gute Übereinstimmung gefunden, so dass das Viskositätsmodell in der Lage ist die Viskosität einer PA12 Kunststoffschmelze als Folge eines kombinierten Lasersinter Zeit- und Temperatureinflusses vorherzusagen. Das Modell wurde im zweiten Schritt angewandt, um die Viskosität während des Lasersinter-Prozesses in Abhängigkeit von der Energiedichte zu berechnen. Hierzu wurden Prozessdaten, wie Schmelzetemperatur und Belichtungszeit benutzt, die mithilfe einer High-Speed Thermografiekamera on-line gemessen wurden. Abschließend wurde der Einfluss der Strukturänderung auf das Viskositätsniveau im Prozess aufgezeigt.

Climate spectrum estimation in the presence of timescale errors

We introduce an algorithm (called REDFITmc2) for spectrum estimation in the presence of timescale errors. It is based on the Lomb-Scargle periodogram for unevenly spaced time series, in combination with the Welch's Overlapped Segment Averaging procedure, bootstrap bias correction and persistence estimation. The timescale errors are modelled parametrically and included in the simulations for determining (1) the upper levels of the spectrum of the red-noise AR(1) alternative and (2) the uncertainty of the frequency of a spectral peak. Application of REDFITmc2 to ice core and stalagmite records of palaeoclimate allowed a more realistic evaluation of spectral peaks than when ignoring this source of uncertainty. The results support qualitatively the intuition that stronger effects on the spectrum estimate (decreased detectability and increased frequency uncertainty) occur for higher frequencies. The surplus information brought by algorithm REDFITmc2 is that those effects are quantified. Regarding timescale construction, not only the fixpoints, dating errors and the functional form of the age-depth model play a role. Also the joint distribution of all time points (serial correlation, stratigraphic order) determines spectrum estimation.

On trend estimation under monotone Gaussian subordination with long-memory: application to fossil pollen series

Fossil pollen data from stratigraphic cores are irregularly spaced in time due to non-linear age-depth relations. Moreover, their marginal distributions may vary over time. We address these features in a nonparametric regression model with errors that are monotone transformations of a latent continuous-time Gaussian process Z(T). Although Z(T) is unobserved, due to monotonicity, under suitable regularity conditions, it can be recovered facilitating further computations such as estimation of the long-memory parameter and the Hermite coefficients. The estimation of Z(T) itself involves estimation of the marginal distribution function of the regression errors. These issues are considered in proposing a plug-in algorithm for optimal bandwidth selection and construction of confidence bands for the trend function. Some high-resolution time series of pollen records from Lago di Origlio in Switzerland, which go back ca. 20,000 years are used to illustrate the methods.

A Nonstationary Space-Time Gaussian Process Model for Partially Converged Simulations

In the context of expensive numerical experiments, a promising solution for alleviating the computational costs consists of using partially converged simulations instead of exact solutions. The gain in computational time is at the price of precision in the response. This work addresses the issue of fitting a Gaussian process model to partially converged simulation data for further use in prediction. The main challenge consists of the adequate approximation of the error due to partial convergence, which is correlated in both design variables and time directions. Here, we propose fitting a Gaussian process in the joint space of design parameters and computational time. The model is constructed by building a nonstationary covariance kernel that reflects accurately the actual structure of the error. Practical solutions are proposed for solving parameter estimation issues associated with the proposed model. The method is applied to a computational fluid dynamics test case and shows significant improvement in prediction compared to a classical kriging model.

Surveillance of bovine tuberculosis and risk estimation of a future reservoir formation in wildlife in Switzerland and Liechtenstein.

Bovine tuberculosis (bTB) caused by Mycobacterium bovis or M. caprae has recently (re-) emerged in livestock and wildlife in all countries bordering Switzerland (CH) and the Principality of Liechtenstein (FL). Comprehensive data for Swiss and Liechtenstein wildlife are not available so far, although two native species, wild boar (Sus scrofa) and red deer (Cervus elaphus elaphus), act as bTB reservoirs elsewhere in continental Europe. Our aims were (1) to assess the occurrence of bTB in these wild ungulates in CH/FL and to reinforce scanning surveillance in all wild mammals; (2) to evaluate the risk of a future bTB reservoir formation in wild boar and red deer in CH/FL. Tissue samples collected from 2009 to 2011 from 434 hunted red deer and wild boar and from eight diseased ungulates with tuberculosis-like lesions were tested by direct real-time PCR and culture to detect mycobacteria of the Mycobacterium tuberculosis complex (MTBC). Identification of suspicious colonies was attempted by real-time PCR, genotyping and spoligotyping. Information on risk factors for bTB maintenance within wildlife populations was retrieved from the literature and the situation regarding identified factors was assessed for our study areas. Mycobacteria of the MTBC were detected in six out of 165 wild boar (3.6%; 95% CI: 1.4-7.8) but none of the 269 red deer (0%; 0-1.4). M. microti was identified in two MTBC-positive wild boar, while species identification remained unsuccessful in four cases. Main risk factors for bTB maintenance worldwide, including different causes of aggregation often resulting from intensive wildlife management, are largely absent in CH and FL. In conclusion, M. bovis and M. caprae were not detected but we report for the first time MTBC mycobacteria in Swiss wild boar. Present conditions seem unfavorable for a reservoir emergence, nevertheless increasing population numbers of wild ungulates and offal consumption may represent a risk.

Distinct lake level lowstand in Lake Prespa (SE Europe) at the time of the 74 (75) ka Toba eruption

Abstract. The 74 (75) ka Toba eruption in Sumatra, Indonesia, is considered to be one of the largest volcanic events during the Quaternary. Tephra from the Toba eruption has been found in many terrestrial and marine sedimentary deposits, and acidity peaks related to the eruption have been used to synchronize ice core records from Greenland and Antarctica. Seismic profiles and sedimentological data from Lake Prespa on the Balkan Peninsula, SE Europe, indicate a lake level lowstand at 73.6±7.7 ka based on ESR dating of shells. Tephrostratigraphy, radiocarbon dating and tuning of the total organic carbon content with the NGRIP isotope record, corroborate that the lake level lowstand was a short-term event superimposed on the general cooling trend at the end of MIS 5, most likely at the onset of the Greenland Stadial GS- 20. Acknowledging that tectonic events or karst processes could have triggered this lake level lowstand, the chronological correspondence between the lowstand and the Toba eruption is intriguing. Therefore a Toba-driven short-term shift to aridity in the Balkan region, leading to lake level changes and triggering spatial expansion events in one of the lake’s most abundant benthic species, the carino mussel Dreissena presbensis, cannot be excluded.

Mitotic figure counts are significantly overestimated in resection specimens of invasive breast carcinomas

Several authors have demonstrated an increased number of mitotic figures in breast cancer resection specimen when compared with biopsy material. This has been ascribed to a sampling artifact where biopsies are (i) either too small to allow formal mitotic figure counting or (ii) not necessarily taken form the proliferating tumor periphery. Herein, we propose a different explanation for this phenomenon. Biopsy and resection material of 52 invasive ductal carcinomas was studied. We counted mitotic figures in 10 representative high power fields and quantified MIB-1 immunohistochemistry by visual estimation, counting and image analysis. We found that mitotic figures were elevated by more than three-fold on average in resection specimen over biopsy material from the same tumors (20±6 vs 6±2 mitoses per 10 high power fields, P=0.008), and that this resulted in a relative diminution of post-metaphase figures (anaphase/telophase), which made up 7% of all mitotic figures in biopsies but only 3% in resection specimen (P<0.005). At the same time, the percentages of MIB-1 immunostained tumor cells among total tumor cells were comparable in biopsy and resection material, irrespective of the mode of MIB-1 quantification. Finally, we found no association between the size of the biopsy material and the relative increase of mitotic figures in resection specimen. We propose that the increase in mitotic figures in resection specimen and the significant shift towards metaphase figures is not due to a sampling artifact, but reflects ongoing cell cycle activity in the resected tumor tissue due to fixation delay. The dwindling energy supply will eventually arrest tumor cells in metaphase, where they are readily identified by the diagnostic pathologist. Taken together, we suggest that the rapidly fixed biopsy material better represents true tumor biology and should be privileged as predictive marker of putative response to cytotoxic chemotherapy.

5-ALA complete resections go beyond MR contrast enhancement: shift corrected volumetric analysis of the extent of resection in surgery for glioblastoma

BACKGROUND The technique of 5-aminolevulinic acid (5-ALA) tumor fluorescence is increasingly used to improve visualization of tumor tissue and thereby to increase the rate of patients with gross total resections. In this study, we measured the resection volumes in patients who underwent 5-ALA-guided surgery for non-eloquent glioblastoma and compared them with the preoperative tumor volume. METHODS We selected 13 patients who had received a complete resection according to intraoperative 5-ALA induced fluorescence and CRET according to post-operative T1 contrast-enhanced MRI. The volumes of pre-operative contrast enhancing tissue, post-operative resection cavity and resected tissue were determined through shift-corrected volumetric analysis. RESULTS The mean resection cavity (29 cm(3)) was marginally smaller than the pre-operative contrast-enhancing tumor (39 cm(3), p = 0.32). However, the mean overall resection volume (84 cm(3)) was significantly larger than the pre-operative contrast-enhancing tumor (39 cm(3), p = 0.0087). This yields a mean volume of resected 5-ALA positive, but radiological non-enhancing tissue of 45 cm(3). The mean calculated rim of resected tissue surpassed pre-operative tumor diameter by 6 mm (range 0-10 mm). CONCLUSIONS Results of the current study imply that (i) the resection cavity underestimates the volume of resected tissue and (ii) 5-ALA complete resections go significantly beyond the volume of pre-operative contrast-enhancing tumor bulk on MRI, indicating that 5-ALA also stains MRI non-enhancing tumor tissue. Use of 5-ALA may thus enable extension of coalescent tumor resection beyond radiologically evident tumor. The impact of this more extended resection method on time to progression and overall survival has not been determined, and potentially puts adjacent and functionally intact tissue at risk.

Estimation of male-to-female ratios of mutation rate in primates and rodents

Haldane (1935) developed a method for estimating the male-to-female ratio of mutation rate ($\alpha$) by using sex-linked recessive genetic disease, but in six different studies using hemophilia A data the estimates of $\alpha$ varied from 1.2 to 29.3. Direct genomic sequencing is a better approach, but it is laborious and not readily applicable to non-human organisms. To study the sex ratios of mutation rate in various mammals, I used an indirect method proposed by Miyata et al. (1987). This method takes advantage of the fact that different chromosomes segregate differently between males and females, and uses the ratios of mutation rate in sequences on different chromosomes to estimate the male-to-female ratio of mutation rate. I sequenced the last intron of ZFX and ZFY genes in 6 species of primates and 2 species of rodents; I also sequenced the partial genomic sequence of the Ube1x and Ube1y genes of mice and rats. The purposes of my study in addition to estimation of $\alpha$'s in different mammalian species, are to test the hypothesis that most mutations are replication dependent and to examine the generation-time effect on $\alpha$. The $\alpha$ value estimated from the ZFX and ZFY introns of the six primate specise is ${\sim}$6. This estimate is the same as an earlier estimate using only 4 species of primates, but the 95% confidence interval has been reduced from (2, 84) to (2, 33). The estimate of $\alpha$ in the rodents obtained from Zfx and Zfy introns is ${\sim}$1.9, and that deriving from Ube1x and Ube1y introns is ${\sim}$2. Both estimates have a 95% confidence interval from 1 to 3. These two estimates are very close to each other, but are only one-third of that of the primates, suggesting a generation-time effect on $\alpha$. An $\alpha$ of 6 in primates and 2 in rodents are close to the estimates of the male-to-female ratio of the number of germ-cell divisions per generation in humans and mice, which are 6 and 2, respectively, assuming the generation time in humans is 20 years and that in mice is 5 months. These findings suggest that errors during germ-cell DNA replication are the primary source of mutation and that $\alpha$ decreases with decreasing length of generation time. ^

Neuromuscular disuse and atrophy: {\it In vivo\/} magnetic resonance studies of intramuscular lipid metabolism and fluid shift changes

We postulated that neuromuscular disuse results in deleteriously affected tissue-vascular fluid exchange processes and subsequently damages the important oxidative bioenergetic process of intramuscular lipid metabolism. The in-depth research reported in the literature is somewhat limited by the ex vivo nature and sporadic time-course characterization of disuse atrophy and recovery. Thus, an in vivo controlled, localized animal model of disuse atrophy was developed in one of the hindlimbs of laboratory rabbits (employing surgically implanted tetrodotoxin (TTX)-filled mini-osmotic pump-sciatic nerve superfusion system) and tested repeatedly with magnetic resonance (MR) throughout the 2-week period of temporarily induced disuse and during the recovery period (following explantation of the TTX-filled pump) for a period of 3 weeks. Controls consisted of saline/"sham"-implanted rabbit hindlimbs. The validity of this model was established with repeated electrophysiologic nerve conduction testing using a clinically appropriate protocol and percutaneously inserted small needle stimulating and recording electrodes. Evoked responses recorded from proximal (P) and distal (D) sites to the sciatic nerve cuff in the TTX-implanted group revealed significantly decreased (p $<$ 0.001) proximal-to-distal (P/D) amplitude ratios (as much as 50-70% below Baseline/pre-implanted and sham-implanted group values) and significantly increased (p $<$ 0.01) differential latency (PL-DL) values (as much as 1.5 times the pre- and sham-implanted groups). By Day 21 of recovery, observed P/D and PL-DL levels matched Baseline/sham-implemented levels. MRI-determined cross-sectional area (CSA) values of Baseline/pre-implanted, sham- or TTX-implanted, and recovering/explanted and the corresponding contralateral hindlimb tibialis anterior (TA) muscles normalized to tibial bone (TB) CSA (in TA/TB ratios) revealed that there was a significant decline (indicative of atrophic response) from pre- and sham-implanted controls by as much as 20% (p $<$ 0.01) at Day 7 and 50-55% (p $<$ 0.001) at Day 13 of TTX-implantation. In the non-implanted contralaterals, a significant increase (indicative of hypertrophic response) by as much as 10% (p $<$ 0.025) at Day 7 and 27% (p $<$ 0.001) at Day 13 + TTX was found. The induced atrophic/hypertrophic TA muscles were observed to be fully recovered by Day 21 post-explantation as evidenced by image TA/TB ratios. End-point biopsy results from a small group of rabbits revealed comprehensive atrophy of both Type I and Type II fibers, although the heterogeneity of the response supports the use of image-guided, volume-localized proton magnetic resonance spectroscopy (MRS) to noninvasively assess tissue-level metabolic changes. MRS-determined results of a 0.25cc volume of tissue within implanted limb TA muscles under resting/pre-ischemic, ischemic-stressed, and post-ischemic conditions at timepoints during and following disuse atrophy/recovery revealed significantly increased intramuscular spectral lipid levels, as much as 2-3 times (p $<$ 0.01) the Baseline/pre-implanted values at Day 7 and 6-7 times (p $<$ 0.001) at Day 13 + TTX, which approached normal levels (compared to pre- and sham-implanted groups) by Day 21 of post-explanation recovery. (Abstract shortened by UMI.) ^

Estimation of aortic pressure waveforms from 4D phase-contrast MRI

Several approaches for the non-invasive MRI-based measurement of the aortic pressure waveform over the heart cycle have been proposed in the last years. These methods are normally based on time-resolved, two-dimensional phase-contrast sequences with uni-directionally encoded velocities (2D PC-MRI). In contrast, three-dimensional acquisitions with tridirectional velocity encoding (4D PC-MRI) have been shown to be a suitable data source for detailed investigations of blood flow and spatial blood pressure maps. In order to avoid additional MR acquisitions, it would be advantageous if the aortic pressure waveform could also be computed from this particular form of MRI. Therefore, we propose an approach for the computation of the aortic pressure waveform which can be completely performed using 4D PC-MRI. After the application of a segmentation algorithm, the approach automatically computes the aortic pressure waveform without any manual steps. We show that our method agrees well with catheter measurements in an experimental phantom setup and produces physiologically realistic results in three healthy volunteers.

Non-invasive and non-occlusive blood pressure estimation via a chest sensor

The clinical demand for a device to monitor Blood Pressure (BP) in ambulatory scenarios with minimal use of inflation cuffs is increasing. Based on the so-called Pulse Wave Velocity (PWV) principle, this paper introduces and evaluates a novel concept of BP monitor that can be fully integrated within a chest sensor. After a preliminary calibration, the sensor provides non-occlusive beat-by-beat estimations of Mean Arterial Pressure (MAP) by measuring the Pulse Transit Time (PTT) of arterial pressure pulses travelling from the ascending aorta towards the subcutaneous vasculature of the chest. In a cohort of 15 healthy male subjects, a total of 462 simultaneous readings consisting of reference MAP and chest PTT were acquired. Each subject was recorded at three different days: D, D+3 and D+14. Overall, the implemented protocol induced MAP values to range from 80 ± 6 mmHg in baseline, to 107 ± 9 mmHg during isometric handgrip maneuvers. Agreement between reference and chest-sensor MAP values was tested by using intraclass correlation coefficient (ICC = 0.78) and Bland-Altman analysis (mean error = 0.7 mmHg, standard deviation = 5.1 mmHg). The cumulative percentage of MAP values provided by the chest sensor falling within a range of ±5 mmHg compared to reference MAP readings was of 70%, within ±10 mmHg was of 91%, and within ±15mmHg was of 98%. These results point at the fact that the chest sensor complies with the British Hypertension Society (BHS) requirements of Grade A BP monitors, when applied to MAP readings. Grade A performance was maintained even two weeks after having performed the initial subject-dependent calibration. In conclusion, this paper introduces a sensor and a calibration strategy to perform MAP measurements at the chest. The encouraging performance of the presented technique paves the way towards an ambulatory-compliant, continuous and non-occlusive BP monitoring system.