Methods of visualization and analysis of cardiac depolarization in the three dimensional space

The master thesis presents methods for intellectual analysis and visualization 3D EKG in order to increase the efficiency of ECG analysis by extracting additional data. Visualization is presented as part of the signal analysis tasks considered imaging techniques and their mathematical description. Have been developed algorithms for calculating and visualizing the signal attributes are described using mathematical methods and tools for mining signal. The model of patterns searching for comparison purposes of accuracy of methods was constructed, problems of a clustering and classification of data are solved, the program of visualization of data is also developed. This approach gives the largest accuracy in a task of the intellectual analysis that is confirmed in this work. Considered visualization and analysis techniques are also applicable to the multi-dimensional signals of a different kind.

Anemia screening in potential female blood donors: comparison of two different quantitative methods

Anemia screening before blood donation requires an accurate, quick, practical, and easy method with minimal discomfort for the donors. The aim of this study was to compare the accuracy of two quantitative methods of anemia screening: the HemoCue 201(+) (Aktiebolaget Leo Diagnostics) hemoglobin (Hb) and microhematocrit (micro-Hct) tests. Two blood samples of a single fingerstick were obtained from 969 unselected potential female donors to determine the Hb by HemoCue 201(+) and micro-Hct using HemataSTAT II (Separation Technology, Inc.), in alternating order. From each participant, a venous blood sample was drawn and run in an automatic hematology analyzer (ABX Pentra 60, ABX Diagnostics). Considering results of ABX Pentra 60 as true values, the sensitivity and specificity of HemoCue 201(+) and micro-Hct as screening methods were compared, using a venous Hb level of 12.0 g per dL as cutoff for anemia. The sensitivities of the HemoCue 201(+) and HemataSTAT II in detecting anemia were 56 percent (95% confidence interval [CI], 46.1%-65.5%) and 39.5 percent (95% CI, 30.2%-49.3%), respectively (p < 0.001). Analyzing only candidates with a venous Hb level lower than 11.0 g per dL, the deferral rate was 100 percent by HemoCue 201(+) and 77 percent by HemataSTAT II. The specificities of the methods were 93.5 and 93.2 percent, respectively. The HemoCue 201(+) showed greater discriminating power for detecting anemia in prospective blood donors than the micro-Hct method. Both presented equivalent deferral error rates of nonanemic potential donors. Compared to the micro-Hct, HemoCue 201(+) reduces the risk of anemic female donors giving blood, specially for those with lower Hb levels, without increasing the deferral of nonanemic potential donors.

Molecular dynamics computer simulation of scratch resistance testing of polymers: visualization

Experimental scratch resistance testing provides two numbers: the penetration depth Rp and the healing depth Rh. In molecular dynamics computer simulations, we create a material consisting of N statistical chain segments by polymerization; a reinforcing phase can be included. Then we simulate the movement of an indenter and response of the segments during X time steps. Each segment at each time step has three Cartesian coordinates of position and three of momentum. We describe methods of visualization of results based on a record of 6NX coordinates. We obtain a continuous dependence on time t of positions of each of the segments on the path of the indenter. Scratch resistance at a given location can be connected to spatial structures of individual polymeric chains.

Promising blood-derived biomarkers for estimation of the postmortem interval

A precise estimation of the postmortem interval (PMI) is one of the most important topics in forensic pathology. However, the PMI estimation is based mainly on the visual observation of cadaverous pheno- mena (e.g. algor, livor and rigor mortis) and on alternative methods such as thanatochemistry that remain relatively imprecise. The aim of this in vitro study was to evaluate the kinetic alterations of several bio- chemical parameters (i.e. proteins, enzymes, substrates, electrolytes and lipids) during putrefaction of human blood. For this purpose, we performed kinetic biochemical analysis during a 264 hour period. The results showed a significant linear correlation between total and direct bilirubin, urea, uric acid, transferrin, immunoglobulin M (IgM), creatine kinase (CK), aspartate transaminase (AST), calcium and iron with the time of blood putrefaction. These parameters allowed us to develop two mathematical models that may have predictive values and become important complementary tools of traditional methods to achieve a more accurate PMI estimation

