Estimation of Peaked Densities Over the Interval [0,1] Using Two-Sided Power Distribution: Application to Lottery Experiments

Size: px
Start display at page:

Download "Estimation of Peaked Densities Over the Interval [0,1] Using Two-Sided Power Distribution: Application to Lottery Experiments"

Transcription

1 MPRA Munich Peronal RePEc Archive Etimation of Peaed Denitie Over the Interval [0] Uing Two-Sided Power Ditribution: Application to Lottery Experiment Krzyztof Konte Artal Invetment 8. April 00 Online at MPRA Paper No. 378 poted 9. April 00 00:43 UTC

2 Etimation of Peaed Denitie Over the Interval [0] Uing Two-Sided Power Ditribution: Application to Lottery Experiment Krzyztof Konte Artal Invetment Waraw 3 Abtract Thi paper deal with etimating peaed denitie over the interval [0] uing twoided power ditribution (Kotz van Dorp 004). Such data were encountered in experiment determining certainty equivalent of lotterie (Konte 00). Thi paper ummarize the baic propertie of the two-ided power ditribution (TP) and it generalized form (GTP). The GTP maximum lielihood etimator a reult not derived by Kotz and van Dorp i preented. The TP and GTP are ued to etimate certainty equivalent denitie in two data et from lottery experiment. The obtained reult how that even a two-parametric TP ditribution provide more accurate etimate than the mooth three-parametric generalized beta ditribution GBT (Libby Novic 98) in one of the conidered data et. The three-parametric GTP ditribution outperform GBT for thee data. The reult are however the very oppoite for the econd data et in which the data are greatly cattered. The paper demontrate that the TP and GTP ditribution may be extremely ueful in etimating peaed denitie over the interval [0] and in tudying the relative utility function. JEL claification: C0 C3 C4 C6 C C5 C8 C9 D03 D8 D87 Keyword: Denity Ditribution; Maximum Lielihood Etimation; Lottery experiment; Relative Utility Function. In memoriam to Profeor Samuel Kotz who adly paed away on March 6 00 at the age of 80. I never had the honor of nowing or meeting Profeor Kotz but it wa he together with Rene van Dorp who invented the two-ided power ditribution. Quite coincidentally I began woring with thi ditribution on March Over the next couple of day I obtained the reult which form the ubject of thi and two other upcoming paper. I had planned to end thi paper to Prof. Kotz for hi comment and only learnt of hi paing while maing my final reviion. The author i grateful to Stefan Traub from the Department of Economic Univerity of Bremen Ulrich Schmidt from the Department of Economic Univerity of Kiel and Katarzyna Idziowa from the Center of Economic Pychology and Deciion Science Kozmini Univerity Waraw for maing the reult of their experiment available. 3 Contact: ul. Chrościciego 93/ Waraw Poland onte@artal.com.pl onte000@yahoo.com tel.:

3 . Introduction Peaed data are becoming more and more frequently oberved epecially in financial application. Such data may be decribed uing the Laplace ditribution which i an alternative to mooth ditribution function. Le examined are denitie defined over a bounded interval. Thi ha mainly been due to the lac of any ditribution uited to analyze uch cae. It wa only recently that Kotz and van Dorp (004) introduced a few type of two-ided power ditribution defined over the interval [0]. However a imple maximum lielihood etimator wa preent only for it baic form. Every other type required a recurive optimization procedure. Peaed data over the interval [0] were encountered in experiment determining certainty equivalent of lotterie (Konte 00). The denitie were etimated uing beta BT and generalized beta GBT (Libby Novic 98) ditribution. The obtained reult howed that uch mooth function may not be appropriate in all cae and that the ue of a peaed ditribution may give better reult. The aim of thi paper i therefore to preent the application of a two-ided power ditribution TP and it generalized form GTP in the conidered cae. The maximum lielihood etimator for GTP a reult not derived by Kotz and van Dorp i preented beforehand. The baic TP and GTP propertie are firt ummarized in Point and 3. The TP and GTP maximum lielihood etimator for unimodal denitie are preented in Point 4. The derivation of TP and GTP maximum lielihood etimator i demontrated in Point 5 where non-unimodal denitie are alo analyzed. Point 6 demontrate the ue of TP and GTP maximum lielihood etimator for two data et conidered previouly by Konte. Point 7 ummarize the paper and conclude that the TP and GTP ditribution may be extremely ueful in etimating peaed denitie over the interval [0] and in reearching the relative utility function.. Propertie of Two-Sided Power Ditribution.. The two-ided power ditribution TP preented here i the Standard Two-Sided Power (STSP) ditribution a conidered by Kotz and van Dorp (004). It ha two parameter: ( r λ γ) TP ; r γ λ = r γ λ γ γ if if 0 r λ λ< r. (.)

4 The TP implifie to a uniform ditribution for γ = a triangular ditribution for γ = a power ditribution for λ = and a reflected power ditribution for λ = 0. Sample TP hape are preented in Figure.. The function i unimodal for γ > uniform for γ = and U- haped for γ <. Pleae note that there are ome minor error in the figure preented by Kotz and van Dorp (004) a well in the textual decription (pp. 7-7). Figure.. Sample hape of two-ided power ditribution. On the left λ = 0.5 in the middle λ = 0.5 on the right λ = The value of γ parameter i equal to the denity at λ... The TP ha an extremely intereting feature in that it mode i given by the λ parameter: and the value at the mode by the γ parameter ( r λ γ) Mode TP ; = λ (.) ( r ) TP = λ; λ γ = γ. (.3) Thi mae thi ditribution unique in that it parameter may be eaily determined from it hape..3. The mean value of TP i equal to: ( ) λ γ + Mean TP ( r; λ γ) =. γ +.4. The Cumulative TP ditribution i defined a: (.4) ( r λ γ) CTP ; γ r λ if 0 r λ λ = γ r ( λ) if λ < r. λ (.5) which allow after inverting the quantile to be determined: ( r λ γ) Quantile q;tp ; q γ λ if 0 q λ λ = q γ ( λ) if λ< q. λ (.6) 3

5 In the pecial cae where q = 0.5 the median i given by: ( r λ γ) Median TP ;.5. The variance of TP i given by: γ γ γ λ if λ = ( λ)( λ) γ if λ <. ( ) ( ) ( γ + ) ( γ + ) γ γ λ λ Variance TP ( r; λ γ) =. (.7) (.8) 3. Propertie of the Generalized Two-Sided Power Ditribution 3.. The generalized two-ided power ditribution GTP preented here i the Generalized Standard Two-Sided Power (GSTSP) ditribution a conidered by Kotz and van Dorp and i a pecial cae of the Uneven Standard Two-Sided Power (USTSP) ditribution. It ha three parameter and i defined by: ( r λ γ δ) GTP ; r φ λ = r φ λ γ δ if if 0 r λ λ< r. (3.) where γδ φ =. λδ + λ γ ( ) (3.) The function i unimodal for γ > and δ > J-haped for γ > and δ < invere J- haped for γ < and δ > uniform from 0 to λ for γ = uniform from λ to for δ = U- haped with anti-mode at λ for γ < and δ < i power ditribution for λ = and i a reflected power ditribution for λ = 0. Figure 3.. Sample hape of generalized two-ided power GTP ditribution. On the left λ = 0.5 in the middle λ = 0.5 and on the right λ = All but the firt hape are referred to a non-unimodal throughout the paper. Some 4

