CSIS Discussion Paper No.77

需要密度に対する供給効果を最大化する AED の最適配置地点
Optimal allocation of AED’s by maximization of supply effect on demand density

片岡裕介

浅見泰司

浅利靖

郡山一明

Yusuke KATAOKA, Yasushi ASAMI, Yasushi ASARI and Kazuaki KOHRIYAMA東京大学大学院新領域創成科学研究科社会文化環境学専攻
Course of Socio-Cultural and Socio-Physical Environmental Studies,
Graduate School of Frontier Sciences, University of Tokyo
5-1-5, Kashiwanoha, Kashiwa-shi, Chiba 277-8568, Japane-mail: ukataoka@csis.u-tokyo.ac.jp

需要密度に対する供給効果を最大化する AED の最適配置地点
片岡裕介・浅見泰司・浅利靖・郡山一明Optimal allocation of AED’s by maximization of supply effect on demand density
Yusuke KATAOKA, Yasushi ASAMI, Yasushi ASARI and Kazuaki KOHRIYAMA
Abstract: A probability density function of demand for automated external
defibrillators (AED's), which is based on the distribution of demand points by kernel
density estimation method, is estimated from the data on occurrences ofcardiopulmonary arrests. The probability that someone can be saved by AED's isformulated using probability of survival to hospital discharge. “Supply effect” isexpressed as demand density multiplied by this probability of saving, and thelocations of AED's are optimized by maximizing supply effect in the entire region. The
results show that optimal locations tend to cover the center of the city and that the
supply effect is approximately proportional to the number of AED's.
Keywords: AED（automated external defibrillator），密度分布（density distribution），
最大被覆（maximal covering），空間パターン（spatial pattern），生存退院率（probability of survival to hospital discharge）

1. はじめにAED 設置に関する最近の研究をみると，Malcom

AED（自動体外式除細動器）は，突然的な心停止

et al.（2004）は，地域の人口密度と心停止が発生状態への救急措置に有効とされる医療機器である．

した場所との関連についての調査を行った．その結取り扱いが簡便で，既に一般市民による使用が認め果として，人口密度が高い場所であるほど，病院外られている状況のなか，全国の自治体には設置に向で心停止が発生する割合が減少する一方で，自宅外けた早急な対処が求められている．現状の AED ので発生する割合については増加する傾向があったこ

設置場所としてその多くは駅，空港やイベント会場とが報じられている．また，Crocco et al.（2004）などの大規模な集客施設，もしくは医療機関が主とでは，集客施設で必要とされる，AED の設置数になる傾向がみられる．人命に直接的に関わる事態で関する数理モデルの提案がなされた．AED の効果

ある AED の設置に際しては，地域の社会的そしてが注目されつつある現在の社会的状況を鑑みても，経済的な事情も加味されなければならない．そこでその設置地点に関する検討が今後求められていくで殊に強調されるべきは，より多くの装置を用意するあろうし，この要求に応える最も直接的で有効な方

努力もさることながら，装置の配置地点に関する慎法は，理論的に得られた設置地点が示されることで重な検討であろう．

あろう．ならびに GIS の分野としても，危機情報の片岡：〒277-8568 千葉県柏市柏の葉 5-1-5
東京大学大学院新領域創成科学研究科社会文化環境学専攻
Graduate School of Frontier Sciences, University of Tokyo
E-mail ukataoka@csis.u-tokyo.ac.jp

管理および解析等の面から，より一層の貢献が望まれものと考える．
そこで本稿では，AED の最適配置地点に関する

考察をおこなうにあたり，主に以下の点を目的とし点は，自宅，医療・福祉施設等といった一般に日常ている．まず，心停止発生地点（需要点）の空間パ

生活の多くを過ごす場所での発生，および外出先でターンにもとづいた地域全体の装置の需要量を表すの発生からなると考えられる．つまり，地域の人口密度分布に対して，装置との位置関係で決定される分布のみならず他の要因をも加味した，AED の潜救命確率をもって各配置地点（供給点）周辺で重み

在的な需要量を表す密度分布がカーネル法で推定さ付けした「供給効果」を最大化する配置問題を提案れる．

する．さらに，実際の心停止発生地点をデータとし

なお，本稿ではカーネル関数のパラメータとなるて用いた場合の最適配置地点を示すとともに，地域

バンド幅選択に際して，Sheather and Jones（1991）の救命確率および「供給効果」の地理的分布についの方程式解によるプラグイン法を採用する 2)．

