Running Time Analysis of the (1+1)-EA for OneMax and LeadingOnes under Bit-wise Noise
|
|
- Camilla Copeland
- 6 years ago
- Views:
Transcription
1 Ruig Time Aalysis of the +-EA for OeMax ad LeadigOes uder Bit-wise Noise Chao Qia Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia Wu Jiag Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia ABSTRACT Previous ruig time aalyses of evolutioary algorithms EAs i oisy eviromets ofte studied the oe-bit oise model, which flips a radomly chose bit of a solutio before evaluatio. I this paper, we study a atural extesio of oe-bit oise, the bit-wise oise model, which idepedetly flips each bit of a solutio with some probability. We aalyze the ruig time of the +-EA solvig OeMax ad LeadigOes uder bit-wise oise for the first time, ad derive the rages of the oise level for polyomial ad superpolyomial ruig time bouds. The aalysis o LeadigOes uder bit-wise oise ca be easily trasferred to oebit oise, ad improves the previously kow results. CCS CONCEPTS Theory of computatio Theory of radomized search heuristics; KEYWORDS Noisy optimizatio, evolutioary algorithms, ruig time aalysis, computatioal complexity ACM Referece format: Chao Qia, Chao Bia, Wu Jiag, ad Ke Tag. 07. Ruig Time Aalysis of the +-EA for OeMax ad LeadigOes uder Bitwise Noise. I Proceedigs of GECCO 7, Berli, Germay, July 5-9, 07, 8 pages. DOI: INTRODUCTION Evolutioary algorithms EAs have bee widely applied to solve real-world optimizatio tasks, where the exact objective i.e., fitess evaluatio of cadidate solutios is ofte Permissio to make digital or hard copies of all or part of this work for persoal or classroom use is grated without fee provided that copies are ot made or distributed for profit or commercial advatage ad that copies bear this otice ad the full citatio o the first page. Copyrights for compoets of this work owed by others tha ACM must be hoored. Abstractig with credit is permitted. To copy otherwise, or republish, to post o servers or to redistribute to lists, requires prior specific permissio ad/or a fee. Request permissios from permissios@acm.org. GECCO 7, Berli, Germay 07 ACM /7/07... $5.00 DOI: Chao Bia Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia biacht@mail.ustc.edu.c Ke Tag Uiversity of Sciece ad Techology of Chia Hefei 3007, Chia ketag@ustc.edu.c impossible, while we ca obtai oly a oisy oe [3, 3]. However, previous theoretical aalyses of EAs maily focused o oise-free optimizatio, where the fitess evaluatio is exact. For ruig time aalysis, a leadig theoretical aspect [, 4], oly a few pieces of work o oisy evolutioary optimizatio have bee reported. Droste [6] first aalyzed the +-EA o the OeMax problem i the presece of oe-bit oise ad showed the maximal oise probability p = log / allowig a polyomial ruig time, where is the problem size. Gieße ad Kötzig [] recetly studied the LeadigOes problem, ad proved that the expected ruig time is polyomial if p /6e ad expoetial if p = /. For iefficiet optimizatio of the +-EA uder high oise levels, some implicit mechaisms of EAs were proved to be robust to oise. I [], it was show that the µ+- EA with a small populatio of size Θlog ca solve Oe- Max i polyomial time eve if the probability of oe-bit oise reaches. The robustess of populatios to oise was also proved i the settig of o-elitist EAs [4, 7]. However, Friedrich et al. [8] showed the limitatio of populatios by provig that the µ+-ea eeds super-polyomial time for solvig OeMax uder additive Gaussia oise N 0, σ with σ 3. This difficulty ca be overcome by the compact geetic algorithm cga [8] ad a simple At Coloy Optimizatio ACO algorithm [9], both of which fid the optimal solutio i polyomial time with a high probability. The ability of explicit oise hadlig strategies was also theoretically studied. Qia et al. [9] proved that the threshold selectio strategy is robust to oise: the expected ruig time of the +-EA usig threshold selectio o OeMax uder oe-bit oise is always polyomial regardless of the oise level. For the +-EA solvig OeMax ad LeadigOes uder oe-bit or additive Gaussia oise, the resamplig strategy was show able to reduce the ruig time from expoetial to polyomial i high oise levels [8]. Akimoto et al. [] also proved that resamplig with a large sample size ca make optimizatio uder additive ubiased oise behave as optimizatio i a oise-free eviromet. The iterplay betwee resamplig ad implicit oise-hadlig mechaisms has bee statistically studied i [0]. The studies metioed above maily cosidered the oebit oise model, which flips a radom bit of a solutio before
2 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. Table : For the ruig time of the +-EA solvig OeMax ad LeadigOes uder prior oise models, the rages of oise parameters for a polyomial upper boud ad a super-polyomial lower boud are show below. +-EA bit-wise oise p, bit-wise oise, q oe-bit oise OeMax Olog /, ωlog / Olog /, ωlog / [] Olog /, ωlog / [6] LeadigOes Olog /, ωlog / Olog / 3, ωlog / /6e, = / []; Olog /, ωlog / evaluatio with probability p. However, the oise model, which ca chage several bits of a solutio simultaeously, may be more realistic ad eeds to be studied, as metioed i the first oisy theoretical work [6]. I this paper, we study the bit-wise oise model, which is characterized by a pair p, q of parameters. It happes with probability p, ad idepedetly flips each bit of a solutio with probability q before evaluatio. We aalyze the ruig time of the +-EA solvig OeMax ad LeadigOes uder bit-wise oise with two specific parameter settigs p, ad, q. The rages of p ad q for a polyomial upper boud ad a super-polyomial lower boud are derived, as show i the middle two colums of Table. For the +-EA o LeadigOes, we also trasfer the ruig time bouds from bit-wise oise p, to oe-bit oise by usig the same proof procedure. As show i the bottom right of Table, our results improve the previously kow oes []. The rest of this paper is orgaized as follows. Sectio itroduces some prelimiaries. The ruig time aalysis o OeMax ad LeadigOes is preseted i Sectios 3 ad 4, respectively. Sectio 5 cocludes the paper. PRELIMINARIES I this sectio, we first itroduce the optimizatio problems, evolutioary algorithms ad oise models studied i this paper, respectively, ad the preset the aalysis tools that we use throughout this paper.. OeMax ad LeadigOes I this paper, we use two well-kow pseudo-boolea fuctios OeMax ad LeadigOes. The OeMax problem as preseted i Defiitio. aims to maximize the umber of -bits of a solutio. The LeadigOes problem as preseted i Defiitio. aims to maximize the umber of cosecutive -bits coutig from the left of a solutio. Their optimal solutio is... briefly deoted as. It has bee show that the expected ruig time of the +-EA o OeMax ad LeadigOes is Θ log ad Θ, respectively [7]. Defiitio. OeMax. The OeMax Problem of size is to fid a bits biary strig x such that x = arg max x {0,} fx = i= xi. Defiitio. LeadigOes. The LeadigOes Problem of size is to fid a bits biary strig x such that x = arg max x {0,} fx = i. i= j= xj. Bit-wise Noise There are maily two kids of oise models: prior ad posterior [, 3]. The prior oise comes from the variatio o a solutio, while the posterior oise comes from the variatio o the fitess of a solutio. Previous theoretical aalyses ofte focused o a specific prior oise model, oe-bit oise. As preseted i Defiitio.3, it flips a radom bit of a solutio before evaluatio with probability p. However, i may realistic applicatios, oise ca chage several bits of a solutio simultaeously rather tha oly oe bit. We thus cosider the bit-wise oise model. As preseted i Defiitio.4, it happes with probability p, ad idepedetly flips each bit of a solutio with probability q before evaluatio. To the best of our kowledge, oly bit-wise oise with p = ad q [0, ] has bee recetly studied. Gieße ad Kötzig [] proved that for the +-EA o OeMax, the expected ruig time is polyomial if q = Olog / ad super-polyomial if q = ωlog /. I this paper, we will study two specific bit-wise oise models: p [0, ] q = ad p = q [0, ], which are briefly deoted as bit-wise oise p, ad bit-wise oise, q, respectively. Defiitio.3 Oe-bit Noise. Give a parameter p [0, ], let f x ad fx deote the oisy ad true fitess of a biary solutio x {0, }, respectively, the { f fx with probability p, x = fx with probability p, where x is geerated by flippig a uiformly radomly chose bit of x. Defiitio.4 Bit-wise Noise. Give parameters p, q [0, ], let f x ad fx deote the oisy ad true fitess of a biary solutio x {0, }, respectively, the { f fx with probability p, x = fx with probability p, where x is geerated by idepedetly flippig each bit of x with probability q..3 +-EA The +-EA as described i Algorithm is studied i this paper. For oisy optimizatio, oly a oisy fitess value f x istead of the exact oe fx ca be accessed, ad thus step 4 of Algorithm chages to be if f x f x. Note that the reevaluatio strategy is used as i [5, 6, ]. That is, besides evaluatig f x, f x will be reevaluated i each iteratio of the +-EA. The ruig time is usually defied as the umber of fitess evaluatios eeded to fid
3 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay a optimal solutio w.r.t. the true fitess fuctio f for the first time [, 6, ]. ALGORITHM +-EA. Give a fuctio f over {0,} to be maximized, it cosists of the followig steps:. x := uiformly radomly selected from {0, }.. Repeat util the termiatio coditio is met 3. x :=flip each bit of x idepedetly with prob. /. 4. if fx fx 5. x := x..4 Aalysis Tools The process of the +-EA solvig OeMax or LeadigOes ca be directly modeled as a Markov chai {ξ t} + t=0. We oly eed to take the solutio space {0, } as the chai s state space i.e., ξ t X = {0, }, ad take the optimal solutio as the chai s optimal state i.e., X = { }. Give a Markov chai {ξ t} + t=0 ad ξˆt = x, we defie its first hittig time FHT as τ = mi{t ξˆt+t X, t 0}. The mathematical expectatio of τ, E[[τ ξˆt = x]] = + i=0 ip τ = i, is called the expected first hittig time EFHT startig from ξˆt = x. If ξ 0 is draw from a distributio π 0, E[[τ ξ 0 π 0]] = x X π0xe[[τ ξ0 = x]] is called the EFHT of the Markov chai over the iitial distributio π 0. Thus, the expected ruig time of the +-EA startig from ξ 0 π 0 is equal to + E[[τ ξ 0 π 0]], where the term correspods to evaluatig the iitial solutio, ad the factor correspods to evaluatig the offsprig solutio x ad reevaluatig the paret solutio x i each iteratio. Note that we cosider the expected ruig time of the +-EA startig from a uiform iitial distributio i this paper. I the followig, we give three drift theorems that will be used to derive the EFHT of Markov chais i the paper. LEMMA.5 ADDITIVE DRIFT []. Give a Markov chai {ξ t} + t=0 ad a distace fuctio V x, if for ay t 0 ad ay ξ t with V ξ t > 0, there exists a real umber c > 0 such that E[[V ξ t V ξ t+ ξ t]] c, the the EFHT satisfies that E[[τ ξ 0]] V ξ 0/c. LEMMA.6 SIMPLIFIED DRIFT [5, 6]. Let X t, t 0, be real-valued radom variables describig a stochastic process. Suppose there exists a iterval [a, b] R, two costats δ, ɛ > 0 ad, possibly depedig o l := b a, a fuctio rl satisfyig rl = ol/ logl such that for all t 0 the followig two coditios hold:. E[[X t X t+ a < X t < b]] ɛ,. P X t+ X t j X t > a rl + δ j for j N 0. The there is a costat c > 0 such that for T := mi{t 0 : X t a X 0 b} it holds P T cl/rl = Ωl/rl. LEMMA.7 SIMPLIFIED DRIFT WITH SELF-LOOPS [0]. Let X t, t 0, be real-valued radom variables describig a stochastic process. Suppose there exists a iterval [a, b] R, two costats δ, ɛ > 0 ad, possibly depedig o l := b a, a fuctio rl satisfyig rl = ol/ logl such that for all t 0 the followig two coditios hold:. a<i<b : E[[X t X t+ X t =i]] ɛ P X t+ i X t =i,. i > a, j N 0 : P X t+ X t j X t =i rl P Xt+ i Xt =i. +δ j The there is a costat c > 0 such that for T := mi{t 0 : X t a X 0 b} it holds P T cl/rl = Ωl/rl. 3 THE ONEMAX PROBLEM I this sectio, we aalyze the ruig time of the +- EA o OeMax uder bit-wise oise. Note that for bit-wise oise, q, it has bee proved that the maximal value of q allowig a polyomial ruig time is log /, as show i Theorem 3.. THEOREM 3.. [] For the +-EA o OeMax uder bitwise oise, q, the expected ruig time is polyomial if q = Olog / ad super-polyomial if q = ωlog /. For bit-wise oise p,, we prove i Theorems 3.4 ad 3.5 that the maximum value of p allowig a polyomial ruig time is log /. Istead of usig the origial drift theorems, we apply the upper ad lower bouds for the +-EA o stochastic OeMax i []. Let x k deote ay solutio with k umber of -bits, ad f x k deote its oisy objective value, which is a radom variable. Lemma 3. ituitively meas that if the probability of recogizig the true better solutio by oisy evaluatio is large, the ruig time ca be polyomially upper bouded. O the cotrary, Lemma 3.3 shows that if the probability of makig a right compariso is small, the ruig time ca be expoetially lower bouded. Both of them are proved by applyig stadard drift theorems, ad ca be used to simplify our aalysis. LEMMA 3.. [] Suppose there is a positive costat c /5 ad some < l / such that k < : P f x k < f x k+ l ; k < l : P f x k < f x k+ c k, the the +-EA optimizes f i expectatio i O log + Ol iteratios. LEMMA 3.3. [] Suppose there is some l /4 ad a costat c 6 such that l k < : P f x k < f x k+ c k, the the +-EA optimizes f i Ωl iteratios with a high probability. THEOREM 3.4. For the +-EA o OeMax uder bit-wise oise p,, the expected ruig time is polyomial if p = Olog /. PROOF. We prove it by usig Lemma 3.. For ay positive costat b, suppose that p b log /. We set the two parameters i Lemma 3. as c = mi{ b log, b} ad l =, ]. 5 c
4 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. For ay k <, f x k f x k+ implies that f x k k + or f x k+ k, either of which happes with probability at most p. By the uio boud, we get P f x k f x k+ p For ay k < l, we easily get b log = lc l. P f x k f x k+ lc < c k. By Lemma 3., we kow that the expected ruig time is O log + Ob log /c, i.e., polyomial. THEOREM 3.5. For the +-EA o OeMax uder bit-wise oise p,, the expected ruig time is super-polyomial if p = ωlog / ωlog / ad expoetial if p = Olog /. PROOF. We use Lemma 3.3 to prove it. Let c = 6. The case p = ωlog / ωlog / is first aalyzed. For ay positive costat b, let l = b log. For ay k l, we get P f x k f x k+ P f x k = k P f x k+ k. To make f x k = k, it is sufficiet that the oise does ot happe, i.e., P f x k = k p. To make f x k+ k, it is sufficiet to flip oe -bit ad keep other bits uchaged by oise, i.e., P f x k+ k p k+. Thus, P f x k f x k+ p p k + e Sice c k c l = ωlog /. = cb log, the coditio of Lemma 3.3 holds. Thus, the expected ruig time is Ωb log where b is ay costat, i.e., super-polyomial. For the case p = Olog /, let l =. We use aother lower boud p for P f x k = k, sice it is sufficiet that o bit flips by oise. Thus, we have P f x k f x k+ p p k + e = Ω. Sice c k c, the coditio of Lemma 3.3 holds. Thus, the expected ruig time is Ω, i.e., expoetial. 4 THE LEADINGONES PROBLEM I this sectio, we first aalyze the ruig time of the +- EA o LeadigOes uder bit-wise oise p, ad bit-wise oise, q, respectively. The, we trasfer the aalysis from bit-wise oise p, to oe-bit oise; the results are complemetary to the kow oes recetly derived i []. 4. Bit-wise Noise p, For bit-wise oise p,, we prove i the followig three theorems that the expected ruig time is polyomial if p = Olog / ad super-polyomial if p = ωlog /. Their proofs are accomplished by applyig additive drift aalysis, the simplified drift theorem with self-loops ad the simplified drift theorem, respectively. THEOREM 4.. For the +-EA o LeadigOes uder bitwise oise p,, the expected ruig time is polyomial if p = Olog /. PROOF. We use Lemma.5 to prove it. For ay positive costat b, suppose that p b log /. Let LOx deote the true umber of leadig -bits of a solutio x. We first costruct a distace fuctio V x as, for ay x with LOx = i, V x = + c + c i, where c = 4b log +. It is easy to verify that V x X = { } = 0 ad V x / X > 0. The, we ivestigate E[[V ξ t V ξ t+ ξ t = x]] for ay x with LOx < i.e., x / X. Assume that curretly LOx = i, where 0 i. We divide the drift ito two parts: positive E + ad egative E. That is, E[[V ξ t V ξ t+ ξ t = x]] = E + E. For the positive drift, we eed to cosider that the umber of leadig -bits is icreased. By mutatio, we have P LOx i + = i, sice it eeds to flip the i + -th bit which must be 0 of x ad keep the i leadig -bits uchaged. For ay x with LOx i +, f x < f x implies that f x i or f x i +. Note that, P f x i = p i, sice at least oe of the first i leadig -bits of x eeds to be flipped by oise; P f x i + = p i, 3 sice it eeds to flip the first 0-bit of x ad keep the leadig -bits uchaged by oise. By the uio boud, we get P f x f x = P f x < f x p i+ p i +, 4 where the last iequality is by p = Olog /. Furthermore, for ay x with V x i +, V x V x + c i+ + c i= c + c i. 5 By combiig Eqs., 4 ad 5, we have E + i c + c i c + c i, 6 where the last iequality is by i. e 3 For the egative drift, we eed to cosider that the umber of leadig -bits is decreased. By mutatio, we have P LOx i = i, 6 sice it eeds to flip at least oe leadig -bit of x. For ay x with LOx i where i, f x f x implies
5 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay that f x i or f x i. Note that, P f x i p i, 7 sice for the first i bits of x, it eeds to flip the 0-bits whose umber is at least ad keep the -bits uchaged by oise; P f x i = p i, 8 sice at least oe leadig -bit of x eeds to be flipped by oise. By the uio boud, we get P f x f x p p i p i+.9 Furthermore, for ay x with LOx i, V x V x + c i. 0 By combiig Eqs. 6, 9 ad 0, we have E i p i + p + c i p e 3 + c i + c i. Thus, by subtractig E from E +, we have E[[V ξ t V ξ t+ ξ t = x]] + c i c 6 p 3 + c i 4b log + b log 6 3 6, where the secod iequality is by c = 4b log + ad p b log /. Note that V x + c e c = e 4b log + = e 4b. By Lemma.5, we get E[[τ ξ 0]] 6 e 4b = O 4b+, i.e., the expected ruig time is polyomial. THEOREM 4.. For the +-EA o LeadigOes uder bit-wise oise p,, if p = ωlog / o, the expected ruig time is super-polyomial. PROOF. We use Lemma.7 to prove it. Let X t = x 0 be the umber of 0-bits of the solutio x after t iteratios of the +-EA. Let c be ay positive costat. We cosider the iterval [0, c log ], i.e., the parameters a = 0 i.e., the global optimum ad b = c log i Lemma.7. The, we aalyze the drift E[[X t X t+ X t = i]] for i < c log. As i the proof of Theorem 4., we divide the drift ito two parts: positive E + ad egative E. That is, E[[X t X t+ X t = i]] = E + E. For the positive drift, we eed to cosider that the umber of 0-bits is decreased. For mutatio o x where x 0 = i, let X ad Y deote the umber of flipped 0-bits ad -bits, respectively. The, X Bi, ad Y B i,, where B, is the biomial distributio. Let P mutx, x deote the probability of geeratig x by mutatig x. To estimate a upper boud o E +, we assume that the offsprig solutio x with x 0 < i is always accepted. Thus, we have E + i P mutx, x i x 0 = k P X Y = k x : x 0 <i = i = i i k= k i j= j k= j=k k= P X = j P Y = j k k P X = j P Y = j k j= j P X = j = i. For the egative drift, we eed to cosider that the umber of 0-bits is icreased. We aalyze the i cases where oly oe -bit is flipped i.e., x 0 = i +, which happes with probability. Assume that LOx = k i. If the j-th where j k leadig -bit is flipped, the offsprig solutio x will be accepted i.e., f x f x if f x j ad f x j. Note that, P f x j = p+p j p j, where the equality is sice it eeds to keep the j leadig -bits of x uchaged, ad the last iequality is by p = o; P f x j = p j 3 = p j + j p e j pj 3, where the equality is sice at least oe of the first j leadig -bits of x eeds to be flipped by oise. Thus, we get P f x f x pj 6. 4 If oe of the i k o-leadig -bits is flipped, LOx = LOx = k. We ca use the same aalysis procedure as Eq. 4 i the proof of Theorem 4. to derive that P f x f x p k +, 5 where the secod iequality is by p = o. Combiig all the i cases, we get E k pj 6 + i k i+ i 6 j= pkk + + i k pk e 36 + i k. 6 By subtractig E from E +, we get E[[X t X t+ X t = i]] i pk 36 i k. 6 To ivestigate coditio of Lemma.7, we also eed to aalyze the probability P X t+ i X t = i. For X t+ i, it is ecessary that at least oe bit of x is flipped ad the offsprig x is accepted. We cosider two cases: at least oe of the k leadig -bits of x is flipped; the k leadig - bits of x are ot flipped ad at least oe of the last k bits is flipped. For case, the mutatio probability is k
6 GECCO 7, July 5-9, 07, Berli, Germay C. Qia et. al. ad the acceptace probability is at most p k+ by Eq. 9. For case, the mutatio probability is k k k ad the acceptace probability is at most. Thus, P X t+ i X t = i p + k. 7 Whe k < p, we have E[[X t X t+ X t = i]] i i k 6 k p/ 7c log k 6 4 p+ k 8, where the secod iequality is by k > p ad i < c log, the third iequality is by p = ωlog / ad the last is by k > p. Whe k p, we have E[[X t X t+ X t = i]] i pk 9 36 c log p 44 p 88 p + k, 576 where the secod iequality is by p = o ad i < c log, the third is by p = ωlog / ad the last is by k p. Combiig Eqs. 7, 8 ad 9, we get that coditio of Lemma.7 holds with ɛ =. 576 For coditio of Lemma.7, we eed to show P X t+ X t j X t = i rl P X +δ j t+ i X t = i for i. For P X t+ i X t = i, we aalyze the cases where oly oe bit is flipped. Usig the similar aalysis procedure as E, except that flippig ay bit rather tha oly -bit is cosidered here, we easily get P X t+ i X t = i pkk k 6. 0 For X t+ X t j, it is ecessary that at least j bits of x are flipped ad the offsprig solutio x is accepted. We cosider two cases: at least oe of the k leadig -bits is flipped; the k leadig -bits are ot flipped. For case, the mutatio j probability is at most k ad the acceptace probability is at most p k+ by Eq. 9. For case, the mutatio j probability is at most k k ad the acceptace j j probability is at most. Thus, we have P X t+ X t j X t = i k j p k + j + k k j j pkk+ 4 + k j pkk+ + k 44 j j By combiig Eq. 0 with Eq., we get that coditio of Lemma.7 holds with δ = ad rl = 44. Note that l = b a = c log. By Lemma.7, the expected ruig time is Ωc log, where c is ay positive costat. Thus, the expected ruig time is super-polyomial. THEOREM 4.3. For the +-EA o LeadigOes uder bitwise oise p,, the expected ruig time is expoetial if p = Ω. PROOF. We use Lemma.6 to prove it. Let X t = i be the umber of 0-bits of the solutio x after t iteratios of the +-EA. We cosider the iterval i [0, / ]. To aalyze the drift E[[X t X t+ X t = i]] = E + E, we use the same aalysis procedure as Theorem 4.. For the positive drift, we have E + i = o. For the egative drift, we re-aalyze Eqs. 4 ad 5. P f x j p j Thus, Eq. 4 becomes From Eqs. ad 3, we get that ad P f x j pj. 3 P f x f x p j 3 j. For Eq. 5, we eed to aalyze the acceptace probability for LOx = LOx = k. Sice it is sufficiet to keep the first k + bits of x ad x uchaged i oise, Eq. 5 becomes P f x f x p k+ p k+. 3 By applyig the above two iequalities to Eq. 6, we have k E p j j + + i k k =Ω, e 3 j= where the equality is by p = Ω. Thus, E + E = Ω. That is, coditio of Lemma.6 holds. Sice it is ecessary to flip at least j bits of x, we have P X t+ X t j X t j j j!, j which implies that coditio of Lemma.6 holds with δ = ad rl =. Note that l = /. Thus, by Lemma.6, the expected ruig time is expoetial. 4. Bit-wise Noise, q For bit-wise oise, q, we prove i the followig three theorems that the expected ruig time is polyomial if q = Olog / 3 ad super-polyomial if q = ωlog /. The proof idea is similar to that for bit-wise oise p,. The mai differece led by the chage of oise is the probability of acceptig the offsprig solutio, i.e., P f x f x. THEOREM 4.4. For the +-EA o LeadigOes uder bitwise oise, q, the expected ruig time is polyomial if q = Olog / 3. PROOF. The proof is very similar to that of Theorem 4.. The chage of oise oly affects the probability of acceptig the offsprig solutio i the aalysis. For ay positive costat b, suppose that q b log / 3. For the positive drift E +, we eed to re-aalyze P f x f x i.e., Eq. 4 i the proof of Theorem 4. for the paret x with LOx = i ad the offsprig x with LOx i +. By bit-wise oise, q, Eqs. ad 3 chage to P f x i = q i ; P f x i+= q i q. Thus, by the uio boud, Eq. 4 becomes P f x f x q i + q i q 4 = q i+ qi + /,
7 Ruig Time Aalysis of the +-EA uder Bit-wise Noise GECCO 7, July 5-9, 07, Berli, Germay where the last iequality is by q = Olog / 3. For the egative drift E, we eed to re-aalyze P f x f x i.e., Eq. 9 i the proof of Theorem 4. for the paret x with LOx = i where i ad the offsprig x with LOx i. By bit-wise oise, q, Eqs. 7 ad 8 chage to P f x i q q i, P f x i = q i. Thus, by the uio boud, Eq. 9 becomes P f x f x q q i + q i 5 = q i q i q q i+q, where the secod iequality is by q i i q ad q > 0 for q = Olog / 3. By applyig Eq. 4 ad Eq. 5 to E + ad E, respectively, Eq. chages to E[[V ξ t V ξ t+ ξ t = x]] + c i 4b log c i c 6 qi+ 3 b log That is, the coditio of Lemma.5 still holds with. Thus, 6 the expected ruig time is polyomial. THEOREM 4.5. For the +-EA o LeadigOes uder bit-wise oise, q, if q = ωlog / o/, the expected ruig time is super-polyomial. PROOF. We use the same aalysis process as Theorem 4.. The oly differece is the probability of acceptig the offsprig solutio due to the chage of oise. For the positive drift, we still have E + i, sice we assume that x is always accepted i the proof of Theorem 4.. For the egative drift, we eed to re-aalyze P f x f x for the paret solutio x with LOx = k ad the offsprig solutio x with LOx = j where j k +. For j k, to derive a lower boud o P f x f x, we cosider the j cases where f x = l ad f x l for 0 l j. Sice P f x = l = q l q ad P f x l = q l, Eq. 4 chages to j P f x f x q l q q l qj l=0 = qj + q q 6 j q j qj q qj, where the last iequality is by q j qj / sice q = o/. For j = k + i.e., LOx = LOx = k, we ca use the same aalysis as Eq. 4 to derive a lower boud /, sice the last iequality of Eq. 4 still holds with p = o/. Thus, Eq. 5 also holds here, i.e., P f x f x. 7 By applyig Eqs. 6 ad 7 to E, Eq. 6 chages to E qk + i k. 6 Thus, we have E[[X t X t+ X t = i]] = E + E i qk i k. 6 For the upper boud aalysis of P X t+ i X t = i i the proof of Theorem 4., we oly eed to replace the acceptace probability p k+ i the case of LOx < LOx with k + q i.e., Eq. 5. Thus, Eq. 7 chages to P X t+ i X t = i k + q + k q + k. To compare E[[X t X t+ X t = i]] with P X t+ i X t = i, we cosider two cases: k < q ad k q. By usig q = ωlog / ad applyig the same aalysis procedure as Eqs. 8 ad 9, we ca derive that coditio of Lemma.7 holds with ɛ =. 9 For the lower boud aalysis of P X t+ i X t = i, by applyig Eqs. 6 ad 7, Eq. 0 chages to qkk + P X t+ i X t = i + k 6. For the aalysis of X t+ X t j, by replacig the acceptace probability p k+ i the case of LOx < LOx with k + q, Eq. chages to P X t+ X t j X t = i qkk+ 4 + k j j qkk + + k j That is, coditio of Lemma.7 holds with δ =, rl = 48. Thus, the expected ruig time is super-polyomial. THEOREM 4.6. For the +-EA o LeadigOes uder bitwise oise, q, the expected ruig time is expoetial if q = Ω/. PROOF. We use Lemma.6 to prove it. Let X t = i be the umber of 0-bits of the solutio x after t iteratios of the +-EA. We cosider the iterval i [0, / ]. To aalyze the drift E[[X t X t+ X t = i]], we use the same aalysis procedure as the proof of Theorem 4.. We first cosider q = Ω/ o. We eed to aalyze the probability P f x f x, where the offsprig solutio x is geerated by flippig oly oe -bit of x. Let LOx = k. For the case where the j-th where j k leadig -bit is flipped, as the aalysis of Eq. 6, we get P f x f x If q j <, qj 4 qj q j qj q. ; otherwise, qj qj q qj. Thus, we have P f x f x mi{/4, qj/}. For the case that flips oe o-leadig -bit i.e., LOx = LOx = k, to derive a lower boud o P f x f x, we cosider f x = l ad f x l for 0 l k. Thus, k P f x f x q l q q l + q k+ q k qk l=0 + q k+ = + qk q,
8 GECCO 7, July 5-9, 07, Berli, Germay where the last iequality is by q = o. By applyig the above two iequalities to Eq. 6, we get k E { mi e 4, qj } + i k. j= If k, k j= mi{, qj } = Ω sice q = Ω/. If 4 k <, i k = Ω sice i. Thus, E = Ω. For q = Ω, we use the trivial lower boud q for the probability of acceptig the offsprig solutio x, sice it is sufficiet to flip the first leadig -bit of x by oise. The, E iq kq + i kq = = Ω. e e Thus, for q = Ω/, we have E[[X t X t+ X t = i]] = E + E i Ω = Ω. That is, coditio of Lemma.6 holds. Its coditio trivially holds with δ = ad rl =. Thus, the expected ruig time is expoetial. 4.3 Oe-bit Noise For the +-EA o LeadigOes uder oe-bit oise, it has bee kow that the ruig time is polyomial if p 6e ad expoetial if p = []. We exted this result by provig i Theorem 4.7 that the ruig time is polyomial if p = Olog / ad super-polyomial if p = ωlog /. The proof ca be accomplished as same as that of Theorems 4., 4. ad 4.3 for bit-wise oise p,. This is because although the probabilities P f x f x of acceptig the offsprig solutio are differet, their bouds used i the proofs for bit-wise oise p, still hold for oe-bit oise. THEOREM 4.7. For the +-EA o LeadigOes uder oe-bit oise, the expected ruig time is polyomial if p = Olog /, super-polyomial if p = ωlog / o ad expoetial if p = Ω. PROOF. We re-aalyze P f x f x for oe-bit oise, ad show that the bouds o P f x f x used i the proofs for bit-wise oise p, still hold for oe-bit oise. For the proof of Theorem 4., Eqs. ad 3 chage to P f x i = p i, P f x i + = p, ad thus Eq. 4 still holds; Eqs. 7 ad 8 chage to P f x i p, P f x i = p i, ad thus Eq. 9 still holds. For the proof of Theorem 4., Eqs. ad 3 chage to P f x j = p j, P f x j = p j, ad thus Eq. 4 still holds. For the proof of Theorem 4.3, Eq. still holds by the above two equalities; Eq. 3 still holds sice the probability of keepig the first k + bits of a solutio uchaged i oe-bit oise is p k+ p k+. 5 CONCLUSION C. Qia et. al. I this paper, we theoretically study the +-EA solvig OeMax ad LeadigOes uder bit-wise oise. We derive the rages of oise parameters for the ruig time beig polyomial ad super-polyomial, respectively. The previously kow rages for the +-EA solvig LeadigOes uder oe-bit oise are also improved. I the future, we shall improve the curretly derived bouds o LeadigOes, as they do ot cover the whole rage of oise parameters. ACKNOWLEDGMENTS We wat to thak the reviewers for their valuable commets. This work was supported by the NSFC , , the Fudametal Research Fuds for the Cetral Uiversities ad the CCF-Tecet Ope Research Fud. REFERENCES [] Y. Akimoto, S. Astete-Morales, ad O. Teytaud. 05. Aalysis of rutime of optimizatio algorithms for oisy fuctios over discrete codomais. Theoretical Computer Sciece , [] A. Auger ad B. Doerr. 0. Theory of Radomized Search Heuristics: Foudatios ad Recet Developmets. World Scietific, Sigapore. [3] L. Biachi, M. Dorigo, L. Gambardella, ad W. Gutjahr A survey o metaheuristics for stochastic combiatorial optimizatio. Natural Computig 8, 009, [4] D.-C. Dag ad P. K. Lehre. 05. Efficiet optimisatio of oisy fitess fuctios with populatio-based evolutioary algorithms. I Proceedigs of FOGA 5. Aberystwyth, UK, [5] B. Doerr, A. Hota, ad T. Kötzig. 0. Ats easily solve stochastic shortest path problems. I Proceedigs of GECCO. Philadelphia, PA, 7 4. [6] S. Droste Aalysis of the + EA for a oisy OeMax. I Proceedigs of GECCO 04. Seattle, WA, [7] S. Droste, T. Jase, ad I. Wegeer. 00. O the aalysis of the + evolutioary algorithm. Theoretical Computer Sciece 76, - 00, 5 8. [8] T. Friedrich, T. Kötzig, M. Krejca, ad A. Sutto. 06. The compact geetic algorithm is efficiet uder extreme Gaussia oise. IEEE Trasactios o Evolutioary Computatio 06, i press. [9] T. Friedrich, T. Kötzig, M. Krejca, ad A. Sutto. 06. Robustess of at coloy optimizatio to oise. Evolutioary Computatio 4, 06, [0] T. Friedrich, T. Kötzig, F. Quiza, ad A. Sutto. 07. Resamplig vs recombiatio: A statistical ru time estimatio. I Proceedigs of FOGA 7. Copehage, Demark. [] C. Gieße ad T. Kötzig. 06. Robustess of populatios i stochastic eviromets. Algorithmica 75, 3 06, [] J. He ad X. Yao. 00. Drift aalysis ad average time complexity of evolutioary algorithms. Artificial Itelligece 7, 00, [3] Y. Ji ad J. Brake Evolutioary optimizatio i ucertai eviromets-a survey. IEEE Trasactios o Evolutioary Computatio 9, 3 005, [4] F. Neuma ad C. Witt. 00. Bioispired Computatio i Combiatorial Optimizatio: Algorithms ad Their Computatioal Complexity. Spriger-Verlag, Berli, Germay. [5] P. Oliveto ad C. Witt. 0. Simplified drift aalysis for provig lower bouds i evolutioary computatio. Algorithmica 59, 3 0, [6] P. Oliveto ad C. Witt. 0. Erratum: Simplified drift aalysis for provig lower bouds i evolutioary computatio. CORR abs/ [7] A. Prugel-Beett, J. Rowe, ad J. Shapiro. 05. Ru-time aalysis of populatio-based evolutioary algorithm i oisy eviromets. I Proceedigs of FOGA 5. Aberystwyth, UK, [8] C. Qia, Y. Yu, K. Tag, Y. Ji, X. Yao, ad Z.-H. Zhou. 06. O the effectiveess of samplig for evolutioary optimizatio i oisy eviromets. Evolutioary Computatio 06, i press. [9] C. Qia, Y. Yu, ad Z.-H. Zhou. 05. Aalyzig evolutioary optimizatio i oisy eviromets. Evolutioary Computatio 05, i press. [0] J. E. Rowe ad D. Sudholt. 04. The choice of the offsprig populatio size i the,λ evolutioary algorithm. Theoretical Computer Sciece , 0 38.
arxiv: v1 [cs.ne] 2 Nov 2017
Noame mauscript No. will be iserted by the editor) Ruig Time Aalysis of the +)-EA for OeMax ad LeadigOes uder Bit-wise Noise Chao Qia Chao Bia Wu Jiag Ke Tag Received: date / Accepted: date arxiv:7.00956v
More informationOn the Effectiveness of Sampling for Evolutionary Optimization in Noisy Environments
O the Effectiveess of Samplig for Evolutioary Optimizatio i Noisy Eviromets Chao Qia,2 chaoqia@ustc.edu.c Yag Yu 2 yuy@ju.edu.c Ke Tag ketag@ustc.edu.c Yaochu Ji 3 yaochu.ji@surrey.ac.uk Xi Yao,4 x.yao@cs.bham.ac.uk
More informationarxiv: v1 [cs.ai] 20 Nov 2013
Aalyzig Evolutioary Optimizatio i Noisy Eviromets Chao Qia, Yag Yu, Zhi-Hua Zhou Natioal Key Laboratory for Novel Software Techology Najig Uiversity, Najig 20023, Chia arxiv:3.4987v [cs.ai] 20 Nov 203
More informationGeneral Lower Bounds for the Running Time of Evolutionary Algorithms
Geeral Lower Bouds for the Ruig Time of Evolutioary Algorithms Dirk Sudholt Iteratioal Computer Sciece Istitute, Berkeley, CA 94704, USA Abstract. We preset a ew method for provig lower bouds i evolutioary
More informationAnalyzing Evolutionary Optimization in Noisy Environments
Evolutioary Computatio /EVCO_a_0070-Qia Jauary 5, 206 6:7 Aalyzig Evolutioary Optimizatio i Noisy Eviromets Chao Qia qiac@lamda.ju.edu.c Natioal Key Laboratory for Novel Software Techology, Najig Uiversity,
More informationSelection Hyper-heuristics Can Provably be Helpful in Evolutionary Multi-objective Optimization
Selectio Hyper-heuristics Ca Provably be Helpful i Evolutioary Multi-objective Optimizatio Chao Qia 1,2, Ke Tag 1, ad Zhi-Hua Zhou 2 1 UBRI, School of Computer Sciece ad Techology, Uiversity of Sciece
More informationA Rigorous View On Neutrality
A Rigorous View O Neutrality Bejami Doerr Michael Gewuch Nils Hebbighaus Frak Neuma Algorithms ad Complexity Group Max-Plack-Istitut für Iformatik Saarbrücke, Germay Departmet of Computer Sciece Christia-Albrechts-Uiversity
More informationReinforcement Learning Based Dynamic Selection of Auxiliary Objectives with Preserving of the Best Found Solution
Reiforcemet Learig Based Dyamic Selectio of Auxiliary Objectives with Preservig of the Best Foud Solutio arxiv:1704.07187v1 [cs.ne] 24 Apr 2017 Abstract Efficiecy of sigle-objective optimizatio ca be improved
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationarxiv: v1 [cs.ne] 4 Sep 2017
Theoretical Aalysis of Stochastic Search Algorithms Per Kristia Lehre School of Computer Sciece, Uiversity of Birmigham, Birmigham, UK Pietro S. Oliveto Departmet of Computer Sciece, Uiversity of Sheffield,
More informationProblem Set 2 Solutions
CS271 Radomess & Computatio, Sprig 2018 Problem Set 2 Solutios Poit totals are i the margi; the maximum total umber of poits was 52. 1. Probabilistic method for domiatig sets 6pts Pick a radom subset S
More informationOn a Smarandache problem concerning the prime gaps
O a Smaradache problem cocerig the prime gaps Felice Russo Via A. Ifate 7 6705 Avezzao (Aq) Italy felice.russo@katamail.com Abstract I this paper, a problem posed i [] by Smaradache cocerig the prime gaps
More informationUpper and Lower Bounds on Unrestricted Black-Box Complexity of Jump n,l
Upper ad Lower Bouds o Urestricted Black-Box Complexity of Jump,l Maxim Buzdalov 1, Mikhail Kever 1, ad Bejami Doerr 1 ITMO Uiversity, 49 Kroverkskiy av., Sait-Petersburg, Russia, 197101 mbuzdalov@gmail.com,
More informationDrift analysis and average time complexity of evolutionary algorithms
Artificial Itelligece 127 (2001) 57 85 Drift aalysis ad average time complexity of evolutioary algorithms Ju He a,xiyao b, a Departmet of Computer Sciece, Norther Jiaotog Uiversity, Beijig 100044, PR Chia
More informationApplication to Random Graphs
A Applicatio to Radom Graphs Brachig processes have a umber of iterestig ad importat applicatios. We shall cosider oe of the most famous of them, the Erdős-Réyi radom graph theory. 1 Defiitio A.1. Let
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationA Novel Genetic Algorithm using Helper Objectives for the 0-1 Knapsack Problem
A Novel Geetic Algorithm usig Helper Objectives for the 0-1 Kapsack Problem Ju He, Feidu He ad Hogbi Dog 1 arxiv:1404.0868v1 [cs.ne] 3 Apr 2014 Abstract The 0-1 kapsack problem is a well-kow combiatorial
More informationLecture 11: Pseudorandom functions
COM S 6830 Cryptography Oct 1, 2009 Istructor: Rafael Pass 1 Recap Lecture 11: Pseudoradom fuctios Scribe: Stefao Ermo Defiitio 1 (Ge, Ec, Dec) is a sigle message secure ecryptio scheme if for all uppt
More informationRuntime Analysis of the (µ+1) EA on Simple Pseudo-Boolean Functions
Rutime Aalysis of the (µ+1) EA o Simple Pseudo-Boolea Fuctios Carste Witt FB Iformatik, LS 2, Uiversität Dortmud, 44221 Dortmud, Germay carste.witt@cs.ui-dortmud.de November 14, 2005 Abstract Although
More informationElement sampling: Part 2
Chapter 4 Elemet samplig: Part 2 4.1 Itroductio We ow cosider uequal probability samplig desigs which is very popular i practice. I the uequal probability samplig, we ca improve the efficiecy of the resultig
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationLecture 3: August 31
36-705: Itermediate Statistics Fall 018 Lecturer: Siva Balakrisha Lecture 3: August 31 This lecture will be mostly a summary of other useful expoetial tail bouds We will ot prove ay of these i lecture,
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationThe Growth of Functions. Theoretical Supplement
The Growth of Fuctios Theoretical Supplemet The Triagle Iequality The triagle iequality is a algebraic tool that is ofte useful i maipulatig absolute values of fuctios. The triagle iequality says that
More informationOptimally Sparse SVMs
A. Proof of Lemma 3. We here prove a lower boud o the umber of support vectors to achieve geeralizatio bouds of the form which we cosider. Importatly, this result holds ot oly for liear classifiers, but
More informationAdvanced Stochastic Processes.
Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.
More informationPopulations can be essential in tracking dynamic optima
Populatios ca be essetial i trackig dyamic optima Duc-Cuog Dag 1, Thomas Jase 2, ad Per Kristia Lehre 1 arxiv:1607.03317v1 [cs.ne] 12 Jul 2016 1 School of Computer Sciece, Uiversity of Nottigham Jubilee
More informationPosted-Price, Sealed-Bid Auctions
Posted-Price, Sealed-Bid Auctios Professors Greewald ad Oyakawa 207-02-08 We itroduce the posted-price, sealed-bid auctio. This auctio format itroduces the idea of approximatios. We describe how well this
More informationLet us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.
Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,
More informationACO Comprehensive Exam 9 October 2007 Student code A. 1. Graph Theory
1. Graph Theory Prove that there exist o simple plaar triagulatio T ad two distict adjacet vertices x, y V (T ) such that x ad y are the oly vertices of T of odd degree. Do ot use the Four-Color Theorem.
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationLecture 2 February 8, 2016
MIT 6.854/8.45: Advaced Algorithms Sprig 206 Prof. Akur Moitra Lecture 2 February 8, 206 Scribe: Calvi Huag, Lih V. Nguye I this lecture, we aalyze the problem of schedulig equal size tasks arrivig olie
More informationRandomized Algorithms I, Spring 2018, Department of Computer Science, University of Helsinki Homework 1: Solutions (Discussed January 25, 2018)
Radomized Algorithms I, Sprig 08, Departmet of Computer Sciece, Uiversity of Helsiki Homework : Solutios Discussed Jauary 5, 08). Exercise.: Cosider the followig balls-ad-bi game. We start with oe black
More informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
More informationBasics of Probability Theory (for Theory of Computation courses)
Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.
More informationON POINTWISE BINOMIAL APPROXIMATION
Iteratioal Joural of Pure ad Applied Mathematics Volume 71 No. 1 2011, 57-66 ON POINTWISE BINOMIAL APPROXIMATION BY w-functions K. Teerapabolar 1, P. Wogkasem 2 Departmet of Mathematics Faculty of Sciece
More informationMonte Carlo Integration
Mote Carlo Itegratio I these otes we first review basic umerical itegratio methods (usig Riema approximatio ad the trapezoidal rule) ad their limitatios for evaluatig multidimesioal itegrals. Next we itroduce
More informationDETERMINATION OF MECHANICAL PROPERTIES OF A NON- UNIFORM BEAM USING THE MEASUREMENT OF THE EXCITED LONGITUDINAL ELASTIC VIBRATIONS.
ICSV4 Cairs Australia 9- July 7 DTRMINATION OF MCHANICAL PROPRTIS OF A NON- UNIFORM BAM USING TH MASURMNT OF TH XCITD LONGITUDINAL LASTIC VIBRATIONS Pavel Aokhi ad Vladimir Gordo Departmet of the mathematics
More informationarxiv: v1 [cs.ne] 26 Mar 2014
How Crossover Speeds Up Buildig-Block Assembly i Geetic Algorithms arxiv:403.6600v [cs.ne] 26 Mar 204 Dirk Sudholt Departmet of Computer Sciece, Uiversity of Sheffield, Uited Kigdom Abstract We re-ivestigate
More informationExpectation and Variance of a random variable
Chapter 11 Expectatio ad Variace of a radom variable The aim of this lecture is to defie ad itroduce mathematical Expectatio ad variace of a fuctio of discrete & cotiuous radom variables ad the distributio
More informationarxiv: v2 [cs.ne] 26 Nov 2014
How Crossover Speeds Up Buildig-Block Assembly i Geetic Algorithms arxiv:403.6600v2 [cs.ne] 26 Nov 204 Dirk Sudholt Departmet of Computer Sciece, Uiversity of Sheffield, Uited Kigdom Abstract We re-ivestigate
More informationIP Reference guide for integer programming formulations.
IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more
More informationRank Modulation with Multiplicity
Rak Modulatio with Multiplicity Axiao (Adrew) Jiag Computer Sciece ad Eg. Dept. Texas A&M Uiversity College Statio, TX 778 ajiag@cse.tamu.edu Abstract Rak modulatio is a scheme that uses the relative order
More informationComputational Intelligence Winter Term 2018/19
Computatioal Itelligece Witer Term 28/9 Prof. Dr. Güter Rudolph Lehrstuhl für Algorithm Egieerig (LS ) Fakultät für Iformatik TU Dortmud Pla for Today Lecture Evolutioary Algorithms (EA) Optimizatio Basics
More informationSince X n /n P p, we know that X n (n. Xn (n X n ) Using the asymptotic result above to obtain an approximation for fixed n, we obtain
Assigmet 9 Exercise 5.5 Let X biomial, p, where p 0, 1 is ukow. Obtai cofidece itervals for p i two differet ways: a Sice X / p d N0, p1 p], the variace of the limitig distributio depeds oly o p. Use the
More informationStat 421-SP2012 Interval Estimation Section
Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible
More informationMATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4
MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.
More informationVector Quantization: a Limiting Case of EM
. Itroductio & defiitios Assume that you are give a data set X = { x j }, j { 2,,, }, of d -dimesioal vectors. The vector quatizatio (VQ) problem requires that we fid a set of prototype vectors Z = { z
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals
More informationDesign and Analysis of Algorithms
Desig ad Aalysis of Algorithms Probabilistic aalysis ad Radomized algorithms Referece: CLRS Chapter 5 Topics: Hirig problem Idicatio radom variables Radomized algorithms Huo Hogwei 1 The hirig problem
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationAverage Convergence Rate of Evolutionary Algorithms
1 Average Covergece Rate of Evolutioary Algorithms Ju He ad Guagmig Li Abstract arxiv:154.8117v3 [cs.ne] 2 Ju 215 I evolutioary optimizatio, it is importat to uderstad how fast evolutioary algorithms coverge
More informationHypothesis Testing. Evaluation of Performance of Learned h. Issues. Trade-off Between Bias and Variance
Hypothesis Testig Empirically evaluatig accuracy of hypotheses: importat activity i ML. Three questios: Give observed accuracy over a sample set, how well does this estimate apply over additioal samples?
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationLecture 7: October 18, 2017
Iformatio ad Codig Theory Autum 207 Lecturer: Madhur Tulsiai Lecture 7: October 8, 207 Biary hypothesis testig I this lecture, we apply the tools developed i the past few lectures to uderstad the problem
More informationLecture 14: Graph Entropy
15-859: Iformatio Theory ad Applicatios i TCS Sprig 2013 Lecture 14: Graph Etropy March 19, 2013 Lecturer: Mahdi Cheraghchi Scribe: Euiwoog Lee 1 Recap Bergma s boud o the permaet Shearer s Lemma Number
More informationCS / MCS 401 Homework 3 grader solutions
CS / MCS 401 Homework 3 grader solutios assigmet due July 6, 016 writte by Jāis Lazovskis maximum poits: 33 Some questios from CLRS. Questios marked with a asterisk were ot graded. 1 Use the defiitio of
More informationw (1) ˆx w (1) x (1) /ρ and w (2) ˆx w (2) x (2) /ρ.
2 5. Weighted umber of late jobs 5.1. Release dates ad due dates: maximimizig the weight of o-time jobs Oce we add release dates, miimizig the umber of late jobs becomes a sigificatly harder problem. For
More informationEcon 325/327 Notes on Sample Mean, Sample Proportion, Central Limit Theorem, Chi-square Distribution, Student s t distribution 1.
Eco 325/327 Notes o Sample Mea, Sample Proportio, Cetral Limit Theorem, Chi-square Distributio, Studet s t distributio 1 Sample Mea By Hiro Kasahara We cosider a radom sample from a populatio. Defiitio
More informationEscaping Local Optima Using Crossover with Emergent Diversity
This article has bee accepted for publicatio i a future issue of this joural, but has ot bee fully edited Cotet may chage prior to fial publicatio Citatio iformatio: DOI 009/TEVC07740, IEEE Trasactios
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationThe coalescent coalescence theory
The coalescet coalescece theory Peter Beerli September 1, 009 Historical ote Up to 198 most developmet i populatio geetics was prospective ad developed expectatios based o situatios of today. Most work
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationLecture 7: Properties of Random Samples
Lecture 7: Properties of Radom Samples 1 Cotiued From Last Class Theorem 1.1. Let X 1, X,...X be a radom sample from a populatio with mea µ ad variace σ
More informationThe log-behavior of n p(n) and n p(n)/n
Ramauja J. 44 017, 81-99 The log-behavior of p ad p/ William Y.C. Che 1 ad Ke Y. Zheg 1 Ceter for Applied Mathematics Tiaji Uiversity Tiaji 0007, P. R. Chia Ceter for Combiatorics, LPMC Nakai Uivercity
More informationOn forward improvement iteration for stopping problems
O forward improvemet iteratio for stoppig problems Mathematical Istitute, Uiversity of Kiel, Ludewig-Mey-Str. 4, D-24098 Kiel, Germay irle@math.ui-iel.de Albrecht Irle Abstract. We cosider the optimal
More informationLast time, we talked about how Equation (1) can simulate Equation (2). We asserted that Equation (2) can also simulate Equation (1).
6896 Quatum Complexity Theory Sept 23, 2008 Lecturer: Scott Aaroso Lecture 6 Last Time: Quatum Error-Correctio Quatum Query Model Deutsch-Jozsa Algorithm (Computes x y i oe query) Today: Berstei-Vazirii
More information1 Inferential Methods for Correlation and Regression Analysis
1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet
More informationHashing and Amortization
Lecture Hashig ad Amortizatio Supplemetal readig i CLRS: Chapter ; Chapter 7 itro; Sectio 7.. Arrays ad Hashig Arrays are very useful. The items i a array are statically addressed, so that isertig, deletig,
More informationECE 901 Lecture 12: Complexity Regularization and the Squared Loss
ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality
More informationEcon 325 Notes on Point Estimator and Confidence Interval 1 By Hiro Kasahara
Poit Estimator Eco 325 Notes o Poit Estimator ad Cofidece Iterval 1 By Hiro Kasahara Parameter, Estimator, ad Estimate The ormal probability desity fuctio is fully characterized by two costats: populatio
More informationCS284A: Representations and Algorithms in Molecular Biology
CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by
More informationLecture 5: April 17, 2013
TTIC/CMSC 350 Mathematical Toolkit Sprig 203 Madhur Tulsiai Lecture 5: April 7, 203 Scribe: Somaye Hashemifar Cheroff bouds recap We recall the Cheroff/Hoeffdig bouds we derived i the last lecture idepedet
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationSimilarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall
Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig
More informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationCS 330 Discussion - Probability
CS 330 Discussio - Probability March 24 2017 1 Fudametals of Probability 11 Radom Variables ad Evets A radom variable X is oe whose value is o-determiistic For example, suppose we flip a coi ad set X =
More information-ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION
NEW NEWTON-TYPE METHOD WITH k -ORDER CONVERGENCE FOR FINDING SIMPLE ROOT OF A POLYNOMIAL EQUATION R. Thukral Padé Research Cetre, 39 Deaswood Hill, Leeds West Yorkshire, LS7 JS, ENGLAND ABSTRACT The objective
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationLecture 2. The Lovász Local Lemma
Staford Uiversity Sprig 208 Math 233A: No-costructive methods i combiatorics Istructor: Ja Vodrák Lecture date: Jauary 0, 208 Origial scribe: Apoorva Khare Lecture 2. The Lovász Local Lemma 2. Itroductio
More informationOn comparison of different approaches to the stability radius calculation. Olga Karelkina
O compariso of differet approaches to the stability radius calculatio Olga Karelkia Uiversity of Turku 2011 Outlie Prelimiaries Problem statemet Exact method for calculatio stability radius proposed by
More informationMachine Learning Brett Bernstein
Machie Learig Brett Berstei Week 2 Lecture: Cocept Check Exercises Starred problems are optioal. Excess Risk Decompositio 1. Let X = Y = {1, 2,..., 10}, A = {1,..., 10, 11} ad suppose the data distributio
More informationLecture 12: November 13, 2018
Mathematical Toolkit Autum 2018 Lecturer: Madhur Tulsiai Lecture 12: November 13, 2018 1 Radomized polyomial idetity testig We will use our kowledge of coditioal probability to prove the followig lemma,
More informationRandom assignment with integer costs
Radom assigmet with iteger costs Robert Parviaie Departmet of Mathematics, Uppsala Uiversity P.O. Box 480, SE-7506 Uppsala, Swede robert.parviaie@math.uu.se Jue 4, 200 Abstract The radom assigmet problem
More informationRecursive Algorithm for Generating Partitions of an Integer. 1 Preliminary
Recursive Algorithm for Geeratig Partitios of a Iteger Sug-Hyuk Cha Computer Sciece Departmet, Pace Uiversity 1 Pace Plaza, New York, NY 10038 USA scha@pace.edu Abstract. This article first reviews the
More informationSequences and Limits
Chapter Sequeces ad Limits Let { a } be a sequece of real or complex umbers A ecessary ad sufficiet coditio for the sequece to coverge is that for ay ɛ > 0 there exists a iteger N > 0 such that a p a q
More informationBerry-Esseen bounds for self-normalized martingales
Berry-Essee bouds for self-ormalized martigales Xiequa Fa a, Qi-Ma Shao b a Ceter for Applied Mathematics, Tiaji Uiversity, Tiaji 30007, Chia b Departmet of Statistics, The Chiese Uiversity of Hog Kog,
More informationThe Rand and block distances of pairs of set partitions
The Rad ad block distaces of pairs of set partitios Frak Ruskey 1 ad Jeifer Woodcock 1 Dept. of Computer Sciece, Uiversity of Victoria, CANADA Abstract. The Rad distaces of two set partitios is the umber
More informationA Note on Effi cient Conditional Simulation of Gaussian Distributions. April 2010
A Note o Effi ciet Coditioal Simulatio of Gaussia Distributios A D D C S S, U B C, V, BC, C April 2010 A Cosider a multivariate Gaussia radom vector which ca be partitioed ito observed ad uobserved compoetswe
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationECE 901 Lecture 14: Maximum Likelihood Estimation and Complexity Regularization
ECE 90 Lecture 4: Maximum Likelihood Estimatio ad Complexity Regularizatio R Nowak 5/7/009 Review : Maximum Likelihood Estimatio We have iid observatios draw from a ukow distributio Y i iid p θ, i,, where
More information1 Statement of the Game
ANALYSIS OF THE CHOW-ROBBINS GAME JON LU May 10, 2016 Abstract Flip a coi repeatedly ad stop wheever you wat. Your payoff is the proportio of heads ad you wish to maximize this payoff i expectatio. I this
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More information7.1 Convergence of sequences of random variables
Chapter 7 Limit theorems Throughout this sectio we will assume a probability space (Ω, F, P), i which is defied a ifiite sequece of radom variables (X ) ad a radom variable X. The fact that for every ifiite
More informationThe picture in figure 1.1 helps us to see that the area represents the distance traveled. Figure 1: Area represents distance travelled
1 Lecture : Area Area ad distace traveled Approximatig area by rectagles Summatio The area uder a parabola 1.1 Area ad distace Suppose we have the followig iformatio about the velocity of a particle, how
More informationRandom Walks on Discrete and Continuous Circles. by Jeffrey S. Rosenthal School of Mathematics, University of Minnesota, Minneapolis, MN, U.S.A.
Radom Walks o Discrete ad Cotiuous Circles by Jeffrey S. Rosethal School of Mathematics, Uiversity of Miesota, Mieapolis, MN, U.S.A. 55455 (Appeared i Joural of Applied Probability 30 (1993), 780 789.)
More informationEstimation of Backward Perturbation Bounds For Linear Least Squares Problem
dvaced Sciece ad Techology Letters Vol.53 (ITS 4), pp.47-476 http://dx.doi.org/.457/astl.4.53.96 Estimatio of Bacward Perturbatio Bouds For Liear Least Squares Problem Xixiu Li School of Natural Scieces,
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More informationLaw of the sum of Bernoulli random variables
Law of the sum of Beroulli radom variables Nicolas Chevallier Uiversité de Haute Alsace, 4, rue des frères Lumière 68093 Mulhouse icolas.chevallier@uha.fr December 006 Abstract Let be the set of all possible
More informationStochastic Simulation
Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso
More information