6 ample hape are given in Figure 3.. Only unimodal curve are preented a they are of mot interet in thi paper and in practical application. An important cae (a will be dicued in 6.3) i when either γ or δ aume an infinite value. The GTP then reduce to a bounded TP ditribution defined over the interval [λ ] or [0 λ] repectively. 3.. The GTP mode i given by it λ parameter: and the value at the mode i given by: 3.3. The mean value of GTP i equal to: ( r λ γ δ) Mode GTP ; = λ (3.3) ( r ) GTP = λ; λ γ δ = φ. (3.4) γ λ Mean GTP ( r; λ γ δ) = + λδ +. ( γ + )( δ + ) λδ + ( λ) γ (3.5) 3.4. The Cumulative GTP i defined a: ( r λ γ δ) CGTP ; γ λφ r if 0 r λ γ λ = δ ( λ) φ r if λ< r. δ λ (3.6) which allow after inverting the ditribution quantile to be calculated: ( p λ γ δ) Quantile q;gtp ; q γ γ λφ λ if 0 q λφ γ = ( ) ( ) δ q δ λφ λ if < q. ( λ) φ γ (3.7) In the pecial cae where q = 0.5 the median i given by: ( r λ γ δ) Median GTP ; γ γ λφ λ if λφ γ = δ λφ δ ( λ) if <. ( λ) φ γ (3.8) 3.5. The GTP variance i given by: 5

7 ( ) ( ) ( γ) ( γ)( γ δ) ( γ)( γ)( γ γδ δ ) ( )( ) ( ) ( γ δ) γ( γ) ( δ) ( ) ( )( ) ( ) ( ) 3 4 ( + λ+ λ + λ + λ ) F A B C D E Variance GTP ( r; λ γ δ) = G where A= γ + γ + γ B= + + C= D= γ γ δ γ + γ γδ + δ 3 E= + + F = γδ G= + γ + γ + δ + δ γ λ δλ. (3.9) Equation (3.9) ha an unpleaant form but i given here a variance i an important meaure of ri in many financial application. Beide it i required in the denity-baed mean regreion approach (Konte 00). 4. Maximum Lielihood Etimator for TP and GTP Ditribution 4.. Although Kotz and van Dorp provided a maximum lielihood etimator for the TP ditribution they did not preent a imilar olution for GTP. Intead they propoed a recurive maximum lielihood procedure for the more general Uneven Standard Two-Sided (USTSP) ditribution. It appear however that uch a olution doe exit for GTP. The only retriction i that thi imple and traightforward etimator only concern unimodal denitie which in any cae eem to be the only intereting one in mot practical application. The olution for non-unimodal denitie i alo preent but require olving a nonlinear equation and more detailed conideration. At thi point the maximum lielihood etimator for TP and GTP in the cae of unimodal denitie are demontrated without giving any detail of how they were derived. Thee are provided in Point 5 together with the non-unimodal denity analyi. 4.. An important feature of both TP and GTP lielihood function i that they may have multiple maxima. That thee local maxima only appear at the ample point ait in finding the global maximum. Importantly the lielihood value at thee point can be calculated quite imply by uing a formula and do not require any optimization algorithm. Thi produce a very different approach than that commonly ued for mooth ditribution. The etimation procedure i baed on checing the value of the lielihood function at 6

8 the ample point and electing the greatet value. Knowing the point at which the lielihood function achieve it maximum allow the ought parameter to be derived once again by uing a formula without recoure to any optimization algorithm The maximum lielihood etimator in the cae of TP i preented in a lightly different manner than that adopted by Kotz and van Dorp. The etimator of the λ parameter i: λ= ˆ r (4.) where r denote the value of the th point from the ordered ample at which the log-lielihood function defined a: LogL = w ln w (4.) achieve it maximum 4. In (4.) denote the ample count and w = w + w (4.3) w w i i= = ln r r ( r) = ( r ) i + i= + ln (4.4) (4.5) where r i denote the value of the i th point from the ordered ample. Once the point i determined the maximum lielihood etimator of γ i given by: ˆ γ = (4.6) w where w i calculated at point In the cae of GTP the maximum lielihood etimator of the λ parameter i alo given by (4.) i.e. it i the value r of the th point from the ordered ample at which the loglielihood function defined a: ( ( ) ) LogL = w + ln ln r w + r w (4.7) achieve it maximum. In (4.7) w w and w + are defined a above. Once the point i 4 w NB: Kotz and van Dorp do not provide thi form and a different expreion viz. e i maximized. 7

9 determined the maximum lielihood etimator of γ and δ are given by: ˆ γ = ( r ) w w w + r (4.8) and ˆ δ = + r w w w + r (4.9) where the value of w GTP wa not provided by Kotz and van Dorp. w + and r are calculated at the point. A mentioned thi reult for 5. Derivation of Maximum Lielihood Etimator for TP and GTP 5.. Point 5 i devoted olely to the derivation of the maximum lielihood etimator for TP and GTP together with an analyi of non-unimodal denitie. Thi i a quite technical and detailed ubject and may therefore be ipped by reader more intereted in practical etimation reult. 5.. The derivation of the maximum lielihood etimator for TP will be preented firt. Thi i a different derivation than that preented by Kotz and van Dorp but it follow their line of reaoning. Let u conider the lielihood function in the interval between point and + from the ordered ample. According to (.) thi may be expreed a: L γ + P P = γ λ ( λ) γ (5.) where P i= i= + = ri and + P ( ) = ri. The interval between the lower bound of the ditribution (i.e. the value 0) and the firt ample point may be alo conidered. In thi cae = 0 P = and P = ( r) 0 0 i. Similarly for the interval between the lat ample point and the i= upper bound of the ditribution (i.e. the value ) = logarithm of (5.) reult in the log-lielihood function: = i i= P r ( )( ) ( ) ( ) ( ) and P + =. Taing the LogL = lnγ + γ ln P lnλ + γ ln P ln λ. (5.) The econd derivative of (5.) with repect to λ i: 8

10 dλ λ ( λ) d LogL = ( γ ) + (5.3) which i alway poitive for γ >. It follow that in the cae of unimodal denitie the loglielihood function i alway convex between the point and + and it reache it maximum value at one of thee point. Moving from one interval to another however change the log-lielihood function (5.) a it depend on. The corollary of thi i that the log-lielihood function i not differentiable at the ample point and local maxima may appear there. It i therefore ufficient to chec the value of the log-lielihood function at the ample point to find it global maximum. Thi can be done a follow. The maximum value of the log-lielihood function at point can be found by determining it firt derivative with repect to γ and comparing the reult to 0: d LogL dγ = + ln P ln r + ln P ( ) ln( r ) = 0 (5.4) γ where λ ha been ubtituted with r. Equation (5.4) may be repreented in a impler form uing (4.3) (4.5): d LogL dγ = w = 0. (5.5) γ Solving (5.5) lead to the maximum lielihood etimator of γ (4.6). Subtituting thi reult into the log-lielihood function (5.) and rearranging yield the form (4.) which doe not depend on the ditribution parameter and i therefore eay to maximize The procedure for GTP i imilar. According to (3.) and (3.) the lielihood function may be expreed a: L γ + γδ P P = λδ ( λ) γ + λ ( λ) δ. (5.6) Taing the logarithm of (5.6) yield the log-lielihood function: LogL { lnγ lnδ ln λδ ( λ) γ } = ( γ )( P λ) + ( δ ) P ( ) ( λ) ln ln ln ln. (5.7) It econd derivative with repect to λ i: ( ) ( )( ) ( ) ( ) + ( ) d LogL γ δ γ δ = + + dλ λ λ λδ λ γ (5.8) which i alway poitive for γ > and δ >. Applying the reaoning preented in point 5. for 9

