String-averaging expectation-maximization for maximum likelihood estimation in emission tomography

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

Three-dimensional ordered-subset expectation maximization iterative protocol for evaluation of left ventricular volumes and function by quantitative gated SPECT: a dynamic phantom study.

The purposes of this study were to characterize the performance of a 3-dimensional (3D) ordered-subset expectation maximization (OSEM) algorithm in the quantification of left ventricular (LV) function with (99m)Tc-labeled agent gated SPECT (G-SPECT), the QGS program, and a beating-heart phantom and to optimize the reconstruction parameters for clinical applications. METHODS: A G-SPECT image of a dynamic heart phantom simulating the beating left ventricle was acquired. The exact volumes of the phantom were known and were as follows: end-diastolic volume (EDV) of 112 mL, end-systolic volume (ESV) of 37 mL, and stroke volume (SV) of 75 mL; these volumes produced an LV ejection fraction (LVEF) of 67%. Tomographic reconstructions were obtained after 10-20 iterations (I) with 4, 8, and 16 subsets (S) at full width at half maximum (FWHM) gaussian postprocessing filter cutoff values of 8-15 mm. The QGS program was used for quantitative measurements. RESULTS: Measured values ranged from 72 to 92 mL for EDV, from 18 to 32 mL for ESV, and from 54 to 63 mL for SV, and the calculated LVEF ranged from 65% to 76%. Overall, the combination of 10 I, 8 S, and a cutoff filter value of 10 mm produced the most accurate results. The plot of the measures with respect to the expectation maximization-equivalent iterations (I x S product) revealed a bell-shaped curve for the LV volumes and a reverse distribution for the LVEF, with the best results in the intermediate range. In particular, FWHM cutoff values exceeding 10 mm affected the estimation of the LV volumes. CONCLUSION: The QGS program is able to correctly calculate the LVEF when used in association with an optimized 3D OSEM algorithm (8 S, 10 I, and FWHM of 10 mm) but underestimates the LV volumes. However, various combinations of technical parameters, including a limited range of I and S (80-160 expectation maximization-equivalent iterations) and low cutoff values (< or =10 mm) for the gaussian postprocessing filter, produced results with similar accuracies and without clinically relevant differences in the LV volumes and the estimated LVEF.

Maximum Likelihood Estimation of Modal Parameters in Structures Using the Expectation Maximization Algorithm

This paper presents a time-domain stochastic system identification method based on maximum likelihood estimation (MLE) with the expectation maximization (EM) algorithm. The effectiveness of this structural identification method is evaluated through numerical simulation in the context of the ASCE benchmark problem on structural health monitoring. The benchmark structure is a four-story, two-bay by two-bay steel-frame scale model structure built in the Earthquake Engineering Research Laboratory at the University of British Columbia, Canada. This paper focuses on Phase I of the analytical benchmark studies. A MATLAB-based finite element analysis code obtained from the IASC-ASCE SHM Task Group web site is used to calculate the dynamic response of the prototype structure. A number of 100 simulations have been made using this MATLAB-based finite element analysis code in order to evaluate the proposed identification method. There are several techniques to realize system identification. In this work, stochastic subspace identification (SSI)method has been used for comparison. SSI identification method is a well known method and computes accurate estimates of the modal parameters. The principles of the SSI identification method has been introduced in the paper and next the proposed MLE with EM algorithm has been explained in detail. The advantages of the proposed structural identification method can be summarized as follows: (i) the method is based on maximum likelihood, that implies minimum variance estimates; (ii) EM is a computational simpler estimation procedure than other optimization algorithms; (iii) estimate more parameters than SSI, and these estimates are accurate. On the contrary, the main disadvantages of the method are: (i) EM algorithm is an iterative procedure and it consumes time until convergence is reached; and (ii) this method needs starting values for the parameters. Modal parameters (eigenfrequencies, damping ratios and mode shapes) of the benchmark structure have been estimated using both the SSI method and the proposed MLE + EM method. The numerical results show that the proposed method identifies eigenfrequencies, damping ratios and mode shapes reasonably well even in the presence of 10% measurement noises. These modal parameters are more accurate than the SSI estimated modal parameters.

Network lifetime maximization with cross-layer design in wireless sensor networks