Ultrasound-guided suprascapular nerve block, description of a novel supraclavicular approach.

BACKGROUND AND OBJECTIVES: The suprascapular nerve (SSN) block is frequently performed for different shoulder pain conditions and for perioperative and postoperative pain control after shoulder surgery. Blind and image-guided techniques have been described, all of which target the nerve within the supraspinous fossa or at the suprascapular notch. This classic target point is not always ideal when ultrasound (US) is used because it is located deep under the muscles, and hence the nerve is not always visible. Blocking the nerve in the supraclavicular region, where it passes underneath the omohyoid muscle, could be an attractive alternative. METHODS: In the first step, 60 volunteers were scanned with US, both in the supraclavicular and the classic target area. The visibility of the SSN in both regions was compared. In the second step, 20 needles were placed into or immediately next to the SSN in the supraclavicular region of 10 cadavers. The accuracy of needle placement was determined by injection of dye and following dissection. RESULTS: In the supraclavicular region of volunteers, the nerve was identified in 81% of examinations (95% confidence interval [CI], 74%-88%) and located at a median depth of 8 mm (interquartile range, 6-9 mm). Near the suprascapular notch (supraspinous fossa), the nerve was unambiguously identified in 36% of examinations (95% CI, 28%-44%) (P < 0.001) and located at a median depth of 35 mm (interquartile range, 31-38 mm; P < 0.001). In the cadaver investigation, the rate of correct needle placement of the supraclavicular approach was 95% (95% CI, 86%-100%). CONCLUSIONS: Visualization of the SSN with US is better in the supraclavicular region as compared with the supraspinous fossa. The anatomic dissections confirmed that our novel supraclavicular SSN block technique is accurate.

Therapeutic drug monitoring of imatinib: Bayesian and alternative methods to predict trough levels

Background: The imatinib trough plasma concentration (C(min)) correlates with clinical response in cancer patients. Therapeutic drug monitoring (TDM) of plasma C(min) is therefore suggested. In practice, however, blood sampling for TDM is often not performed at trough. The corresponding measurement is thus only remotely informative about C(min) exposure. Objectives: The objectives of this study were to improve the interpretation of randomly measured concentrations by using a Bayesian approach for the prediction of C(min), incorporating correlation between pharmacokinetic parameters, and to compare the predictive performance of this method with alternative approaches, by comparing predictions with actual measured trough levels, and with predictions obtained by a reference method, respectively. Methods: A Bayesian maximum a posteriori (MAP) estimation method accounting for correlation (MAP-ρ) between pharmacokinetic parameters was developed on the basis of a population pharmacokinetic model, which was validated on external data. Thirty-one paired random and trough levels, observed in gastrointestinal stromal tumour patients, were then used for the evaluation of the Bayesian MAP-ρ method: individual C(min) predictions, derived from single random observations, were compared with actual measured trough levels for assessment of predictive performance (accuracy and precision). The method was also compared with alternative approaches: classical Bayesian MAP estimation assuming uncorrelated pharmacokinetic parameters, linear extrapolation along the typical elimination constant of imatinib, and non-linear mixed-effects modelling (NONMEM) first-order conditional estimation (FOCE) with interaction. Predictions of all methods were finally compared with 'best-possible' predictions obtained by a reference method (NONMEM FOCE, using both random and trough observations for individual C(min) prediction). Results: The developed Bayesian MAP-ρ method accounting for correlation between pharmacokinetic parameters allowed non-biased prediction of imatinib C(min) with a precision of ±30.7%. This predictive performance was similar for the alternative methods that were applied. The range of relative prediction errors was, however, smallest for the Bayesian MAP-ρ method and largest for the linear extrapolation method. When compared with the reference method, predictive performance was comparable for all methods. The time interval between random and trough sampling did not influence the precision of Bayesian MAP-ρ predictions. Conclusion: Clinical interpretation of randomly measured imatinib plasma concentrations can be assisted by Bayesian TDM. Classical Bayesian MAP estimation can be applied even without consideration of the correlation between pharmacokinetic parameters. Individual C(min) predictions are expected to vary less through Bayesian TDM than linear extrapolation. Bayesian TDM could be developed in the future for other targeted anticancer drugs and for the prediction of other pharmacokinetic parameters that have been correlated with clinical outcomes.

