Jurnal Matematika Vol. 2 No. 2, Desember 2012. ISSN : 1693-1394

Pendekatan-pendekatan Titik Sadel: Suatu Ka jian Pustaka

I Wayan Sumarjaya

Jurusan Matematika FMIPA, Universitas Udayana Kampus Bukit Jimbaran Badung, Bali e-mail: sumarjaya@unud.ac.id

Abstract: Dalam banyak area riset seperti biologi dan ilmu kedokteran, sering dijumpai jumlah amatan terlalu sedikit atau informasi yang diperoleh sangat sedikit atau terbatas. Inferensi dengan menggunakan teorema limit pusat untuk menyimpulkan statistik yang diteliti mendekati normal tentu tidaklah tepat. Salah satu cara untuk mengatasi masalah ini adalah dengan menggunakan pendekatan titik sadel (saddlepoint approximation). Artikel ini mengulas konsep dan beberapa perkembangan penting pendekatan titik sadel seperti pendekatan titik sadel untuk rata-rata, fungsi distribusi kumulatif, distribusi bersyarat, statistik umum, dan distribusi posterior.

Keywords: pendekatan titik sadel, ekspansi Edgeworth, asimptotik sampel kecil, asimtotik tingkat tinggi

  • 1.    Pendahuluan

Dalam banyak area riset seperti biologi dan ilmu kedokteran, sering dijumpai jumlah amatan terlalu sedikit atau informasi yang diperoleh sangat sedikit atau terbatas. Salah satu cara untuk mengatasi masalah ini adalah dengan menggunakan teorema limit pusat untuk menyimpulkan statistik yang sedang diteliti mendekati normal. Pendekatan ini disebut asimptotik tingkat pertama (first order asymptotics). Namun, pendekatan ini tidaklah tepat karena tidak memberikan hasil yang akurat (lihat [11]).

Metode lain untuk menyelesaikan permasalahan di atas adalah dengan menggunakan ekspansi asimptotik (asymptotic expansion) atau asimptotik kemungkinan modern (modern likelihood asymptotics). Konsep dasar asimptotik kemungkinan modern adalah dengan memodifikasi teorema limit pusat ”naif” dengan ekspansi tingkat tinggi. Secara garis besar terdapat dua ekspansi asimptotik penting dan saling berhubungan yaitu ekspansi Edgeworth dan pendekatan titik sadel (saddlepoint approximation).

Ekspansi Edgeworth yang berbentuk

k

F˜n(x) = Φ(x) +

i=1


qi(x) ni


+ o(n k/2)


(1)


memberikan pendekatan yang baik pada pusat distribusi. Namun, ekspansi Edgeworth tidak bagus pada daerah ekor (tail) distribusi. Hal ini dapat dilihat dengan munculnya kenegatifan pada ekor yang disebabkan oleh suku polinomial (lihat [15]). Salah satu upaya untuk mengatasi masalah kenegatifan ini adalah menggunakan pendekatan titik sadel.

Artikel ini diatur sebagai berikut. Bagian 2 akan membahas konsep dasar pendekatan titik sadel dalam aplikasi statistika. Pada bagian 3, akan dibahas beberapa perkembangan dalam pendekatan titik sadel seperti pendekatan titik sadel untuk rata-rata, fungsi distribusi kumulatif, distribusi bersyarat, statistik umum, dan distribusi posterior. Selanjutnya, aplikasi pendekatan titik sadel yang saat ini aktif dikembangkan dan digunakan akan dibahas pada bagian 4. Isu-isu yang berhubungan dengan pendekatan titik sadel akan dibahas pada bagian 5. Akhirnya, pada bagian 6 artikel ini akan memberikan simpulan.

  • 2.    Pendekatan Titik Sadel

Pada bagian ini akan dibahas konsep dasar pendekatan titik sadel. Secara umum pendekatan titik sadel adalah bagian dari analisis asimptotik. Aplikasi titik sadel dalam konteks statistika diusulkan oleh [15]. Pengkajian tentang titik sadel dan aplikasinya dalam bidang statistika telah diulas secara mendalam oleh [8] dan [32]. Selain itu terdapat pula monograf dan buku yang membahas pendekatan titik sadel dari tingkat dasar sampai lanjut, misalnya [21], [25], [26], dan [14].

