Cost and benefit including value of life and limb measured in time units

Size: px
Start display at page:

Download "Cost and benefit including value of life and limb measured in time units"

Transcription

1 Secial Worksho on Risk Accetance and Risk Communication March 26-27, 2007, Stanford University Cost and benefit including value of life and limb measured in time units Ove Ditlevsen and Peter Friis-Hansen Technical University of Denmark Nils Koels Allé, Building 403, DK 2800 Kgs. Lyngby, Denmark, Abstract: Key elements of the authors ork on money equivalent time allocation to costs and benefits in risk analysis are ut together as an entity. This includes the data suorted dimensionless analysis of an equilibrium relation beteen total oulation ork time and gross domestic roduct leading to the definition of the life quality time allocation index (LQTAI). A ostulate of invariance of the LQTAI allocates societal value in terms of time to avoid life shortening fatalities as ell as serious injuries that shorten the life in good health. The observed large difference beteen the actual costs used by the oner for economically otimizing an activity motivates a simle risk accet criterion suited to be imosed on the oner by the ublic. The aer ends ith an examle concerning allocation of economical means for mitigation of loss of life and limb on a ferry in fire. 1. Introduction An engineering activity is lanned to be realized ithin a statistically homogeneous economical region. By dividing all costs and benefits in the exected net cost equation by the average age er time unit over the oulation of the region the equation gets hysical dimension as time, and the equation becomes indeendent of local inflation an urchase oer. The equation may contain terms that reresent loss of life and limb, and ossibly also environmental damages. Such terms have in the ast been considered as intangibles causing them to be excluded from the cost-benefit analysis. Hoever, ideas based on macro economical reasoning have in the last decade oened ossibilities of making rational evaluations (e.g. by use of the life quality index (LQI) defined in (Nathani, Lind, & Pandey 1997, Pandey, Nathani, & Lind 2006) and extensively studied and alied in (Rackitz 2002, Rackitz, Lentz, & Faber 2005), or, for environmental damage, the Nature reservation illingness index defined in (Friis-Hansen & Ditlevsen 2003)). In the time formulation the life and limb losses may at a first glance simly be ritten as the increment of the exected life time in good health caused by the loss giving accident. Hoever, this ould be an oversimlification because a art of the loss is ork time of larger societal value than the free time. The correction can be made by use of the criterion of invariance of the LQI, or, directly in time units, by use of invariance of the life quality time allocation index (LQTAI) defined in (Ditlevsen & Friis-Hansen 2005, Ditlevsen & Friis-Hansen 2007, Ditlevsen 2007). The results may be slightly different because the LQTAI is an extended more general version of the LQI. This aer gives a short recaitulation of the authors thoughts behind the LQTAI including its emirical suort. To the authors surrise these thoughts and the suorting emirical findings have turned out to be controversial. For large rojects the decision making by the oner is restricted by ublic requirements. The oner is rimarily focused on maximizing the rofit only including the direct costs to be sent on insurance remiums, uncovered losses and damage comensations. These costs are often much less than the societal value of life and limb as obtained from the invariance of the LQI or the LQTAI. Therefore the society must consider the ossibility of this larger loss and rotect itself under the consideration that the society has a ositive 1

2 interest in the realization of the roject. Rational reasoning leads to a ublic accet criterion formulated in (Friis-Hansen & Ditlevsen 2003, Ditlevsen 2003) and is recaitulated herein. Finally the aer gives an examle of using the LQTAI to assess the exected societal time value loss of life and limb due to a fire on a ferry. 2. Standard exected loss function examle The standard examle of technical decision roblem is related to an oeration (shi transort, factory, structure, etc.) that in the mean roduces a net gain g er time unit facing the fact that n different indeendent categories of cost generating accidents occur in time at random time oints that for the ith category is in accordance ith a homogeneous Poisson rocess of intensity λ i, i = 1,..., n. Denoting the mean of the random cost of an accident in category i by µ i, the exected loss function related to the oeration becomes asymtotically as t, (Ditlevsen 2003) L(λ 1,..., λ n, µ 1,..., µ n ) = c(λ 1,..., λ n, µ 1,..., µ n ) + 1 γ n λ i µ i g γ i=1 (1) here c(λ 1,..., λ n, µ 1,..., µ n ) is the initial cost invested at time zero and γ is the interest rate. Of dimensional reasons it follos that the initial cost function must be linear in µ 1,..., µ n. The unit of (1) may be the currency of the relevant country, and as such the amounts are time deendent on inflation and current urchase oer of the currency in the actual country. Thus the µ-values must be corrected to equivalent values at time zero. This correction is made automatic if the current values are divided by the average current age (salary) er time unit over the entire oulation of the country. Thus the mean losses µ i and the gain g become exressed as the time it takes to earn the amounts by the average age. To the first aroximation these amounts measured in time units become indeendent of time and to some degree of the same size from country to country. The oner of the oeration then acts otimally if the oeration is designed so that λ 1,..., λ n, µ 1,..., µ n get values that minimize the loss function (1). Even though some of the µ-values may reresent values of life and limb these values are not intangibles from the oners oint of vie. These values are simly those that are imlied by common insurance remiums and damage comensation ractice. These oner s values are usually considerably less than the societal values obtained from life quality index considerations (LQI) as first advocated by Nathani, Lind, and Pandey a decade ago (Nathani, Lind, & Pandey 1997), and since then further develoed by the same authors (Pandey, Nathani, & Lind 2006) and by Rackitz (Rackitz 2002, Rackitz, Lentz, & Faber 2005). The observation that the societal loss value is much larger than the oners loss value related to some accident causes the society to set u certain restrictions on the oner s activity to ensure that the societal loss is fully comensated by the benefit the society obtains from the oner s activity. A rational ublic accetance rule based on the concet of oner s value versus societal value is formulated in (Ditlevsen 2003, Friis-Hansen & Ditlevsen 2003), and is recaitulated in Section 6. Having a ish to better understand the first rinciles behind the LQI formulation, and to see to hich degree the LQI may be emirically justifiable and to hich degree it is based on a ostulate of consistent behavior, the authors came to a slightly different and more general societal life quality model based on dimensionless formulation of reference equilibrium beteen ork time and free time (leisure time). The formulation of this model is in short form recaitulated in the next section based on the authors aers (Ditlevsen & Friis-Hansen 2005, Ditlevsen & Friis-Hansen 2007, Ditlevsen 2007). 2