This paper investigates a cross-layer design approach for minimizing energy consumption and maximizing network lifetime (NL) of a multiple-source and single-sink (MSSS) WSN with energy constraints. The optimization problem for MSSS WSN can be formulated as a mixed integer convex optimization problem with the adoption of time division multiple access (TDMA) in medium access control (MAC) layer, and it becomes a convex problem by relaxing the integer constraint on time slots. Impacts of data rate, link access and routing are jointly taken into account in the optimization problem formulation. Both linear and planar network topologies are considered for NL maximization (NLM). With linear MSSS and planar single-source and single-sink (SSSS) topologies, we successfully use Karush-Kuhn-Tucker (KKT) optimality conditions to derive analytical expressions of the optimal NL when all nodes are exhausted simultaneously. The problem for planar MSSS topology is more complicated, and a decomposition and combination (D&C) approach is proposed to compute suboptimal solutions. An analytical expression of the suboptimal NL is derived for a small scale planar network. To deal with larger scale planar network, an iterative algorithm is proposed for the D&C approach. Numerical results show that the upper-bounds of the network lifetime obtained by our proposed optimization models are tight. Important insights into the NL and benefits of cross-layer design for WSN NLM are obtained.

Within- and between-person variations as determinant factors to calculate the number of observations to estimate usual dietary intake of adolescents

This cross-sectional study describes the variation in the dietary intake of energy and macronutrients based on estimated coefficients of within- and between-person variation and intends to calculate the number of days required to evaluate the dietary intake of adolescents in Piracicaba, Brazil. The sample was made up by ninety-two adolescents aged between eleven and sixteen years. Interviews were performed to assess demographic, anthropometric and dietary variables, which were evaluated by the answers obtained through two 24-hour recalls. Descriptive analyses were performed and variances and Coefficients of variation within- and between-person were obtained by results of one-way ANOVA. The mean energy intake observed was 2,326.2 kilocalories; mean macronutrient intake came to 89.0 grams of fatty acids, 305.3 grams of carbohydrates and 82.2 grams of proteins. Coefficients of within-person variation ranged from 36% to 45%, while between-person coefficients varied from 26% to 31%. Variance ratios ranged from 1.35 (carbohydrate) to 2.62 (protein). The lowest number of observations required to correctly evaluate the usual diet calculated (r = 0.90) was six days, for carbohydrates while the highest was eleven, for proteins. Coefficients of variation were similar to those observed in international studies. We conclude that, in this study, two observations were sufficient to obtain reasonable correlations. However, six evaluations are necessary for these adolescents in order to obtain excellent correlations.

Pharmacogenetics of OATP Transporters Reveals That SLCO1B1 c.388A > G Variant Is Determinant of Increased Atorvastatin Response

Aims: The relationship between variants in SLCO1B1 and SLCO2B1 genes and lipid-lowering response to atorvastatin was investigated. Material and Methods: One-hundred-thirty-six unrelated individuals with hypercholesterolemia were selected and treated with atorvastatin (10 mg/day/4 weeks). They were genotyped with a panel of ancestry informative markers for individual African component of ancestry (ACA) estimation by SNaPshot (R) and SLCO1B1 (c.388A>G, c.463C>A and c.521T>C) and SLCO2B1 (-71T>C) gene polymorphisms were identified by TaqMan (R) Real-time PCR. Results: Subjects carrying SLCO1B1 c.388GG genotype exhibited significantly high low-density lipoprotein (LDL) cholesterol reduction relative to c.388AA+c.388AG carriers (41 vs. 37%, p = 0.034). Haplotype analysis revealed that homozygous of SLCO1B1*15 (c.521C and c.388G) variant had similar response to statin relative to heterozygous and non-carriers. A multivariate logistic regression analysis confirmed that c.388GG genotype was associated with higher LDL cholesterol reduction in the study population (OR: 3.2, CI95%: 1.3-8.0, p < 0.05). Conclusion: SLCO1B1 c.388A>G polymorphism causes significant increase in atorvastatin response and may be an important marker for predicting efficacy of lipid-lowering therapy.

ZETA DETERMINANT AND OPERATOR DETERMINANTS

We apply techniques of zeta functions and regularized products theory to study the zeta determinant of a class of abstract operators with compact resolvent, and in particular the relation with other spectral functions.

A contribution to the exact modal solution of in-plane beam structures

The exact vibration modes and natural frequencies of planar structures and mechanisms, comprised Euler-Bernoulli beams, are obtained by solving a transcendental. nonlinear, eigenvalue problem stated by the dynamic stiffness matrix (DSM). To solve this kind of problem, the most employed technique is the Wittrick-Williams algorithm, developed in the early seventies. By formulating a new type of eigenvalue problem, which preserves the internal degrees-of-freedom for all members in the model, the present study offers an alternative to the use of this algorithm. The new proposed eigenvalue problem presents no poles, so the roots of the problem can be found by any suitable iterative numerical method. By avoiding a standard formulation for the DSM, the local mode shapes are directly calculated and any extension to the beam theory can be easily incorporated. It is shown that the method here adopted leads to exact solutions, as confirmed by various examples. Extensions of the formulation are also given, where rotary inertia, end release, skewed edges and rigid offsets are all included. (C) 2008 Elsevier Ltd. All rights reserved.