11 TP it i ufficient to chec the value of the log-lielihood function at the ample point in order to find it global maximum. Thee can be found a follow. The firt derivative of the log-lielihood function (5.7) with repect to γ and δ are compared to 0: d LogL r δ = w = 0 dγ γ r δ + ( r ) γ d LogL r + = w = 0 dδ δ r δ + ( r ) γ (5.9) (5.0) where (4.4) and (4.5) were ued to implify the reult. Solving the et of equation (5.9) and (5.0) reult in the maximum lielihood etimator of γ (4.8) and δ (4.9). It i worth noting that the econd pair of olution alo exit with the root in the denominator being negative. Thi pair however lead to very high value of γ of δ and ditribution which do not fit the data. Subtituting the etimated value (4.8) and (4.9) into the log-lielihood function (5.7) reult in the form (4.7) which doe not depend on the ditribution parameter and i ued for maximization It i intereting (albeit from a theoretical rather than a practical point of view) to determine the maximum lielihood etimator for non-unimodal denitie. In thi cae the econd derivative of the log-lielihood function may aume a negative value. Thi would indicate that the log-lielihood function i concave and it maximum may be located between point and +. The TP etimator will be conidered firt In order to find a maximum the firt derivative of the log-lielihood function (5.) with repect to λ i compared to 0: d LogL = ( γ ) = 0 dλ λ λ (5.) which reult in: ˆ λ=. (5.) Equation (5.) determine the value of λ at which the log-lielihood function (5.) achieve it maximum. It can eaily be checed whether thi value i located between point and +. If not then the log-lielihood function achieve the maximum either at point or at point + a in the cae of unimodal denitie. If the value of λ determined by (5.) i located inide the conidered interval then the 0

12 log-lielihood value may be calculated a follow. Firt we define w w and w + more generally than (4.3) (4.4) and (4.5) a function of λ within the interval [ + ] rather than the value at point : ( λ) ( λ) ( λ) w = w + w (5.3) ri = ln λ i= ( ) w λ ( ) w λ ( r) = ( λ) i + i= + ln. (5.4) (5.5) Clearly (5.3) - (5.5) reult in the ame value a (4.3) - (4.5) for λ= r i.e. at the bound of the conidered interval [ + ]. Repeating the tep preented in point 5.. yield the loglielihood value (5.6) which ha a very imilar form to (4.): LogL = w( λ) ln. w( λ) (5.6) The only difference i that w i now determined for uing (4.3). λ= uing (5.3) rather than for λ= r The procedure preented above mut be repeated for all the interval [ + ] in order to find the global maximum of the log-lielihood function. Very imilarly to unimodal denitie the maximum lielihood etimator of γ i then given by: ˆ γ = (5.7) w ( λ) which i calculated uing (5.3) for the determined value of λ and The firt derivative of the log-lielihood function (5.7) with repect to λ i compared to 0: 5.5. The cae of the GTP i more complex. which reult in two olution: ( ) + ( ) + ( ) ( ) ( ) + d LogL λ γ λδ λ γ λδ = = 0 dλ λ λ γ λδ (5.8)