3 3. Preference equilibrium beteen ork time and free time in good health This section closely follos the text in (Ditlevsen 2007) ith a minor but imortant generalization resent in the earlier aer (Ditlevsen & Friis-Hansen 2005) but not exlicitly addressed in the cited aers on the LQI. The model formulation in Section 2 of (Ditlevsen & Friis-Hansen 2005) starts ith random variable modeling of the relevant quantities for the individual member of the oulation hereuon exectations are taken over the oulation. Herein as in (Ditlevsen & Friis-Hansen 2007, Ditlevsen 2007) the modeling starts directly on the level of mean values over the entire oulation. The slight generalization is that only a fraction r of the time unit reresents time in good health hile the fraction 1 r reresents time ith illness or bad health, that is, time that cannot be used as (moneyroducing) ork time. Thus r here is the fraction of the time unit used for ork, and the ratio /r 1 is the fraction of the time in good health used for ork Time value modeling An increment d has a monetary value directly roortional to the monetary value of the roduction in the time increment, of course. Hoever, hat herein is called exerienced orth is reasonably defined as the dimensionless ratio of this increment value to the value of the roduction in the entire ork time fraction. This definition of exerienced orth of ork time is simly based on the sychological feeling that a small increment of much has less orth than a small increment of little. The definition is consistent ith the emirical Weber-Fechner la that states that there is a logarithmic relation beteen the hysical magnitudes of stimuli and the erceived intensity of the stimuli (source htt://en.ikiedia.org/iki/weber- Fechner la). Thus the exerienced orth of the increment d is measured relatively as d( )/, here is a factor defined as the ratio beteen the total yearly GDP (Gross Domestic Product) G ithin a geograhical region such as a country and the total yearly ork based salary S aid out to the inhabitants of the region. The factor may be interreted as an enhancement factor that increases (or decreases) the societal value of the ork time fraction from to. The corresonding money value is then G = S = (S/), here S/ is the salary er ork time unit. It is seen that aears as an overall (or equivalent ) roductivity that enhances S to G. This enhancement includes all ossible sources to G, that is, also those sources not directly related to the used manoer. Basically G is a measure of the economical ealth of the society and S is a measure of ho much of this economical ealth is distributed directly to the orking oulation. In the folloing ill be denoted as the time equivalent roductivity realizing that the ord roductivity is used differently in roduction economics. More details are discussed in (Ditlevsen 2007). Similarly it is assumed that the exerienced orth of the free time in good health is measured relatively as d(r )/(r ) hich for a constant r is the same as d(1 /r)/(1 /r). Note that d(1 /r)/(1 /r) is negative for a ositive ork fraction increment d. Therefore a balance beteen d()/ = d(/r)/(/r) and d(1 /r)/(1 /r) may exist for hich constant orth is maintained. This consideration leads to the balance equation c d( /r ) + (1 c) d( 1 /r ) = 0 (2) /r 1 /r for some constant c beteen 0 and 1. More correctly, c should be ritten as a function c(r) of r, but for simlicity of riting (r) is suressed herever it is not needed exlicitly in the folloing. For exlanation of (2) it is noted that any relation beteen differentials er definition is linear in the differentials. That the necessarily ositive coefficients sum to one is just the result of a convenient normalization obtained by dividing the equation by the sum of the coefficients. Assume no that the ork market tries to maintain this balance hatever the values of and and assume that locally may be a function of (or, mathematically equivalent, that locally may be a function of ). 3

4 Then (2) reduces to the simle differential equation ( c /r 1 c ) d(/r) = c d (3) 1 /r beteen /r and. It is seen that d = 0 for /r = c and an elementary analysis of the variation of the sign of d ith /r in the neighborhood of /r = c shos that is minimal for /r = c. The general solution to (3) is = min (1 /r) 1 1/c (/r) 1 (1 c) 1 1/c c 1 (4) here min = min Data The relation (4) is suorted by data from different OECD countries. The numerical data for Denmark are given in full ith detailed exlanations in (Ditlevsen & Friis-Hansen 2005), and for the other OECD countries in diagrammatic form in (Ditlevsen & Friis-Hansen 2007) all based on the data in (OECD 2004) Denmark ADAMS Denmark ADAMS 0.16 Denmark ADAMS; c =0.084; K = = ex[!0.044 (y!1948)] year Germany 0) year Germany Germany; c =0.080; K = )15 0)14-0)080.0)12& ex23!0)0&4 45!1" )13 0)12 0)11 0) )0" ) "50 1"60 1"&0 1"80 1"" ear Canada year Canada Canada; c =0.095; K = = ex[!0.041 (y!1948)] year year Figure 1. To: Data oints for Denmark in the years Recorded ork time fraction and time equivalent roductivity (ratio of salary to GDP er same time unit). The smooth curves are jointly least square fits of (5) (left diagram) and (4) (for r = 1) ith (5) substituted (middle diagram) and are consistent ith the curve in the right diagram for (, ) ( min = 1.81, c = 0.084). Middle: Data for Germany in the years ( min = 1.80, c = 0.080). Bottom: Data for Canada in the years ( min = 1.89, c = 0.095). 4

5 The to first to diagrams in Fig. 1 sho the data for and as functions of the year for Denmark in the eriod from 1948 to The data oints for the air (, ) are shon in the third diagram. The largest value of corresonds to year 1948, and the connecting lines beteen the data oints sho the succession of the oints ith the year. The second and third ro of diagrams sho the corresonding data time series for Germany and Canada only reaching back to 1960, hoever, and therefore giving less strong suort than the Danish data. The smooth curves in the left diagrams are fits of the decreasing exonential function = c + a e b(t 1948) c as t (5) adoting the hyothesis that the asymtotic value as t is equal to the value of at hich () = (/1, 1) takes its minimal value min. The unknon arameters a, b, c, min are estimated by a simultaneous least square fit of (4) and (5) to the to data sets for (t, ) and (t, ). The squared deviations are divided by the average squared deviation for each of the to data sets to bring the squared deviations into the same size scale. In (Ditlevsen 2007) the time ordered Danish data samles of the ork time ratio and the time equivalent roductivity are analyzed don to details that make it ossible to formulate a hite noise driven stochastic rocess model for the joint time develoment of the to dimensionless quantities. It is remarkable that the data residuals ith resect to the redicted relation as function of time sho stationary random behavior in all essentials. Both rocesses are ell described as autoregressive rocesses in discrete time. Simulations ith the model sho samle curves of great similarity to the data samle curves. Thus all significant systematic behavior of the Danish data is catured by the relation (4) for r = 1 and the exonential time decay (5). It is seen in all three ros of diagrams in Fig. 1 that the Danish and the German data have reached a statistically stationary state after about 1975 hile a stationary state is reached before 1960 for the Canadian data. As judged by the eye, the diagrams in (Ditlevsen & Friis-Hansen 2007) for 18 OECD countries reveal that only the data for Iceland, Ireland, Jaan and Ne Zealand seem to be deviate significantly from the redictions of the model The time reduction factor r and the ealth index I r As it is seen from (4), the time reduction factor r simly defines an affinity maing that moves all data oints (, ) to (, /r). The grah of (, ) defined by (4) for r = 1 mas to a curve that deviates from the grah of (4). Hoever, a simle test calculation as done in (Ditlevsen & Friis-Hansen 2007) shos that the deviation in beteen the to grahs for fixed ithin the relevant value domain of is at most about 0.5% for r > 2/3. Since the oint (c, min ) mas to the oint (c/r, min ) it follos that c(r 0 ) = c/r 0 for the actual value r 0 of r. In rincile the value of r 0 can be estimated through a oulation health investigation. A value of about 0.9 may be a reasonable assessment. The deviations beteen the maed curve and the curve defined by (4), both corresonding to r and coinciding in the oint of minimum of, are negligible as comared to the random fluctuations of the data oints along the curves. Therefore it is not ossible to estimate the value of r from the considered data. Hoever, the data fits indicate that it ould be in conflict ith the model to claim that r has been varying substantially during the observation eriod. If the differential on the left side of (2) takes a value different from zero it is a sign of unbalance beteen the ork effort and the available free time. For a ositive value of the left side the total societal orth (or societal living standard) is increasing, and for a negative value the orth is decreasing. A basic non-dimensional ealth index I r may therefore be defined relatively by di r I r = V (r) [c d( /r ) + (1 c) d( 1 /r ) ] /r 1 /r (6) 5