The murine cytomegalovirus chemokine homolog, m131/129, is a determinant of viral pathogenicity

Chemokines are important mediators of the early inflammatory response to infection and modify a wide range of host immune responses. Functional homologs of cellular chemokines have been identified in a number of herpesviruses, suggesting that the subversion of the host chemokine response contributes to the pathogenesis of these viruses. Transcriptional and reverse transcription-PCR analyses demonstrated that the murine cytomegalovirus (MCMV) chemokine homolog, m131, was spliced at the 3' end to the adjacent downstream open reading frame, m129, resulting in a predicted product of 31 kDa, which is significantly larger than most known chemokines. The in vivo impact of m131/129 was investigated by comparing the replication of MCMV mutants having m131/129 deleted (Delta m131/129) with that of wild-type (wt) MCMV. Our studies demonstrate that both wt and Delta m131/129 viruses replicated to equivalent levels during the first 2 to 3 days following in vivo infection. However, histological studies demonstrated that the early inflammatory response elicited by Delta m131/129 was reduced compared with that of wt MCMV. Furthermore, the Delta m131/129 mutants failed to establish a high-titer infection in the salivary glands, These results suggest that m131/129 possesses proinflammatory properties in vivo and is important for the dissemination of MCMV to or infection of the salivary gland. Notably, the Delta m131/129 mutants were cleared more rapidly from the spleen and liver during acute infection compared with wt MCMV. The accelerated clearance of the mutants was dependent on NK cells and cells of the CD4(+) CD8(+) phenotype. These data suggest that m131/129 may also contribute to virus mechanisms of immune system evasion during early infection, possibly through the interference of NK cells and T cells.

Fast-convergent iterative scheme for filtering velocity signals and finding Kolmogorov scales

The present fundamental knowledge of fluid turbulence has been established primarily from hot- and cold-wire measurements. Unfortunately, however, these measurements necessarily suffer from contamination by noise since no certain method has previously been available to optimally filter noise from the measured signals. This limitation has impeded our progress of understanding turbulence profoundly. We address this limitation by presenting a simple, fast-convergent iterative scheme to digitally filter signals optimally and find Kolmogorov scales definitely. The great efficacy of the scheme is demonstrated by its application to the instantaneous velocity measured in a turbulent jet.

Frequency of thrombophilia determinant factors in patients with livedoid vasculopathy and treatment with anticoagulant drugs - a prospective study

Background Livedoid vasculopathy (LV) is a chronic idiopathic disease characterized by painful purpuric macules on lower extremities. Its exact aetiology remains uncertain, but thrombotic and microcirculatory phenomena have been implicated as possible pathogenic factors. Objectives To assess prospectively the frequency of thrombophilia and to verify the effectiveness of anticoagulant therapy among LV patients. Methods Thirty-four LV patients were tested for prothrombin time, activated partial thromboplastin time, antithrombin activity, protein C and S activity, anticardiolipin antibodies, lupus anticoagulant, prothrombin gene mutation, factor V Leiden mutation, methylenetetrahydrofolate reductase mutation, plasma homocysteine and fibrinogen. Thirteen of these patients were treated with anticoagulant drugs (either warfarin or heparin). Results Of 34 patients, 18 (52%) presented laboratory abnormalities of procoagulant conditions. Positive treatment response to anticoagulant therapy was observed in 11 patients. Improvement of pain was obtained in 1-3 weeks, an average of 1.8 week. Complete healing of the lesions was observed in about 2.3 months. Remission was sustained even after treatment interruption and lasted an average 7.8 months. No severe adverse effects were noticed. Conclusion The authors suggest all patients with diagnosis of LV to be investigated for thrombophilic status. Anticoagulant drugs were well tolerated and seemed to be effective in treating not only LV symptoms but also its ulcerations.

Is Full Postpleurodesis Lung Expansion a Determinant of a Successful Outcome After Talc Pleurodesis?