Analysis of matched matrices

We consider the joint visualization of two matrices which have common rowsand columns, for example multivariate data observed at two time pointsor split accord-ing to a dichotomous variable. Methods of interest includeprincipal components analysis for interval-scaled data, or correspondenceanalysis for frequency data or ratio-scaled variables on commensuratescales. A simple result in matrix algebra shows that by setting up thematrices in a particular block format, matrix sum and difference componentscan be visualized. The case when we have more than two matrices is alsodiscussed and the methodology is applied to data from the InternationalSocial Survey Program.

Intermittent atrial tachycardia facilitates atrial fibrillation by a shortening of activation recovery interval

Introduction: We recently observed in a chronic ovine model that a shortening of action potential duration (APD) as assessed by the activation recovery interval (ARI) may be a mechanism whereby pacing-induced atrial tachycardia (PIAT) facilitates atrial fibrillation (AF), mediated by a return to 1:1 atrial capture after the effective refractory period has been reached. The aim of the present study is to evaluate the effect of long term intermittent burst pacing on ARI before induction of AF.Methods: We specifically developed a chronic ovine model of PIAT using two pacemakers (PM) each with a right atrial (RA) lead separated by ∼2cm. The 1st PM (Vitatron T70) was used to record a broadband unipolar RA EGM (800 Hz, 0.4 Hz high pass filter). The 2nd was used to deliver PIAT during electrophysiological protocols at decremental pacing CL (400 beats, from 400 to 110ms) and long term intermittent RA burst pacing to promote electrical remodeling (5s of burst followed by 2s of sinus rhythm) until onset of sustained AF. ARI was defined as the time difference between the peak of the atrial repolarization wave and the first atrial depolarization. The mean ARIs of paired sequences (before and after remodeling), each consisting of 20 beats were compared.Results: As shown in the figure, ARIs (n=4 sheep, 46 recordings) decreased post remodeling compared to baseline (86±19 vs 103±12 ms, p<0.05). There was no difference in atrial structure as assessed by light microscopy between control and remodeled sheep.Conclusions: Using standard pacemaker technology, atrial ARIs as a surrogate of APDs were successfully measured in vivo during the electrical remodeling process leading to AF. The facilitation of AF by PIAT mimicking salvos from pulmonary veins is heralded by a significant shortening of ARI.

Study on period search methods of variable stars: Application to ASAS and CRTS databases