ても表すことで，本稿で提案される方法の有効性を実証する．
2. 需要量の空間パターンと供給効果の定式化
本節では，まず地域における需要量の空間パターンを把握する．さらに，装置の設置に伴って装置周辺に生じると考えられる救命確率を定義し，需要量に救命確率を加味した「供給効果」を定式化する．
2.1. 需要点にもとづく密度分布の推定
図 1 点分布と密度分布

通常，需要量を把握する際には，対象となるものの特性を考慮し，それに応じて適切な方法が用いられる．ここでは，実際に装置が必要とされた地点，

カーネル法によると，平面上の点(x,y)における確つまり過去に心停止が発生した地点を AED の「需率密度関数 lf ( x, y )は，点分布の各点 i の座標(xi,yi)を要点（各点で均質）」とし，需要点の分布パターンに用いて以下のように表せる（Bowman and Azzelini,

より地域の需要量が決定されるものと考えた．とは

1997）
．

いえ，心停止が発生した地点のみに需要があるとするのは，同一地点で再発する可能性を考えてみても妥当性を欠く．よって，過去に数多く発生した地域

 x − xi   y − yi lf ( x, y ) = 1
K K 
∑in hx   hy (1)

（ K ( ⋅)：カーネル関数，n：点の総数，
hx, hy：各座標値に関するバンド幅）を潜在的に需要が高い地域であると見なすこととする．従って各需要点周辺の需要量は，需要点で最大であり，周囲に拡がるにつれて減衰する確率密度分

また式(1)で，カーネル関数として正規分布を用い布で表現することができよう．ただし，ごく短期間て，各需要点の座標が(x1,i,y1,i)とすると，点分布にものうちでは，心停止の発生が偶然的であることを免とづく地域全体の確率密度関数は，以下の式(2)のよれない．そのため需要点数については，需要を代表うにできる．

1)

するのに十分な数である必要がある ．
需要量の地域全体の確率密度分布についてはカーネル法を用いて推定する（Silverman, 1986）．カ

lf ( x, y ) = ( x − x )2 ( y − y )2 
n1
1
1,i
1,i
∑ exp  − 2h 2 − 2h 2 
2π hx hy n1 i =1
x
y



(2)

（n1：需要点数）

ーネル法とは，点分布を形成する各点の中心に，カーネルと呼ばれる密度分布の山を置き，それら全て

式(2)は，需要点から推定された確率密度関数であのカーネルの合計をもって全体の確率密度関数を求ると同時に，需要量の地理的分布を表す空間パター

める密度推定法である（図 1 参照）．心停止発生地ンとして捉えられる．

2.2. 供給効果の定式化

刻を争う事態では行動が直感に大きく左右され，状次に，設置される装置周辺の需要量に対して，装況として通常の施設選択の場合とは異なるともいえ置による効果の程度を表す「供給効果」について考る．つまり，救助者は必ずしも現場から最も近くにえる．何よりも装置が有効に使用されるためには，ある装置を常に把握しているわけではなく，ある地

心停止発生時点から限られた時間内に処置が施され点において装置が心停止にもたらす効果は，心停止る必要がある．その際，各施設立地点から一定距離

発生地点から最も近くにある装置との関係性のみに内にある領域に含まれる需要量を最大化する配置問