13 γ λ= γ + ( ) δ (5.9) γ λ=. (5.0) γ δ Here determining whether the obtained value of λ are located within the conidered interval [ + ] i not poible a it wa in the cae of TP becaue γ and δ are not nown. One poibility i therefore to maximize the log-lielihood function (5.7) and to chec the obtained value of λ. Such a procedure however i not o convenient becaue (5.7) ha 3 parameter The ta can be implified a follow. Calculating the firt derivative of the log-lielihood function (5.7) with repect to γ and δ reult in: d LogL λδ = w = dγ γ λδ + ( λ) γ ( λ) 0 d LogL λ + = w( λ) = 0 dδ δ λδ + ( λ) γ (5.) (5.) where (5.4) and (5.5) are ued to implify the reult (cf. (5.9) and (5.0)). Solving the et of equation (5.) and (5.) yield etimator of γ and δ a function of λ (cf.(4.8) and (4.9)) ˆ γ λ = ( ) w ( λ) + ( λ) w( λ) w( λ) λ (5.3) and ˆ δ( λ) =. + λ w( λ) w( λ) w( λ) + λ (5.4) (4.7)) Subtituting (5.3) and (5.4) to (5.7) yield the log-lielihood function (cf. ( ) LogL ( λ) = w( λ) + ln ln λ w( λ) ( λ) w( λ) (5.5) which i a function of only one parameter λ and i therefore much eaier to maximize than (5.7). A maximizing a function i equivalent to olving it firt derivative with repect to 0 the preented way of proceeding may be further invetigated Returning to the full form of the log-lielihood function by ubtituting (5.3) - (5.5) to (5.5) calculating it firt derivative with repect to λ then putting bac (5.3) -

14 (5.5) in order to implify the reult and comparing it to 0 yield: + w( λ) m+ w( λ ) ( + m) d LogL w( ) ( ) w( ) m λ λ λ λ = + = 0. dλ λ λ λw λ λ w λ ( ) + ( ) ( ) (5.6) Solving (5.6) with repect to λ and leaving all the w(λ) expreion on the right-hand ide give the following three nonlinear equation: λ= m w + w( λ) ( λ) + w( λ) + ( λ) ( λ) 4w w (5.7) λ= (5.8) ( λ) ( λ) + 4w w λ= + (5.9) which can be olved numerically with repect to λ. Each equation may have 0 or more olution The maximum of the log-lielihood function (5.7) in the conidered interval [ + ] may be any of λ value determined by olving (5.7) (5.8) and (5.9). The whole maximum lielihood procedure for non-unimodal denitie would therefore require the following tep: Solve (5.7) (5.8) and (5.9) numerically with repect to λ Chec whether any of the reulting λ value are located within the interval [ + ] Calculate γ and δ for the λ value() uing (5.3) and (5.4) Calculate the econd derivative of the log-lielihood function uing (5.8) and chec whether thi value i negative (a poitive value would indicate a minimum of the loglielihood function) If the above condition are atified then calculate the value of the log-lielihood function uing (5.5). If any one of them i not atified then the maximum i either at point or + a i the cae with unimodal denitie. Calculate the log-lielihood value at the bound of the interval uing (5.5) or (4.7). Choe the greatet value of the log-lielihood function. The decribed procedure mut be repeated for all interval [ + ] in order to find 3

15 the global maximum The procedure decribed in may be implified if olving (5.7) - (5.9) i replaced by maximizing (5.5). Such a procedure hould be more efficient a fewer equation and poible λ value have to be examined A preented in thi Point both the TP and GTP maximum lielihood etimator are fairly eay to calculate for unimodal denitie. The procedure for non-unimodal denitie i not o traightforward but can be eaily implemented by a oftware pacage. 6. Etimation of Lottery Reult 6.. The TP and GTP maximum lielihood etimator are here ued in an application conidered in a previou paper by Konte (00). Two data et are examined. Set - the experimental data preented by Traub and Schmidt (009) whoe reearch concerned the relationhip between WTP (Willingne to Pay) and WTA (Willingne to Accept). Set - the experimental data of Idziowa (009) whoe reearch concern the quetion of whether the form in which probability i preented ha any impact on the hape of the probability weighting function. 6.. Konte (00) propoed a two-tage regreion procedure to decribe the experimental reult. In the firt tage the relative certainty equivalent are determined uing the tranformation: ce A r= (6.) P A where r denote the relative certainty equivalent ce denote the certainty equivalent P = Max(x) = the maximum lottery outcome and A = Min(x) = the minimum lottery outcome. The relationhip (6.) enure that r aume value in the range [0 ]. The denitie of the relative certainty equivalent for given value of probability are then etimated. The remaining part of thi paper concentrate on thee etimation reult. In the econd tage the obtained denitie are ued to etimate the relationhip between the relative certainty equivalent and the probability of winning the main prize: ( ) r Q p = (6.) where p denote the probability of winning the prize and Q denote a relative utility function which ha the form of a cumulative denity function defined over the range [0]. Further detail are available in Konte (00). 4

16 6.3. Although the maximum lielihood etimator for TP and GTP preented in Point 4 are fairly imple and traightforward a few practical comment are in order. Firt the TP and GTP maximum lielihood etimator appear to be computationally extremely fat when compared with tandard etimator of mooth ditribution which require running optimization oftware. The difference may be a much a two order of magnitude when compared with the GBT etimation The calculation time can be even further reduced whenever there i a large ample whoe data aume a limited number of value. Thi i becaue it i not neceary to calculate the value of the log-lielihood function (4.) or (4.7) for all the ample point but only for ditinct point. In financial application however the ample value are uually all different and any attempt to reduce the number of point will only lengthen the procedure. Second the procedure preented here only hold for unimodal denitie. In any cae when γ or δ aume a value le than the obtained reult may not be the maximum lielihood etimator. Thee cae were imply excluded from further conideration in the practical tudy preented in thi Point. Thi i a ound approach when the reearcher ha reaon to believe that the ought ditribution i unimodal. A tated γ and δ value le than appeared for ample point located away from the mode and the reulting ditribution exhibited a U- hape pointing to the anti-mode of the denity. The log-lielihood function value were low in thee cae which made it natural to eliminate uch olution. Third olution where either γ or δ aumed a value of infinity were alo encountered. Thi i a corollary of w or w + auming a value of 0 (cf. (4.8) and (4.9)). Such olution were alo excluded from thi reearch. The motivation for excluding uch olution even if they provide the highet lielihood value may not be all that obviou. However an infinite value for either γ or δ reduce the GTP to a bounded TP defined over part of the interval [0] and the lielihood function value appear to be incommenurable with ditribution defined over the entire interval The TP maximum lielihood etimation reult for Set are preented in Figure 6.. Each box preent an etimation reult for a ingle probability. A can be een the TP which ha parameter perform better than the GBT which ha 3 parameter in 7 out of cae. Thi i indicated by the maximized value of the log-lielihood function. Even without checing thee value the TP curve appear to be better uited for uch peaed data than a mooth ditribution. 5

17 Figure 6.. Denitie of r for repective probabilitie in Set. The TP etimation reult are mared in red and the GBT etimation reult are mared with dahed red line. The tp and gbt name indicate the maximized value of the log-lielihood function. The parenthee contain the mode median and mean of the ditribution calculated uing the formula given in Point (for TP) and in Konte (00) for GBT. The preented etimation how that the data are poitively ewed for low probabilitie and negatively ewed for high probabilitie 6

18 6.5. Figure 6. preent the etimation reult for Set uing the GTP maximum lielihood etimator. Figure 6.. Denitie of r for repective probabilitie in Set. The GTP etimation reult are mared in red and the GBT etimation reult are mared with dahed red line. The gtp and gbt name indicate the maximized value of the log-lielihood function. The parenthee contain the mode median and mean of the ditribution calculated uing the formula given in Point 3 (for GTP) and in Konte (00) for GBT. 7

19 Here the GTP etimator outperform GBT in all cae. Thi how that GTP i much better uited to decribe peaed denitie than the commonly ued mooth ditribution. The trong correpondence of the mode with the probability value i epecially worth noting. Thi how that the mot liely relative certainty equivalent i roughly equal to the probability Figure 6.3 preent the etimation reult for Set uing the GTP maximum lielihood etimator. In thi cae however the GTP perform poorly compared with GBT with only one cae out of nine yielding a better reult than the GBT. Thi how that GTP i not a ditribution to be ued when the data are cattered. Figure 6.3. Denitie of r for repective probabilitie in Set. The GTP etimation reult are mared in red and the GBT etimation reult are mared in dahed red. The gtp and gbt value indicate the maximized log-lielihood function. The parenthee contain the mode median and mean of the ditribution. 8

20 6.7. In the econd phae of the regreion procedure decribed in 6.. the relationhip between relative certainty equivalent and probabilitie i etimated. Thi i done uing the denity function hown in Figure and 6.3. The reult are only briefly preented here a the regreion procedure digree from the main ubject of thi paper. Figure 6.4 preent the etimation of the relative utility function Q for Set uing the GTP denitie hown in Figure 6.. The mean quantile (including median) and mode (maximum lielihood) regreion etimation are hown on the one graph. A can be een the lottery valuation are only nonlinear with probability when median and mean are conidered. Such nonlinearity i not confirmed for mode. The parameter α and β are both equal to in thi cae which mean the relationhip i linear. It follow that the mot liely lottery valuation are cloe to their expected value and that the mot liely behavior of a group i fully rational. Figure 6.4. Etimation reult of the relative utility function Q for Set. MN Mean (Leat Square); Q Lower Quartile; MED Median; Q3 Upper Quartile; ML Maximum Lielihood (Mode). The orange area mar data between the lower and upper quartile. The functional form of Q i a cumulative beta ditribution with parameter α and β. 7. Concluion Thi paper preented etimation procedure for peaed denitie over the interval [0]. The two-ided power ditribution TP and the generalized ditribution GTP (Kotz and van Dorp 004) were conidered for thi purpoe. The paper preented the TP and GTP maximum lielihood etimator. The latter wa not derived by Kotz and van Dorp. The obtained etimation demontrated that GTP outperform generalized beta ditribution GBT (Libby Novic 98) in one of the conidered data et. However the reult for much cat- 9

21 tered data a preent in the econd et are oppoite. The obtained denitie are then ued to etimate the relative utility function (Konte 00). The previouly preented reult that lottery valuation are only nonlinear with probability when mean and median are conidered i confirmed. Such nonlinearity diappear for mode. The paper demontrate that TP and GTP ditribution may be extremely ueful in etimating peaed denitie over the interval [0] and in reearching the relative utility function. Reference:. Idziowa K. (009). Determinant of the probability weighting function preented under name Katarzyna Domurat at SPUSM Conference Rovereto Italy Augut 009 paper Id = 8 Konte K. (00). Denity Baed Regreion for Inhomogeneou Data; Application to Lottery Experiment MPRA Paper Available at SSRN: 3. Kotz S. van Dorp J. R. (004). Beyond Beta; Other Continuou Familie of Ditribution with Bounded Support and Application World Scientific Publihing Singapore. 4. Libby D. L. Novic M. R. (98). Multivariate generalized beta ditribution with application to utility aement Journal of Educational Statitic 7 pp Traub S. Schmidt U. (009). An Experimental Invetigation of the Diparity between WTA and WTP for Lotterie Theory & Deciion 66 (009) pp

Social Studies 201 Notes for November 14, 2003

Social Studies 201 Notes for November 14, 2003 1 Social Studie 201 Note for November 14, 2003 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

Social Studies 201 Notes for March 18, 2005

Social Studies 201 Notes for March 18, 2005 1 Social Studie 201 Note for March 18, 2005 Etimation of a mean, mall ample ize Section 8.4, p. 501. When a reearcher ha only a mall ample ize available, the central limit theorem doe not apply to the

More information

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang

ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION. Xiaoqun Wang Proceeding of the 2008 Winter Simulation Conference S. J. Maon, R. R. Hill, L. Mönch, O. Roe, T. Jefferon, J. W. Fowler ed. ON THE APPROXIMATION ERROR IN HIGH DIMENSIONAL MODEL REPRESENTATION Xiaoqun Wang

More information

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis

Source slideplayer.com/fundamentals of Analytical Chemistry, F.J. Holler, S.R.Crouch. Chapter 6: Random Errors in Chemical Analysis Source lideplayer.com/fundamental of Analytical Chemitry, F.J. Holler, S.R.Crouch Chapter 6: Random Error in Chemical Analyi Random error are preent in every meaurement no matter how careful the experimenter.

More information

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281

7.2 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 281 72 INVERSE TRANSFORMS AND TRANSFORMS OF DERIVATIVES 28 and i 2 Show how Euler formula (page 33) can then be ued to deduce the reult a ( a) 2 b 2 {e at co bt} {e at in bt} b ( a) 2 b 2 5 Under what condition

More information

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems

A Constraint Propagation Algorithm for Determining the Stability Margin. The paper addresses the stability margin assessment for linear systems A Contraint Propagation Algorithm for Determining the Stability Margin of Linear Parameter Circuit and Sytem Lubomir Kolev and Simona Filipova-Petrakieva Abtract The paper addree the tability margin aement

More information

Bogoliubov Transformation in Classical Mechanics

Bogoliubov Transformation in Classical Mechanics Bogoliubov Tranformation in Claical Mechanic Canonical Tranformation Suppoe we have a et of complex canonical variable, {a j }, and would like to conider another et of variable, {b }, b b ({a j }). How

More information

Optimal Coordination of Samples in Business Surveys

Optimal Coordination of Samples in Business Surveys Paper preented at the ICES-III, June 8-, 007, Montreal, Quebec, Canada Optimal Coordination of Sample in Buine Survey enka Mach, Ioana Şchiopu-Kratina, Philip T Rei, Jean-Marc Fillion Statitic Canada New

More information

Comparing Means: t-tests for Two Independent Samples

Comparing Means: t-tests for Two Independent Samples Comparing ean: t-tet for Two Independent Sample Independent-eaure Deign t-tet for Two Independent Sample Allow reearcher to evaluate the mean difference between two population uing data from two eparate

More information

Clustering Methods without Given Number of Clusters

Clustering Methods without Given Number of Clusters Clutering Method without Given Number of Cluter Peng Xu, Fei Liu Introduction A we now, mean method i a very effective algorithm of clutering. It mot powerful feature i the calability and implicity. However,

More information

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas) Lecture 7: Analytic Function and Integral (See Chapter 4 in Boa) Thi i a good point to take a brief detour and expand on our previou dicuion of complex variable and complex function of complex variable.