Pendekatan titik sadel dapat dilihat dari beberapa perspektif: metode Laplace, exponential tilting, dan ekspansi Edgeworth. Pada bagian ini akan dibahas penurunan pendekatan titik sadel menggunakan metode Laplace. Pembahasan berikut diadaptasi dari artikel [24]. Ide pendekatan titik sadel tersebut adalah sebagai berikut. Pertama, lakukan pendekatan terhadap fungsi positif f (x) menggunakan deret Taylor dan expan-sikan integran dan kenali bahwa kernel dari integran adalah sama dengan kernel dari distribusi normal. Langkah ini disebut dengan pendekatan Laplace. Langkah berikutnya adalah mengombinasikan pendekatan Taylor dan pendekatan Laplace.

Misalkan akan didekati fungsi positif f (x) yang tidak mesti fungsi densitas peluang. Misalkan g(x) = log f(x). Tulis f(x) = exp(g(x)) dan gunakan ekspansi Taylor pada x0 sehingga diperoleh

f(x) exp


g(x0) + (x - x0)g(x0)


+


(x - χ0)2 „    )

2      g′′(x0) .


(2)


Untuk menyederhanakan pendekatan pilih x0 = xˆ dengan g(xˆ) = 0, sehingga ekspansi pada sisi kanan (2) menjadi

f(x) exp g(xˆ) + (x -2x) g′′(xˆ) .

(3)


Dapat dilihat pendekatan pada (3) adalah pasti jika g(x) adalah fungsi kuadratik. Namun, jika g(x) bukan fungsi kuadratik tentu saja suku-suku ekpansi tingkat yang lebih tingi seperti (x - xˆ)3 atau (x - xˆ)4 akan berpengaruh terhadap pendekatan dan hasilnya tidak akan bagus [24]. Lakukan ekspansi integran pada (3) sehingga diperoleh

f (x) dx


exp


g(xˆ) +


(x -2 xˆ)2 g′′(xˆ) dx.


(4)


Jika xˆ pada (4) maksimum dan g′′(xˆ) adalah negatif, dapat dilihat bahwa kernel dari integral tersebut adalah sama dengan kernel densitas normal N (xˆ, -1/g′′(xˆ)). Selanjutnya dari sifat densitas normal dengan rata-rata μ dan varians σ2 diperoleh

y exp(—(x — μ)2/2σ2) dx = (2πσ2)1 /2.                      (5)

Lakukan ekspansi pada sisi kanan persamaan (4) sehingga diperoleh

f(x) dx


exp g(x)( -⅜ )  

g (x)


(6)


Teknik yang digunakan pada (6) disebut pendekatan Laplace. Pada langkah berikut akan dikombinasikan ekspansi Taylor dan pendekatan Laplace. Tulis f sebagai

f(x) =


h(x, t) dt


(7)


untuk beberapa h(x, t) positif. Definisikan k(x, t) = log h(x, t) dan dengan menggunakan (3) serta (4) diperoleh

(t — £(x))2 2k(x,t)

f(x)   exp k(x,tˆ(x)) +                           dt

2 dt 2     (χ x)J