影響するとは限らない 4)．すると装置の「認知度」

題 と し て 知 ら れ る ， 最 大 被 覆 問 題 （ Maximalの拡がりが，装置からの距離に伴う救命確率の減衰

Covering Location Problem）があげられる．最大に代替されているとしたときに，近くに設置されて被覆問題は，Church and ReVelle（1974）により提いる装置が多いほど，設置に伴う効果は上昇するも案された後，救急施設などを対象とした離散平面にのとも考えられる．以上より，ある地点におけるおける配置問題として発展するなかで，様々な解法

AED 設置によって与えられる影響の程度は，各装がこれまでに提案されてきた（Brotcorne et al.,

置の救命確率の和として説明できるものと仮定する．

2003）．ただ先述のとおり，地域の需要量は点分布

よって式(3)を用いると，地点(x,y)での救命確率のであるよりも，各地点の潜在的な需要を表す確率密

和は，各装置の設置地点の座標(x2,j,y2,j)を用いて以下

度分布で表現される方が適切であろう．
さらに AEDのようにできる．

の場合では，心停止に始まり処置に至るまでの時間に，救命確率が大きく依存するとされるため

n2

( (

g ( x, y ) = ∑ exp − a ( x − x2, j ) + ( y − y2, j )j =1（American Heart Association, 2000），被覆される

2

2

))

(4)

（n2：装置の設置数）

領域内においても場所によっては効果の程度は著しく異なると考えられる．同様に Alsalloum et al.

以上をふまえ，装置の効果が期待される装置周辺（2006）は，救急車の配置モデルにおいて，通常用において，救命確率を需要量に加味することで得らいられる一様に被覆される領域を，定めた制限時間

れる，
需要と装置の効果に依存する量を本稿では「供内に被覆される確率に置き換えた方法を提案してい給効果」と定義する．なお，理解を容易にするためるが，AED の場合においても類似した状況にあるに一次元で考えると，需要量と「供給効果」の関係と考えてよい．

は以下の図 2 のようになる．

そこで，装置から距離 d の地点における救命確率
は，経過時間に依存する救命確率を，距離的な問題に置き換えた以下の関数型で近似できるとする 3)．
g ( d ) = exp ( − ad 2 )

需要密度
供給効果

(3)
（a：パラメータ）

本稿では，心停止が発生した現場に居合わせた「救助者」が装置を取りに行くことを前提としてい装置の配置地点

る．救助者の置かれた状況を考えると，設置場所がわかっている装置については，装置に向かって取り

図 2 需要密度に対する供給効果に行くことが可能である．ところが，当然ながら救助者は装置の全設置点を把握しているわけではない．

これより，ある地点における需要に対して救命可

事前知識も不確実であることに加え，このように一能な量となる「供給効果」は，その地点の需要量に

救命確率の和を乗じたもので表されることから，「供置を複数台設置するという結果を得ることを意味す

給効果」の地域全体の総和（図 2 の網掛け部分）は，る．ところが，ある装置のごく近辺で，複数の心停

以下の式(5)のように定式化される．

止患者がほぼ同時に発生する状況は想定しづらく，
いわば「容量の制限がない」と条件付けられる AED

L1 ( x, y ) = ∫∫ lf ( x, y ) ⋅ g ( x, y ) dxdyの問題に限っては，同一地点に解が重複して得られる状態は好ましくない．そのため，ここでは装置間 ( x − x )2 ( y − y )2 
n1
1

1,i
1,i

= ∫∫ 
exp  −−∑
2hx 2
2hy 2  2π hx hy n1 i =1



( (

n2

× ∑ exp −a ( x − x2, j ) + ( y − y2, j )
j =1

=

2

2

))

の距離についての制約を設けることにより，救命確率の和が1 を超える地点が存在しないものとする 5)．

 dxdy

つまり，装置の効果が期待できる各装置の設置地点
からの限界の距離を dmin とおき，ある装置から dmin