More information

1. The F-test for Equality of Two Variances

1. The F-test for Equality of Two Variances . The F-tet for Equality of Two Variance Previouly we've learned how to tet whether two population mean are equal, uing data from two independent ample. We can alo tet whether two population variance are

More information

Estimation of Current Population Variance in Two Successive Occasions

Estimation of Current Population Variance in Two Successive Occasions ISSN 684-8403 Journal of Statitic Volume 7, 00, pp. 54-65 Etimation of Current Population Variance in Two Succeive Occaion Abtract Muhammad Azam, Qamruz Zaman, Salahuddin 3 and Javed Shabbir 4 The problem

More information

Stratified Analysis of Probabilities of Causation

Stratified Analysis of Probabilities of Causation Stratified Analyi of Probabilitie of Cauation Manabu Kuroki Sytem Innovation Dept. Oaka Univerity Toyonaka, Oaka, Japan mkuroki@igmath.e.oaka-u.ac.jp Zhihong Cai Biotatitic Dept. Kyoto Univerity Sakyo-ku,

More information

Chapter 4. The Laplace Transform Method

Chapter 4. The Laplace Transform Method Chapter 4. The Laplace Tranform Method The Laplace Tranform i a tranformation, meaning that it change a function into a new function. Actually, it i a linear tranformation, becaue it convert a linear combination

More information

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...4

More information

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay

Reliability Analysis of Embedded System with Different Modes of Failure Emphasizing Reboot Delay International Journal of Applied Science and Engineering 3., 4: 449-47 Reliability Analyi of Embedded Sytem with Different Mode of Failure Emphaizing Reboot Delay Deepak Kumar* and S. B. Singh Department

More information

Control Systems Analysis and Design by the Root-Locus Method

Control Systems Analysis and Design by the Root-Locus Method 6 Control Sytem Analyi and Deign by the Root-Locu Method 6 1 INTRODUCTION The baic characteritic of the tranient repone of a cloed-loop ytem i cloely related to the location of the cloed-loop pole. If

More information

Beta Burr XII OR Five Parameter Beta Lomax Distribution: Remarks and Characterizations

Beta Burr XII OR Five Parameter Beta Lomax Distribution: Remarks and Characterizations Marquette Univerity e-publication@marquette Mathematic, Statitic and Computer Science Faculty Reearch and Publication Mathematic, Statitic and Computer Science, Department of 6-1-2014 Beta Burr XII OR

More information

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall

Suggested Answers To Exercises. estimates variability in a sampling distribution of random means. About 68% of means fall Beyond Significance Teting ( nd Edition), Rex B. Kline Suggeted Anwer To Exercie Chapter. The tatitic meaure variability among core at the cae level. In a normal ditribution, about 68% of the core fall

More information

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011

NCAAPMT Calculus Challenge Challenge #3 Due: October 26, 2011 NCAAPMT Calculu Challenge 011 01 Challenge #3 Due: October 6, 011 A Model of Traffic Flow Everyone ha at ome time been on a multi-lane highway and encountered road contruction that required the traffic

More information

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions

Molecular Dynamics Simulations of Nonequilibrium Effects Associated with Thermally Activated Exothermic Reactions Original Paper orma, 5, 9 7, Molecular Dynamic Simulation of Nonequilibrium Effect ociated with Thermally ctivated Exothermic Reaction Jerzy GORECKI and Joanna Natalia GORECK Intitute of Phyical Chemitry,

More information

Design By Emulation (Indirect Method)