= exp[k (x,l( χ ))](- ∂ 2 k ( Xt)/∂t 2


1 / 2


*(x)


(8)


untuk setiap x, t(x) memenuhi ∂k(x,t)/∂t = 0 dan 2k(x,t)/∂t2 < 0 yang memaksimalkan k(x, t). Teknik ini selanjutnya diterapkan pada densitas f (x). Misalkan ϕX (t) adalah fungsi pembangkit momen dari f (x). Formula inversi (lihat [20]) memberikan

1

f(x) = 2π -∞ ϕX (it) exp(itx) dt

=         exp(ψX (it) - itx) dt,

(9)


2π -∞

dengan i = -1 dan ψX (it) adalah fungsi pembangkit kumulan (cumulant generating function). Misalkan t = it; selanjutnya persamaan (9) dapat ditulis sebagai

f(x) = 1

2πi τ



exp(ψX (t) - tx) dt.


(10)


Sekarang misalkan k(x, t) = ψX (it) - tx dan cari titik tˆ(x) yang memenuhi

ψX (tˆ(x)) = x.                                    (11)

Lakukan ekspansi ψX (t) - tx di sekitar tˆ(x) sehingga diperoleh

ψX (t) - tx ≈ ψX (tˆ(x)) - tˆ(x)x + (t - t2(x)) ψX′′ (tˆ(x)).              (12)

Fig 1. Plot densitas N (0,1) dan pendekatan titik sadel.


Substitusikan (12) ke dalam (9) dan mengintegralkan f(x) diperoleh

∕     1     V / 2

f (x) o                 exP[ψx(^(x)) - ^(x)x]

(13)


2 ∏ΨX (t (x))√

Titik t(x) pada (11) disebut titik sadel (saddlepoint) dari ψχ (t) — tx dan f (x) disebut densitas titik sadel.

Contoh 1. Misalkan peubah acak X berdistribusi N (0, 1). Kumulan distribusi ini adalah ψx(t) = t2/2 untuk t (-∞, ∞). Titik sadel dalam hal ini adalah £ = x. Dengan demikian diperoleh densitas titik sadel

f(x) = φ(x) =


1

2π


exp -x22


(14)


yang dalam hal ini tepat sama dengan densitas N (0, 1). Gambar 1 memperlihatkan plot densitas N(0, 1) dan densitas titik sadel.

Contoh 2. Misalkan perubah acak X berdistribusi Gamma(θ, 1). Kumulan X adalah ψx (t) = log(1 — t), t (-∞, 1) dan titik sadel tunggalnya adalah it = 1 — θ∕x. Selanjutnya diperoleh pendekatan titik sadel untuk Gamma(θ, 1) adalah

f (x) = (√2πθθ 1 /2 exp()) 1 xθ 1 exp(—x), x> 0.             (15)

Gambar 2 memperlihatkan plot densitas Gamma(2, 1) dan densitas titik sadel. Terlihat bahwa bentuk pendekatan titik sadel f sama dengan f , hanya berbeda dalam konstanta penormalan. Pendekatan Stirling untuk Γ (θ) (lihat [2]) memberikan

f (θ) = √2πθθ-1 /2 exp() ≈ Γ(θ),

(16)


sehingga dapat dilihat


fˆ(x) =   ( ) f (x).

Γ(θ)


(17)


densitas Gamma(2,1) densitas titik sadel

024


10


Fig 2. Plot densitas Gamma(2,1) dan pendekatan titik sadel.

  • 3.    Beberapa Perkembangan Pendekatan Titik Sadel

    • 3.1.    Pendekatan titik sadel untuk rata-rata sampel

Artikel [15] memperkenalkan penggunaan pendekatan titik sadel pada densitas rata-rata sampel untuk peubah acak bebas dan berdistribusi identik Xι,... ,Xn. Misalkan X adalah rata-rata sampel, φX (t) adalah fungsi pembangkit momen sampel, dan ψχ (t) adalah fungsi pembangkit kumulan sampel. Dengan menerapkan teknik titik sadel pada bagian sebelumnya diperoleh densitas rata-rata sampel sebagai berikut

n /2

fXX(x) ≈ θ ffif-w     exp{n[ψX(^(x)) - ^(x)x]}

(18)


2πψX′′ (t(x))

dengan ψX (tˆ(x)) = x¯.

Contoh 3. Misalkan Xι,...,Xn berdistribusi bebas dan identik N(0, 1). Titik sadel tunggal yang menyelesaikan persamaan t = x¯ memberikan pendekatan

fX      -∏ exp(- nτ)


(19)


yang tepat sama dengan densitas eksak dari X yaitu N(0, 1/n).

Contoh 4. Misalkan X1,..., Xn berdistribusi bebas dan identik Poisson(θ). Fungsi pembangkit momennya adalah φ(t) = exp[θ(exp(t) 1)] dan fungsi pembangkit kumulan-nya adalah ψ(t) = θ(exp(t) 1). Titik sadel tunggal yang menyelesaikan persamaan θ exp(tˆ) = x¯ adalah tˆ = log(x¯). Dengan demikian diperoleh pendekatan titik sadel

fX(X) = ^2^ exp(-n)


θnx

.

- nx + o

x +2


(20)


  • 3.2.    Pendekatan titik sadel untuk ekor distribusi

Pada banyak aplikasi statistika biasanya yang menjadi daya tarik dalam melakukan inferensi adalah menghitung peluang ekor distribusi (tail probability). Sebagai contoh bidang keuangan, ilmu aktuaria, atau kejadian langka (rare event) biasanya yang menjadi objek penelitian adalah bagaimana menghitung peluang ekor atau fungsi distribusi kumulatif.

Beberapa pendekatan untuk fungsi distribusi telah dilakukan. Sebagai contoh pendekatan simpangan besar (large deviation) dari [6], transformasi Laplace dari deret Edgeworth [33], pendekatan integrasi ruang sampel [36], metode penyesuaian [7], dan metode penurunan tercuram (method of steepest descent) [28].

Pada bagian ini akan dibicarakan salah satu pendekatan untuk fungsi distribusi yaitu metode yang diusulkan oleh [28]. Metode ini sering disebut formula Lugannani-Rice atau formula ajaib dalam literatur statistika asimptotik. Formula ini memberikan hasil yang sangat akurat dalam mendekati fungsi distribusi. Ide pendekatan fungsi distribusi yang diusulkan [28] ini adalah sebagai berikut. Pertama, nyatakan pendekatan fungsi distribusi kumulatif sebagai integral kumulatif berhingga dari densitas titik sadel

1          1/2

f(x)                      exp[ψX (t(x)) - t(x)x].

(21)


2πψX′′ (t(x))

Selanjutnya, gunakan pendekatan Temme dengan bentuk

w0


h(w)φ(w) dw


h(0)Φ(w0) +φ(w0) h(0) - h(w0) , w0


(22)


dengan Φ dan φ masing-masing menyatakan fungsi distribusi kumulatif dan fungsi den-sitas normal. Pendekatan terhadap fungsi distribusi kumulatif menjadi

f(x) dx

-∞


y            1          1/2

=            ′′            exp[ψX (t(x)) - t(x)x] dx

-∞ 2πψX′′ (t(x))

y1

(23)


= -∞ψX′′(tˆ(x))φ(w)dx,

dengan ketergantungan terhadap tˆ dan wˆ dikendalikan (lihat [14]). Jadi diperoleh

dengan


y 1                 w      wˆ

′′ ˆ   φ(wˆ)dx =                 (wˆ)dwˆ

-∞ ψX (t(x))              -∞ tˆ ψX′′ (tˆ(x))

= Φ(wˆy) + φ(wˆy) 1 - h(wy), wy


h(0) = lim tˆ0


wˆ

t√ψ(t)


= 1.


(24)

(25)


Dengan demikian, pendekatan titik sadel untuk fungsi distribusi yang diusulkan oleh

[28] adalah


φ(w) + φ(w)(1 - Γ), wu

F (x) =   1 +      ψX′′′ (0)

κ 2   6√2π(ψχ(0))3/2 ,


jika x ̸= µ;

jika x = µ;


(26)


dengan wˆ =


sign(tˆ) 2(tˆx - ψX (tˆ))


uˆ = tˆ  ψX′′ (tˆ)


dan µ = EX . Untuk proses penu-


runan untuk pendekatan fungsi distribusi secara lebih rinci dapat dilihat pada [28], [25],


[26], dan [14].


  • 3.3.    Pendekatan titik sadel untuk distribusi bersyarat

Salah satu pengembangan penting dalam pendekatan titik sadel adalah pendekatan un-

tuk distribusi bersyarat. [34] merampatkan formula Lugannani-Rice untuk distribusi bersyarat berbentuk Pr(Y ≥ y|X = x) dengan (X, Y ) adalah rata-rata dari n vek-

tor acak bivariat bebas. Ide dari pendekatan ini adalah sebagai berikut. Misalkan (X1 , Y1 ), . . . , (Xn , Yn ) adalah peubah-peubah acak bivariat yang berdistribusi identik dan bebas dengan denstias f (x, y) dan fungsi pembangkit kumulan κ yang diberikan

oleh

κ(s, u) = ln


f(x, y) exp(sx + uy) du ds,


s, u C,


(27)


yang diasumsikan ada pada persekitaran (0, 0). Misalkan X = n-1 Xi dan Y = n-1 Yi dan tetapkan titik (x˜, y˜). Misalkan turunan kumulan κ dituliskan sebagai berikut

κs(s, u) =

sκ(s,u),

(28)

κu(s, u) =

∂uκ(s,u),

(29)

2

κss(s, u) =

∂s2 κ(s, u).

(30)

Selanjutnya, [34] memperoleh pendekatan titik sadel untuk distribusi bersyarat di atas

Pr(Y¯ ≥ y|X¯ = x) 1 - Φ(nw˜) + φ(nw˜)


[κss(s˜0, 0)]1/2

nu˜(s˜, u˜)|1/2


- ( nw˜)-1


(31)


dengan s˜0 adalah solusi dari persamaan titik sadel

κs(s˜0, 0) = x˜, s˜0 R,

(32)


Φ adalah fungsi distribusi kumulatif normal, φ adalah fungsi densitas normal, (s˜, u˜)| adalah determinan dari κ(s˜, u˜), dan

w˜ = sign(u˜){2[s˜x˜ + u˜y˜ - κ(s˜u˜)] - 2[s˜0x˜ - κ(s˜0, 0)]}1/2.             (33)

Beberapa kondisi yang diperlukan untuk pendekatan di atas adalah sebagai berikut (lihat [34]). Pertama, beberapa pangkat positif fungsi karakteristik untuk (Xi, Yi) adalah terintegralkan. Kedua, fungsi pembangkit kumulan κ untuk (Xi, Yi) ada pada perseki-taran (0, 0). Ketiga, matriks varians κ(0, 0) adalah taksingular.

  • 3.4.    Pendekatan titik sadel umum

Pada subbagian sebelumnya telah dibicarakan konsep pendekatan titik sadel dan aplikasinya, terutama untuk rata-rata sampel dan fungsi distribusi. Untuk statistik tertentu, katakanlah Vn, [19] mengusulkan pendekatan titik sadel umum (general saddlepoint approximation). Pendekatan ini menggunakan keempat kumulant dari statistik Vn.

Konsep pendekatan titik sadel umum yang diusulkan [19] adalah sebagai berikut. Misalkan X1 , . . . , Xn adalah peubah acak real bebas dan berdistribusi identik dengan fungsi densitas peluang f dan Vn(X1, . . . , Xn) adalah statistik bernilai real dengan den-sitas fn, ϕn(t) adalah fungsi pembangkit momen, dan ψn(t) adalah fungsi pembangkit kumulan. Anggap fungsi pembangkit momen ada untuk bilangan real t di dalam selang yang taklenyap (non-vanishing) yang mengandung titik asal. Misalkan

Rn(t)=ψn(nt).(34)

n

Selanjutnya akan didekati (34) dengan teknik titik sadel. Anggap ekpansi Edgeworth fn untuk fn ada dan meliputi tingkat suku (order term) n-1 , yaitu,

fn(x) = f˜n(x) + o(n-1),(35)

dengan

fn(x) = φ(x) [1 + κ3n (x3 - 3x) + κ4n (x4 - 6x2 + 3) + κ3n (x6 - 15x4 + 45x2 - 15)] (36) 6             2472

serta κ3n dan κ4n merupakan bilangan-bilangan yang diketahui. Tujuannya adalah untuk mendekati Rn(t) dengan empat kumulan pertama dengan bentuk

Rn(t) = µnt+


n2 t2 2


+ κ3nσn3n2t3 +6


κ4nσn4 n3t4 24


(37)


dengan µn adalah rata-rata dari statistik Vn , σn2 adalah varians dari statistik Vn , κ3n dan κ4n masing-masing adalah kumulan ketiga dan keempat dari statistik Vn . Secara umum µn dan σn tidak diketahui dengan pasti. Mensubstitusi µn dan σn dengan

µn = µ + a1 + o(n-1), n

dan

σn = n1σ/2 + nb31/2 + o(n-3/2)

(38)


(39)


tidak mengubah tingkat (order) pendekatan [19].

Dengan menggunakan formula inversi pada [20], yaitu,


n  τ++

fn(x) = 2πi τ-i exp[n(Rn(t) - tx)] dt.


(40)


dan menerapkan teknik titik sadel pada (40), [19] memperoleh pendekatan titik sadel


untuk fn dengan galat seragam dengan tingkat n



1 yang berbentuk


gˆn(x) =       ˜n′′ ˆ exp[n(R˜n(tˆ) - tˆx)],

2πRn (t)


(41)


dengan t adalah titik sadel yang dihitung sebagai akar dari persamaan


A

Rn (t) = x,


(42)


˜


˜


˜


dengan Rn(t) dan Rn (t) menyatakan turunan pertama dan kedua dari Rn(t).


3.5. Pendekatan titik sadel untuk momen posterior


Dalam inferensi Bayes masalah utama yang biasanya dihadapi adalah menghitung distribusi posterior. Untuk analisis Bayes univariat jika prior yang dipilih adalah prior sekawan (conjugate prior), distribusi posteriornya juga biasanya sekawan. Namun, pada banyak aplikasi penggunaan prior sekawan seringkali tidaklah tepat atau terlalu menyederhanakan masalah. Salah satu cara untuk mendekati distribusi posterior adalah dengan pendekatan titik sadel. [37] mengusulkan pendekatan untuk rata-rata dan varians posterior dari fungsi positif dari suatu parameter real atau parameter bernilai vektor. Untuk pendekatan ini kondisi utama yang diperlukan adalah perkalian kemungkinan (likelihood) dengan prior bermodus tunggal (unimodal).

Ide pendekatan yang diusulkan [37] adalah sebagai berikut. Misalkan h adalah fungsi mulus positif pada ruang parameter. Nilai tengah posterior h(θ) dapat ditulis sebagai


E n [ h ] =E[ h (θ) X(n)] =


h(θ) exp((θ))p(θ) dθ exp((θ))p(θ) dθ


(43)


dengan (θ) adalah fungsi log-kemungkinan dan p adalah densitas prior. Sekarang misalkan L = logp + ℓ/n dan L = log h + logp + ℓ/n. Selanjutnya persamaan (43) dapat


ditulis sebagai


exp(nL)dθ

En[h] = ∫ exp(nL) dθ .


(44)


Kemudian dengan pendekatan Laplace, [37] memperoleh pendekatan untuk (44)


En[h] = (σ) exp[n(L(θ) - L(θ))]


(45)


dengan θ adalah modus dari L, σ2 = 1 /L(θ), θ* adalah modus dari L*, dan σ*2 = A

-1/L′′(θ). Lebih lanjut [37] mengembangkan pendekatan di atas untuk kasus multiparameter, distribusi peramalan (predictive distribution), dan densitas posterior marginal.


  • 4.    Aplikasi Pendekatan Titik Sadel

Penelitian dan aplikasi pendekatan titik sadel beragam, baik dalam bidang statistika teoretis maupun terapannya. Berikut ini beberapa contoh pengembangan dan aplikasi pendekatan titik sadel sejak tahun 2000. Untuk ulasan aplikasi titik sadel sampai tahun 2000 silakan lihat [35]. Daftar berikut tidaklah lengkap.

Pada bidang statistika teoretis misalnya dikembangkan pendekatan titik sadel untuk pendekatan konstanta penormalan distribusi Fisher–Bingham [27]; pendekatan ekor atas statistik Bagai [29]; uji kebebasan pada kelas statistik rank [1]; peluang ekor rasio Student [3]; distribusi penduga fungsi eigen matriks kovarians dan korelasi [10]; dan pendekatan distribusi takterpusat t [12].

Pada bidang terapan misalnya pendekatan titik sadel digunakan pada metode estimasi peluang dalam analisis reliabilitas,[17], [40], [39]; bagan kendali nonparametrik berdasarkan statistik Mood [30]; distribusi jumlah klaim total diskonto [22]; pendugaan peluang distribusi peubah demografi dalam penelitian hidupan liar (wild life) [16]; opti-misasi jaringan sensor [4],[5]; dan penghitungan risiko portofolio [13].

  • 5.    Beberapa Isu Penting dan Solusi untuk Pendekatan Titik Sadel

Pada data yang berisi banyak nol (0) pendekatan titik sadel mungkin tidak akan memberikan hasil yang memuaskan (lihat [18]). Salah satu solusinya adalah dengan melakukan penyesuaian pada fungsi varians pangkat (power variance function).

Pendekatan titik sadel umum sebagai pengembangan dari pendekatan titik sadel biasa telah memungkinkan untuk mendekati densitas peluang untuk sembarang statistik Vn . Hal ini tentu saja didukung dengan asumsi hanya menggunakan keempat kumu-lan dari statistik dan ekspansi Edgeworth untuk statistik tersebut ada. Sebagaimana diklaim oleh [19] tingkat besaran (order of magnitude) pada pendekatan titik sadel umum terjaga (preserved). Namun, beberapa masalah muncul. Pertama, beberapa hasil komputasi numerik dalam artikel [19] justru memperlihatkan hasil numerik yang tidak meyakinkan (lihat juga [9]). Masalah kedua yang muncul adalah tidak terdapatnya solusi tunggal pada persamaan titik sadel pada (42). Untuk memperoleh solusi tunggal [38] menyarankan mengganti (37) dengan

Hn(t; b) = κ 1 nt + κ2nt2 + ( κ3nt3 + κ4nt4) gb(t) (46) 2! 3! 4!

dengan gb(t) = exp(2nb2t2/2n) dan b > 0 adalah konstanta yang dipilih dengan wajar. Masalah ketiga adalah tidak selalu muncul full support [23]. Masalah terakhir yang muncul adalah kinerja pendekatan titik sadel umum yang hanya bagus pada shrinking neighbourhood nilai tengah apabila ukuran sampel dinaikkan (lihat [31]). Untuk mengatasi masalah tersebut [31] menyarankan untuk selalu menggunakan fungsi pembangkit momen sebenarnya (jika ada) untuk statistik yang ingin didekati.

  • 6.    Kesimpulan

Secara umum pendekatan titik sadel akan memberikan hasil yang sangat akurat karena pendekatan ini memiliki ”tingkah laku” seperti densitas peluang. Penggunaan pen-

dekatan titik sadel umum akan lebih bagus apabila diketahui kumulan dari statistik yang ingin diteliti.

Daftar Pustaka

  • [1]    Ehab F. Abd-Elfattah. 2009. Testing for Independence: Saddlepoint Approximation to Associated Permutation Distributions. Electronic Journal of Statistics. 3: 625– 632.

  • [2]    Milton Abramowitz dan Irene A. Stegun (Editors). 1964. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. National Bureau of Standards. Washington.

  • [3]    Kingsley Agho, Wen Dai, dan John Robinson. 2005. Empirical Saddlepoint Approximations of the Studentized Ratio and Regression Estimates for Finite Populations. Statistics and Probability Letters. 71: 237–247.

  • [4]    Saeed A. Aldosari dan Jos´e M. F. Moura. 2005. Saddlepoint Approximation for Sensor Network Optimization. ICASSP 2005. 741–744.

  • [5]    Saeed A. Aldosari dan Jos´e M. F. Moura. 2007. Detection in Sensor Networks: The Saddlepoint Approximation. IEEE Transactions on Signal Processing. 55: 327–340.

  • [6]    R. R. Bahadur dan R. Ranga Rao. 1960. On deviations of the sample mean. The Annals of Mathematical Statistics. 31: 1015–1027.

  • [7]    O. Barndorff-Nielsen. 1986. Inference on Full or Partial Parameters Based on the Standardized Signed Log Likelihood Ratio. Biometrika. 73: 307–322.

  • [8]    O. Barndorff-Nielsen dan D. R. Cox. 1979. Edgeworth and Saddle-point Approximations with Statistical Applications. Journal of the Royal Statistical Society. Series B( Methodological). 41: 279–312.

  • [9]    M. Bogdan. 1994. Asymptotic Distributions of Linear Combinations of Order Statistics. Applicationes Mathematicae. 22: 201–225.

  • [10]    Robert J. Boik. 2005. Second-order Accurate Inference on Eigenvalues of Covariance and Correlation Matrices. Journal of Multivariate Analysis. 96: 136–171.

  • [11]    Allessandra Brazzale, Anthony C. Davison, dan Nancy Reid.2007. Applied Asymptotics: Case Studies in Small-sample Statistics. Cambridge University Press. Cambridge.

  • [12]    Simon Broda dan Marc S. Paolella. 2007. Saddlepoint Approximations for the Doubly Noncentral t Distribution. Computational Statistics and Data Analysis. 51: 2907–2918.

  • [13]    Simon A. Broda dan Marc S. Paolella. 2009. CHICAGO: A Fast and Accurate Method for Portfolio Risk Calculation. Journal of Financial Econometrics. 7: 412– 436.

  • [14]    Ronald W. Butler. 2007. Saddlepoint Approximations with Applications. Cambridge University Press. Cambridge.

  • [15]    H. E. Daniels. 1954. Saddlepoint Approximations in Statistics. The Annals of Mathematical Statistics. 25: 631–650.

  • [16]    Stephen J. Demaso, Fidel Herna´andez, Leonard A. Brennan, dan Ralph L. Bingham. 2011. Application of the Simple Saddlepoint Approximation to Estimate

Probability Distributions in Wildlife Research. Journal of Wild life Management. 75: 740–746.

  • [17]    Xiaopin Du. 2010. System Reliability Analysis with Saddlepoint Approximation. Struct. Multidisc. Optim. 42: 193–208.

  • [18]    Peter K. Dunn dan Gordon K. Smyth. 2005. Series Evaluation of Tweedie Exponential Dispersion Model Densities. Statistics and Computing. 15: 267–280.

  • [19]    George S. Easton dan Elvezio Ronchetti. 1986. General Saddlepoint Approximations with Application to L Statistics. Journal of the American Statistical Association. 81: 420–430.

  • [20]    William Feller. 1966. An Introduction to Probability Theory and Its Applications. Volume II. John Wiley and Sons. New York.

  • [21]    Christopher Field and Elvezio Ronchetti. 1990. Smal l Sample Asymptotics. Institute of Mathematical Statistics. California.

  • [22]    Riccardo Gatto. 2010. A Saddlepoint Approximation to the Distribution of Inho-mogenous Discounted Compounded Poisson Processes. Methodol. Comput. Appl. Probab. 12: 533–551.

  • [23]    Colin S. Gillespie dan Eric Renshaw. 2007. An Improved Saddlepoint Approximation. Mathematical Biosciences. 208: 359–374.

  • [24]    Constantino Goutis dan George Casella. 1999. Explaining the Saddlepoint Approximation. The American Statistician. 53: 216–224.

  • [25]    Jens Ledet Jensen. 1995. Sadd lepoint Approximations. Oxford University Press. Oxford.

  • [26]    John E. Kollassa. 2006. Series Approximation Method in Statistics. Edisi ketiga. Springer. New York.

  • [27]    A. Kume dan S. G. Walker. 2009. On the Fisher-Bingham Distribution. Stat. Com-put.. 19: 167–172.

  • [28]    Robert Lugannani dan Stephen Rice. 1980. Saddle Point Approximation for the Distribution of the Sum of Independent Random Variables. Advances in Applied Probability. 12: 475–490.

  • [29]    Hidetoshi Murakami. 2010. A Numerical Comparison of the Normal and Some Saddlepoint Approximations to a Distribution-Free Test for Stochastic Ordering in the Competing Risks Model. Comput. Stat.. 25: 633–643.

  • [30]    Hidetoshi Murakami dan Takashi Matsuki. 2009. A Nonparametric Control Chart based on the Mood Statistic for Dispersion. Int. J. Adv. Manuf. Technol.. 49: 757– 763.

  • [31]    Spiridon Penev dan John Robinson. 2012. On the Validity of the General Saddlepoint Method. Naskah tidak diterbitkan. The University of New South Wales. Australia.

  • [32]    N. Reid. 1988. Saddlepoint Methods and Statistical Inference. Statistical Science. 3: 213–238.

  • [33]    J. Robinson. 1982. Saddlepoint Approximations for Permutation Tests and Confidence Intervals. Journal of the Royal Statistical Society. Series B (Methodological). 44: 91–101.

  • [34]    Ib M. Skovgaard. 1987. Saddlepoint Expansions for Conditional Distributions. Journal of Applied Probability. 24: 875–887.

  • [35]    Robert L. Strawderman. 2000. Higher-order Asymptotic Approximation: Laplace, Saddlepoint, and Related Methods. Journal of the American Statistical Association. 95: 1358–1364.

  • [36]    N. M. Temme. 1982. The Uniform Asymptotic Expansion of a Class of Integrals Related to Cumulative Distribution Functions. SIAM Journal on Mathematical Analysis. 13: 239–252.

  • [37]    Luke Tierney dan Joseph B. Kadane. 1986. Accurate Approximations for Posterior Moments and Marginal Densities. Journal of the American Statistical Association. 81: 82–86.

  • [38]    Suojin Wang. 1992. General Saddlepoint Approximations in the Bootstrap. Statistics and Probability Letters. 13: 61–66.

  • [39]    Ning-Cong Xiao, Hong-Zhong Huang, Zhonglai Wang, Yu Liu, Xiao-Ling Zhang. 2012. Unified Uncertainty Analysis by the Mean Value First Order Saddlepoint Approximation. Struct. Multidisc. Optim. 46: 803–812.

  • [40]    Zhimin Zi, Chao Hu, dan Byeng D. Youn. 2012. A Comparative Study of Probability Estimation Methods for Reliability Analysis. Struct. Multidisc Optim. 45: 33–52.