1の距離内には他の装置が設置されることがないもの

4n1 1 + 2ahx 2 1 + 2ahy 2とする．よって，各装置間の距離を ds,t で表し，制約

 a ( x − x )2 a ( y − y )2 
n1 n2
1,i
2, j
1,i
2, j

× ∑∑ exp  −−
1 + 2ahx 2
1 + 2ahy 2 
i =1 j =1

条件として式(9)とおくことで，最適化されるときの各配置地点の座標値は，式(7)が最大化されるときの x − x1,i + 2ahx 2 ( x − x2, j ) × erf 
hx 2 + 4ahx 2

解集合として与えられる．つまり，本稿で定式化さ y − y + 2ah 2 ( x − y ) 
y
1,i
2, j

× erf 
2


+
h
2
4
ah
y
y



れた配置問題は，以下の多変数最大化問題に帰着す(5)

式(5)において，各項における積分定数についてはる．
Maximize

(

)

(8)

∀ s, t(9)

L2 x2,1 , x2,2 ," , x2,n2 , y2,1 , y2,2 ," , y2,n2

省略した．また，erf(⋅)は誤差関数をあらわし，以下subject toのように定義される．
2

erf ( z ) =zd s ,t =

e dt
π ∫
−t 2

(x − x ) +( y − y ) ≥ d
2

2, s

2,t

2

2, s

2,t

min

(6)

0

4. 適用事例
前節で提案された AED の配置問題の適用事例を扱う．需要点として弘前市の心停止発生地点のデー

3. 供給効果を最大化する配置問題
本節では，2 節で定式化された「供給効果」を最

4.1. 心停止発生地点とその密度分布大化する，AED の最適配置問題を提案する．
まず，式(5)で積分区間を無限領域としてよいこと(=

で，原因疾患が心疾患と推定された 411 例を分析に

)

用いた．

1

始めに，当データ内の属性として含まれている住n1 1 + 2ahx 2 1 + 2ahy 2
 a(x − x ) a( y − y ) 
1,i
2, j
1,i
2, j

× ∑∑ exp  −
−
2
2


1
2
1
2
ah
ah
+
+
i =1 j =1
x
y


n1

n2

需要点となる心停止発生地点としては，弘前地区
消防事務組合が過去 3 年間に扱った内因性心肺停止から，以下の式(7)を得る．
L2 x2,1 , x2,2 ," , x2,n2 , y2,1 , y2,2 ," , y2,n2

タを用いた場合の，AED の最適配置地点を示す．

2

所表記から位置情報を取得した 6)．各需要点の位置

2(7)と，式(2)より得られる各地点の確率密度を表す等値
線（間隔は 1.0×10-9）を図 3 に示す．なお，カーネ式(7)で定式化された「供給効果」を，このようにル関数のバンド幅選択は，2.1.で述べたように，設置場所の候補地が全くない状況で，単純に最大化

Sheather and Jones（1991）によるプラグイン法をしようとするならば，需要量の多い領域を幾度も被採用した．

覆してしまう可能性がある．それは，同一地点に装図 3 の需要点にもとづく密度分布について見てみると，弘前駅や弘前城の南側に隣接する市役所から．なお，実際のパララメータを求めた（a=4.32×10-5）

少し離れた地域，および土手町と呼ばれるメインス

メータの決定においては，本稿で設定された仮定にトリート付近を中心として，心停止発生密度が非常

もとづき，分速 80m で装置に向かって直進するとしに高いとされる地域が拡がっていることが確認されたうえで，往復にかかる時間を考え 5 分経過した

た．また，市街地部周縁である地域帯で密度分布が400m の時点で救命確率が 1.0×10-3 となるとした 8)．

大きく変化しており，概して人口分布の変化に対応
1

救命確率

しているものと理解される．

0.8
0.6
0.4
0.2
0

(m)
0