6 in hich V (r) is some differentiable function of r. The fraction di r /I r is taken to be a measure of the total exerienced orth of ork time and free time together. Integration of (6) gives the general solution I r = J(r) [ ( ) c ( 1 ) ] 1 c V (r) r r { [ ( = J(r) ex V (r) c log ) ( + (1 c) log 1 )]} r r (7) in hich J(r) is another differentiable function of r. In the final definition of the ealth index I r the to functions V (r) and J(r) ill both be set to 1, relying on an argument of slo variation. 4. Invariance rincile Assume that some increment d(r E) of the exected life in good health at birth is obtained by some risk reducing action. The art (/r) d(r E) is used for ork and has the exerienced orth (/r) d(r E)/ (E) = d(r E)/(r E), and the rest is free time that has the exerienced orth (1 /r) d(r E)/[(1 /r)r E] = d(r E)/(r E). According to the combination rule (6) the total exerienced orth of the time increment d(r E) is then exerienced orth of d(r E) = V (r) d(r E) r E [ dr = V (r) r + de ] E The invariance rincile ostulates that the amount d taken from to ay for the risk reducing activity that roduces the value in (8) should be such that di r /I r under constant ork effort is balanced by the exerienced orth of d(r E). By differentiation of (7) under constant it follos that di r I r = V (r) c d + V (r) { ( ) c log r 1 c r + r(r ) + V (r) [ ( c log ) ( + (1 c) log 1 )] + J } (r) dr (9) V (r) r r J(r) here means derivative ith resect to r. Making the reasonable assumtion that the to functions V (r) and J(r) are sufficiently sloly varying that the terms ith V (r)/v (r) and J (r)/j(r) can be neglected, it follos by adding (8) and (9) and setting to zero that the amount d taken from given by (8) d = 1 c de E + [ 1 c r + c log ( r )] dr c (10) balances the gain by obtaining the increments de of exected life and dr of healthy life. For the actual value r 0 of r and the stationary state here = c(r 0 ) = c(1)/r 0 and = min, (10) becomes d = r [ 0 de 1 min c(1) E + c(1) + r ( )] 0 min c(1) c(1) c (r 0 ) log dr (11) r 0 c(1) The sum of the to equations (8) and (9) leading to (10) can be obtained directly as d(i r r E)/(I r r E) = 0 under constant. The roduct I r r E (setting V (r) = 1 and J(r) = 1) is called the Life Quality Time Allocation Index (LQTAI) (here revised relative to (Ditlevsen & Friis-Hansen 2005) by the factor r). To aly (11) for value assignment to dr an assessment of the value of the derivative c (r 0 ) is needed. Data suort is hardly obtainable. Thus it is necessary to make some normative choice based on some areciable 6

7 assumtion. Such an assumtion is used in (Ditlevsen & Friis-Hansen 2007) setting c(r) = c/r and obtained as seen above by ostulating that is not influenced by variations of r (imlying that the GDP and the total salary changes in the same ace ith the variation of r). For the urose to obtain c (r 0 ) the ostulate may be relaxed to be valid only for small variations of r in the vicinity of r 0. Acceting this ostulate the factor to the logarithm in (11) becomes r 0 c(1) c (r 0 ) = 1 r 0 (12) The to terms on the right side of (11) are the fractions of min that the oulation should agree that the society invests er erson for removing the accident source that causes the exected life to decrease by de and the exected life in good health to decrease by E dr, resectively. Dividing by the mean accident rate er erson gives the equivalent time values. After multilication by the salary er time the amounts are obtained in the actual money unit. The first amount is the so-called Imlied Cost of Averting a Fatality (ICAF). The second may be called the Imlied Cost of Averting an Injury (ICAI). To aly these results in a decision roblem for a secific roject, the increments de and dr must be assessed in relation to the roject. This assessment is made by use of engineering modeling of the consequences of the ossible accidents that are taken into account in the cost-benefit analysis. Examles are given in (Ditlevsen 2003, Ditlevsen 2004). Even though it is based on an idea of logically consistent evaluation of exerienced orth, the invariance rincile alied to the ealth index I r to obtain ICAF and ICAI is a normative rincile because it can only be suorted by subjective oinion. Of logical and ethical reasons it is argued that it is the right rincile to aly, but the illingness to ay and current ractice may be different. The ay out is to interret the ICAF and the ICAI as societal values that should enter the ublic accet rules as formulated in Section Public accetance rules To rotect the elfare of its citizens against undue damage and exloitation the ublic sets certain restrictions on otherise free money making activities. The decisions related to the design and oeration of the technical installations needed for these activities are in a free society left to the oners ho aim for otimizing their gains er time unit. Their cost- benefit analyses tyically only include the direct costs to be sent on insurance remiums and damage comensations. These costs are often much less than the societal value of life and limb as obtained in Section 4. Therefore the society must consider the ossibility of this larger loss and rotect itself under the consideration that the society has a ositive interest in the realization of the roject. A ay to ensure this is to require that the comany tax yield to the society is sufficiently large to cover the loss of the society in excess of the oners direct comensation after the occurrence of the damage. The concet of society is indeendent of country borders imlying that in the modeling it is not imortant hether the tax is aid in the one or the other country ithin a region that may embrace several countries. The comany tax rate as ell as other arameters of the model should be considered as local or regional assessment arameters that may vary from region to region around the orld. Accordingly the model defines a rocedure that in the average over all kinds of risky oerations is a guide to formulate ublic accetance criteria. Referring to (1), the net gain rate after loss due to the adverse events is n i=1 (g i λ i µ oi ) here µ oi is the oners mean cost value of the adverse event and g i is the gain er time unit, both associated to the ith adverse event category. The gain g i is reasonably defined so that g 1 /µ o1 = g 2 /µ o2 =... = g n /µ on and n i=1 g i = g. With a comany tax rate of ρ, the requirement of full comensation for damage in any single category then imlies that the oeration is only acceted by the ublic if (g i λ i µ oi )ρ λ i µ i here µ i is 7

