Variational Physiologically Informed Solution to Hemodynamic and Perfusion Response Estimation from ASL fmri Data
|
|
- Phillip Singleton
- 5 years ago
- Views:
Transcription
1 Variational Pysioloically Informed Solution to Hemodynamic and Perfusion Response Estimation from ASL fmri Data Aina Frau-Pascual, Florence Forbes, Pilippe Ciuciu To cite tis version: Aina Frau-Pascual, Florence Forbes, Pilippe Ciuciu. Variational Pysioloically Informed Solution to Hemodynamic and Perfusion Response Estimation from ASL fmri Data International Worksop on Pattern Reconition in NeuroImain, Jun 2015, Stanford, CA, United States. IEEE, pp.57-60, 2015, 2015 International Worksop on Pattern Reconition in NeuroImain, Stanford, CA, USA, June 10-12, < /PRNI >. <al > HAL Id: al ttps://al.arcives-ouvertes.fr/al Submitted on 7 Jan 2016 HAL is a multi-disciplinary open access arcive for te deposit and dissemination of scientific researc documents, weter tey are publised or not. Te documents may come from teacin and researc institutions in France or abroad, or from public or private researc centers. L arcive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recerce, publiés ou non, émanant des établissements d enseinement et de recerce français ou étraners, des laboratoires publics ou privés.
2 Variational pysioloically informed solution to emodynamic and perfusion response estimation from ASL fmri data Aina Frau-Pascual 1,3, Florence Forbes 1, Pilippe Ciuciu 2,3 1 INRIA, MISTIS, Grenoble University, LJK, Grenoble, France 2 CEA/DSV/IBM NeuroSpin center, Bât. 145, F Gif-sur-Yvette, France 3 INRIA, Parietal, F Orsay, France Abstract Functional Arterial Spin Labelin (fasl) MRI can provide a quantitative measurement of cerebral blood flow. A joint detection-estimation (JDE) framework as been considered to extract task-related perfusion and emodynamic responses not restricted to canonical response function sapes. In tis work, we provide a variational expectation-maximization (VEM) aloritm for emodynamic and perfusion responses estimation. Tis approac provides a lower computational load compared to previous attempts, and facilitates te incorporation of prior knowlede and constraints in te estimation. Validation on simulated and real data sets as been performed. I. INTRODUCTION Functional Arterial Spin Labelin (fasl) is a functional MRI modality tat is able to provide a quantitative measurement of cerebral blood flow (CBF) and its variations elicited by specific tasks. Altou ASL [1] is mainly used as a static probe of cerebral perfusion, it as also been used in functional MRI as an alternative to te standard blood-oxyenlevel-dependent (BOLD) [2] modality, as it can ive more specific information about brain function. FASL data consists of alternatin pairs of control and manetically taed ( ta ) imaes. Local CBF or perfusion canes can be measured by doin te control-ta difference, so as to et rid of te emodynamic component (BOLD effect) of te sinal contained in bot control and ta imaes. Many pairs (ą 50) of control-ta imaes need to be acquired to compensate for te low Sinal-to-Noise Ratio (SNR) of tis difference. Te standard approac for fasl analysis is eneral linear model (GLM) [3], [4] usin te stimulus-induced canonical emodynamic response function (HRF) as reressor for bot components of te sinal: emodynamic and perfusion components. Te latter one encodes control-ta differences. Altou te canonical HRF is enerally used, it as been just calibrated on BOLD experiments and it reflects simultaneous variations of CBF, cerebral blood volume (CBV) and cerebral oxyen consumption (CMRO 2 ). In contrast, te perfusion response function (PRF) only reflects te variation in CBF and tus departs from te canonical sape. Moreover, tere as been evidence tat tis response can vary between reions and across subjects so tat assumin a constant response sape mit be misleadin. Some flexibility in te response sape as already been considered in BOLD GLM analysis by usin time/dispersion derivatives [5] or finite impulse response models. In a joint detection-estimation (JDE) framework bot taskrelated perfusion and emodynamic responses can be estimated as well as perfusion-related and BOLD-related maps of evoked activity. JDE as been oriinally developed for BOLD data analysis [6], [7] and as been extended in [8] for fasl data analysis, considerin Monte Carlo Markov Cains (MCMC) metods. In tis work, followin te spirit of [7], we provide an alternative solution based on te variational expectation-maximization (VEM) aloritm. Tis framework is more convenient to deal wit constraints (e, normalization or positivity) in te M-step. JDE-VEM provides comparable results to JDE-MCMC for a muc lower computational load. One of te novelties is te introduction of prior knowlede as in [9] trou te relationsip between perfusion and emodynamic responses derived from pysioloical models [10], [11]. Tis relationsip allows us to inform te PRF estimation from te HRF one, as te emodynamic component as a ier SNR tan te perfusion one due to te acquisition procedure. Experiments on simulated and real data sow te ood performance of te metod, wit a fast converence of te parameter estimates. II. JOINT DETECTION ESTIMATION MODEL FOR FASL In fasl te manetically taed imae reflects a perfusion effect besides te emodynamic effect, tat can be subtracted by comparison wit te control imae. Te control imae is a noisy BOLD sinal because of te effect of te manetic tain. Te different ASL sinal components can be explained wit te ASL JDE model [8], [12], a reion-based approac tat considers functionally omoeneous reions. In a reion P comprisin J voxels, te enerative model for j P J: Mÿ y j a looomooon m j X m ` looooomooooon c m j W X m ` lomon P l j ` lomon α j w ` lomon b j m1 paq pbq pcq pdq peq (1) in te presence of M experimental conditions. Te data y j P R N can be decomposed into (a) task-related emodynamic
3 and (b) perfusion components; (c) a drift component P l j ; (d) a perfusion baseline term α j w 1 wic completes te modellin of te perfusion component and (e) a noise term, assumed wite Gaussian of variance σj 2. Te control/ta effect is implicit in te ASL JDE model iven te use of W diapwq. Vectors and represent te D-dimensional (D ă N) unknown HRF and PRF sapes, constant witin P. Te manitudes of activation or response levels for emodynamic and perfusion components are a ( ( a m j and c c m j and denoted as BRLs and PRLs ereafter. X P R NˆD is a binary matrix tat encodes te laed onset stimuli. Te response levels are assumed to follow different spatial Gaussian mixture models but overned by M common binary idden Markov random fields q m wit q m tq mj, j P Pu encodin voxels activation states for eac experimental condition m. BRLs and PRLs are assumed independent conditionally to q tq m, m 1 : Mu. For furter detail, please refer to [8]. Te perfusion component in te ASL sinal as a very low SNR owin to its small size captured by te controlta subtraction. To address tis issue, a conditional prior as been used to inform te PRF estimation from te less noisy HRF recovery durin inference. To link tese two responses, a relationsip Ω was derived in [9] from pysioloical models. In our model, we consider tat HRF and PRF sapes follow prior Gaussian distributions N p0, v Σ q and N pω, v Σ q, wit covariance matrices Σ and Σ encodin a constraint on te second order derivative so as to account for temporal smootness. We also consider constraints on te response functions to enforce teir L 2 -norm to be 1, i.e., P S 2 pr D`1 q, S 2 pr D`1 q bein te L 2 unit ball of R D`1. III. VARIATIONAL EM In our fasl model, tere are missin parameters tat need to be estimated: a P A, P H, c P C, P G, q P Q. In tis work, an Expectation-Maximization (EM)framework is proposed to perform te estimation. EM can be viewed [13] as an alternatin maximization procedure of a function F suc tat for any p P D, bein D te set of all probability distributions on A ˆ H ˆ C ˆ G ˆ Q, F p p, θq E p lo ppy, a,, c,, q ; θq ` Ir ps. Here Ir ps E p lo ppa,, c,, qq is te entropy of p, and E p denotes te expectation wit respect to p. Maximizin function F is equivalent to minimizin te Kullback-Leibler diverence KLp p, ppa,, c,, q yqq, tat ives a measure of te distance between two distributions, in tis case te approximation p and te true ppa,, c,, q yq distributions. Tis view of EM as led to a number of variants in wic te E-step is solved over a restricted class of probability distributions, D. Te variational approac corresponds to D cosen as te set of distributions tat factorize as ppa,, c,, qq p a paq p pq p c pcq p pq p q pqq were p a P D A, p P D H, p c P D C, p P D G and p q P D Q, te sets of probability distributions on A, H, C, G, Q respectively. Note 1 Vector w is N-dimensional suc tat w tn 1{2 if t n is even (control) and w tn 1{2 oterwise (taed). tat te dependence between random variables, as in MCMC, is translated into dependence between statistical moments in VEM. Te E-step becomes an approximate E-step tat can be furter decomposed into five staes updatin te different variables in turn. At iteration prq, wit current estimates denoted by p pr 1q a, p pr 1q, p pr 1q c, p pr 1q, p pr 1q q and θ prq, te updatin formulae are of te form: E-H-step: ar max p prq E-G-step: p prq p PD H F p p pr 1q a p p pr 1q c p pr 1q p pr 1q q ; θ prq q (2) ar max F p p a pr 1q p prq p c pr 1q p p pr 1q q ; θ prq q (3) p PD G wit similar expressions for te oter steps obtained by permutin te roles of te variables. Hereafter, for te ease of presentation, te prq and pr 1q superscripts are omitted. We also consider normalization constraints on and so tat our solution as unit L 2 -norm. For tat matter, we modify te seek variational approximation to p p a δ p c δ p q, were te probabilities on and are replaced by Dirac functions. Tis reduces te searc to pointwise estimates and. Te E-H and E-G steps in Eqs. (2)-(3) yield ten maximizations problems wic are easily constrained to account for normalization: E-H: ar max E p a p c p q lo pp y, a, c,, q; θq (4) E-G: ar max E pa p c p q lo pp y, a,, c, q; θq (5) It is straitforward to see tat (4) and (5) amount to minimizin a quadratic function under a quadratic constraint, namely }} 2 1 and }} 2 1 respectively. Te oter E-steps can be derived from standard expressions replacin expectations over and by and, e..: E-Q: p q pqq9 exp E pa p c lo ppq y, a,, c, ; θq, (6) wit similar expressions for te E-A and E-C steps obtained by permutations of te variables. BRLs and PRLs ave Gaussian distributions as te likeliood and te priors are Gaussian too. Reardin te labels, a furter factorization is needed: p qm pq m q ś p qmj pq mj q. Te labels depend on te BRLs jpj and PRLs as well as on te binary activation states q mj. Te correspondin M-step is iven by: M: θ ar max E pa p c lo ppy a,, c, ; α, l, σ 2 q θpθ ` lo pp ; v q ` lo pp ; v q ` E pa p q lo ppa q; µa, σ a q ` E pc p q lo ppc q; µc, σ c q ` E pq lo ppq; βq j were θ α, l, σ 2, µ a,c, σ a,c, v, v, β (. Given te separability of te prior probability density functions, te M-step can be divided into separate M-steps, as in [7]. Note ere
4 tat a mean-field approximation is used to compute te MRF parameter β. For furter details on tese aspects, te reader is invited to refer to [7]. IV. RESULTS Different data sets ave been analysed to test te performance of tis aloritm. First, some artificial data as been enerated wit te ASL JDE enerative model. Ten, real data acquired on different individuals from te AINSI initiative 2 ave been analysed to validate te proposed approac. A. Artificial data N 288 ASL artificial imaes (i.e. 144 control/ta pairs) ave been simulated usin a realistic low SNR accordin to te observation model in Eq. (1). Different levels of SNR ave been used, in order to sow te performance of te metod dependin on te noise level. To emulate te slow samplin rate of ASL imaes, Eq. (1) was syntesized at t 0.5s and ten down-sampled at a certain time of repetition (TR), wic means tat te temporal resolution of rows (TR) and columns ( t) of X m is different. Here, we considered a fast event-related paradim comprisin two conditions (M 2), wit mean ISI 5s. TR 3s. is considered as a realistic ASL experiment, compared to te TR 1s. tat could be used for a realistic experiment wen usin BOLD sinal. In te experiments, and are enerated as depicted in Fi.2(a)-(b) by dased lines. P is a polynomial basis of order O 4. Drift coefficients and noise realizations were drawn accordin to l j N p0, 10.I O q and b j N p0, 2.I N q, respectively. BRLs were sampled from pa m j qm j 1q N p2.2, 0.3q (for activatin voxels) and from pa m j qm j 0q N p0, 0.3q (for inactivatin voxels). PRLs were enerated wit a lower contrast tan BRLs: pc m j qm j 1q N p1.6, 0.3q and pc m j qm j 0q N p0, 0.3q. PRLs and BRLs were cosen so as to make tis syntetic settin realistic: PRLs lower tan BRLs, and activatin/non-activatin voxels distribution means close. Activation states (assinment variables Q) were set by a and-drawn map. perfusion sinal perfusion sinal (a) BOLD sinal BOLD sinal (b) Fi. 2. Artificial data: Ground-trut response curves (black dased lines) and estimated perfusion (a) and emodynamic (b) response functions. Te first and second rows correspond to te response functions for TR 3s and noise variances 1 (i.e. SNR 2.4 db, top) and 6 (SNR 0.5 db, bottom). simulated T R 1 s. T R 3 s. (a) (b) (c) Fi. 3. Results on artificial data for response levels. Top row: roundtrut maps. Center and bottom rows: estimated maps for T R 1 s and T R 3 s. (a) Hemodynamic response levels. (b) Perfusion response levels. (c) Assinment variables (activation states). noise variance Fi. 1. RMSE of te response functions BRF and PRF estimated in experiments enerated wit different SNRs. Here te difference of te errors in te case of TR 1s and TR 3s is sown. 2 ttp://talie.ujf-renoble.fr/ainsi Fi. 1 sows te root mean squared error (RMSE) of te estimated PRFs and HRFs for experiments enerated wit noise variances from 1 to 6, wic corresponds to decreasin te SNR from 2 down to 0.5 db. Curves are depicted for TR 1s and TR 3s, so as to sow te impact of increasin TR on te performance of te metod and tus mimicin part of te increased difficulty wen movin from BOLD to ASL fmri data. For TR 1s, we observed as expected ood HRF and PRF estimates wit an accurate recovery of bot peaks (Fi. 2 top), as well as ood BRLs and PRLs map estimates (Fi. 3 center). In terms of RMSE (Fi. 1), te values are pretty close to eac oter altou slitly larer for te PRF. However, in te case of TR 3s, we recovered a muc better PRF estimate owin to a lare undersoot in te HRF sapes for very low SNR scenarios (Fi. 2 bottom rit). For te worse SNR scenario, we obtained over-smooted sapes, wit peaks displaced compared to te round trut.
5 sinal (a) (c) Fi. 4. Results on real fasl data for a subject of te AINSI database, for te auditory condition (radioloical convention: Left is rit). (a) BRLs, (b) PRLs, (c) response functions (d) reion of interest of te response functions. In (c), te red and blue curves represent te PRF and HRF estimates. As a reference, we depicted te canonical HRF in black dased line. Tis impacts te recovery of te emodynamic and perfusion activation maps too (see Fi. 3 bottom vs center), muc worse in tis case. Fis. 1-2 sow tat despite te dependence on te SNR level, our sape estimates are pretty reliable. B. Real data Real ASL data were recorded durin an experiment desined to map auditory and visual brain functions, wic consisted of N 291 scans lastin T R 3 s, wit T E 18 ms, FoV 192 mm, eac yieldin a 3-D volume composed of 64 ˆ 64 ˆ 22 voxels (resolution of 3 ˆ 3 ˆ 7 mm 3 ). Te tain sceme used was PICORE Q2T, wit T I ms, T I ms. Te paradim was a fast event-related desin (mean ISI 5.1 s) comprisin sixty auditory and visual stimuli. Ward parcellations were computed in auditory and visual cortices for reionally analysin te evoked activity elicited by auditory and visual stimuli, respectively. Te estimation results are depicted in Fi. 4. First, it is wort notin tat as expected we recovered BOLD response levels wit larer manitude compared to tat of te perfusion response levels (Fi. 4(a)-(b)). Interestinly, in response to auditory stimuli, we retrieved bilateral evoked activity in te auditory cortices for te BOLD and perfusion components, altou not exactly in te same parcels. Subsequent tresoldin of individual posterior probability maps sould allow us to derive reliable statistical maps. Besides, we also recovered plausible sapes for BRF and PRF estimates, as illustrated in Fi. 4(c), were we selected a parcel correspondin to i BRL values. As expected, te HRF sape in te primary auditory cortex (left emispere) is well captured by te canonical sape. Te pysioloical prior elps te parameter estimation, in te perfusion component and enforces te temporal precedence of te PRF estimate over te HRF one. (b) (d) V. DISCUSSION We proposed a variational Expectation-Maximization aloritm to address te issue of jointly detectin activity and estimatin emodynamic and perfusion responses from functional ASL data. Te variational approac as te advantae to provide estimations in analytic form for eac variable of interest. It facilitates te inclusion of additional information from pysioloical models and te incorporation of constraints on te responses to favor stability in te estimations. In particular, we considered a pysioloically informed link between normalized emodynamic and perfusion responses so as to compensate te low sinal-to-noise ratio of te perfusion component. Te results obtained on simulated and real data indicate a ood performance of te approac and suest te variational approac as a potential robust, fast and pramatic metod to analyse callenin ASL data. Addressin perfusion quantification is ten te main perspective of tis work as providin a ood estimation of perfusion modulation effects mit translate into more accurate quantification of perfusion. REFERENCES [1] D. Williams, J. Detre, J. Lei, and A. Koretsky, Manetic resonance imain of perfusion usin spin inversion of arterial water, Proceedins of te National Academy of Sciences, vol. 89, no. 1, pp , [2] S. Oawa, D.W. Tank, R. Menon, J.M. Ellermann, S-G. Kim, H. Merkle, and K. Uurbil, Intrinsic sinal canes accompanyin sensory stimulation: functional brain mappin wit manetic resonance imain, Proceedins of te National Academy of Sciences, vol. 89, pp , [3] L. Hernandez-Garcia, H. Jaanian, and D. B. Rowe, Quantitative analysis of arterial spin labelin fmri data usin a eneral linear model, Manetic resonance imain, vol. 28, no. 7, pp , [4] J. A. Mumford, L. Hernandez-Garcia, G. R. Lee, and T. E. Nicols, Estimation efficiency and statistical power in arterial spin labelin fmri, Neuroimae, vol. 33, no. 1, pp , [5] Karl J Friston, P Fletcer, Oliver Joseps, A Holmes, MD Ru, and Robert Turner, Event-related fmri: caracterizin differential responses, Neuroimae, vol. 7, no. 1, pp , [6] T. Vincent, L. Risser, and P. Ciuciu, Spatially adaptive mixture modelin for analysis of witin-subject fmri time series, IEEE Trans. on Medical Imain, vol. 29, no. 4, pp , Apr [7] L. Caari, T. Vincent, F. Forbes, M. Dojat, and P. Ciuciu, Fast joint detection-estimation of evoked brain activity in event-related fmri usin a variational approac, IEEE Trans. on Medical Imain, vol. 32, no. 5, pp , May [8] T. Vincent, J. Warnkin, M. Villien, A. Krainik, P. Ciuciu, and F. Forbes, Bayesian Joint Detection-Estimation of cerebral vasoreactivity from ASL fmri data, in 16t Proc. MICCAI, LNCS Spriner Verla, Naoya, Japan, Sept. 2013, vol. 2, pp [9] A. Frau-Pascual, T. Vincent, J. Sloboda, P. Ciuciu, and F. Forbes, Pysioloically informed Bayesian analysis of ASL fmri data, in Bayesian and rapical Models for Biomedical Imain, pp Spriner, [10] R. B. Buxton, E. C. Won, and Frank. L. R., Dynamics of blood flow and oxyenation canes durin brain activation: te balloon model, Manetic Resonance in Medicine, vol. 39, pp , June [11] K. J. Friston, A. Mecelli, R. Turner, and C. J. Price, Nonlinear responses in fmri: te balloon model, Volterra kernels, and oter emodynamics, Neuroimae, vol. 12, pp , June [12] T. Vincent, F. Forbes, and P. Ciuciu, Bayesian BOLD and perfusion source separation and deconvolution from functional ASL imain, in 38t Proc. IEEE ICASSP, Vancouver, Canada, May 2013, pp [13] R.M. Neal and G. E. Hinton, A view of te EM aloritm tat justifies incremental, sparse, and oter variants, in Learnin in rapical models, pp Spriner, 1998.
Comparison of stochastic and variational solutions to ASL fmri data analysis
Comparison of stochastic and variational solutions to ASL fmri data analysis Aina Frau-Pascual 1,3, Florence Forbes 1, and Philippe Ciuciu 2,3 p1q INRIA, Univ. Grenoble Alpes, LJK, Grenoble, France p2q
More informationComparison of Stochastic and Variational Solutions to ASL fmri Data Analysis
Comparison of Stochastic and Variational Solutions to ASL fmri Data Analysis Aina Frau-Pascual, Florence Forbes, Philippe Ciuciu To cite this version: Aina Frau-Pascual, Florence Forbes, Philippe Ciuciu.
More informationVariational solution to hemodynamic and perfusion response estimation from ASL fmri data
Variational solution to hemodynamic and perfusion response estimation from ASL fmri data Aina Frau-Pascual, Florence Forbes, Philippe Ciuciu June, 2015 1 / 18 BOLD: Qualitative functional MRI Blood Oxygen
More informationarxiv: v1 [stat.ap] 6 Jan 2015
Physiologically informed Bayesian analysis of ASL fmri data Aina Frau-Pascual 1,3, Thomas Vincent 1, Jennifer Sloboda 1, Philippe Ciuciu 2,3, and Florence Forbes 1 p1q INRIA, MISTIS, Grenoble University,
More informationPHYSIOLOGICAL MODELS COMPARISON FOR THE ANALYSIS OF ASL FMRI DATA
PHYSIOLOGICAL MODELS COMPARISON FOR THE ANALYSIS OF ASL FMRI DATA Aina Frau-Pascual p1,2q Florence Forbes p1q Philippe Ciuciu p2,3q p1q INRIA, MISTIS, Grenoble University, LJK, Grenoble, France p2q INRIA,
More informationSecurity Constrained Optimal Power Flow
Security Constrained Optimal Power Flow 1. Introduction and notation Fiure 1 below compares te optimal power flow (OPF wit te security-constrained optimal power flow (SCOPF. Fi. 1 Some comments about tese
More informationBayesian Joint Detection-Estimation of cerebral vasoreactivity from ASL fmri data
Bayesian Joint Detection-Estimation of cerebral vasoreactivity from ASL fmri data Thomas Vincent, Jan Warnking, Marjorie Villien, Alexandre Krainik, Philippe Ciuciu, Florence Forbes To cite this version:
More informationHee Jung * and Mark L. Psiaki. Cornell University, Ithaca, N.Y Abstract
Tests of Manetometer/Sun-Sensor Orbit Determination Usin Flit Data Hee Jun * and Mark L. Psiaki Cornell University, Itaca, N.Y. 4853-750 Abstract A manetometer-based orbit determination batc filter as
More informationA priori error indicator in the transformation method for problems with geometric uncertainties
A priori error indicator in te transformation metod for problems wit geometric uncertainties Duc Hung Mac, Stépane Clenet, Jean-Claude Mipo, Igor Tsukerman To cite tis version: Duc Hung Mac, Stépane Clenet,
More informationRegularized Regression
Regularized Regression David M. Blei Columbia University December 5, 205 Modern regression problems are ig dimensional, wic means tat te number of covariates p is large. In practice statisticians regularize
More informationNon Linear Observation Equation For Motion Estimation
Non Linear Observation Equation For Motion Estimation Dominique Béréziat, Isabelle Herlin To cite this version: Dominique Béréziat, Isabelle Herlin. Non Linear Observation Equation For Motion Estimation.
More informationMathematical Theory of the h- and g-index in Case of Fractional Counting of Authorship
Matematical Teory of te - and -Index in Case of Fractional Countin of Autorsip Leo Ee Universiteit Hasselt (UHasselt), Campus Diepenbeek, Aoralaan, B-3590 Diepenbeek, Belium and Universiteit Antwerpen
More informationIsotropic diffeomorphisms: solutions to a differential system for a deformed random fields study
Isotropic diffeomorphisms: solutions to a differential system for a deformed random fields study Marc Briant, Julie Fournier To cite this version: Marc Briant, Julie Fournier. Isotropic diffeomorphisms:
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationMinimizing D(Q,P) def = Q(h)
Inference Lecture 20: Variational Metods Kevin Murpy 29 November 2004 Inference means computing P( i v), were are te idden variables v are te visible variables. For discrete (eg binary) idden nodes, exact
More informationThe Measurement of the Gravitational Constant g with Kater s Pendulum
e Measurement of te Gravitational Constant wit Kater s Pendulum Abstract A Kater s pendulum is set up to measure te period of oscillation usin a lamppotocell module and a ektronix oscilloscope. Usin repeated
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationUltra low frequency pressure transducer calibration
Ultra low frequency pressure transducer calibration A. Semenov To cite this version: A. Semenov. Ultra low frequency pressure transducer calibration. Journal de Physique IV Colloque, 1994, 04 (C5), pp.c7-251-c7-254.
More informationCMU Fall VLSI CAD
CMU Fall 00 8-760 VLSI CAD [5 pts] HW 5. Out: Tue Nov 7, Due: Tu. Dec 0, in class. (V). Quadratic Placement [5 pts] Consider tis simple netlist wit fixed pins, wic as placeable objects. All te -point wires
More informationQuasi-geostrophic motion
Capter 8 Quasi-eostropic motion Scale analysis for synoptic-scale motions Simplification of te basic equations can be obtained for synoptic scale motions. Consider te Boussinesq system ρ is assumed to
More informationTHE BAYESIAN ABEL BOUND ON THE MEAN SQUARE ERROR
THE BAYESIAN ABEL BOUND ON THE MEAN SQUARE ERROR Alexandre Renaux, Philippe Forster, Pascal Larzabal, Christ Richmond To cite this version: Alexandre Renaux, Philippe Forster, Pascal Larzabal, Christ Richmond
More informationThe ASL signal. Parenchy mal signal. Venous signal. Arterial signal. Input Function (Label) Dispersion: (t e -kt ) Relaxation: (e -t/t1a )
Lecture Goals Other non-bold techniques (T2 weighted, Mn contrast agents, SSFP, Dynamic Diffusion, ASL) Understand Basic Principles in Spin labeling : spin inversion, flow vs. perfusion ASL variations
More informationLearning based super-resolution land cover mapping
earning based super-resolution land cover mapping Feng ing, Yiang Zang, Giles M. Foody IEEE Fellow, Xiaodong Xiuua Zang, Siming Fang, Wenbo Yun Du is work was supported in part by te National Basic Researc
More informationDispersion relation results for VCS at JLab
Dispersion relation results for VCS at JLab G. Laveissiere To cite this version: G. Laveissiere. Dispersion relation results for VCS at JLab. Compton Scattering from Low to High Momentum Transfer, Mar
More informationTwo-step centered spatio-temporal auto-logistic regression model
Two-step centered spatio-temporal auto-logistic regression model Anne Gégout-Petit, Shuxian Li To cite this version: Anne Gégout-Petit, Shuxian Li. Two-step centered spatio-temporal auto-logistic regression
More informationProbabilistic Graphical Models Homework 1: Due January 29, 2014 at 4 pm
Probabilistic Grapical Models 10-708 Homework 1: Due January 29, 2014 at 4 pm Directions. Tis omework assignment covers te material presented in Lectures 1-3. You must complete all four problems to obtain
More informationEaster bracelets for years
Easter bracelets for 5700000 years Denis Roegel To cite this version: Denis Roegel. Easter bracelets for 5700000 years. [Research Report] 2014. HAL Id: hal-01009457 https://hal.inria.fr/hal-01009457
More informationExogenous input estimation in Electronic Power Steering (EPS) systems
Exogenous input estimation in Electronic Power Steering (EPS) systems Valentina Ciarla, Carlos Canudas de Wit, Franck Quaine, Violaine Cahouet To cite this version: Valentina Ciarla, Carlos Canudas de
More informationVibro-acoustic simulation of a car window
Vibro-acoustic simulation of a car window Christophe Barras To cite this version: Christophe Barras. Vibro-acoustic simulation of a car window. Société Française d Acoustique. Acoustics 12, Apr 12, Nantes,
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More informationWidely Linear Estimation with Complex Data
Widely Linear Estimation with Complex Data Bernard Picinbono, Pascal Chevalier To cite this version: Bernard Picinbono, Pascal Chevalier. Widely Linear Estimation with Complex Data. IEEE Transactions on
More informationLarge brain effective network from EEG/MEG data and dmr information
Large brain effective network from EEG/MEG data and dmr information Brahim Belaoucha, Théodore Papadopoulo To cite this version: Brahim Belaoucha, Théodore Papadopoulo Large brain effective network from
More informationQuantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry
Quantum efficiency and metastable lifetime measurements in ruby ( Cr 3+ : Al2O3) via lock-in rate-window photothermal radiometry A. Mandelis, Z. Chen, R. Bleiss To cite this version: A. Mandelis, Z. Chen,
More informationSolution to Sylvester equation associated to linear descriptor systems
Solution to Sylvester equation associated to linear descriptor systems Mohamed Darouach To cite this version: Mohamed Darouach. Solution to Sylvester equation associated to linear descriptor systems. Systems
More informationNonparametric estimation of the average growth curve with general nonstationary error process
Nonparametric estimation of te average growt curve wit general nonstationary error process Karim Benenni, Mustapa Racdi To cite tis version: Karim Benenni, Mustapa Racdi. Nonparametric estimation of te
More informationMethylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma.
Methylation-associated PHOX2B gene silencing is a rare event in human neuroblastoma. Loïc De Pontual, Delphine Trochet, Franck Bourdeaut, Sophie Thomas, Heather Etchevers, Agnes Chompret, Véronique Minard,
More informationFinite Volume for Fusion Simulations
Finite Volume for Fusion Simulations Elise Estibals, Hervé Guillard, Afeintou Sangam To cite this version: Elise Estibals, Hervé Guillard, Afeintou Sangam. Finite Volume for Fusion Simulations. Jorek Meeting
More informationAnalyzing large-scale spike trains data with spatio-temporal constraints
Analyzing large-scale spike trains data with spatio-temporal constraints Hassan Nasser, Olivier Marre, Selim Kraria, Thierry Viéville, Bruno Cessac To cite this version: Hassan Nasser, Olivier Marre, Selim
More informationImpedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing
Impedance Transmission Conditions for the Electric Potential across a Highly Conductive Casing Hélène Barucq, Aralar Erdozain, David Pardo, Victor Péron To cite this version: Hélène Barucq, Aralar Erdozain,
More informationSparse filter models for solving permutation indeterminacy in convolutive blind source separation
Sparse filter models for solving permutation indeterminacy in convolutive blind source separation Prasad Sudhakar, Rémi Gribonval To cite this version: Prasad Sudhakar, Rémi Gribonval. Sparse filter models
More informationOn Newton-Raphson iteration for multiplicative inverses modulo prime powers
On Newton-Raphson iteration for multiplicative inverses modulo prime powers Jean-Guillaume Dumas To cite this version: Jean-Guillaume Dumas. On Newton-Raphson iteration for multiplicative inverses modulo
More informationSensitivity of hybrid filter banks A/D converters to analog realization errors and finite word length
Sensitivity of hybrid filter banks A/D converters to analog realization errors and finite word length Tudor Petrescu, Jacques Oksman To cite this version: Tudor Petrescu, Jacques Oksman. Sensitivity of
More informationNOTES WEEK 14 DAY 2 SCOT ADAMS
NOTES WEEK 14 DAY 2 SCOT ADAMS We igligt tat it s possible to ave two topological spaces and a continuous bijection from te one to te oter wose inverse is not continuous: Let I : r0, 2πq and let C : tpx,
More informationBin Wang, Bernard Brogliato, Vincent Acary, Ahcene Boubakir, Franck Plestan. HAL Id: hal https://hal.archives-ouvertes.
Experimental comparisons between implicit and explicit implementations of discrete-time sliding mode controllers: Towards cattering suppression in output and input signals Bin Wang, Bernard Brogliato,
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationNONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura
NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., georgeta.budura@etc.utt.ro Abstract:
More informationNumerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral. Equation Approach
Numerical Modeling of Eddy Current Nondestructive Evaluation of Ferromagnetic Tubes via an Integral Equation Approach Anastassios Skarlatos, Grégoire Pichenot, Dominique Lesselier, Marc Lambert, Bernard
More informationAdaptive Neural Filters with Fixed Weights
Adaptive Neural Filters wit Fixed Weigts James T. Lo and Justin Nave Department of Matematics and Statistics University of Maryland Baltimore County Baltimore, MD 150, U.S.A. e-mail: jameslo@umbc.edu Abstract
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationTropical Graph Signal Processing
Tropical Graph Signal Processing Vincent Gripon To cite this version: Vincent Gripon. Tropical Graph Signal Processing. 2017. HAL Id: hal-01527695 https://hal.archives-ouvertes.fr/hal-01527695v2
More informationMIMO decorrelation for visible light communication based on angle optimization
MIMO decorrelation for visible ligt communication based on angle optimization Haiyong Zang, and Yijun Zu Citation: AIP Conference Proceedings 80, 09005 (07); View online: ttps://doi.org/0.03/.4977399 View
More informationHardware Operator for Simultaneous Sine and Cosine Evaluation
Hardware Operator for Simultaneous Sine and Cosine Evaluation Arnaud Tisserand To cite this version: Arnaud Tisserand. Hardware Operator for Simultaneous Sine and Cosine Evaluation. ICASSP 6: International
More informationBlock-Processing Approach to Fractional Sample Rate Conversion with Adjustable Timing
Block-Processin Approac to ractional Sample Rate Conversion wit Adjustable Timin Alexandra Grot, Heinz G. Göckler Diital Sinal Processin Group Rur-Universität Bocum Universitätsstr. 5, D-448 Bocum,Germany
More informationNear-Earth Asteroids Orbit Propagation with Gaia Observations
Near-Earth Asteroids Orbit Propagation with Gaia Observations David Bancelin, Daniel Hestroffer, William Thuillot To cite this version: David Bancelin, Daniel Hestroffer, William Thuillot. Near-Earth Asteroids
More informationOverdispersed Variational Autoencoders
Overdispersed Variational Autoencoders Harsil Sa, David Barber and Aleksandar Botev Department of Computer Science, University College London Alan Turing Institute arsil.sa.15@ucl.ac.uk, david.barber@ucl.ac.uk,
More informationLogarithmic functions
Roberto s Notes on Differential Calculus Capter 5: Derivatives of transcendental functions Section Derivatives of Logaritmic functions Wat ou need to know alread: Definition of derivative and all basic
More informationLearning an Adaptive Dictionary Structure for Efficient Image Sparse Coding
Learning an Adaptive Dictionary Structure for Efficient Image Sparse Coding Jérémy Aghaei Mazaheri, Christine Guillemot, Claude Labit To cite this version: Jérémy Aghaei Mazaheri, Christine Guillemot,
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationRESIDUAL 11-YEAR SIGNAL IN COEFFICIENTS OF MAIN FIELD MODELS. Several recently developed main geomagnetic field models, based on both observatory and
RESIDUAL -YEAR SIGNAL IN COEFFICIENTS OF MAIN FIELD MODELS CRISTIANA ŞTEFAN, VENERA DOBRICĂ, CRIŞAN DEMETRESCU Sabba S. Ştefănescu Institute of Geodynamics of te Romanian Academy, 9-, Jean-Louis Calderon
More informationEstimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
International Journal of Clinical Medicine Researc 2016; 3(5): 76-80 ttp://www.aascit.org/journal/ijcmr ISSN: 2375-3838 Estimating Peak Bone Mineral Density in Osteoporosis Diagnosis by Maximum Distribution
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
More informationA new simple recursive algorithm for finding prime numbers using Rosser s theorem
A new simple recursive algorithm for finding prime numbers using Rosser s theorem Rédoane Daoudi To cite this version: Rédoane Daoudi. A new simple recursive algorithm for finding prime numbers using Rosser
More informationFinancial Econometrics Prof. Massimo Guidolin
CLEFIN A.A. 2010/2011 Financial Econometrics Prof. Massimo Guidolin A Quick Review of Basic Estimation Metods 1. Were te OLS World Ends... Consider two time series 1: = { 1 2 } and 1: = { 1 2 }. At tis
More informationAnalysis of Boyer and Moore s MJRTY algorithm
Analysis of Boyer and Moore s MJRTY algorithm Laurent Alonso, Edward M. Reingold To cite this version: Laurent Alonso, Edward M. Reingold. Analysis of Boyer and Moore s MJRTY algorithm. Information Processing
More information1.5 Function Arithmetic
76 Relations and Functions.5 Function Aritmetic In te previous section we used te newly deined unction notation to make sense o epressions suc as ) + 2 and 2) or a iven unction. It would seem natural,
More informationThe sound power output of a monopole source in a cylindrical pipe containing area discontinuities
The sound power output of a monopole source in a cylindrical pipe containing area discontinuities Wenbo Duan, Ray Kirby To cite this version: Wenbo Duan, Ray Kirby. The sound power output of a monopole
More informationModelling temporal structure (in noise and signal)
Modelling temporal structure (in noise and signal) Mark Woolrich, Christian Beckmann*, Salima Makni & Steve Smith FMRIB, Oxford *Imperial/FMRIB temporal noise: modelling temporal autocorrelation temporal
More informationTrench IGBT failure mechanisms evolution with temperature and gate resistance under various short-circuit conditions
Trench IGBT failure mechanisms evolution with temperature and gate resistance under various short-circuit conditions Adel Benmansour, Stephane Azzopardi, Jean-Christophe Martin, Eric Woirgard To cite this
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationEarly detection of thermal contrast in pulsed stimulated thermography
Early detection of thermal contrast in pulsed stimulated thermography J.-C. Krapez, F. Lepoutre, D. Balageas To cite this version: J.-C. Krapez, F. Lepoutre, D. Balageas. Early detection of thermal contrast
More informationOnline Appendix to Estimating Fiscal Limits: The Case of Greece
Online Appendix to Estimating Fiscal Limits: Te Case of Greece Huixin Bi and Nora Traum May 5, Tis appendix includes details of nonlinear numerical solutions, estimation diagnostics and additional results
More informationLecture XVII. Abstract We introduce the concept of directional derivative of a scalar function and discuss its relation with the gradient operator.
Lecture XVII Abstract We introduce te concept of directional derivative of a scalar function and discuss its relation wit te gradient operator. Directional derivative and gradient Te directional derivative
More informationSolution for the Homework 4
Solution for te Homework 4 Problem 6.5: In tis section we computed te single-particle translational partition function, tr, by summing over all definite-energy wavefunctions. An alternative approac, owever,
More informationPosterior Covariance vs. Analysis Error Covariance in Data Assimilation
Posterior Covariance vs. Analysis Error Covariance in Data Assimilation François-Xavier Le Dimet, Victor Shutyaev, Igor Gejadze To cite this version: François-Xavier Le Dimet, Victor Shutyaev, Igor Gejadze.
More informationSmart Bolometer: Toward Monolithic Bolometer with Smart Functions
Smart Bolometer: Toward Monolithic Bolometer with Smart Functions Matthieu Denoual, Gilles Allègre, Patrick Attia, Olivier De Sagazan To cite this version: Matthieu Denoual, Gilles Allègre, Patrick Attia,
More informationUnbiased minimum variance estimation for systems with unknown exogenous inputs
Unbiased minimum variance estimation for systems with unknown exogenous inputs Mohamed Darouach, Michel Zasadzinski To cite this version: Mohamed Darouach, Michel Zasadzinski. Unbiased minimum variance
More informationContinuous Stochastic Processes
Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling
More informationSome approaches to modeling of the effective properties for thermoelastic composites
Some approaches to modeling of the ective properties for thermoelastic composites Anna Nasedkina Andrey Nasedkin Vladimir Remizov To cite this version: Anna Nasedkina Andrey Nasedkin Vladimir Remizov.
More informationMAT244 - Ordinary Di erential Equations - Summer 2016 Assignment 2 Due: July 20, 2016
MAT244 - Ordinary Di erential Equations - Summer 206 Assignment 2 Due: July 20, 206 Full Name: Student #: Last First Indicate wic Tutorial Section you attend by filling in te appropriate circle: Tut 0
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationQuantum Numbers and Rules
OpenStax-CNX module: m42614 1 Quantum Numbers and Rules OpenStax College Tis work is produced by OpenStax-CNX and licensed under te Creative Commons Attribution License 3.0 Abstract Dene quantum number.
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationBest linear unbiased prediction when error vector is correlated with other random vectors in the model
Best linear unbiased prediction when error vector is correlated with other random vectors in the model L.R. Schaeffer, C.R. Henderson To cite this version: L.R. Schaeffer, C.R. Henderson. Best linear unbiased
More informationEEG/MEG Inverse Solution Driven by fmri
EEG/MEG Inverse Solution Driven by fmri Yaroslav Halchenko CS @ NJIT 1 Functional Brain Imaging EEG ElectroEncephaloGram MEG MagnetoEncephaloGram fmri Functional Magnetic Resonance Imaging others 2 Functional
More informationCompleteness of the Tree System for Propositional Classical Logic
Completeness of the Tree System for Propositional Classical Logic Shahid Rahman To cite this version: Shahid Rahman. Completeness of the Tree System for Propositional Classical Logic. Licence. France.
More informationDerivatives of Exponentials
mat 0 more on derivatives: day 0 Derivatives of Eponentials Recall tat DEFINITION... An eponential function as te form f () =a, were te base is a real number a > 0. Te domain of an eponential function
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More information3.1 Extreme Values of a Function
.1 Etreme Values of a Function Section.1 Notes Page 1 One application of te derivative is finding minimum and maimum values off a grap. In precalculus we were only able to do tis wit quadratics by find
More informationAccounting for Missing Data in Sparse Wavelet Representation of Observation Error Correlations
Accounting for Missing Data in Sparse Wavelet Representation of Observation Error Correlations Vincent Chabot, Maëlle Nodet, Arthur Vidard To cite this version: Vincent Chabot, Maëlle Nodet, Arthur Vidard.
More informationA non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications
A non-commutative algorithm for multiplying (7 7) matrices using 250 multiplications Alexandre Sedoglavic To cite this version: Alexandre Sedoglavic. A non-commutative algorithm for multiplying (7 7) matrices
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationComplete Quadratic Lyapunov functionals using Bessel-Legendre inequality
Complete Quadratic Lyapunov functionals using Bessel-Legendre inequality Alexandre Seuret Frédéric Gouaisbaut To cite tis version: Alexandre Seuret Frédéric Gouaisbaut Complete Quadratic Lyapunov functionals
More informationDecentralized Interference Channels with Noisy Feedback Possess Pareto Optimal Nash Equilibria
Decentralized Interference Channels with Noisy Feedback Possess Pareto Optimal Nash Equilibria Samir M. Perlaza Ravi Tandon H. Vincent Poor To cite this version: Samir M. Perlaza Ravi Tandon H. Vincent
More informationFast Computation of Moore-Penrose Inverse Matrices
Fast Computation of Moore-Penrose Inverse Matrices Pierre Courrieu To cite this version: Pierre Courrieu. Fast Computation of Moore-Penrose Inverse Matrices. Neural Information Processing - Letters and
More information