100

200

300

400

装置からの距離

図 4 各装置からの救命確率
4.3. AED の最適配置地点
次に式(8)を用いて，AED の設置数を 2 点から 30点までの間で所与とし最適地を求める．前項での考察をふまえて，式(9)において dmin=400m とし，式(8)を多変数の制約条件付き最大化問題として数値計算をおこなった結果 9)，分析の対象とした設置数の全ての場合で各最適地が対象地域内に得られた．そのなかで，設置数が「20 点」，「30 点」のときの各最
0

1

2

(km)

適地，および任意の地点における「供給効果」の密度分布を示す等値線（間隔は 5.0×10-8）を図 5 に，（図中の表記については以下の通り．●：需要点，および需要量密度を示す等値線．背景地図として，
河川，鉄道，弘前城を表記．
）

図 3 心停止発生地点とその密度分布

また，装置による救命確率の密度分布を示す等値線（間隔は 2.0×10-1）を図 6 に示す．参考のために，
図 3 の需要点の密度分布についても各図に併せて示した．
得られた最適配置地点の分布の特徴として，4.1.

4.2. 処置までの時間にもとづく救命確率式

で装置の需要密度が高いとされる地域において，設

本稿では，式(3)で定式化された装置からの距離に置数が増加するに伴い，最適地が密度分布の「山」

もとづく救命確率式を求めるために，処置までの時の頂点付近からその周辺に向かって，満遍なくカバ間に依存する「生存退院率」を用いることとする．

ーするように得られるという傾向が表れている．

「生存退院率」とは処置による蘇生の成功率を意味

また，図 6 の救命確率の密度分布は，需要分布にし，倒れてから除細動の開始までに要する時間が 1

依らない各地点での装置の効果を示しており，各等分経過するごとに約10%低下するとされることから

値線が救命確率で表された装置までの一定の到達距

7)，AED のような問題においては国内外を問わず重離の分布を意味している．従って，一刻を争う状況

要な知見とされている．これにもとづき，10 分経過

下では，視覚化された判断基準として特に有効な手

した時点で救命確率がほぼ0 になるように式(3)のパがかりとなることも期待できるであろう．

20 点配置
図5

20 点配置
図630 点配置AED の最適配置地点と供給効果の密度分布

30 点配置AED の最適配置地点と救命確率の密度分布

0

1

2

(km)

（図 5,6 中の表記については以下の通り．□：AED の最適配置地点，および各々の密度（濃線），
需要量密度（淡線）を示す等値線．背景地図として，河川，鉄道，弘前城を表記．）

なお，図 5 と図 6 における密度分布が非常によく況にあるなか，適切な配置地点についての議論は一似た様相を呈している理由は，最適地が得られる範層重要性を増すことが予想される．現実には，すべ囲内で需要量に大きな変化がないことによる．

ての最適地に AED を配置することは非常に困難ではあるが，問題の重要性からも現状の配置状況の適4.4. 設置数に伴う効果

切性は常に問われるものである．そのためには，何図 7 に，各設置数と式(8)の目的関数値との関係を示す．例えば，「供給効果」を表す目的関数値が 0.01

よりも理論的な最適配置地点，およびその効果については厳密に示されている必要があろう．

となるときでは，最適地に設置された装置によって地域全体の需要量の 1%がカバーされることを示している．

5. おわりに
本稿では，AED の配置地点に対する検討が必要であるとの認識にたち，実際の心停止発生地点から

目的関数値

0.08

なる需要点より推定された確率密度分布に，生存退0.07

院率にもとづく救命確率を加味した AED の最適配

0.06

置問題を提案した．なかでも，地域全体の需要量に

0.05
0.04

対する「供給効果」を定式化するとともに，配置地

0.03

点の座標値に関する多変数最適化問題としてその解0.02を数値的に得た．AED の最適配置問題では，装置0.01がカバーできる範囲が限られているとともに，装置