Design By Emulation (Indirect Method) Deign By Emulation (Indirect Method he baic trategy here i, that Given a continuou tranfer function, it i required to find the bet dicrete equivalent uch that the ignal produced by paing an input ignal

More information

By Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago

By Xiaoquan Wen and Matthew Stephens University of Michigan and University of Chicago Submitted to the Annal of Applied Statitic SUPPLEMENTARY APPENDIX TO BAYESIAN METHODS FOR GENETIC ASSOCIATION ANALYSIS WITH HETEROGENEOUS SUBGROUPS: FROM META-ANALYSES TO GENE-ENVIRONMENT INTERACTIONS

More information

Factor Analysis with Poisson Output

Factor Analysis with Poisson Output Factor Analyi with Poion Output Gopal Santhanam Byron Yu Krihna V. Shenoy, Department of Electrical Engineering, Neurocience Program Stanford Univerity Stanford, CA 94305, USA {gopal,byronyu,henoy}@tanford.edu

More information

Jan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size

Jan Purczyński, Kamila Bednarz-Okrzyńska Estimation of the shape parameter of GED distribution for a small sample size Jan Purczyńki, Kamila Bednarz-Okrzyńka Etimation of the hape parameter of GED ditribution for a mall ample ize Folia Oeconomica Stetinenia 4()/, 35-46 04 Folia Oeconomica Stetinenia DOI: 0.478/foli-04-003

More information

Convex Hulls of Curves Sam Burton

Convex Hulls of Curves Sam Burton Convex Hull of Curve Sam Burton 1 Introduction Thi paper will primarily be concerned with determining the face of convex hull of curve of the form C = {(t, t a, t b ) t [ 1, 1]}, a < b N in R 3. We hall

More information

Nonlinear Single-Particle Dynamics in High Energy Accelerators

Nonlinear Single-Particle Dynamics in High Energy Accelerators Nonlinear Single-Particle Dynamic in High Energy Accelerator Part 6: Canonical Perturbation Theory Nonlinear Single-Particle Dynamic in High Energy Accelerator Thi coure conit of eight lecture: 1. Introduction

More information

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is EE 4G Note: Chapter 6 Intructor: Cheung More about ZSR and ZIR. Finding unknown initial condition: Given the following circuit with unknown initial capacitor voltage v0: F v0/ / Input xt 0Ω Output yt -

More information

One Class of Splitting Iterative Schemes

One Class of Splitting Iterative Schemes One Cla of Splitting Iterative Scheme v Ciegi and V. Pakalnytė Vilniu Gedimina Technical Univerity Saulėtekio al. 11, 2054, Vilniu, Lithuania rc@fm.vtu.lt Abtract. Thi paper deal with the tability analyi

More information

Problem Set 8 Solutions

Problem Set 8 Solutions Deign and Analyi of Algorithm April 29, 2015 Maachuett Intitute of Technology 6.046J/18.410J Prof. Erik Demaine, Srini Devada, and Nancy Lynch Problem Set 8 Solution Problem Set 8 Solution Thi problem

More information

Avoiding Forbidden Submatrices by Row Deletions

Avoiding Forbidden Submatrices by Row Deletions Avoiding Forbidden Submatrice by Row Deletion Sebatian Wernicke, Jochen Alber, Jen Gramm, Jiong Guo, and Rolf Niedermeier Wilhelm-Schickard-Intitut für Informatik, niverität Tübingen, Sand 13, D-72076

More information

Lecture 8: Period Finding: Simon s Problem over Z N

Lecture 8: Period Finding: Simon s Problem over Z N Quantum Computation (CMU 8-859BB, Fall 205) Lecture 8: Period Finding: Simon Problem over Z October 5, 205 Lecturer: John Wright Scribe: icola Rech Problem A mentioned previouly, period finding i a rephraing

More information

DIFFERENTIAL EQUATIONS Laplace Transforms. Paul Dawkins

DIFFERENTIAL EQUATIONS Laplace Transforms. Paul Dawkins DIFFERENTIAL EQUATIONS Laplace Tranform Paul Dawkin Table of Content Preface... Laplace Tranform... Introduction... The Definition... 5 Laplace Tranform... 9 Invere Laplace Tranform... Step Function...

More information

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis

Advanced Digital Signal Processing. Stationary/nonstationary signals. Time-Frequency Analysis... Some nonstationary signals. Time-Frequency Analysis Advanced Digital ignal Proceing Prof. Nizamettin AYDIN naydin@yildiz.edu.tr Time-Frequency Analyi http://www.yildiz.edu.tr/~naydin 2 tationary/nontationary ignal Time-Frequency Analyi Fourier Tranform

More information

Calculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient

Calculation of the temperature of boundary layer beside wall with time-dependent heat transfer coefficient Ŕ periodica polytechnica Mechanical Engineering 54/1 21 15 2 doi: 1.3311/pp.me.21-1.3 web: http:// www.pp.bme.hu/ me c Periodica Polytechnica 21 RESERCH RTICLE Calculation of the temperature of boundary

More information

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004

Lecture 21. The Lovasz splitting-off lemma Topics in Combinatorial Optimization April 29th, 2004 18.997 Topic in Combinatorial Optimization April 29th, 2004 Lecture 21 Lecturer: Michel X. Goeman Scribe: Mohammad Mahdian 1 The Lovaz plitting-off lemma Lovaz plitting-off lemma tate the following. Theorem

More information

Stochastic Neoclassical Growth Model

Stochastic Neoclassical Growth Model Stochatic Neoclaical Growth Model Michael Bar May 22, 28 Content Introduction 2 2 Stochatic NGM 2 3 Productivity Proce 4 3. Mean........................................ 5 3.2 Variance......................................

More information

Laplace Transformation

Laplace Transformation Univerity of Technology Electromechanical Department Energy Branch Advance Mathematic Laplace Tranformation nd Cla Lecture 6 Page of 7 Laplace Tranformation Definition Suppoe that f(t) i a piecewie continuou

More information

Standard Guide for Conducting Ruggedness Tests 1

Standard Guide for Conducting Ruggedness Tests 1 Deignation: E 69 89 (Reapproved 996) Standard Guide for Conducting Ruggedne Tet AMERICA SOCIETY FOR TESTIG AD MATERIALS 00 Barr Harbor Dr., Wet Conhohocken, PA 948 Reprinted from the Annual Book of ASTM

More information

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH

A BATCH-ARRIVAL QUEUE WITH MULTIPLE SERVERS AND FUZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Mathematical and Computational Application Vol. 11 No. pp. 181-191 006. Aociation for Scientific Reearch A BATCH-ARRIVA QEE WITH MTIPE SERVERS AND FZZY PARAMETERS: PARAMETRIC PROGRAMMING APPROACH Jau-Chuan

More information

[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY

[Saxena, 2(9): September, 2013] ISSN: Impact Factor: INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY [Saena, (9): September, 0] ISSN: 77-9655 Impact Factor:.85 IJESRT INTERNATIONAL JOURNAL OF ENGINEERING SCIENCES & RESEARCH TECHNOLOGY Contant Stre Accelerated Life Teting Uing Rayleigh Geometric Proce

More information

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays

Gain and Phase Margins Based Delay Dependent Stability Analysis of Two- Area LFC System with Communication Delays Gain and Phae Margin Baed Delay Dependent Stability Analyi of Two- Area LFC Sytem with Communication Delay Şahin Sönmez and Saffet Ayaun Department of Electrical Engineering, Niğde Ömer Halidemir Univerity,

More information

Lecture 10 Filtering: Applied Concepts

Lecture 10 Filtering: Applied Concepts Lecture Filtering: Applied Concept In the previou two lecture, you have learned about finite-impule-repone (FIR) and infinite-impule-repone (IIR) filter. In thee lecture, we introduced the concept of filtering

More information

in a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim

in a circular cylindrical cavity K. Kakazu Department of Physics, University of the Ryukyus, Okinawa , Japan Y. S. Kim Quantization of electromagnetic eld in a circular cylindrical cavity K. Kakazu Department of Phyic, Univerity of the Ryukyu, Okinawa 903-0, Japan Y. S. Kim Department of Phyic, Univerity of Maryland, College

More information

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014 Phyic 7 Graduate Quantum Mechanic Solution to inal Eam all 0 Each quetion i worth 5 point with point for each part marked eparately Some poibly ueful formula appear at the end of the tet In four dimenion

More information

Math Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK

Math Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK ppendix 5 Scientific Notation It i difficult to work with very large or very mall number when they are written in common decimal notation. Uually it i poible to accommodate uch number by changing the SI

More information

Theoretical Computer Science. Optimal algorithms for online scheduling with bounded rearrangement at the end

Theoretical Computer Science. Optimal algorithms for online scheduling with bounded rearrangement at the end Theoretical Computer Science 4 (0) 669 678 Content lit available at SciVere ScienceDirect Theoretical Computer Science journal homepage: www.elevier.com/locate/tc Optimal algorithm for online cheduling

More information

HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES

HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES 15 TH INTERNATIONAL CONFERENCE ON GEOMETRY AND GRAPHICS 0 ISGG 1-5 AUGUST, 0, MONTREAL, CANADA HELICAL TUBES TOUCHING ONE ANOTHER OR THEMSELVES Peter MAYRHOFER and Dominic WALTER The Univerity of Innbruck,

More information

Assessing the Discriminatory Power of Credit Scores under Censoring

Assessing the Discriminatory Power of Credit Scores under Censoring Aeing the Dicriminatory Power of Credit Score under Cenoring Holger Kraft, Gerald Kroiandt, Marlene Müller Fraunhofer Intitut für Techno- und Wirtchaftmathematik (ITWM) Thi verion: Augut 27, 2003 Abtract:

More information

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog

Chapter 2 Sampling and Quantization. In order to investigate sampling and quantization, the difference between analog Chapter Sampling and Quantization.1 Analog and Digital Signal In order to invetigate ampling and quantization, the difference between analog and digital ignal mut be undertood. Analog ignal conit of continuou

More information

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems

A Simplified Methodology for the Synthesis of Adaptive Flight Control Systems A Simplified Methodology for the Synthei of Adaptive Flight Control Sytem J.ROUSHANIAN, F.NADJAFI Department of Mechanical Engineering KNT Univerity of Technology 3Mirdamad St. Tehran IRAN Abtract- A implified

More information

Introduction to Laplace Transform Techniques in Circuit Analysis

Introduction to Laplace Transform Techniques in Circuit Analysis Unit 6 Introduction to Laplace Tranform Technique in Circuit Analyi In thi unit we conider the application of Laplace Tranform to circuit analyi. A relevant dicuion of the one-ided Laplace tranform i found

More information

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment

Advanced D-Partitioning Analysis and its Comparison with the Kharitonov s Theorem Assessment Journal of Multidiciplinary Engineering Science and Technology (JMEST) ISSN: 59- Vol. Iue, January - 5 Advanced D-Partitioning Analyi and it Comparion with the haritonov Theorem Aement amen M. Yanev Profeor,

More information

A Bluffer s Guide to... Sphericity

A Bluffer s Guide to... Sphericity A Bluffer Guide to Sphericity Andy Field Univerity of Suex The ue of repeated meaure, where the ame ubject are teted under a number of condition, ha numerou practical and tatitical benefit. For one thing

More information

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem

An Inequality for Nonnegative Matrices and the Inverse Eigenvalue Problem An Inequality for Nonnegative Matrice and the Invere Eigenvalue Problem Robert Ream Program in Mathematical Science The Univerity of Texa at Dalla Box 83688, Richardon, Texa 7583-688 Abtract We preent

More information

Alternate Dispersion Measures in Replicated Factorial Experiments

Alternate Dispersion Measures in Replicated Factorial Experiments Alternate Diperion Meaure in Replicated Factorial Experiment Neal A. Mackertich The Raytheon Company, Sudbury MA 02421 Jame C. Benneyan Northeatern Univerity, Boton MA 02115 Peter D. Krau The Raytheon

More information

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL =

μ + = σ = D 4 σ = D 3 σ = σ = All units in parts (a) and (b) are in V. (1) x chart: Center = μ = 0.75 UCL = Our online Tutor are available 4*7 to provide Help with Proce control ytem Homework/Aignment or a long term Graduate/Undergraduate Proce control ytem Project. Our Tutor being experienced and proficient

More information

Preemptive scheduling on a small number of hierarchical machines

Preemptive scheduling on a small number of hierarchical machines Available online at www.ciencedirect.com Information and Computation 06 (008) 60 619 www.elevier.com/locate/ic Preemptive cheduling on a mall number of hierarchical machine György Dóa a, Leah Eptein b,

More information

Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions

Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions Stochatic Optimization with Inequality Contraint Uing Simultaneou Perturbation and Penalty Function I-Jeng Wang* and Jame C. Spall** The John Hopkin Univerity Applied Phyic Laboratory 11100 John Hopkin

More information

online learning Unit Workbook 4 RLC Transients

online learning Unit Workbook 4 RLC Transients online learning Pearon BTC Higher National in lectrical and lectronic ngineering (QCF) Unit 5: lectrical & lectronic Principle Unit Workbook 4 in a erie of 4 for thi unit Learning Outcome: RLC Tranient

More information

March 18, 2014 Academic Year 2013/14

March 18, 2014 Academic Year 2013/14 POLITONG - SHANGHAI BASIC AUTOMATIC CONTROL Exam grade March 8, 4 Academic Year 3/4 NAME (Pinyin/Italian)... STUDENT ID Ue only thee page (including the back) for anwer. Do not ue additional heet. Ue of

More information

Chapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem

Chapter 5 Consistency, Zero Stability, and the Dahlquist Equivalence Theorem Chapter 5 Conitency, Zero Stability, and the Dahlquit Equivalence Theorem In Chapter 2 we dicued convergence of numerical method and gave an experimental method for finding the rate of convergence (aka,

More information

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002 Correction for Simple Sytem Example and Note on Laplace Tranform / Deviation Variable ECHE 55 Fall 22 Conider a tank draining from an initial height of h o at time t =. With no flow into the tank (F in

More information

Compact finite-difference approximations for anisotropic image smoothing and painting

Compact finite-difference approximations for anisotropic image smoothing and painting CWP-593 Compact finite-difference approximation for aniotropic image moothing and painting Dave Hale Center for Wave Phenomena, Colorado School of Mine, Golden CO 80401, USA ABSTRACT Finite-difference

More information

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries

Coupling of Three-Phase Sequence Circuits Due to Line and Load Asymmetries Coupling of Three-Phae Sequence Circuit Due to Line and Load Aymmetrie DEGO BELLAN Department of Electronic nformation and Bioengineering Politecnico di Milano Piazza Leonardo da inci 01 Milano TALY diego.ellan@polimi.it

More information

Week 3 Statistics for bioinformatics and escience

Week 3 Statistics for bioinformatics and escience Week 3 Statitic for bioinformatic and escience Line Skotte 28. november 2008 2.9.3-4) In thi eercie we conider microrna data from Human and Moue. The data et repreent 685 independent realiation of the

More information

The Use of MDL to Select among Computational Models of Cognition

The Use of MDL to Select among Computational Models of Cognition The Ue of DL to Select among Computational odel of Cognition In J. yung, ark A. Pitt & Shaobo Zhang Vijay Balaubramanian Department of Pychology David Rittenhoue Laboratorie Ohio State Univerity Univerity

More information

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE

MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE Proceeding 3rd Annual Conference IEEE/EMBS Oct.5-8, 1, Itanbul, TURKEY MULTI-LAYERED LOSSY FINITE LENGTH DIELECTRIC CYLINDIRICAL MODEL OF MAN AT OBLIQUE INCIDENCE S.S. Şeker, B. Yeldiren Boğaziçi Univerity,

More information

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD

SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD SMALL-SIGNAL STABILITY ASSESSMENT OF THE EUROPEAN POWER SYSTEM BASED ON ADVANCED NEURAL NETWORK METHOD S.P. Teeuwen, I. Erlich U. Bachmann Univerity of Duiburg, Germany Department of Electrical Power Sytem

More information

Chapter 12 Simple Linear Regression

Chapter 12 Simple Linear Regression Chapter 1 Simple Linear Regreion Introduction Exam Score v. Hour Studied Scenario Regreion Analyi ued to quantify the relation between (or more) variable o you can predict the value of one variable baed

More information

Unified Correlation between SPT-N and Shear Wave Velocity for all Soil Types

Unified Correlation between SPT-N and Shear Wave Velocity for all Soil Types 6 th International Conference on Earthquake Geotechnical Engineering 1-4 ovember 15 Chritchurch, ew Zealand Unified Correlation between SPT- and Shear Wave Velocity for all Soil Type C.-C. Tai 1 and T.

More information

Dynamic Behaviour of Timber Footbridges

Dynamic Behaviour of Timber Footbridges Contança RIGUEIRO MSc PhD Student EST-IPCB contança@et.ipcb.pt Dynamic Behaviour of Timber Footbridge Graduated in Civil Engineering in Univ. of Coimbra (1992). MSc, Univ. of Coimbra (1997). João NEGRÃO

More information

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach

Control of Delayed Integrating Processes Using Two Feedback Controllers R MS Approach Proceeding of the 7th WSEAS International Conference on SYSTEM SCIENCE and SIMULATION in ENGINEERING (ICOSSSE '8) Control of Delayed Integrating Procee Uing Two Feedback Controller R MS Approach LIBOR

More information

If Y is normally Distributed, then and 2 Y Y 10. σ σ

If Y is normally Distributed, then and 2 Y Y 10. σ σ ull Hypothei Significance Teting V. APS 50 Lecture ote. B. Dudek. ot for General Ditribution. Cla Member Uage Only. Chi-Square and F-Ditribution, and Diperion Tet Recall from Chapter 4 material on: ( )

More information

Rao Transforms: A New Approach to Integral and Differential Equations

Rao Transforms: A New Approach to Integral and Differential Equations Rao Tranform: A New Approach to Integral and Differential Equation Dr. Muralidhara SubbaRao (Rao) SUNY at Stony Brook, murali@ece.unyb.edu,, rao@integralreearch.net Rao Tranform (RT) provide a brand new

More information

USPAS Course on Recirculated and Energy Recovered Linear Accelerators

USPAS Course on Recirculated and Energy Recovered Linear Accelerators USPAS Coure on Recirculated and Energy Recovered Linear Accelerator G. A. Krafft and L. Merminga Jefferon Lab I. Bazarov Cornell Lecture 6 7 March 005 Lecture Outline. Invariant Ellipe Generated by a Unimodular

More information

List coloring hypergraphs

List coloring hypergraphs Lit coloring hypergraph Penny Haxell Jacque Vertraete Department of Combinatoric and Optimization Univerity of Waterloo Waterloo, Ontario, Canada pehaxell@uwaterloo.ca Department of Mathematic Univerity

More information

A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT

A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: CORRESPONDENCE: ABSTRACT A FUNCTIONAL BAYESIAN METHOD FOR THE SOLUTION OF INVERSE PROBLEMS WITH SPATIO-TEMPORAL PARAMETERS AUTHORS: Zenon Medina-Cetina International Centre for Geohazard / Norwegian Geotechnical Intitute Roger

More information

Singular perturbation theory

Singular perturbation theory Singular perturbation theory Marc R. Rouel June 21, 2004 1 Introduction When we apply the teady-tate approximation (SSA) in chemical kinetic, we typically argue that ome of the intermediate are highly

More information

Bayesian-Based Decision Making for Object Search and Characterization

Bayesian-Based Decision Making for Object Search and Characterization 9 American Control Conference Hyatt Regency Riverfront, St. Loui, MO, USA June -, 9 WeC9. Bayeian-Baed Deciion Making for Object Search and Characterization Y. Wang and I. I. Huein Abtract Thi paper focue

More information

Thermodynamics of cation ordering in karrooite (MgTi 2 O 5 )

Thermodynamics of cation ordering in karrooite (MgTi 2 O 5 ) American Mineralogit, Volume 8, page 37 37, 999 Thermodynamic of cation ordering in karrooite (Mg O 5 ) MARK S. GHIORSO,, * HEIONG YANG, AND ROBERT M. HAZEN Department of Geological Science, Box 353, Univerity

More information

An estimation approach for autotuning of event-based PI control systems

An estimation approach for autotuning of event-based PI control systems Acta de la XXXIX Jornada de Automática, Badajoz, 5-7 de Septiembre de 08 An etimation approach for autotuning of event-baed PI control ytem Joé Sánchez Moreno, María Guinaldo Loada, Sebatián Dormido Departamento

More information

PARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES

PARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES PARAMETERS OF DISPERSION FOR ON-TIME PERFORMANCE OF POSTAL ITEMS WITHIN TRANSIT TIMES MEASUREMENT SYSTEM FOR POSTAL SERVICES Daniel Salava Kateřina Pojkarová Libor Švadlenka Abtract The paper i focued

More information

To appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm

To appear in International Journal of Numerical Methods in Fluids in Stability analysis of numerical interface conditions in uid-structure therm To appear in International Journal of Numerical Method in Fluid in 997. Stability analyi of numerical interface condition in uid-tructure thermal analyi M. B. Gile Oxford Univerity Computing Laboratory

More information

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty

Simple Observer Based Synchronization of Lorenz System with Parametric Uncertainty IOSR Journal of Electrical and Electronic Engineering (IOSR-JEEE) ISSN: 78-676Volume, Iue 6 (Nov. - Dec. 0), PP 4-0 Simple Oberver Baed Synchronization of Lorenz Sytem with Parametric Uncertainty Manih

More information

DYNAMIC MODELS FOR CONTROLLER DESIGN

DYNAMIC MODELS FOR CONTROLLER DESIGN DYNAMIC MODELS FOR CONTROLLER DESIGN M.T. Tham (996,999) Dept. of Chemical and Proce Engineering Newcatle upon Tyne, NE 7RU, UK.. INTRODUCTION The problem of deigning a good control ytem i baically that

More information

Acceptance sampling uses sampling procedure to determine whether to

Acceptance sampling uses sampling procedure to determine whether to DOI: 0.545/mji.203.20 Bayeian Repetitive Deferred Sampling Plan Indexed Through Relative Slope K.K. Sureh, S. Umamahewari and K. Pradeepa Veerakumari Department of Statitic, Bharathiar Univerity, Coimbatore,

More information

Interaction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load

Interaction of Pile-Soil-Pile in Battered Pile Groups under Statically Lateral Load Interaction of Pile-Soil-Pile in Battered Pile Group under Statically Lateral Load H. Ghaemadeh 1*, M. Alibeikloo 2 1- Aitant Profeor, K. N. Tooi Univerity of Technology 2- M.Sc. Student, K. N. Tooi Univerity

More information

Critical Height of Slopes in Homogeneous Soil: the Variational Solution

Critical Height of Slopes in Homogeneous Soil: the Variational Solution Critical Height of Slope in Homogeneou Soil: the Variational Solution Chen, Rong State Key Laboratory of Coatal and Offhore Engineering & Intitute of Geotechnical Engineering, Dalian Univerity of Technology,

More information

The Laplace Transform

The Laplace Transform Chapter 7 The Laplace Tranform 85 In thi chapter we will explore a method for olving linear differential equation with contant coefficient that i widely ued in electrical engineering. It involve the tranformation

More information

arxiv: v1 [cond-mat.stat-mech] 22 Sep 2009

arxiv: v1 [cond-mat.stat-mech] 22 Sep 2009 Exact and limit ditribution of the larget fitne on correlated fitne landcape arxiv:99.3933v [cond-mat.tat-mech] 22 Sep 29 Kavita Jain, Abhihek Dagupta 2 and Gayatri Da 3 Theoretical Science Unit and Evolutionary

More information

Computers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order

Computers and Mathematics with Applications. Sharp algebraic periodicity conditions for linear higher order Computer and Mathematic with Application 64 (2012) 2262 2274 Content lit available at SciVere ScienceDirect Computer and Mathematic with Application journal homepage: wwweleviercom/locate/camwa Sharp algebraic

More information

Lecture 7: Testing Distributions

Lecture 7: Testing Distributions CSE 5: Sublinear (and Streaming) Algorithm Spring 014 Lecture 7: Teting Ditribution April 1, 014 Lecturer: Paul Beame Scribe: Paul Beame 1 Teting Uniformity of Ditribution We return today to property teting

More information

Convergence criteria and optimization techniques for beam moments

Convergence criteria and optimization techniques for beam moments Pure Appl. Opt. 7 (1998) 1221 1230. Printed in the UK PII: S0963-9659(98)90684-5 Convergence criteria and optimization technique for beam moment G Gbur and P S Carney Department of Phyic and Atronomy and

More information

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation

IEOR 3106: Fall 2013, Professor Whitt Topics for Discussion: Tuesday, November 19 Alternating Renewal Processes and The Renewal Equation IEOR 316: Fall 213, Profeor Whitt Topic for Dicuion: Tueday, November 19 Alternating Renewal Procee and The Renewal Equation 1 Alternating Renewal Procee An alternating renewal proce alternate between

More information

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS

CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS CHAPTER 8 OBSERVER BASED REDUCED ORDER CONTROLLER DESIGN FOR LARGE SCALE LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.1 INTRODUCTION 8.2 REDUCED ORDER MODEL DESIGN FOR LINEAR DISCRETE-TIME CONTROL SYSTEMS 8.3

More information

A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho

A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES. Sanghyun Cho A SIMPLE NASH-MOSER IMPLICIT FUNCTION THEOREM IN WEIGHTED BANACH SPACES Sanghyun Cho Abtract. We prove a implified verion of the Nah-Moer implicit function theorem in weighted Banach pace. We relax the

More information

Optimization model in Input output analysis and computable general. equilibrium by using multiple criteria non-linear programming.

Optimization model in Input output analysis and computable general. equilibrium by using multiple criteria non-linear programming. Optimization model in Input output analyi and computable general equilibrium by uing multiple criteria non-linear programming Jing He * Intitute of ytem cience, cademy of Mathematic and ytem cience Chinee

More information

Homework #7 Solution. Solutions: ΔP L Δω. Fig. 1

Homework #7 Solution. Solutions: ΔP L Δω. Fig. 1 Homework #7 Solution Aignment:. through.6 Bergen & Vittal. M Solution: Modified Equation.6 becaue gen. peed not fed back * M (.0rad / MW ec)(00mw) rad /ec peed ( ) (60) 9.55r. p. m. 3600 ( 9.55) 3590.45r.

More information