8 the societal cost value in excess of the oner s exected cost value µ oi. The accetance rule can be ritten on the dimensionless form ( λ i µ oi ) µ 1 i (13) g i ρ µ oi It is interesting to note that accetance by requiring that the LQI (or the LQTAI) does not decrease by the lanned activity is a milder requirement than (13). In fact, it is shon in (Ditlevsen 2003) that the accetance criterion becomes (13) ith ρ = A model for determining de and dr generated by an accident As a reasonable simlification it is assumed that the robability that any erson exeriences at most one serious injury due to an accident during the life time is close enough to 1 that the ossibility of more than one serious injury er erson can be neglected. Moreover, it is assumed that the time from birth to the injury or the death haens can be modeled as the random variable min{x, Y } here X is an exonentially distributed random variable of mean κ 1 and Y is the random time to natural death (that is, the life time in case the source of the considered fatal event is removed). This model as considered in (Ditlevsen 2003) and generalized slightly in (Ditlevsen 2004) for the time to death. Here it is used for the time to the occurrence of serious injury also. To distinguish beteen death and injury the life time L and the life time in good health L h are modeled as { { min{x, Y } R min{x, Y } both ith robability P f L = and L h = (14) Y RY T 1 X<Y both ith robability 1 P f resectively, here R is the random reduction factor, T RY is the time it takes to recover from the injury, P f is the robability that the accident is fatal given that it occurs, and 1 X<Y is the indicator function for the event X < Y. Without deeer investigations it is assumed that R is uncorrelated ith both max{x, Y } and Y (not excluding stochastic deendence due to a ossible nonlinear relation beteen R and Y ) and that X and Y are indeendent. Then the exectations become, (Ditlevsen 2003), E[L] = P f 1 κ 0 (1 e κt ) f Y (t) dt + (1 P f )E[Y ] (15) E[L h ] = re[l] (1 P f )E[T ]E[1 X<Y ] = (r + dr)e (16) here r = E[R], E = E[L] and here T is assumed to be indeendent of X and Y. The last equality is according to the definition of the mean life in good health. It follos from (15) that E[Y ] E[L] 1 P κe[l] 2 f 2 (1 + V Y 2 ) or de E P 1 f 2 (1 + V Y 2 ) κ E (17) as κ 0. V Y ( 0.2) is the coefficient of variation of Y. This is a simle modification of the result in (Ditlevsen 2003) corresonding to P f = 1 and roved in the same ay. Moreover it follos from (16) that dr = (1 P f )κe[t ] (18) using that E[1 X<Y ] = P(X < Y ) = κ E, (Ditlevsen 2003). 8. Distribution on accident categories The results (17) for de and (18) for dr are for a single accident category. Assume that it is ossible to categorize all essential accidents into n different categories and that the oulation of size N can be assigned to these categories so that N i ersons are assigned to category i = 1,..., n. By assumtion 8

9 N N n = N since no erson can belong to to or more categories. The total -allocation er time unit to accidents is then N d = N 1 ( d) N n ( d) n here ( d) i = r [ κ i c(1) 2 (1 + V Y 2 ) E min P fi + c(1) + 1 ( )] r0 c(1) log E[T i ] min (1 P fi ) (19) r 0 min c(1) according to (10), (17) and (18). Since κ i is the mean fraction er time unit of the oulation associated to category i that suffers from an accident, each erson associated to the ith category gets allocated a time loss equal to the right side of (19). 9. Examle: Ferry on fire This examle illustrates the socio-economic assessment of life and limb by alication to a risk analysis of a fire on a RoRo ferry carrying 150 assengers inclusive cre. Although fires reresent a serious hazard for shis, many studies have shon that fire is the second largest hazards for cre and assengers on shis. Foundering due to collision, grounding or hull structural failure is generally considered being more dominant hazards. Nonetheless, for illustrative uroses the focus is here on assessing the number of injuries, the related duration of injury, and the number of fatalities. The most critical hase, as regards fatalities and injuries of cre and assengers, is not measured directly after the occurrence of the incident. In the hase that follos, the catain may try to revent the incident from evolving into a serious accident by, for examle, intentionally beaching a shi that takes in ater and thereby try to kee it from sinking. Hoever, if such risk mitigating measures fail, evacuation rovides the last oortunity to minimize injuries and fatalities. When an evacuation is initiated the evacuation erformance is very imortant. Clearly an orderly and timely evacuation is crucial. To model the critical situation, a Bayesian Netork construct of the fire scenario is very useful. The netork constructed for this examle considers initiation of fires in the accommodation areas and in the engine room. Also the escalation of the fires is included as a model element. Simulated evacuation times reorted in (Vanem & Skjong 2006) are used. The number of fatalities in the fire scenario is assessed in deendence of the available evacuation time. Similarly, the robability of a erson being injured during the evacuation is modeled artly as a function of the seriousness of the incident and of the available evacuation time. The seriousness of the occurring injuries in deendence of the available evacuation time is assessed on a subjective basis although suorted by the findings in (Hansen & Vinter 2003). Results: The annual frequency of fire in Scandinavian aters is estimated to κ = y 1. Given the occurrence of a fire estimates of the robability of a fatal accident and of injury is P f = and P i = 0.214, resectively. Thus, the robability is that the fire ocurrence does not harm an arbitrary assenger or cre member. Conditional on injury the exected recovery time is estimated to E[T injury] = 0.27 y. Hence, E[T ] = ( )/( ) = y = 22 days. Inserting these values into (19) and using c = 0.084, min = 1.81, r 0 = 0.95, E = 80 y, and V Y = 0.2, the to terms in (19) become d f /κ = 6.68 y and d i /κ = 1.34 y. Multilying these by S = G D P/ min = 33, 340/1.81 = 18, 420 /y their monetary equivalents are mill and mill, resectively. Thus the societal cost of an injury is about 20% of the cost of a fatality. For N = 150 assengers the total exected socio-economic loss of life and limb er fire becomes 150 ( ) = 22.2 mill. The oner loss exectations are estimated to / fatality and / injury. The exected oner loss due to fatalities and injuries caused by a fire occurrence is thus µ o = ( ) 150 = mill. The total annual gain before tax for the RoRo-ferry is mill / y of hich 30% is allocated to cover the fire risk. The to sides of the accetance inequality (13) then become left side: 1.2 ( = , right side: ) ρ ρ=0.3 = (20) Thus the ublic accetance criterion (13) is satisfied. For ρ = 1 the criterion is the so-called LQI (or LQTAI) criterion hich is less restrictive than if ρ < 1. For ρ = 1 the right side of (13) becomes

10 Conclusions A recaitulation is given of the reviously ublished dimensionless analysis of the equilibrium relation beteen the gross domestic roduct, the salary, and the ork time all on average over the entire oulation of a macro economical homogeneous geograhical region. The analysis suggests to measure all monetary quantities in the unit of time and thereby to eliminate the influence of inflation, the currency unit, and the local urchase oer. The theoretically obtained ork and leisure time distributing equilibrium rincile has an emirical back-u that makes it reasonable also to use this equilibrium rincile to define a normative rule for allocating time to life and limb saving measures in risky activities. The allocation rule is generated from invariance of the life quality time allocation index (LQTAI) derived conveniently from the dimensionless analysis behind the equilibrium rincile. The insiration to this construct comes decisively from a study of the literature on the life quality index (LQI) mainly authored by Nathani, Lind, Pandey, and Rackitz, and a concern to get a deeer understanding of hat the basis is for the LQI. The result of the study is a more general mathematical formulation than that of the LQI (in different versions). If the macro economy including the ork time fraction is long time stationary the to indices are almost identical even though their background hilosohies are strongly different. One imortant difference is that the LQTAI is defined so that it is alicable also for a non-stationary macro economy ith varying ork time ratio. A result of the more general formulation is that not only a rincile is obtained for socio-econnomical allocation of time (i.e. money) to fatality mitigation, but also for allocation of reasonable and ethically defendable time for injury mitigation. The concluding examle ith a ferry in fire illustrates the fact that there is a large difference beteen the damage comensation required from the oner of a risky oeration and the societal value of life and limb. This difference motivates the formulation of a ublic accetance rule (reviously ublished) hich is considerably more restrictive than the LQI (or LQTAI) criterion hich states that the gain er time unit must be at least as large as the societal loss value er time unit as obtained by the rincile of LQI invariance. Acknoledgment The first author thanks Otto Mønsteds Fond for conference travel suort. References Ditlevsen, O. (2003). Decision modeling and accetance criteria. Structural Safety 25: Ditlevsen, O. (2004). Life quality index revisited. Structural Safety 26: Ditlevsen, O. (2007). Model of observed stochastic balance beteen ork and free time suorting the LQTAI definition. Structural Safety. In ress. Prerint on htt://.mek.dtu.dk/staff/od/aers.htm until ublication. Ditlevsen, O. & P. Friis-Hansen (2005). Life quality allocation index an equilibrium economy consistent version of the current life quality index. Structural Safety 27: Ditlevsen, O. & P. Friis-Hansen (2007). Life quality index an emirical or a normative concet? International Journal of Risk Assessment and Management, in ress. Friis-Hansen, P. & O. Ditlevsen (2003). Nature reservation accetance model alied to tanker oil sill simulations. Structural Safety 25:1 34. Hansen, H. L. & M. Vinter (2003). Occuational accidents and shi design: Imlications for revention. In World Maritime Conference, San Francisco. Nathani, J. S., N. C. Lind, & M. D. Pandey (1997). Affordable Safety By Choice: The Life Quality Method. Waterloo, Ontario, Canada: Institute for risk Research, University of Waterloo. OECD, S. (2004). OECD Economic Outlook Database, ne.sourceoecd.org/database/oecdeconomicoutlook. Pandey, M. D., J. S. Nathani, & N. C. Lind (2006). The derivation and calibration of the life quality index (LQI) from economical rinciles. Structural Safety 28: Rackitz, R. (2002). Otimization and risk accetability based on the life quality index. Structural Safety 24: Rackitz, R., A. Lentz, & M. Faber (2005). Socio-economically sustainable civil engineering infrastructures by otimization. Structural Safety 27: Vanem, E. & R. Skjong (2006). Designing for safety in assenger shis utilizing advanced evacuation analysis - a risk based aroach. Safety Science 44:

Research on Evaluation Method of Organization s Performance Based on Comparative Advantage Characteristics

Research on Evaluation Method of Organization s Performance Based on Comparative Advantage Characteristics Vol.1, No.10, Ar 01,.67-7 Research on Evaluation Method of Organization s Performance Based on Comarative Advantage Characteristics WEN Xin 1, JIA Jianfeng and ZHAO Xi nan 3 Abstract It as under the guidance

More information

Feedback-error control

Feedback-error control Chater 4 Feedback-error control 4.1 Introduction This chater exlains the feedback-error (FBE) control scheme originally described by Kawato [, 87, 8]. FBE is a widely used neural network based controller

More information

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO) Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment

More information

Finding Shortest Hamiltonian Path is in P. Abstract

Finding Shortest Hamiltonian Path is in P. Abstract Finding Shortest Hamiltonian Path is in P Dhananay P. Mehendale Sir Parashurambhau College, Tilak Road, Pune, India bstract The roblem of finding shortest Hamiltonian ath in a eighted comlete grah belongs

More information

Statics and dynamics: some elementary concepts

Statics and dynamics: some elementary concepts 1 Statics and dynamics: some elementary concets Dynamics is the study of the movement through time of variables such as heartbeat, temerature, secies oulation, voltage, roduction, emloyment, rices and

More information

Notes on Instrumental Variables Methods

Notes on Instrumental Variables Methods Notes on Instrumental Variables Methods Michele Pellizzari IGIER-Bocconi, IZA and frdb 1 The Instrumental Variable Estimator Instrumental variable estimation is the classical solution to the roblem of

More information

Tempo effect in first marriage table: Japan and China. Kiyosi Hirosima (Shimane University, Japan)

Tempo effect in first marriage table: Japan and China. Kiyosi Hirosima (Shimane University, Japan) Temo effect in first marriage table: Jaan and China Kiyosi Hirosima (Shimane University, Jaan) Introduction Recently many East Asian countries and regions are undergoing the raid increase in age at marriage.

More information

arxiv: v1 [physics.data-an] 26 Oct 2012

arxiv: v1 [physics.data-an] 26 Oct 2012 Constraints on Yield Parameters in Extended Maximum Likelihood Fits Till Moritz Karbach a, Maximilian Schlu b a TU Dortmund, Germany, moritz.karbach@cern.ch b TU Dortmund, Germany, maximilian.schlu@cern.ch

More information

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression

A Comparison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Journal of Modern Alied Statistical Methods Volume Issue Article 7 --03 A Comarison between Biased and Unbiased Estimators in Ordinary Least Squares Regression Ghadban Khalaf King Khalid University, Saudi

More information

DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS METHODS

DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS METHODS Eleventh International IBPSA Conference Glasgo, Scotland July 27-30, 2009 DEVELOPMENT AND VALIDATION OF A VERSATILE METHOD FOR THE CALCULATION OF HEAT TRANSFER IN WATER-BASED RADIANT SYSTEMS Massimiliano

More information

Developing A Deterioration Probabilistic Model for Rail Wear

Developing A Deterioration Probabilistic Model for Rail Wear International Journal of Traffic and Transortation Engineering 2012, 1(2): 13-18 DOI: 10.5923/j.ijtte.20120102.02 Develoing A Deterioration Probabilistic Model for Rail Wear Jabbar-Ali Zakeri *, Shahrbanoo

More information

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation Paer C Exact Volume Balance Versus Exact Mass Balance in Comositional Reservoir Simulation Submitted to Comutational Geosciences, December 2005. Exact Volume Balance Versus Exact Mass Balance in Comositional

More information

The Binomial Approach for Probability of Detection

The Binomial Approach for Probability of Detection Vol. No. (Mar 5) - The e-journal of Nondestructive Testing - ISSN 45-494 www.ndt.net/?id=7498 The Binomial Aroach for of Detection Carlos Correia Gruo Endalloy C.A. - Caracas - Venezuela www.endalloy.net

More information

An Improved Calibration Method for a Chopped Pyrgeometer

An Improved Calibration Method for a Chopped Pyrgeometer 96 JOURNAL OF ATMOSPHERIC AND OCEANIC TECHNOLOGY VOLUME 17 An Imroved Calibration Method for a Choed Pyrgeometer FRIEDRICH FERGG OtoLab, Ingenieurbüro, Munich, Germany PETER WENDLING Deutsches Forschungszentrum

More information

Distributed Rule-Based Inference in the Presence of Redundant Information

Distributed Rule-Based Inference in the Presence of Redundant Information istribution Statement : roved for ublic release; distribution is unlimited. istributed Rule-ased Inference in the Presence of Redundant Information June 8, 004 William J. Farrell III Lockheed Martin dvanced

More information

An Analysis of Reliable Classifiers through ROC Isometrics

An Analysis of Reliable Classifiers through ROC Isometrics An Analysis of Reliable Classifiers through ROC Isometrics Stijn Vanderlooy s.vanderlooy@cs.unimaas.nl Ida G. Srinkhuizen-Kuyer kuyer@cs.unimaas.nl Evgueni N. Smirnov smirnov@cs.unimaas.nl MICC-IKAT, Universiteit

More information

CHAPTER 5 STATISTICAL INFERENCE. 1.0 Hypothesis Testing. 2.0 Decision Errors. 3.0 How a Hypothesis is Tested. 4.0 Test for Goodness of Fit

CHAPTER 5 STATISTICAL INFERENCE. 1.0 Hypothesis Testing. 2.0 Decision Errors. 3.0 How a Hypothesis is Tested. 4.0 Test for Goodness of Fit Chater 5 Statistical Inference 69 CHAPTER 5 STATISTICAL INFERENCE.0 Hyothesis Testing.0 Decision Errors 3.0 How a Hyothesis is Tested 4.0 Test for Goodness of Fit 5.0 Inferences about Two Means It ain't

More information

Economics 101. Lecture 7 - Monopoly and Oligopoly

Economics 101. Lecture 7 - Monopoly and Oligopoly Economics 0 Lecture 7 - Monooly and Oligooly Production Equilibrium After having exlored Walrasian equilibria with roduction in the Robinson Crusoe economy, we will now ste in to a more general setting.

More information

State Estimation with ARMarkov Models

State Estimation with ARMarkov Models Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,

More information

1 Gambler s Ruin Problem

1 Gambler s Ruin Problem Coyright c 2017 by Karl Sigman 1 Gambler s Ruin Problem Let N 2 be an integer and let 1 i N 1. Consider a gambler who starts with an initial fortune of $i and then on each successive gamble either wins

More information

Morten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014

Morten Frydenberg Section for Biostatistics Version :Friday, 05 September 2014 Morten Frydenberg Section for Biostatistics Version :Friday, 05 Setember 204 All models are aroximations! The best model does not exist! Comlicated models needs a lot of data. lower your ambitions or get

More information

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition A Qualitative Event-based Aroach to Multile Fault Diagnosis in Continuous Systems using Structural Model Decomosition Matthew J. Daigle a,,, Anibal Bregon b,, Xenofon Koutsoukos c, Gautam Biswas c, Belarmino

More information

Factor Analysis of Convective Heat Transfer for a Horizontal Tube in the Turbulent Flow Region Using Artificial Neural Network

Factor Analysis of Convective Heat Transfer for a Horizontal Tube in the Turbulent Flow Region Using Artificial Neural Network COMPUTATIONAL METHODS IN ENGINEERING AND SCIENCE EPMESC X, Aug. -3, 6, Sanya, Hainan,China 6 Tsinghua University ess & Sringer-Verlag Factor Analysis of Convective Heat Transfer for a Horizontal Tube in

More information

8.7 Associated and Non-associated Flow Rules

8.7 Associated and Non-associated Flow Rules 8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more

More information

Road Traffic Accidents in Saudi Arabia: An ARDL Approach and Multivariate Granger Causality

Road Traffic Accidents in Saudi Arabia: An ARDL Approach and Multivariate Granger Causality MPRA Munich Personal RePEc Archive Road Traffic Accidents in Saudi Arabia: An ARDL Aroach and Multivariate Granger Causality Mohammed Moosa Ageli King Saud University, RCC, Riyadh, Saudi Arabia 24. Aril

More information

Information collection on a graph

Information collection on a graph Information collection on a grah Ilya O. Ryzhov Warren Powell February 10, 2010 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements

More information

Estimation of Separable Representations in Psychophysical Experiments

Estimation of Separable Representations in Psychophysical Experiments Estimation of Searable Reresentations in Psychohysical Exeriments Michele Bernasconi (mbernasconi@eco.uninsubria.it) Christine Choirat (cchoirat@eco.uninsubria.it) Raffaello Seri (rseri@eco.uninsubria.it)

More information

Robust Beamforming via Matrix Completion

Robust Beamforming via Matrix Completion Robust Beamforming via Matrix Comletion Shunqiao Sun and Athina P. Petroulu Deartment of Electrical and Comuter Engineering Rutgers, the State University of Ne Jersey Piscataay, Ne Jersey 8854-858 Email:

More information

General Linear Model Introduction, Classes of Linear models and Estimation

General Linear Model Introduction, Classes of Linear models and Estimation Stat 740 General Linear Model Introduction, Classes of Linear models and Estimation An aim of scientific enquiry: To describe or to discover relationshis among events (variables) in the controlled (laboratory)

More information

Estimation of component redundancy in optimal age maintenance

Estimation of component redundancy in optimal age maintenance EURO MAINTENANCE 2012, Belgrade 14-16 May 2012 Proceedings of the 21 st International Congress on Maintenance and Asset Management Estimation of comonent redundancy in otimal age maintenance Jorge ioa

More information

Aggregate Prediction With. the Aggregation Bias

Aggregate Prediction With. the Aggregation Bias 100 Aggregate Prediction With Disaggregate Models: Behavior of the Aggregation Bias Uzi Landau, Transortation Research nstitute, Technion-srael nstitute of Technology, Haifa Disaggregate travel demand

More information

Information collection on a graph

Information collection on a graph Information collection on a grah Ilya O. Ryzhov Warren Powell October 25, 2009 Abstract We derive a knowledge gradient olicy for an otimal learning roblem on a grah, in which we use sequential measurements

More information

Linear diophantine equations for discrete tomography

Linear diophantine equations for discrete tomography Journal of X-Ray Science and Technology 10 001 59 66 59 IOS Press Linear diohantine euations for discrete tomograhy Yangbo Ye a,gewang b and Jiehua Zhu a a Deartment of Mathematics, The University of Iowa,

More information

Proof: We follow thearoach develoed in [4]. We adot a useful but non-intuitive notion of time; a bin with z balls at time t receives its next ball at

Proof: We follow thearoach develoed in [4]. We adot a useful but non-intuitive notion of time; a bin with z balls at time t receives its next ball at A Scaling Result for Exlosive Processes M. Mitzenmacher Λ J. Sencer We consider the following balls and bins model, as described in [, 4]. Balls are sequentially thrown into bins so that the robability

More information

Hotelling s Two- Sample T 2

Hotelling s Two- Sample T 2 Chater 600 Hotelling s Two- Samle T Introduction This module calculates ower for the Hotelling s two-grou, T-squared (T) test statistic. Hotelling s T is an extension of the univariate two-samle t-test

More information

AI*IA 2003 Fusion of Multiple Pattern Classifiers PART III

AI*IA 2003 Fusion of Multiple Pattern Classifiers PART III AI*IA 23 Fusion of Multile Pattern Classifiers PART III AI*IA 23 Tutorial on Fusion of Multile Pattern Classifiers by F. Roli 49 Methods for fusing multile classifiers Methods for fusing multile classifiers

More information

Wolfgang POESSNECKER and Ulrich GROSS*

Wolfgang POESSNECKER and Ulrich GROSS* Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND

More information

4. Score normalization technical details We now discuss the technical details of the score normalization method.

4. Score normalization technical details We now discuss the technical details of the score normalization method. SMT SCORING SYSTEM This document describes the scoring system for the Stanford Math Tournament We begin by giving an overview of the changes to scoring and a non-technical descrition of the scoring rules

More information

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data

Lower Confidence Bound for Process-Yield Index S pk with Autocorrelated Process Data Quality Technology & Quantitative Management Vol. 1, No.,. 51-65, 15 QTQM IAQM 15 Lower onfidence Bound for Process-Yield Index with Autocorrelated Process Data Fu-Kwun Wang * and Yeneneh Tamirat Deartment

More information

Econ 101A Midterm 2 Th 8 April 2009.

Econ 101A Midterm 2 Th 8 April 2009. Econ A Midterm Th 8 Aril 9. You have aroximately hour and minutes to answer the questions in the midterm. I will collect the exams at. shar. Show your work, and good luck! Problem. Production (38 oints).

More information

Estimation of the large covariance matrix with two-step monotone missing data

Estimation of the large covariance matrix with two-step monotone missing data Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo

More information

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management

More information

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek

Use of Transformations and the Repeated Statement in PROC GLM in SAS Ed Stanek Use of Transformations and the Reeated Statement in PROC GLM in SAS Ed Stanek Introduction We describe how the Reeated Statement in PROC GLM in SAS transforms the data to rovide tests of hyotheses of interest.

More information

The Poisson Regression Model

The Poisson Regression Model The Poisson Regression Model The Poisson regression model aims at modeling a counting variable Y, counting the number of times that a certain event occurs during a given time eriod. We observe a samle

More information

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential

dn i where we have used the Gibbs equation for the Gibbs energy and the definition of chemical potential Chem 467 Sulement to Lectures 33 Phase Equilibrium Chemical Potential Revisited We introduced the chemical otential as the conjugate variable to amount. Briefly reviewing, the total Gibbs energy of a system

More information

DEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES

DEPARTMENT OF ECONOMICS ISSN DISCUSSION PAPER 20/07 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES DEPARTMENT OF ECONOMICS ISSN 1441-549 DISCUSSION PAPER /7 TWO NEW EXPONENTIAL FAMILIES OF LORENZ CURVES ZuXiang Wang * & Russell Smyth ABSTRACT We resent two new Lorenz curve families by using the basic

More information

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation

Uniformly best wavenumber approximations by spatial central difference operators: An initial investigation Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations

More information

arxiv: v1 [nlin.cg] 7 Dec 2015

arxiv: v1 [nlin.cg] 7 Dec 2015 When is a bottleneck a bottleneck? Andreas Schadschneider, Johannes Schmidt and Vladislav Pokov arxiv:52.02626v [nlin.cg] 7 Dec 205 Abstract Bottlenecks, i.e. local reductions of caacity, are one of the

More information

A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION

A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST BASED ON THE WEIBULL DISTRIBUTION O P E R A T I O N S R E S E A R C H A N D D E C I S I O N S No. 27 DOI:.5277/ord73 Nasrullah KHAN Muhammad ASLAM 2 Kyung-Jun KIM 3 Chi-Hyuck JUN 4 A MIXED CONTROL CHART ADAPTED TO THE TRUNCATED LIFE TEST

More information

UNDERWATER SOUND RECEIVED FROM SOME DEFENCE ACTIVITIES IN SHALLOW OCEAN REGIONS. Signal and Noise Transmission and Data Processing

UNDERWATER SOUND RECEIVED FROM SOME DEFENCE ACTIVITIES IN SHALLOW OCEAN REGIONS. Signal and Noise Transmission and Data Processing UNDERWATER SOUND RECEIVED FROM SOME DEFENCE ACTIVITIES IN SHALLOW OCEAN REGIONS Adrian D. Jones(1) and Paul A. Clarke() (1) Program Office, DSTO Edinburgh, SA 5111, Australia () Maritime Oerations Division,

More information

Prospect Theory Explains Newsvendor Behavior: The Role of Reference Points

Prospect Theory Explains Newsvendor Behavior: The Role of Reference Points Submitted to Management Science manuscrit (Please, rovide the mansucrit number! Authors are encouraged to submit new aers to INFORMS journals by means of a style file temlate, which includes the journal

More information

Approximate Dynamic Programming for Dynamic Capacity Allocation with Multiple Priority Levels

Approximate Dynamic Programming for Dynamic Capacity Allocation with Multiple Priority Levels Aroximate Dynamic Programming for Dynamic Caacity Allocation with Multile Priority Levels Alexander Erdelyi School of Oerations Research and Information Engineering, Cornell University, Ithaca, NY 14853,

More information

Optimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment

Optimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment Neutrosohic Sets and Systems Vol 14 016 93 University of New Mexico Otimal Design of Truss Structures Using a Neutrosohic Number Otimization Model under an Indeterminate Environment Wenzhong Jiang & Jun

More information

NUMERICAL STUDY OF REVERSAL FLOW IN TUNNEL FIRES

NUMERICAL STUDY OF REVERSAL FLOW IN TUNNEL FIRES , Volume, Number 4,.94-98, NUMERIAL STUD OF REVERSAL FLOW IN TUNNEL FIRES Jiangguo hen, Haixin hen and Song Fu Deartment of Engineering Mechanics, Tsinghua University, Beijing 84, hina ABSTRAT In the fire

More information

Theory of Externalities Partial Equilibrium Analysis

Theory of Externalities Partial Equilibrium Analysis Theory of Externalities Partial Equilibrium Analysis Definition: An externality is resent whenever the well being of a consumer or the roduction ossibilities of a firm are directly affected by the actions

More information

Elementary Analysis in Q p

Elementary Analysis in Q p Elementary Analysis in Q Hannah Hutter, May Szedlák, Phili Wirth November 17, 2011 This reort follows very closely the book of Svetlana Katok 1. 1 Sequences and Series In this section we will see some

More information

On split sample and randomized confidence intervals for binomial proportions

On split sample and randomized confidence intervals for binomial proportions On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have

More information

The directivity of the forced radiation of sound from panels and openings including the shadow zone

The directivity of the forced radiation of sound from panels and openings including the shadow zone The directivity of the forced radiation of sound from anels and oenings including the shadow zone J. Davy RMIT University, Alied Physics, GPO Box 476V, 3001 Melbourne, Victoria, Australia john.davy@rmit.edu.au

More information

arxiv:cond-mat/ v2 25 Sep 2002

arxiv:cond-mat/ v2 25 Sep 2002 Energy fluctuations at the multicritical oint in two-dimensional sin glasses arxiv:cond-mat/0207694 v2 25 Se 2002 1. Introduction Hidetoshi Nishimori, Cyril Falvo and Yukiyasu Ozeki Deartment of Physics,

More information

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK Comuter Modelling and ew Technologies, 5, Vol.9, o., 3-39 Transort and Telecommunication Institute, Lomonosov, LV-9, Riga, Latvia MATHEMATICAL MODELLIG OF THE WIRELESS COMMUICATIO ETWORK M. KOPEETSK Deartment

More information

A continuous review inventory model with the controllable production rate of the manufacturer

A continuous review inventory model with the controllable production rate of the manufacturer Intl. Trans. in O. Res. 12 (2005) 247 258 INTERNATIONAL TRANSACTIONS IN OERATIONAL RESEARCH A continuous review inventory model with the controllable roduction rate of the manufacturer I. K. Moon and B.

More information

Notes on Optical Pumping Procedure & Theory

Notes on Optical Pumping Procedure & Theory Notes on Otical Puming Procedure & Theory Pre-lab 1. Why is the exeriment called otical uming? What is umed? 2. What is the exerimental signature of having cancelled all magnetic fields in the samle cell?

More information

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning

Uncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment

More information

y p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2

y p 2 p 1 y p Flexible System p 2 y p2 p1 y p u, y 1 -u, y 2 Component Breakdown z 1 w 1 1 y p 1 P 1 P 2 w 2 z 2 ROBUSTNESS OF FLEXIBLE SYSTEMS WITH COMPONENT-LEVEL UNCERTAINTIES Peiman G. Maghami Λ NASA Langley Research Center, Hamton, VA 368 Robustness of flexible systems in the resence of model uncertainties at

More information

Chapter 1 Fundamentals

Chapter 1 Fundamentals Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors

More information

Hidden Predictors: A Factor Analysis Primer

Hidden Predictors: A Factor Analysis Primer Hidden Predictors: A Factor Analysis Primer Ryan C Sanchez Western Washington University Factor Analysis is a owerful statistical method in the modern research sychologist s toolbag When used roerly, factor

More information

START Selected Topics in Assurance

START Selected Topics in Assurance START Selected Toics in Assurance Related Technologies Table of Contents Introduction Statistical Models for Simle Systems (U/Down) and Interretation Markov Models for Simle Systems (U/Down) and Interretation

More information

Johan Lyhagen Department of Information Science, Uppsala University. Abstract

Johan Lyhagen Department of Information Science, Uppsala University. Abstract Why not use standard anel unit root test for testing PPP Johan Lyhagen Deartment of Information Science, Usala University Abstract In this aer we show the consequences of alying a anel unit root test that

More information

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK

Towards understanding the Lorenz curve using the Uniform distribution. Chris J. Stephens. Newcastle City Council, Newcastle upon Tyne, UK Towards understanding the Lorenz curve using the Uniform distribution Chris J. Stehens Newcastle City Council, Newcastle uon Tyne, UK (For the Gini-Lorenz Conference, University of Siena, Italy, May 2005)

More information

Elliptic Curves and Cryptography

Elliptic Curves and Cryptography Ellitic Curves and Crytograhy Background in Ellitic Curves We'll now turn to the fascinating theory of ellitic curves. For simlicity, we'll restrict our discussion to ellitic curves over Z, where is a

More information

PROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2)

PROFIT MAXIMIZATION. π = p y Σ n i=1 w i x i (2) PROFIT MAXIMIZATION DEFINITION OF A NEOCLASSICAL FIRM A neoclassical firm is an organization that controls the transformation of inuts (resources it owns or urchases into oututs or roducts (valued roducts

More information

International Journal of Statistics and Mathematics Vol. 1(1), pp , April, ISSN: x

International Journal of Statistics and Mathematics Vol. 1(1), pp , April, ISSN: x International Journal of Statistics and Mathematics Vol. 1(1),. 2-8, Aril, 214. www.remierublishers.org ISSN: 2375-499 x IJSM Research Article The distortion measurement strategy in manufacturing trading

More information

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites

More information

John Weatherwax. Analysis of Parallel Depth First Search Algorithms

John Weatherwax. Analysis of Parallel Depth First Search Algorithms Sulementary Discussions and Solutions to Selected Problems in: Introduction to Parallel Comuting by Viin Kumar, Ananth Grama, Anshul Guta, & George Karyis John Weatherwax Chater 8 Analysis of Parallel

More information

rate~ If no additional source of holes were present, the excess

rate~ If no additional source of holes were present, the excess DIFFUSION OF CARRIERS Diffusion currents are resent in semiconductor devices which generate a satially non-uniform distribution of carriers. The most imortant examles are the -n junction and the biolar

More information

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model Shadow Comuting: An Energy-Aware Fault Tolerant Comuting Model Bryan Mills, Taieb Znati, Rami Melhem Deartment of Comuter Science University of Pittsburgh (bmills, znati, melhem)@cs.itt.edu Index Terms

More information

INTRODUCTION. Please write to us at if you have any comments or ideas. We love to hear from you.

INTRODUCTION. Please write to us at if you have any comments or ideas. We love to hear from you. Casio FX-570ES One-Page Wonder INTRODUCTION Welcome to the world of Casio s Natural Dislay scientific calculators. Our exeriences of working with eole have us understand more about obstacles eole face

More information

The non-stochastic multi-armed bandit problem

The non-stochastic multi-armed bandit problem Submitted for journal ublication. The non-stochastic multi-armed bandit roblem Peter Auer Institute for Theoretical Comuter Science Graz University of Technology A-8010 Graz (Austria) auer@igi.tu-graz.ac.at

More information

Age of Information: Whittle Index for Scheduling Stochastic Arrivals

Age of Information: Whittle Index for Scheduling Stochastic Arrivals Age of Information: Whittle Index for Scheduling Stochastic Arrivals Yu-Pin Hsu Deartment of Communication Engineering National Taiei University yuinhsu@mail.ntu.edu.tw arxiv:80.03422v2 [math.oc] 7 Ar

More information

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.

More information

arxiv: v1 [math-ph] 29 Apr 2016

arxiv: v1 [math-ph] 29 Apr 2016 INFORMATION THEORY AND STATISTICAL MECHANICS REVISITED arxiv:64.8739v [math-h] 29 Ar 26 JIAN ZHOU Abstract. We derive Bose-Einstein statistics and Fermi-Dirac statistics by Princile of Maximum Entroy alied

More information

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density

A Note on Massless Quantum Free Scalar Fields. with Negative Energy Density Adv. Studies Theor. Phys., Vol. 7, 13, no. 1, 549 554 HIKARI Ltd, www.m-hikari.com A Note on Massless Quantum Free Scalar Fields with Negative Energy Density M. A. Grado-Caffaro and M. Grado-Caffaro Scientific

More information

Variable Selection and Model Building

Variable Selection and Model Building LINEAR REGRESSION ANALYSIS MODULE XIII Lecture - 38 Variable Selection and Model Building Dr. Shalabh Deartment of Mathematics and Statistics Indian Institute of Technology Kanur Evaluation of subset regression

More information

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE ANALYTICAL MODEL FOR THE BYPASS ALE IN A LOOP HEAT PIPE Michel Seetjens & Camilo Rindt Laboratory for Energy Technology Mechanical Engineering Deartment Eindhoven University of Technology The Netherlands

More information

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules

CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is

More information

Topic 7: Using identity types

Topic 7: Using identity types Toic 7: Using identity tyes June 10, 2014 Now we would like to learn how to use identity tyes and how to do some actual mathematics with them. By now we have essentially introduced all inference rules

More information

COMMUNICATION BETWEEN SHAREHOLDERS 1

COMMUNICATION BETWEEN SHAREHOLDERS 1 COMMUNICATION BTWN SHARHOLDRS 1 A B. O A : A D Lemma B.1. U to µ Z r 2 σ2 Z + σ2 X 2r ω 2 an additive constant that does not deend on a or θ, the agents ayoffs can be written as: 2r rθa ω2 + θ µ Y rcov

More information

Adaptive estimation with change detection for streaming data

Adaptive estimation with change detection for streaming data Adative estimation with change detection for streaming data A thesis resented for the degree of Doctor of Philosohy of the University of London and the Diloma of Imerial College by Dean Adam Bodenham Deartment

More information

An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators

An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical

More information

Session 5: Review of Classical Astrodynamics

Session 5: Review of Classical Astrodynamics Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve

More information

All-fiber Optical Parametric Oscillator

All-fiber Optical Parametric Oscillator All-fiber Otical Parametric Oscillator Chengao Wang Otical Science and Engineering, Deartment of Physics & Astronomy, University of New Mexico Albuquerque, NM 87131-0001, USA Abstract All-fiber otical

More information

Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process

Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process Using the Divergence Information Criterion for the Determination of the Order of an Autoregressive Process P. Mantalos a1, K. Mattheou b, A. Karagrigoriou b a.deartment of Statistics University of Lund

More information

PHYSICAL REVIEW LETTERS

PHYSICAL REVIEW LETTERS PHYSICAL REVIEW LETTERS VOLUME 81 20 JULY 1998 NUMBER 3 Searated-Path Ramsey Atom Interferometer P. D. Featonby, G. S. Summy, C. L. Webb, R. M. Godun, M. K. Oberthaler, A. C. Wilson, C. J. Foot, and K.

More information

Modeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation

Modeling and Estimation of Full-Chip Leakage Current Considering Within-Die Correlation 6.3 Modeling and Estimation of Full-Chi Leaage Current Considering Within-Die Correlation Khaled R. eloue, Navid Azizi, Farid N. Najm Deartment of ECE, University of Toronto,Toronto, Ontario, Canada {haled,nazizi,najm}@eecg.utoronto.ca

More information

Analysis of some entrance probabilities for killed birth-death processes

Analysis of some entrance probabilities for killed birth-death processes Analysis of some entrance robabilities for killed birth-death rocesses Master s Thesis O.J.G. van der Velde Suervisor: Dr. F.M. Sieksma July 5, 207 Mathematical Institute, Leiden University Contents Introduction

More information

Modelling 'Animal Spirits' and Network Effects in Macroeconomics and Financial Markets Thomas Lux

Modelling 'Animal Spirits' and Network Effects in Macroeconomics and Financial Markets Thomas Lux Modelling 'Animal Sirits' and Network Effects in Macroeconomics and Financial Markets Kiel Institute for the World Economy & Banco de Esaña Chair of Comutational Economics University of Castellón GSDP

More information

SAS for Bayesian Mediation Analysis

SAS for Bayesian Mediation Analysis Paer 1569-2014 SAS for Bayesian Mediation Analysis Miočević Milica, Arizona State University; David P. MacKinnon, Arizona State University ABSTRACT Recent statistical mediation analysis research focuses

More information

8 STOCHASTIC PROCESSES

8 STOCHASTIC PROCESSES 8 STOCHASTIC PROCESSES The word stochastic is derived from the Greek στoχαστικoς, meaning to aim at a target. Stochastic rocesses involve state which changes in a random way. A Markov rocess is a articular

More information

An Analysis of TCP over Random Access Satellite Links

An Analysis of TCP over Random Access Satellite Links An Analysis of over Random Access Satellite Links Chunmei Liu and Eytan Modiano Massachusetts Institute of Technology Cambridge, MA 0239 Email: mayliu, modiano@mit.edu Abstract This aer analyzes the erformance

More information

Finite Mixture EFA in Mplus

Finite Mixture EFA in Mplus Finite Mixture EFA in Mlus November 16, 2007 In this document we describe the Mixture EFA model estimated in Mlus. Four tyes of deendent variables are ossible in this model: normally distributed, ordered

More information