0
0

5

10

15

20

25

（台）

30

設置数

図 7 設置数と目的関数値までの距離が非常に重大な要因である．本稿では，対象地域内で密度分布が比較的高い場所が広範囲に拡がっていたために，装置の最適配置地点は特定の地域に集中して得られる傾向が顕著であった．しか図 7 を見る限りでは，設置数の増加に従って「供し，心停止発生地点に影響を及ぼすとされる，高齢給効果」はほぼ線型に増加することから，少なくと者の人口分布や，人的そして物的な集積度に大きなも設置数が 30 点までにおいては，需要がほぼ同等

偏りがある地域では，地域全体としての需要量の地に高い地域に装置が割り当てられていると考えられ理的な変化が激しくなる．その際には，求まる最適る．
「供給効果」の変化の割合は，結局のところ需要地点もある程度分散するといった，今回とは異なる

量の分布状況に大きく依存するものであるが，一般結果を得ることも予想される．要は，AED の最適

的には，設置数が少ない初期の段階では急激に増加

配置地点は直感的に導き出すことが困難であるだけ

し，
その後緩やかに増加を続けることが予測できる．

に，十分な検討が望まれる問題であると考えてよい

言い換えれば，設置数が増加するにつれ，需要密度であろう．

の高い地域から需要密度の低い地域へと，装置で新たにカバーする領域が移行して行くといえる．

また，AED の場合に，設置に伴う用地取得や周辺環境への配慮などの問題は，通常の施設立地の場

さらにここで得た結果は，より多くの AED を地合ほど重大ではないと考えられる．装置に関する維域全体として有することで，地域医療に貢献できる

持管理などの問題はあるにせよ，明らかに立地不可ことを意味するとともに，供給効果の上昇の程度が能であるような条件を除けば，この点で大きく性質

鈍らないこの範囲の設置数においては，AED の配が異なることが理解される．つまり，現段階におい置の検討が特に重要とされるべきであることを示唆て候補地として装置の設置地点を事前に限定するのしている．

ではなく，連続的な平面上に得ようとする本稿の試

従来の限られた施設以外にも設置が検討される状

みは，現状に対しても有意義で実効性の高い問題と位置付けて良いと言える．
最後に，今後の課題について述べる．AED の設人口が多い居住地域ほど密度が高くなる傾向を示すことなどから，需要量をよく代表しているものと判断し置が希望される施設として，
「駅」
，
「学校」
，
「公共施

た．

設」などの，現段階で積極的に設置が進んでいる施

2) Sheather and Jones のプラグイン法では，式(1)の設を挙げる声が一般には大きいものの，「コンビニエlf の推定値が，理論的な考察によるバンド幅からではンスストア」あるいは「交番」への設置を希望する

なく，データから推定されるという点に特徴がある．

「コンビニエンスストア」や意見も少なくない 10)．

Simonoff（1996）では，一般に最良なバンド幅選択法「交番」は，地域住民に比較的認知されている施設が発見されていないことを踏まえつつも，プラグイン

であり，例えば深夜のような時間を問わない緊急時

法が「最先端」の方法であるとともに，多方面で有効への対応の面でも期待できる．また，市街地のみなな選択法であると紹介されている．

らず居住地域も含めて幅広く分布していることから

3) 処置に至るまでの時間の経過と後出の生存退院率

も，有効な AED の配置地点として，需要量との関との関係に加えて，道路状況などの地理的な条件を考係性を含めた検討をおこなうことも重要であると考慮した結果，距離の増大に伴う効果の減衰を示す関数

える．これらをふまえて今後は，既存施設への AED式として適当であると判断した．

の設置を視野に入れた配置問題についても考察して

4) 例えば，本稿で想定される AED の場合とは異なり，

いきたい．

救急車の出動が前提となるような事例では，発生地点から最も近くにある立地点が選択されるのが妥当であ

