A proposed nonparametric mixture density estimation using B-spline functions
|
|
- Karin Ferguson
- 5 years ago
- Views:
Transcription
1 A proposed nonparametric mixture density estimation using B-spine functions Atizez Hadrich a,b, Mourad Zribi a, Afif Masmoudi b a Laboratoire d Informatique Signa et Image de a Côte d Opae (LISIC-EA 4491), ULCO, 50 Rue Ferdinand Buisson BP 719, Caais Cedex, France. b Laboratory of Probabiity and Statistics, Facuty of Sciences of Sfax, Sfax University. B. P Sfax. Tunisia Abstract In this paper, we suppose that a density of probabiity f is expressed as a finite inear combination of second order B-spine functions. Then, we obtain a finite mixture of B-spine. We extend the Expectation Maximization (EM) agorithm in order to estimate the new mixture density. The experiments show that the proposed estimator using B-spine functions can produce a satisfactory estimation of mixture density than gaussian cassica theory. Key words: function. B-spine functions, EM agorithm, Estimator, Mixture density, Probabiity density 1. Introduction In statistics, concerning the density estimation in the construction, there are severa competitive casses of nonparametric estimators, the two most popuar being the Kerne estimators introduced by Rosenbatt [8] and the Orthogona series estimators [8]. Usuay, a Kerne estimator has a greater efficiency than an estimator based on Orthogona series, athough it may be more compex to cacuate and update. In case the density has its support confined to the positive haf ine and may not be continuous at the origin, the Kerne method wi not be so attractive since the estimator may perform poory in the neighborhood of the origin. In such a situation, an estimator based on Orthogona series may be more appropriate. Athough scarcey used, the Orthogona series density estimator bears a striking resembance to the B-spine density estimator [1]. However, the Orthogona series estimator may sometimes ead to a function that is not a pdf. Again, the ony reason for incuding this estimator is that the B-spine estimator actuay has the same form, except that the basis functions are not Orthogona. In this paper, we propose a nonparametric method to estimate a mixture density using the second order B-spine and the EM agorithm [3]. The combination of the B-spine functions with the EM agorithm aows us to define a new agorithm which we denote as EM Generaized B-spine (EMGB) agorithm. The paper is organized as foows: In section 2, we describe the nonparametric density estimation methods as we as the B-spine functions. Convergence properties of the iterative Maximum Likeihood Estimation (MLE) are given in section 3. The performance comparison based on the Mean Squared Errors (MSE) is presented in section 4 in order to show the accuracy of the proposed estimator. Finay, the concusion appears in section 5. Emai addresses: hadrich atizez@isic.univ-ittora.fr (Atizez Hadrich a,b ), Mourad.Zribi@isic.univ-ittora.fr (Mourad Zribi a ), Afif.Masmoudi@fss.rnu.tn (Afif Masmoudi b ) URL: (Atizez Hadrich a,b ), (Atizez Hadrich a,b ) Preprint submitted to Esevier Juy 13, 2016
2 2. Proposed estimation methods by using B-spine functions. In the mathematica subfied of numerica anaysis, a B-spine is a spine function that has a minimum support with respect to a given degree, smoothness and domain partition. It s we known that every spine density function can be represented as a finite inear combination of B-spine [1]. The term B-spine stands for basis spines according to Isaac Jacob Schoenberg [2]. A B-spine nonparametric density estimator with uniformy spaced knots convenient for arge data sets was discussed by Gehringer [4]. Curry and Schoenberg (1966) [2] have proved that every spine function S of degree d (d 1, 2,...) with m knots (m 1, 2,...) has a unique expansion S(x) m+d b h B d (x), for a < x < b and h B d (x)dx, (2.1) where a, b R and b s are unknown parameters that need to be estimated. Note that b 0 and m+d b 1 is a specia requirement when using B-spines in order to estimate probabiity density functions. The B-spine of d degree are defined recursivey by B d (x) x x B d 1 (x) + x +d+1 x B d 1 +1 (x) (2.2) x +d x x +d+1 x +1 where B 0 (x) { 1, if x [x, x +1 ) 0, esewhere. For genera spines, this specia requirement is not true. The case d 1 corresponds to a piecewise inear approximation which is attractivey simpe but produces a visibe roughness, uness the knots are cose to each other. In more technica terms, a spine function S of degree d with m interior knots, a x 1 < x 2 <... < x b is a (d 1) continuousy differentiabe function, such that S P d [2], (P d is the cass of poynomias of a maximum degree d in each of the intervas (a, x 2 ), (x 2, x 3 ),..., (x m+1, b)). When we consider the approximation by the B-spines of d degree, the density function f is supposed to be (d 1) continuousy differentiabe in each of the intervas (a, x 2 ), (x 2, x 3 ),..., (x m+1, b). Hence, we have chosen the quadratic B-spines functions which are ony continuousy differentiabe, that is to say, having a minima requirement, as compared with the B-spines of degree d > 2. In our work, we have used the second order B-spines functions (B 2)... defined by (x x ) 2 (x +1 x )(x +2 x ), x [x, x +1 ) 1 B 2 [ (x x )(x +2 x) + (x x +1)(x +3 x) ], x [x (x) +1, x +2 ) x +2 x +1 x +2 x x +3 x +1 (x +3 x) 2 (x +3 x +1 )(x +3 x +2 ), x [x +2, x +3 ) 0, esewhere, where is an integer as usua. We notice that 0 B 2 (x) 1 is aways verified. The support of each spine covers three intervas. It is a quadratic poynomia on each support interva. Note that the peak vaue of B 2(x) is ess than 1. If a the B2 spines coud be potted in any interva, there woud be contributions from the three spines. Note that any pdf f can be approximated by the foowing mixture (2.3) of B-spines: f(x) 2 b B 2,h (x), (2.3)
3 + x+3 with b 0, b 1, B,h 2 (x) B2 (x) and h B 2 (x)dx B 2 (x)dx. The estimation h x of the density f reduces the estimation of the finite-dimensiona parameters (b 1,..., b ) that characterize f Estimation of b by the EM agorithm Assuming independence between the observations X 1, X 2,..., X n, we define the Log-ikeihood function by (X 1,..., X n b) Log( b B,h 2 (X i )). To obtain the MLE of (b 1, b 2,..., b ), we appy the EM agorithm [3]. We initiaize the B-spine coefficients by b (0) 1 B,h 2 n (X j ). E-Step:,j b(p) B,h 2 (X j ) b (p) s Bs,h 2 (X j ) s1 ; 1,..., m + 2, j 1,..., n where,j is the conditiona probabiity at the pth iteration. M-Step: The coefficients of B-spine at (p + 1)th iteration are given by b (p+1) 1 n,j ; 1,..., m + 2 After some iterations of the EM agorithm, we obtain the mixture B-spine estimator of f as f(x) b B,h 2 (x) where b is the maximum ikeihood of b. In what foows, we suggest to extend the second order B-spine estimator to the probem of the mixture density estimation [1]. The combination of the B-spine functions with the EM agorithm aows us to define a new agorithm denoted as EM Generaized B-spine (EMGB) agorithm Description of the EMGB agorithm Let s consider the foowing mixture density of mixture of B-spine f(x) π k f(x b.,k ) π k b,k B 2,h (x), where π k 0 and π k 1. Θ { θ (π 1,..., π K, b.,1,..., b.,k ); π k 0, and f(x b.,k ) π k 1; k 1,..., K; b.,k [0, 1] ; b,k 1 } b,k B 2,h (x). Let X 1, X 2,..., X n be n random variabes with an unknown common pdf f. For each X i, we consider a random variabe Z i (Z 1,i,..., Z K,i ) such that Z i foows a 3
4 mutinomia distribution Mu(1, π 1,..., π K ). Let X (X 1,..., X n ) and Z (Z 1,..., Z n ), the og maximum ikeihood is given by (X, Z θ) [ Z Log(π k f(x i b.,k ))]. We define the mean of the og ikeihood Q(θ, θ (p) ) E((X, Z θ) θ (p) ) [ E(Z X, θ (p) )Log(π k f(x i b.,k ))]. If E(Z X, θ (p) ), then Q Q(θ, θ (p) ) [ Log(π kf(x i b.,k ))] K [ Log(π kf(x i b.,k )) + 1,i Log((1 π k )f(x i b.,1 ))]. k2 k2 The EMGB agorithm seeks to find the MLE by appying iterativey the foowing two steps: E-Step: We estimate a posterior probabiity beonging to the cass k at the pth iteration: E(Z X i, θ (p) ) π (p) k f(x i b (p).,k ). π (p) k f(x i b (p).,k ) M-Step: we cacuate the parameter θ (p+1) that maximizes θ (p+1) arg max θ Θ On the one hand, Q π k which impies π k k2 k2 ( π k n 1,i 1 (1 1,i k2 k2 π k ) 0. This means that πk π k ). Therefore, π (p+1) k Q(θ, θ (p) ). (1 1,i 1 n π k ) Q On the second hand, (B,h 2 (X i ) B1,h 2 (X i )) 0; 2,..., m + 2 b,k f(x i b.,k ) and 2 Q n b 2 (B,h 2 (X i ) B1,h 2 (X i )) 2; 2,..., m + 2.,k f(x i b.,k ) In order to sove Q 0, we use the Newton Raphson method [5] or the steepest descent method, b,k k2. 4
5 we obtain b (p+1),k b (p),k + (B2,h (X i ) B1,h 2 (X i ) f(x i b (p).,k ) ) (B2,h (X i ) B1,h 2 (X i ) f(x i b (p).,k ) ) 2 for 2,..., m + 2 andk 1,..., K. (2.4) The cacuation cyces are made from one step to another unti we achieve the convergence. As starting vaues for the EMGB agorithm, we take the initia vaues of π k and b,k which have to be equa respectivey to n k n and 1 n k B,h 2 n (X j ); where n k is the tota number of observations in the k kth cass. 3. Convergence of the proposed estimator In this section, we wish to show first that the EMGB iterations converge to a vaue ( π k, b,k ), second that this vaue ( π k, b,k ) indeed satisfies the ikeihood equations. Proposition 3.1. Let θ (π 1,..., π K, b.,1,..., b.,k ) Θ. If θ (p+1) converges to θ, as p +, then θ satisfies the ikeihood equations. Proof. On the one hand, for a fixed k, the derivative of the Log-ikeihood function with respect to b,k is given by (X 1,..., X n θ) (B,h 2 π (X j ) B1,h 2 (X j )) k ; 2,..., m + 2. where b,k f(x j ) f(x j ) and f k (X j ) b,k B,h 2 (X j ). By etting p + in (2.3), by using the fact that b (p),k b,k and k,j τ k,j π f k k (X j ), one has τ k,j (B 2,h (X j ) B 2 1,h (X j )) f k (X j ) 0; 2,...,. Therefore 0. And, b,k satisfies the derivation of the ikeihood equations. π k (B 2,h (X j ) B 2 1,h (X j )) f(x j ) 5
6 On the other hand, π k 1 n 0 n Let f k (X j ). This impies that π 1 f1 (X j ) + π 2 f2 (X j ) π K 1 fk 1 (X j ) + (1 π 1 π 2... π K 1 ) f K (X j ) f k (X j ) π 1 ( f 1 (X j ) f K (X j )) π K 1 ( f K 1 (X j ) f K (X j )) + f K (X j ) f k (X j ). H (3.5) π 1 ( f 1 (X j ) f K (X j )) π K 1 ( f K 1 (X j ) f K (X j )). (3.6) Then, mutipying (3.4) by π k for k 1,..., K 1 and by summing the resuts of the (K 1) successive mutipications, we obtain ( π 1 + π π K 1 )H H. Since π k 1, then (1 π K )H H and π K H 0. Note that we chose π K arbitrariy in (3.5). We coud have chosen any one of the π k, k 1, 2,..., K. So, we can say that π k H 0, k 1,..., K. Necessariy, H 0. From (3.4), this impies (X 1,..., X n θ) π k f k (X j ) f K (X j ) 0 for k 1,..., K 1. Therefore, π 1, π 2,..., π K of the EMGB satisfy the ikeihood equations. 4. Performance comparison In our study we simuate n 1000 observations according to different distributions (Norma (N), Beta (B) and Gamma (G)) with given true parameters. We cacuate the density estimator by using Histogram, Kerne, Orthogona and B-spine methods. Then, we compute the Mean Squared Errors MSE 1 (f(x j ) n f(x j )) 2 for each method. Among the different methods, Tabe 1 shows that the B-spine method is the best one by giving the owest MSE. This finding is vaid for a unique density as we as a mixture density. In Tabe 2, we have first computed the MSE between the empirica distribution and the estimated mixture density by using the Gaussian EM agorithm. Second, we have computed the MSE between the empirica distribution and the estimated mixture density by using the EMBG agorithm. This work was done for two images (Lena and Boat). We notice that the EMGB method gives a much ower MSE for both images.in Figure 1, in the case of mixture densities, the B-spine method remains the best one in terms of being cose to the rea density. 6
7 Tabe 1: MSE of different estimations methods. True mode Methods Histogram Kerne Orthogona B-spine N(0, 1) B(3, 5) G(2, 3) Mixed of N(-2, 2), G(3, 1) and B(11, 5) Tabe 2: MSE of mixture distribution obtained by EM and EMGB. Methods Images MSE of Lena MSE of Boat Gaussian EM EMGB Figure 1: Estimation of mixture density using both cassica EM agorithm and EMGB agorithm. 5. Concusion In this paper, we have introduced a new nonparametric B-spine estimator for mixture distributions. Many resuts presented show that the estimation of pdf density by using the proposed estimator is better than those of other methods. This comparison is obtained by the computation of the Mean Squared Errors between the estimated density and the empirica distribution. References [1] Atigan, Bozdogan, (1992). Convergence properties of MLE s and asymptotic simutaneous confidence intervas in fitting cardina B-spine for density estimation. Statistics and Probabiity Letters, Vo. 13, pp , North-Hoand. [2] Curry H.B, Schoenberg I.J., (1966). On poy freqency functions IV: The fundamenta spine functions and their imits. J.Ana. Math, Vo. 17, pp [3] Dempster A.P., Laird N.M., Rubin D.B, (1977). Maximum Likeihood from Incompete Data via the E.M. Agorithm, Journa of the Roya Statistica Society, Ser. B, Vo. 39, pp [4] Gehringer, K.R. and Redner, R.A., (1992), Nonparametric density estimation using tensor product spines. Communications in Statistics-Simuation and Computation, Vo. 21, pp
8 [5] Hazewinke, Michie, ed. (2001), Newton method, Encycopedia of Mathematics, Springer, ISBN [6] Kim J. and a., (2005). A Nonparametric Statistica Method for image segmentation using information theory and curve evoution. IEEE transactions on image processing, Vo. 14, pp [7] Peter Ha, (1983). Orthogona series distribution function estimation, with appication. J, R, Statist, Vo 48, pp [8] Rosenbatt M., (1956). Remarks on some nonparametric estimates of a density function. Ann. Math.Statist, Vo. 27, pp
The EM Algorithm applied to determining new limit points of Mahler measures
Contro and Cybernetics vo. 39 (2010) No. 4 The EM Agorithm appied to determining new imit points of Maher measures by Souad E Otmani, Georges Rhin and Jean-Marc Sac-Épée Université Pau Veraine-Metz, LMAM,
More informationA Brief Introduction to Markov Chains and Hidden Markov Models
A Brief Introduction to Markov Chains and Hidden Markov Modes Aen B MacKenzie Notes for December 1, 3, &8, 2015 Discrete-Time Markov Chains You may reca that when we first introduced random processes,
More informationSTA 216 Project: Spline Approach to Discrete Survival Analysis
: Spine Approach to Discrete Surviva Anaysis November 4, 005 1 Introduction Athough continuous surviva anaysis differs much from the discrete surviva anaysis, there is certain ink between the two modeing
More informationCS229 Lecture notes. Andrew Ng
CS229 Lecture notes Andrew Ng Part IX The EM agorithm In the previous set of notes, we taked about the EM agorithm as appied to fitting a mixture of Gaussians. In this set of notes, we give a broader view
More informationMelodic contour estimation with B-spline models using a MDL criterion
Meodic contour estimation with B-spine modes using a MDL criterion Damien Loive, Ney Barbot, Oivier Boeffard IRISA / University of Rennes 1 - ENSSAT 6 rue de Kerampont, B.P. 80518, F-305 Lannion Cedex
More informationFRST Multivariate Statistics. Multivariate Discriminant Analysis (MDA)
1 FRST 531 -- Mutivariate Statistics Mutivariate Discriminant Anaysis (MDA) Purpose: 1. To predict which group (Y) an observation beongs to based on the characteristics of p predictor (X) variabes, using
More information6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17. Solution 7
6.434J/16.391J Statistics for Engineers and Scientists May 4 MIT, Spring 2006 Handout #17 Soution 7 Probem 1: Generating Random Variabes Each part of this probem requires impementation in MATLAB. For the
More information4 1-D Boundary Value Problems Heat Equation
4 -D Boundary Vaue Probems Heat Equation The main purpose of this chapter is to study boundary vaue probems for the heat equation on a finite rod a x b. u t (x, t = ku xx (x, t, a < x < b, t > u(x, = ϕ(x
More informationBayesian Learning. You hear a which which could equally be Thanks or Tanks, which would you go with?
Bayesian Learning A powerfu and growing approach in machine earning We use it in our own decision making a the time You hear a which which coud equay be Thanks or Tanks, which woud you go with? Combine
More information4 Separation of Variables
4 Separation of Variabes In this chapter we describe a cassica technique for constructing forma soutions to inear boundary vaue probems. The soution of three cassica (paraboic, hyperboic and eiptic) PDE
More informationSeparation of Variables and a Spherical Shell with Surface Charge
Separation of Variabes and a Spherica She with Surface Charge In cass we worked out the eectrostatic potentia due to a spherica she of radius R with a surface charge density σθ = σ cos θ. This cacuation
More informationExpectation-Maximization for Estimating Parameters for a Mixture of Poissons
Expectation-Maximization for Estimating Parameters for a Mixture of Poissons Brandon Maone Department of Computer Science University of Hesini February 18, 2014 Abstract This document derives, in excrutiating
More informationConverting Z-number to Fuzzy Number using. Fuzzy Expected Value
ISSN 1746-7659, Engand, UK Journa of Information and Computing Science Vo. 1, No. 4, 017, pp.91-303 Converting Z-number to Fuzzy Number using Fuzzy Expected Vaue Mahdieh Akhbari * Department of Industria
More informationConvergence Property of the Iri-Imai Algorithm for Some Smooth Convex Programming Problems
Convergence Property of the Iri-Imai Agorithm for Some Smooth Convex Programming Probems S. Zhang Communicated by Z.Q. Luo Assistant Professor, Department of Econometrics, University of Groningen, Groningen,
More informationWeek 6 Lectures, Math 6451, Tanveer
Fourier Series Week 6 Lectures, Math 645, Tanveer In the context of separation of variabe to find soutions of PDEs, we encountered or and in other cases f(x = f(x = a 0 + f(x = a 0 + b n sin nπx { a n
More informationDo Schools Matter for High Math Achievement? Evidence from the American Mathematics Competitions Glenn Ellison and Ashley Swanson Online Appendix
VOL. NO. DO SCHOOLS MATTER FOR HIGH MATH ACHIEVEMENT? 43 Do Schoos Matter for High Math Achievement? Evidence from the American Mathematics Competitions Genn Eison and Ashey Swanson Onine Appendix Appendix
More informationIterative Decoding Performance Bounds for LDPC Codes on Noisy Channels
Iterative Decoding Performance Bounds for LDPC Codes on Noisy Channes arxiv:cs/060700v1 [cs.it] 6 Ju 006 Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department University
More informationA GENERALIZED SKEW LOGISTIC DISTRIBUTION
REVSTAT Statistica Journa Voume 11, Number 3, November 013, 317 338 A GENERALIZED SKEW LOGISTIC DISTRIBUTION Authors: A. Asgharzadeh Department of Statistics, University of Mazandaran Babosar, Iran a.asgharzadeh@umz.ac.ir
More informationA. Distribution of the test statistic
A. Distribution of the test statistic In the sequentia test, we first compute the test statistic from a mini-batch of size m. If a decision cannot be made with this statistic, we keep increasing the mini-batch
More informationSmoothers for ecient multigrid methods in IGA
Smoothers for ecient mutigrid methods in IGA Cemens Hofreither, Stefan Takacs, Water Zuehner DD23, Juy 2015 supported by The work was funded by the Austrian Science Fund (FWF): NFN S117 (rst and third
More informationDIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM
DIGITAL FILTER DESIGN OF IIR FILTERS USING REAL VALUED GENETIC ALGORITHM MIKAEL NILSSON, MATTIAS DAHL AND INGVAR CLAESSON Bekinge Institute of Technoogy Department of Teecommunications and Signa Processing
More informationMath 124B January 17, 2012
Math 124B January 17, 212 Viktor Grigoryan 3 Fu Fourier series We saw in previous ectures how the Dirichet and Neumann boundary conditions ead to respectivey sine and cosine Fourier series of the initia
More informationAnalysis of rounded data in mixture normal model
Stat Papers (2012) 53:895 914 DOI 10.1007/s00362-011-0395-0 REGULAR ARTICLE Anaysis of rounded data in mixture norma mode Ningning Zhao Zhidong Bai Received: 13 August 2010 / Revised: 9 June 2011 / Pubished
More informationXSAT of linear CNF formulas
XSAT of inear CN formuas Bernd R. Schuh Dr. Bernd Schuh, D-50968 Kön, Germany; bernd.schuh@netcoogne.de eywords: compexity, XSAT, exact inear formua, -reguarity, -uniformity, NPcompeteness Abstract. Open
More informationSTABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM 1. INTRODUCTION
Journa of Sound and Vibration (996) 98(5), 643 65 STABILITY OF A PARAMETRICALLY EXCITED DAMPED INVERTED PENDULUM G. ERDOS AND T. SINGH Department of Mechanica and Aerospace Engineering, SUNY at Buffao,
More informationAutomobile Prices in Market Equilibrium. Berry, Pakes and Levinsohn
Automobie Prices in Market Equiibrium Berry, Pakes and Levinsohn Empirica Anaysis of demand and suppy in a differentiated products market: equiibrium in the U.S. automobie market. Oigopoistic Differentiated
More informationResearch Article On the Lower Bound for the Number of Real Roots of a Random Algebraic Equation
Appied Mathematics and Stochastic Anaysis Voume 007, Artice ID 74191, 8 pages doi:10.1155/007/74191 Research Artice On the Lower Bound for the Number of Rea Roots of a Random Agebraic Equation Takashi
More informationA GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS
A GENERAL METHOD FOR EVALUATING OUTAGE PROBABILITIES USING PADÉ APPROXIMATIONS Jack W. Stokes, Microsoft Corporation One Microsoft Way, Redmond, WA 9852, jstokes@microsoft.com James A. Ritcey, University
More informationSome Measures for Asymmetry of Distributions
Some Measures for Asymmetry of Distributions Georgi N. Boshnakov First version: 31 January 2006 Research Report No. 5, 2006, Probabiity and Statistics Group Schoo of Mathematics, The University of Manchester
More informationWAVELET LINEAR ESTIMATION FOR DERIVATIVES OF A DENSITY FROM OBSERVATIONS OF MIXTURES WITH VARYING MIXING PROPORTIONS. B. L. S.
Indian J. Pure App. Math., 41(1): 275-291, February 2010 c Indian Nationa Science Academy WAVELET LINEAR ESTIMATION FOR DERIVATIVES OF A DENSITY FROM OBSERVATIONS OF MIXTURES WITH VARYING MIXING PROPORTIONS
More informationTHINKING IN PYRAMIDS
ECS 178 Course Notes THINKING IN PYRAMIDS Kenneth I. Joy Institute for Data Anaysis and Visuaization Department of Computer Science University of Caifornia, Davis Overview It is frequenty usefu to think
More informationAlgorithms to solve massively under-defined systems of multivariate quadratic equations
Agorithms to sove massivey under-defined systems of mutivariate quadratic equations Yasufumi Hashimoto Abstract It is we known that the probem to sove a set of randomy chosen mutivariate quadratic equations
More informationA Novel Learning Method for Elman Neural Network Using Local Search
Neura Information Processing Letters and Reviews Vo. 11, No. 8, August 2007 LETTER A Nove Learning Method for Eman Neura Networ Using Loca Search Facuty of Engineering, Toyama University, Gofuu 3190 Toyama
More informationExplicit overall risk minimization transductive bound
1 Expicit overa risk minimization transductive bound Sergio Decherchi, Paoo Gastado, Sandro Ridea, Rodofo Zunino Dept. of Biophysica and Eectronic Engineering (DIBE), Genoa University Via Opera Pia 11a,
More informationFORECASTING TELECOMMUNICATIONS DATA WITH AUTOREGRESSIVE INTEGRATED MOVING AVERAGE MODELS
FORECASTING TEECOMMUNICATIONS DATA WITH AUTOREGRESSIVE INTEGRATED MOVING AVERAGE MODES Niesh Subhash naawade a, Mrs. Meenakshi Pawar b a SVERI's Coege of Engineering, Pandharpur. nieshsubhash15@gmai.com
More informationSequential Decoding of Polar Codes with Arbitrary Binary Kernel
Sequentia Decoding of Poar Codes with Arbitrary Binary Kerne Vera Miosavskaya, Peter Trifonov Saint-Petersburg State Poytechnic University Emai: veram,petert}@dcn.icc.spbstu.ru Abstract The probem of efficient
More informationFirst-Order Corrections to Gutzwiller s Trace Formula for Systems with Discrete Symmetries
c 26 Noninear Phenomena in Compex Systems First-Order Corrections to Gutzwier s Trace Formua for Systems with Discrete Symmetries Hoger Cartarius, Jörg Main, and Günter Wunner Institut für Theoretische
More informationData Search Algorithms based on Quantum Walk
Data Search Agorithms based on Quantum Wak Masataka Fujisaki, Hiromi Miyajima, oritaka Shigei Abstract For searching any item in an unsorted database with items, a cassica computer takes O() steps but
More informationA Comparison Study of the Test for Right Censored and Grouped Data
Communications for Statistica Appications and Methods 2015, Vo. 22, No. 4, 313 320 DOI: http://dx.doi.org/10.5351/csam.2015.22.4.313 Print ISSN 2287-7843 / Onine ISSN 2383-4757 A Comparison Study of the
More informationAALBORG UNIVERSITY. The distribution of communication cost for a mobile service scenario. Jesper Møller and Man Lung Yiu. R June 2009
AALBORG UNIVERSITY The distribution of communication cost for a mobie service scenario by Jesper Møer and Man Lung Yiu R-29-11 June 29 Department of Mathematica Sciences Aaborg University Fredrik Bajers
More informationStatistical Learning Theory: A Primer
Internationa Journa of Computer Vision 38(), 9 3, 2000 c 2000 uwer Academic Pubishers. Manufactured in The Netherands. Statistica Learning Theory: A Primer THEODOROS EVGENIOU, MASSIMILIANO PONTIL AND TOMASO
More informationTarget Location Estimation in Wireless Sensor Networks Using Binary Data
Target Location stimation in Wireess Sensor Networks Using Binary Data Ruixin Niu and Pramod K. Varshney Department of ectrica ngineering and Computer Science Link Ha Syracuse University Syracuse, NY 344
More informationImproving the Reliability of a Series-Parallel System Using Modified Weibull Distribution
Internationa Mathematica Forum, Vo. 12, 217, no. 6, 257-269 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/imf.217.611155 Improving the Reiabiity of a Series-Parae System Using Modified Weibu Distribution
More informationStatistical Astronomy
Lectures for the 7 th IAU ISYA Ifrane, nd 3 rd Juy 4 p ( x y, I) p( y x, I) p( x, I) p( y, I) Statistica Astronomy Martin Hendry, Dept of Physics and Astronomy University of Gasgow, UK http://www.astro.ga.ac.uk/users/martin/isya/
More informationMARKOV CHAINS AND MARKOV DECISION THEORY. Contents
MARKOV CHAINS AND MARKOV DECISION THEORY ARINDRIMA DATTA Abstract. In this paper, we begin with a forma introduction to probabiity and expain the concept of random variabes and stochastic processes. After
More informationMultilayer Kerceptron
Mutiayer Kerceptron Zotán Szabó, András Lőrincz Department of Information Systems, Facuty of Informatics Eötvös Loránd University Pázmány Péter sétány 1/C H-1117, Budapest, Hungary e-mai: szzoi@csetehu,
More informationA Solution to the 4-bit Parity Problem with a Single Quaternary Neuron
Neura Information Processing - Letters and Reviews Vo. 5, No. 2, November 2004 LETTER A Soution to the 4-bit Parity Probem with a Singe Quaternary Neuron Tohru Nitta Nationa Institute of Advanced Industria
More informationOn generalized quantum Turing machine and its language classes
Proceedings of the 11th WSEAS Internationa onference on APPLIED MATHEMATIS, Daas, Texas, USA, March -4, 007 51 On generaized quantum Turing machine and its anguage casses SATOSHI IRIYAMA Toyo University
More informationA unified framework for Regularization Networks and Support Vector Machines. Theodoros Evgeniou, Massimiliano Pontil, Tomaso Poggio
MASSACHUSETTS INSTITUTE OF TECHNOLOGY ARTIFICIAL INTELLIGENCE LABORATORY and CENTER FOR BIOLOGICAL AND COMPUTATIONAL LEARNING DEPARTMENT OF BRAIN AND COGNITIVE SCIENCES A.I. Memo No. 1654 March23, 1999
More informationA MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES
6 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES A MODEL FOR ESTIMATING THE LATERAL OVERLAP PROBABILITY OF AIRCRAFT WITH RNP ALERTING CAPABILITY IN PARALLEL RNAV ROUTES Sakae NAGAOKA* *Eectronic
More informationBP neural network-based sports performance prediction model applied research
Avaiabe onine www.jocpr.com Journa of Chemica and Pharmaceutica Research, 204, 6(7:93-936 Research Artice ISSN : 0975-7384 CODEN(USA : JCPRC5 BP neura networ-based sports performance prediction mode appied
More informationLearning Fully Observed Undirected Graphical Models
Learning Fuy Observed Undirected Graphica Modes Sides Credit: Matt Gormey (2016) Kayhan Batmangheich 1 Machine Learning The data inspires the structures we want to predict Inference finds {best structure,
More informationT.C. Banwell, S. Galli. {bct, Telcordia Technologies, Inc., 445 South Street, Morristown, NJ 07960, USA
ON THE SYMMETRY OF THE POWER INE CHANNE T.C. Banwe, S. Gai {bct, sgai}@research.tecordia.com Tecordia Technoogies, Inc., 445 South Street, Morristown, NJ 07960, USA Abstract The indoor power ine network
More informationAn explicit Jordan Decomposition of Companion matrices
An expicit Jordan Decomposition of Companion matrices Fermín S V Bazán Departamento de Matemática CFM UFSC 88040-900 Forianópois SC E-mai: fermin@mtmufscbr S Gratton CERFACS 42 Av Gaspard Coriois 31057
More informationORTHOGONAL MULTI-WAVELETS FROM MATRIX FACTORIZATION
J. Korean Math. Soc. 46 2009, No. 2, pp. 281 294 ORHOGONAL MLI-WAVELES FROM MARIX FACORIZAION Hongying Xiao Abstract. Accuracy of the scaing function is very crucia in waveet theory, or correspondingy,
More informationHigh-order approximations to the Mie series for electromagnetic scattering in three dimensions
Proceedings of the 9th WSEAS Internationa Conference on Appied Mathematics Istanbu Turkey May 27-29 2006 (pp199-204) High-order approximations to the Mie series for eectromagnetic scattering in three dimensions
More informationAkaike Information Criterion for ANOVA Model with a Simple Order Restriction
Akaike Information Criterion for ANOVA Mode with a Simpe Order Restriction Yu Inatsu * Department of Mathematics, Graduate Schoo of Science, Hiroshima University ABSTRACT In this paper, we consider Akaike
More information6 Wave Equation on an Interval: Separation of Variables
6 Wave Equation on an Interva: Separation of Variabes 6.1 Dirichet Boundary Conditions Ref: Strauss, Chapter 4 We now use the separation of variabes technique to study the wave equation on a finite interva.
More informationAppendix of the Paper The Role of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Model
Appendix of the Paper The Roe of No-Arbitrage on Forecasting: Lessons from a Parametric Term Structure Mode Caio Ameida cameida@fgv.br José Vicente jose.vaentim@bcb.gov.br June 008 1 Introduction In this
More informationSchedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness
1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea- Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process Management,
More information14 Separation of Variables Method
14 Separation of Variabes Method Consider, for exampe, the Dirichet probem u t = Du xx < x u(x, ) = f(x) < x < u(, t) = = u(, t) t > Let u(x, t) = T (t)φ(x); now substitute into the equation: dt
More informationAn Algorithm for Pruning Redundant Modules in Min-Max Modular Network
An Agorithm for Pruning Redundant Modues in Min-Max Moduar Network Hui-Cheng Lian and Bao-Liang Lu Department of Computer Science and Engineering, Shanghai Jiao Tong University 1954 Hua Shan Rd., Shanghai
More informationCombining reaction kinetics to the multi-phase Gibbs energy calculation
7 th European Symposium on Computer Aided Process Engineering ESCAPE7 V. Pesu and P.S. Agachi (Editors) 2007 Esevier B.V. A rights reserved. Combining reaction inetics to the muti-phase Gibbs energy cacuation
More informationPhysics 127c: Statistical Mechanics. Fermi Liquid Theory: Collective Modes. Boltzmann Equation. The quasiparticle energy including interactions
Physics 27c: Statistica Mechanics Fermi Liquid Theory: Coective Modes Botzmann Equation The quasipartice energy incuding interactions ε p,σ = ε p + f(p, p ; σ, σ )δn p,σ, () p,σ with ε p ε F + v F (p p
More informationSchedulability Analysis of Deferrable Scheduling Algorithms for Maintaining Real-Time Data Freshness
1 Scheduabiity Anaysis of Deferrabe Scheduing Agorithms for Maintaining Rea-Time Data Freshness Song Han, Deji Chen, Ming Xiong, Kam-yiu Lam, Aoysius K. Mok, Krithi Ramamritham UT Austin, Emerson Process
More informationMaximum likelihood decoding of trellis codes in fading channels with no receiver CSI is a polynomial-complexity problem
1 Maximum ikeihood decoding of treis codes in fading channes with no receiver CSI is a poynomia-compexity probem Chun-Hao Hsu and Achieas Anastasopouos Eectrica Engineering and Computer Science Department
More informationarxiv: v2 [cond-mat.stat-mech] 14 Nov 2008
Random Booean Networks Barbara Drosse Institute of Condensed Matter Physics, Darmstadt University of Technoogy, Hochschustraße 6, 64289 Darmstadt, Germany (Dated: June 27) arxiv:76.335v2 [cond-mat.stat-mech]
More informationBALANCING REGULAR MATRIX PENCILS
BALANCING REGULAR MATRIX PENCILS DAMIEN LEMONNIER AND PAUL VAN DOOREN Abstract. In this paper we present a new diagona baancing technique for reguar matrix pencis λb A, which aims at reducing the sensitivity
More informationPartial permutation decoding for MacDonald codes
Partia permutation decoding for MacDonad codes J.D. Key Department of Mathematics and Appied Mathematics University of the Western Cape 7535 Bevie, South Africa P. Seneviratne Department of Mathematics
More informationTwo-sample inference for normal mean vectors based on monotone missing data
Journa of Mutivariate Anaysis 97 (006 6 76 wwweseviercom/ocate/jmva Two-sampe inference for norma mean vectors based on monotone missing data Jianqi Yu a, K Krishnamoorthy a,, Maruthy K Pannaa b a Department
More informationSVM: Terminology 1(6) SVM: Terminology 2(6)
Andrew Kusiak Inteigent Systems Laboratory 39 Seamans Center he University of Iowa Iowa City, IA 54-57 SVM he maxima margin cassifier is simiar to the perceptron: It aso assumes that the data points are
More informationImproved Min-Sum Decoding of LDPC Codes Using 2-Dimensional Normalization
Improved Min-Sum Decoding of LDPC Codes sing -Dimensiona Normaization Juntan Zhang and Marc Fossorier Department of Eectrica Engineering niversity of Hawaii at Manoa Honouu, HI 968 Emai: juntan, marc@spectra.eng.hawaii.edu
More informationTitle: Crack Diagnostics via Fourier Transform: Real and Imaginary Components vs. Power Spectral Density
Tite: rack Diagnostics via ourier Transform: ea and maginary omponents vs. Power pectra Density Authors: Leonid M. Geman Victor Giurgiutiu van V. Petrunin 3 rd nternationa Workshop on tructura Heath Monitoring
More informationDISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE
DISTRIBUTION OF TEMPERATURE IN A SPATIALLY ONE- DIMENSIONAL OBJECT AS A RESULT OF THE ACTIVE POINT SOURCE Yury Iyushin and Anton Mokeev Saint-Petersburg Mining University, Vasiievsky Isand, 1 st ine, Saint-Petersburg,
More informationIntegrating Factor Methods as Exponential Integrators
Integrating Factor Methods as Exponentia Integrators Borisav V. Minchev Department of Mathematica Science, NTNU, 7491 Trondheim, Norway Borko.Minchev@ii.uib.no Abstract. Recenty a ot of effort has been
More informationGauss Law. 2. Gauss s Law: connects charge and field 3. Applications of Gauss s Law
Gauss Law 1. Review on 1) Couomb s Law (charge and force) 2) Eectric Fied (fied and force) 2. Gauss s Law: connects charge and fied 3. Appications of Gauss s Law Couomb s Law and Eectric Fied Couomb s
More informationEfficiently Generating Random Bits from Finite State Markov Chains
1 Efficienty Generating Random Bits from Finite State Markov Chains Hongchao Zhou and Jehoshua Bruck, Feow, IEEE Abstract The probem of random number generation from an uncorreated random source (of unknown
More informationStatistics for Applications. Chapter 7: Regression 1/43
Statistics for Appications Chapter 7: Regression 1/43 Heuristics of the inear regression (1) Consider a coud of i.i.d. random points (X i,y i ),i =1,...,n : 2/43 Heuristics of the inear regression (2)
More informationEvolutionary Product-Unit Neural Networks for Classification 1
Evoutionary Product-Unit Neura Networs for Cassification F.. Martínez-Estudio, C. Hervás-Martínez, P. A. Gutiérrez Peña A. C. Martínez-Estudio and S. Ventura-Soto Department of Management and Quantitative
More informationReliability: Theory & Applications No.3, September 2006
REDUNDANCY AND RENEWAL OF SERVERS IN OPENED QUEUING NETWORKS G. Sh. Tsitsiashvii M.A. Osipova Vadivosto, Russia 1 An opened queuing networ with a redundancy and a renewa of servers is considered. To cacuate
More informationMath 124B January 31, 2012
Math 124B January 31, 212 Viktor Grigoryan 7 Inhomogeneous boundary vaue probems Having studied the theory of Fourier series, with which we successfuy soved boundary vaue probems for the homogeneous heat
More informationMinimizing Total Weighted Completion Time on Uniform Machines with Unbounded Batch
The Eighth Internationa Symposium on Operations Research and Its Appications (ISORA 09) Zhangiaie, China, September 20 22, 2009 Copyright 2009 ORSC & APORC, pp. 402 408 Minimizing Tota Weighted Competion
More informationSupervised i-vector Modeling - Theory and Applications
Supervised i-vector Modeing - Theory and Appications Shreyas Ramoji, Sriram Ganapathy Learning and Extraction of Acoustic Patterns LEAP) Lab, Eectrica Engineering, Indian Institute of Science, Bengauru,
More informationProblem set 6 The Perron Frobenius theorem.
Probem set 6 The Perron Frobenius theorem. Math 22a4 Oct 2 204, Due Oct.28 In a future probem set I want to discuss some criteria which aow us to concude that that the ground state of a sef-adjoint operator
More informationNEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION
NEW DEVELOPMENT OF OPTIMAL COMPUTING BUDGET ALLOCATION FOR DISCRETE EVENT SIMULATION Hsiao-Chang Chen Dept. of Systems Engineering University of Pennsyvania Phiadephia, PA 904-635, U.S.A. Chun-Hung Chen
More informationMinimum Enclosing Circle of a Set of Fixed Points and a Mobile Point
Minimum Encosing Circe of a Set of Fixed Points and a Mobie Point Aritra Banik 1, Bhaswar B. Bhattacharya 2, and Sandip Das 1 1 Advanced Computing and Microeectronics Unit, Indian Statistica Institute,
More informationarxiv: v1 [math.fa] 23 Aug 2018
An Exact Upper Bound on the L p Lebesgue Constant and The -Rényi Entropy Power Inequaity for Integer Vaued Random Variabes arxiv:808.0773v [math.fa] 3 Aug 08 Peng Xu, Mokshay Madiman, James Mebourne Abstract
More informationStructural Control of Probabilistic Boolean Networks and Its Application to Design of Real-Time Pricing Systems
Preprints of the 9th Word Congress The Internationa Federation of Automatic Contro Structura Contro of Probabiistic Booean Networks and Its Appication to Design of Rea-Time Pricing Systems Koichi Kobayashi
More informationSmoothness equivalence properties of univariate subdivision schemes and their projection analogues
Numerische Mathematik manuscript No. (wi be inserted by the editor) Smoothness equivaence properties of univariate subdivision schemes and their projection anaogues Phiipp Grohs TU Graz Institute of Geometry
More informationAssignment 7 Due Tuessday, March 29, 2016
Math 45 / AMCS 55 Dr. DeTurck Assignment 7 Due Tuessday, March 9, 6 Topics for this week Convergence of Fourier series; Lapace s equation and harmonic functions: basic properties, compuations on rectanges
More informationProvisions estimation for portfolio of CDO in Gaussian financial environment
Technica report, IDE1123, October 27, 2011 Provisions estimation for portfoio of CDO in Gaussian financia environment Master s Thesis in Financia Mathematics Oeg Maximchuk and Yury Vokov Schoo of Information
More informationResearch Article Numerical Range of Two Operators in Semi-Inner Product Spaces
Abstract and Appied Anaysis Voume 01, Artice ID 846396, 13 pages doi:10.1155/01/846396 Research Artice Numerica Range of Two Operators in Semi-Inner Product Spaces N. K. Sahu, 1 C. Nahak, 1 and S. Nanda
More informationNumerical solution of one dimensional contaminant transport equation with variable coefficient (temporal) by using Haar wavelet
Goba Journa of Pure and Appied Mathematics. ISSN 973-1768 Voume 1, Number (16), pp. 183-19 Research India Pubications http://www.ripubication.com Numerica soution of one dimensiona contaminant transport
More informationConsistent linguistic fuzzy preference relation with multi-granular uncertain linguistic information for solving decision making problems
Consistent inguistic fuzzy preference reation with muti-granuar uncertain inguistic information for soving decision making probems Siti mnah Binti Mohd Ridzuan, and Daud Mohamad Citation: IP Conference
More informationPHYS 110B - HW #1 Fall 2005, Solutions by David Pace Equations referenced as Eq. # are from Griffiths Problem statements are paraphrased
PHYS 110B - HW #1 Fa 2005, Soutions by David Pace Equations referenced as Eq. # are from Griffiths Probem statements are paraphrased [1.] Probem 6.8 from Griffiths A ong cyinder has radius R and a magnetization
More informationSupport Vector Machine and Its Application to Regression and Classification
BearWorks Institutiona Repository MSU Graduate Theses Spring 2017 Support Vector Machine and Its Appication to Regression and Cassification Xiaotong Hu As with any inteectua project, the content and views
More informationSVM-based Supervised and Unsupervised Classification Schemes
SVM-based Supervised and Unsupervised Cassification Schemes LUMINITA STATE University of Pitesti Facuty of Mathematics and Computer Science 1 Targu din Vae St., Pitesti 110040 ROMANIA state@cicknet.ro
More informationOn the evaluation of saving-consumption plans
On the evauation of saving-consumption pans Steven Vanduffe Jan Dhaene Marc Goovaerts Juy 13, 2004 Abstract Knowedge of the distribution function of the stochasticay compounded vaue of a series of future
More informationPhysics 235 Chapter 8. Chapter 8 Central-Force Motion
Physics 35 Chapter 8 Chapter 8 Centra-Force Motion In this Chapter we wi use the theory we have discussed in Chapter 6 and 7 and appy it to very important probems in physics, in which we study the motion
More informationApproximated MLC shape matrix decomposition with interleaf collision constraint
Agorithmic Operations Research Vo.4 (29) 49 57 Approximated MLC shape matrix decomposition with intereaf coision constraint Antje Kiese and Thomas Kainowski Institut für Mathematik, Universität Rostock,
More information