Study on variable stars is an important topic of modern astrophysics. After the invention of powerful telescopes and high resolving powered CCD’s, the variable star data is accumulating in the order of peta-bytes. The huge amount of data need lot of automated methods as well as human experts. This thesis is devoted to the data analysis on variable star’s astronomical time series data and hence belong to the inter-disciplinary topic, Astrostatistics. For an observer on earth, stars that have a change in apparent brightness over time are called variable stars. The variation in brightness may be regular (periodic), quasi periodic (semi-periodic) or irregular manner (aperiodic) and are caused by various reasons. In some cases, the variation is due to some internal thermo-nuclear processes, which are generally known as intrinsic vari- ables and in some other cases, it is due to some external processes, like eclipse or rotation, which are known as extrinsic variables. Intrinsic variables can be further grouped into pulsating variables, eruptive variables and flare stars. Extrinsic variables are grouped into eclipsing binary stars and chromospheri- cal stars. Pulsating variables can again classified into Cepheid, RR Lyrae, RV Tauri, Delta Scuti, Mira etc. The eruptive or cataclysmic variables are novae, supernovae, etc., which rarely occurs and are not periodic phenomena. Most of the other variations are periodic in nature. Variable stars can be observed through many ways such as photometry, spectrophotometry and spectroscopy. The sequence of photometric observa- xiv tions on variable stars produces time series data, which contains time, magni- tude and error. The plot between variable star’s apparent magnitude and time are known as light curve. If the time series data is folded on a period, the plot between apparent magnitude and phase is known as phased light curve. The unique shape of phased light curve is a characteristic of each type of variable star. One way to identify the type of variable star and to classify them is by visually looking at the phased light curve by an expert. For last several years, automated algorithms are used to classify a group of variable stars, with the help of computers. Research on variable stars can be divided into different stages like observa- tion, data reduction, data analysis, modeling and classification. The modeling on variable stars helps to determine the short-term and long-term behaviour and to construct theoretical models (for eg:- Wilson-Devinney model for eclips- ing binaries) and to derive stellar properties like mass, radius, luminosity, tem- perature, internal and external structure, chemical composition and evolution. The classification requires the determination of the basic parameters like pe- riod, amplitude and phase and also some other derived parameters. Out of these, period is the most important parameter since the wrong periods can lead to sparse light curves and misleading information. Time series analysis is a method of applying mathematical and statistical tests to data, to quantify the variation, understand the nature of time-varying phenomena, to gain physical understanding of the system and to predict future behavior of the system. Astronomical time series usually suffer from unevenly spaced time instants, varying error conditions and possibility of big gaps. This is due to daily varying daylight and the weather conditions for ground based observations and observations from space may suffer from the impact of cosmic ray particles. Many large scale astronomical surveys such as MACHO, OGLE, EROS, xv ROTSE, PLANET, Hipparcos, MISAO, NSVS, ASAS, Pan-STARRS, Ke- pler,ESA, Gaia, LSST, CRTS provide variable star’s time series data, even though their primary intention is not variable star observation. Center for Astrostatistics, Pennsylvania State University is established to help the astro- nomical community with the aid of statistical tools for harvesting and analysing archival data. Most of these surveys releases the data to the public for further analysis. There exist many period search algorithms through astronomical time se- ries analysis, which can be classified into parametric (assume some underlying distribution for data) and non-parametric (do not assume any statistical model like Gaussian etc.,) methods. Many of the parametric methods are based on variations of discrete Fourier transforms like Generalised Lomb-Scargle peri- odogram (GLSP) by Zechmeister(2009), Significant Spectrum (SigSpec) by Reegen(2007) etc. Non-parametric methods include Phase Dispersion Minimi- sation (PDM) by Stellingwerf(1978) and Cubic spline method by Akerlof(1994) etc. Even though most of the methods can be brought under automation, any of the method stated above could not fully recover the true periods. The wrong detection of period can be due to several reasons such as power leakage to other frequencies which is due to finite total interval, finite sampling interval and finite amount of data. Another problem is aliasing, which is due to the influence of regular sampling. Also spurious periods appear due to long gaps and power flow to harmonic frequencies is an inherent problem of Fourier methods. Hence obtaining the exact period of variable star from it’s time series data is still a difficult problem, in case of huge databases, when subjected to automation. As Matthew Templeton, AAVSO, states “Variable star data analysis is not always straightforward; large-scale, automated analysis design is non-trivial”. Derekas et al. 2007, Deb et.al. 2010 states “The processing of xvi huge amount of data in these databases is quite challenging, even when looking at seemingly small issues such as period determination and classification”. It will be beneficial for the variable star astronomical community, if basic parameters, such as period, amplitude and phase are obtained more accurately, when huge time series databases are subjected to automation. In the present thesis work, the theories of four popular period search methods are studied, the strength and weakness of these methods are evaluated by applying it on two survey databases and finally a modified form of cubic spline method is intro- duced to confirm the exact period of variable star. For the classification of new variable stars discovered and entering them in the “General Catalogue of Vari- able Stars” or other databases like “Variable Star Index“, the characteristics of the variability has to be quantified in term of variable star parameters.

Improvements in the ray tracing of implicit surfaces based on interval arithmetic