謝辞

ろう．

本研究をおこなうにあたり，国立保健医療科学院

5) 計算上，被覆される領域外でも 1 を超える地点は存

技術評価部丹後俊郎部長，東京大学空間情報科学研

在するが，本稿でも述べるように影響を無視できるほ究センター丸山祐造助教授および東京大学大学院工

ど微少であることから，装置から一定距離以上の地点

学系研究科都市工学専攻貞広幸雄助教授より貴重なでの効果は無視できることが前提とされている．

コメントを頂いた．また，弘前地区消防事務組合に

6) 位置情報の取得に際しては，東京大学空間情報科学

は救急車による搬送データの提供をはじめ，多大な研究センターが提供するアドレスマッチングサービス

ご協力を頂いた．

を用いた．

本研究は，
厚生労働科学研究費補助金（健康科学総7) American Heart Association (2000)によれば，心停合研究事業）「地理及び社会状況を加味した地域分析

止発生から 1 分経過するにつれ，生存退院率が約

方法の開発に関する研究」（代表：浅見泰司）の支援を

7%~10%で減少するとされている．これに基づいて本受けた．ならびに，東京大学空間情報科学研究セン

稿では，10 分経過した時点で生存退院率がほぼ 0 にな

ターの研究用空間データ利用を伴う共同研究（研究るとした．

番号 75）として，以下のデータを利用した．

8) AED の効果が期待できる範囲がごく限られている（株）ゼンリン提供：ZmapTownII 青森県弘前市

ことに加え，対象地域では道路密度の粗密が分析に大記して謝意を表する．

きな影響を与えないと考えられることから，本稿では直線距離を採用した．

補注

9) 最 適 化 計 算 で は ， MATLAB の Optimization

1) 本稿の適用事例で扱っている心停止発生数に関し

Toolbox を使用した．

ては，季節などの時期による発生数の変動の可能性や，

10) AED（自動体外式除細動器）普及を推進する会に設置する AED の数との関係を吟味し，過去 3 年分のよる「AED の設置場所に関するアンケート」より＜411 例を本稿で分析を行うのに十分なサンプル数とし

http://www.aedjapan.com/vote/qnaire2.cgi＞．

た．実際の密度分布を見る限りでも，市街地のなかで参考文献

Science Association, 32, 101-18.

American Heart Association (2000) Guidelines 2000 for

Crocco, T.J., Sayre, M.R., Liu, T., Davis, S.M., Cannon, C.

emergency

and Potluri, J. (2004) Mathematical determination ofcardiovascular care: International consensus onexternal defibrillators needed at mass gatherings.

science, Circulation.

Prehospital Emergency Care, 8, 292-297.

cardiopulmonary

resuscitation

andAlsalloum, O.I. and Rand, G.K. (2006) Extensions to

Malcom, G.E., Thompson, T.M. and Coule, P.L. (2004) The

emergency vehicle location models. Computers &location and incidence of out-of hospital cardiac Arrest

Operations Research, 33, 2725-2743.

in Georgia. Prehospital Emergency Care, 8, 10-14.

Bowman, A. and Azzelini, A. (1997) Applied Smoothing

Sheather, S.J. and Jones, M.C. (1991) A reliable

Techniques for Data Analysis. Oxford: Oxford

data-based bandwidth selection method for kernel

University Press.

density estimation. Journal of the Royal Statistical

Brotcorne, L., Gilbert Laporte, G. and Semet, F. (2003)
Ambulance location and relocation models. European

Society series B, 53, 683-690.
Silverman, B.W. (1986) Density Estimation for Statistics

Journal of Operational Research, 147, 451-463.

and Data Analysis. London: Chapman and Hall.

Church R.L. and ReVelle C.S. (1974) The maximal

Simonoff, J.S. (1996) Smoothing Methods in Statistics.

covering location problem. Papers of the Regional

New York: Springer-Verlag.