Study objectives: To analyze and compare radiologic lung expansion after tale pleurodesis performed either by videothoracoscopy or chest tube and correlate it with clinical outcome. Secondary end points evaluated were its follows: clinical efficacy; quality of life; safety; and survival. Methods: Prospective randomized study that included 60 patients (45 women, 15 men; mean age, 55.2 years) with recurrent malignant pleural effusion between January, 2005 and January 2008. They were randomized into the following two groups: video-assisted thoracic surgery (VATS) talc poudrage; and tale slurry (TS) administered through a chest tube. Lung expansion was evaluated through chest CT scans obtained 0, 1, 3 and 6 months after pleurodesis. Complications, drainage time, hospital stay,and quality of life (Medical Outcomes Study 36-item short form and World Health Organization quality-of-life questionnaires) were also analyzed. Results: There were no significant differences in preprocedure clinical and pathologic variables between groups. The immediate total (ie, > 90%) lung expansion was observed in 27 patients (45%) and wits more frequent in the VATS group (60% vs 30%, respectively; p = 0.027). During follow-up, 71% of the patients showed unaltered or improved lung expansion and 9 patients (15%) needed new pleural procedures (VATS group, 5 recurrences; TS group, 4 recurrences; p = 0.999). No differences, were found between groups regarding quality of life, complications, drainage time, hospital stay, and survival. Immediate lung expansion (lid not correlate with radiologic recurrence, clinical recurrence, or complications (p = 0.60, 0.15, and 0.20, respectively). Conclusion: Immediate partial lung expansion was a frequent finding and was more frequent after TS. Nonetheless, no correlation between immediate lung expansion and clinical outcome was found in this study. (CHEST 2009; 136:361-368)

5-Lipoxygenase is a key determinant of acute myocardial inflammation and mortality during Trypanosoma cruzi infection

This study provides evidence supporting the idea that although inflammatory cells migration to the cardiac tissue is necessary to control the growth of Trypanosoma cruzi, the excessive influx of such cells during acute myocarditis may be deleterious to the host. Production of lipid mediators of inflammation like leukotrienes (LTs) along with cytokines and chemokines largely influences the severity of inflammatory injury in response to tissue parasitism. T cruzi infection in mice deficient in 5-lipoxygenase (5-LO), the enzyme responsible for the synthesis of LTs and other lipid inflammatory mediators, resulted in transiently increased parasitemia, and improved survival rate compared with WT mice. Myocardia from 5-LO(-/-) mice exhibited reduced inflammation, collagen deposition, and migration of CD4(+), CD8(+), and IFN-gamma-producer cells compared with WT littermates. Moreover, decreased amounts of TNF-alpha, IFN-gamma, and nitric oxide synthase were found in the hearts of 5-LO(-/-) mice. Interestingly, despite of early higher parasitic load, 5-LO(-/-) mice survived, and controlled T cruzi infection. These results show that efficient parasite clearance is possible in a context of moderate inflammatory response, as occurred in 5-LO(-/-) mice, in which reduced myocarditis protects the animals during T cruzi infection. (c) 2010 Elsevier Masson SAS. All rights reserved.

Tyrosine residue 271 of the norepinephrine transporter is an important determinant of its pharmacology

The aim was to examine the functional importance in the norepinephrine transporter (NET) of (i) the phenylalanine residue at position 531 in transmembrane domain (TMD) 11 by mutating it to tyrosine in the rat (rF531Y) and human (hF531Y) NETs and (ii) the highly conserved tyrosine residues at positions 249 in TMD 4 of human NET (hNET) (mutated to alanine: hY249A) and 271 in TMD 5, by mutating to alanine (hY271A), phenylalanine (hY271F) and histidine (hY271H). The effects of the mutations on NET function were for uptake of the substrates, examined by expressing the mutant and wildtype NETs in COS-7 cells and measuring the K-m and V-max for uptake of the substrates, [H-3]norepinephrine, [H-3]MPP+ and [H-3]dopamine, the K-D and B-max for [H-3]nisoxetine binding and the K-i of the inhibitors, nisoxetine, desipramine and cocaine, for inhibition of [H-3]norepinephrine uptake. The K-m values of the substrates were lower for the mutants at amino acid 271 than hNET and unaffected for the other mutants, and each mutant had a significantly lower than NET for substrate uptake. The mutations at position 271 caused an increase in the K-i or K-D values of nisoxetine, desipramine and cocaine, but there were no effects for the other mutations. Hence, the 271 tyrosine residue in TMD 5 is an important determinant of NET function, with the mutants showing an increase in the apparent affinities of substrates and a decrease in the apparent affinities of inhibitors, but the 249 tyrosine and 531 phenylalanine residues do not have a major role in determining NET function. (C) 2001 Elsevier Science B.V. All rights reserved.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.