Las superfícies implícitas son útiles en muchas áreasde los gráficos por ordenador. Una de sus principales ventajas es que pueden ser fácilmente usadas como primitivas para modelado. Aun asi, no son muy usadas porque su visualización toma bastante tiempo. Cuando se necesita una visualización precisa, la mejor opción es usar trazado de rayos. Sin embargo, pequeñas partes de las superficies desaparecen durante la visualización. Esto ocurre por la truncación que se presenta en la representación en punto flotante de los ordenadores; algunos bits se puerden durante las operaciones matemáticas en los algoritmos de intersección. En este tesis se presentan algoritmos para solucionar esos problemas. La investigación se basa en el uso del Análisis Intervalar Modal el cual incluye herramientas para resolver problemas con incertidumbe cuantificada. En esta tesis se proporcionan los fundamentos matemáticos necesarios para el desarrollo de estos algoritmos.

Evaluating interval forecasts of high-frequency financial data

A number of methods of evaluating the validity of interval forecasts of financial data are analysed, and illustrated using intraday FTSE100 index futures returns. Some existing interval forecast evaluation techniques, such as the Markov chain approach of Christoffersen (1998), are shown to be inappropriate in the presence of periodic heteroscedasticity. Instead, we consider a regression-based test, and a modified version of Christoffersen's Markov chain test for independence, and analyse their properties when the financial time series exhibit periodic volatility. These approaches lead to different conclusions when interval forecasts of FTSE100 index futures returns generated by various GARCH(1,1) and periodic GARCH(1,1) models are evaluated.

Toward the Optimization of an e-Tongue System Using Information Visualization: A Case Study with Perylene Tetracarboxylic Derivative Films in the Sensing Units

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Combined TUNEL and TRAP methods suggest that apoptotic bone cells are inside vacuoles of alveolar bone osteoclasts in young rats

Although it is generally accepted that osteoclasts breakdown and resorb bone matrix, the possibility that they may also be able to engulf apoptotic osteoblasts/ lining cells and/or osteocytes remains controversial. Apoptosis of osteoblasts/ lining cells and/or osteocytes and interactions between these cells and osteoclasts are extremely rapid events that are difficult to observe in viva. A suitable in viva model for studying these events is the alveolar bone of young rats because it is continuously. Thus, sections of aldehyde fixed alveolar undergoing intense resorption/remodeling bone of young rats were stained by the combined terminal deoxynucleotidyl transferase-mediated dUTP nick end labeling (TUNEL) method and the tartrate-resistant acid phosphatase (TRAP) method for the simultaneous visualization of apoptotic cells and osteoclasts in the same section. The combined TUNEL and TRAP reactions, in the same section, greatly facilitated visualization of relationship between osteoclasts and apoptotic bone cells during alveolar bone remodeling. Our results showed that several TRAP-positive osteoclasts exhibited large vacuoles containing TUNEL positive apoptotic structures, probably derived from osteoblasts/lining cells and/or osteocytes. These results support the idea that alveolar bone osteoclasts are able to internalize dying apoptotic bone cells.

Components of the acoustic swallowing signal: Preliminary study

Purpose: To analyze the components of the acoustic signal of swallowing using a specific software. Methods: Fourteen healthy subjects ranging in age from 20 to 50 years (mean age 31±10 years), were evaluated. Data collection consisted on the simultaneous capture of the swallowing audio with a microphone and of the swallowing videofluoroscopic image. The bursts of the swallowing acoustic signal were identified and their duration and the interval between them were later analyzed using a specific software, which allowed the simultaneous analyses between the acoustic wave and the videofluoroscopic image. Results: Three burst components were identified in most of the swallows evaluated. The first burst presented mean time of 87.3 milliseconds (ms) for water and 78.2 for the substance. The second burst presented mean time of 112.9 ms for water and 85.5 for the pasty substance. The mean interval between first and second burst was 82.1 ms for water and 95.3 ms for the pasty consistency, and between second and third burst was 339.8 ms for water and 322.0 ms for the pasty consistency. Conclusion: The software allowed the visualization of three bursts during the swallowing of healthy individuals, and showed that the swallowing signal in normal subjects is highly variable.

Evaluation of early changes induced by diuron in the rat urinary bladder using different processing methods for scanning electron microscopy

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)