Phân phối ổn định dương trong R


9

Phân phối ổn định dương được mô tả bởi bốn tham số: tham số độ lệch , tham số tỷ lệ , tham số vị trí và vì vậy tham số chỉ mục -called \ alpha \ in (0,2] . Khi \ beta bằng 0, phân phối đối xứng quanh \ mu , khi nó dương (tương ứng âm) phân phối bị lệch sang phải (tương ứng bên trái) Phân phối ổn định cho phép đuôi chất béo khi \ alpha giảm.β[1,1]σ>0μ(,)α(0,2]βμα

Khi α hoàn toàn nhỏ hơn một và β=1 sự hỗ trợ của phân phối giới hạn ở (μ,) .

Hàm mật độ chỉ có biểu thức dạng đóng cho một số kết hợp giá trị cụ thể cho các tham số. Khi μ=0 , α<1 , β=1σ=α thì đó là (xem công thức (4.4) tại đây ):

f(y)=1πyk=1Γ(kα+1)k!(yα)ksin(αkπ)

Nó có ý nghĩa vô hạn và phương sai.

Câu hỏi

Tôi muốn sử dụng mật độ đó trong R. Tôi sử dụng

> alpha <- ...
> dstable(y, alpha=alpha, beta=1, gamma=alpha, delta=0, pm=1)

nơi dstable chức năng đi kèm với các gói fBasics.

Bạn có thể xác nhận đây là cách đúng để tính mật độ đó trong R không?

Cảm ơn bạn trước!

BIÊN TẬP

Một lý do tại sao tôi nghi ngờ là, trong đầu ra, giá trị của delta khác với giá trị đầu vào. Thí dụ:

> library(fBasics)
> alpha <- 0.4
> dstable(4, alpha=alpha, beta=1, gamma=alpha, delta=0, pm=1)
[1] 0.02700602
attr(,"control")
   dist alpha beta gamma    delta pm
stable   0.4    1   0.4 0.290617  1

Câu trả lời:


6

Câu trả lời ngắn gọn là của bạn vẫn ổn, nhưng của bạn sai. Để có được phân phối ổn định dương được cung cấp bởi công thức của bạn trong R, bạn cần đặt δγ

γ=|1itan(πα/2)|1/α.

Ví dụ sớm nhất tôi có thể tìm thấy công thức bạn đưa ra là (Feller, 1971), nhưng tôi chỉ tìm thấy cuốn sách đó ở dạng vật lý. Tuy nhiên (Hougaard, 1986) đưa ra công thức tương tự, cùng với biến đổi Laplace Từ hướng dẫn sử dụng ( được sử dụng trong ), việc tham số hóa là từ (Samorodnitsky và Taqqu, 1994), một tài nguyên khác mà việc tái tạo trực tuyến đã lảng tránh tôi. Tuy nhiên (Weron, 2001) cung cấp chức năng đặc trưng trong tham số hóa của Samorodnitsky và Taqqu cho là

L(s)=E[exp(sX)]=exp(sα).
stablediststabledistfBasicspm=1α1
φ(t)=E[exp(itX)]=exp[iδtγα|t|α(1iβsign(t)tanπα2)].
Tôi đã đổi tên một số thông số từ giấy của Weron thành coinside với ký hiệu chúng tôi đang sử dụng. Anh ta sử dụng cho và cho . Trong mọi trường hợp, cắm và , chúng tôi nhận được μδσγβ=1δ=0
φ(t)=exp[γα|t|α(1isign(t)tanπα2)].

Lưu ý rằng cho và rằng . Chính thức, , do đó, bằng cách đặt trong chúng tôi nhận được Một điểm thú vị cần lưu ý là tương ứng với cũng là , vì vậy nếu bạn muốn thử hoặc(1itan(πα/2))/|1itan(πα/2)|=exp(iπα/2)α(0,1)iα=exp(iπα/2)L(s)=φ(is)γ=|1itan(πα/2)|1/αφ(t)

φ(is)=exp(sα)=L(s).
γα=1/21/2γ=αγ=1α, đó thực sự không phải là một xấp xỉ xấu, bạn kết thúc chính xác cho .α=1/2

Đây là một ví dụ trong R để kiểm tra tính chính xác:

library(stabledist)

# Series representation of the density
PSf <- function(x, alpha, K) {
  k <- 1:K
  return(
    -1 / (pi * x) * sum(
      gamma(k * alpha + 1) / factorial(k) * 
        (-x ^ (-alpha)) ^ k * sin(alpha * k * pi)
    )
  )
}

# Derived expression for gamma
g <- function(a) {
  iu <- complex(real=0, imaginary=1)
  return(abs(1 - iu * tan(pi * a / 2)) ^ (-1 / a))
}

x=(1:100)/100
plot(0, xlim=c(0, 1), ylim=c(0, 2), pch='', 
     xlab='x', ylab='f(x)', main="Density Comparison")
legend('topright', legend=c('Series', 'gamma=g(alpha)'),
       lty=c(1, 2), col=c('gray', 'black'),
       lwd=c(5, 2))
text(x=c(0.1, 0.25, 0.7), y=c(1.4, 1.1, 0.7), 
     labels=c(expression(paste(alpha, " = 0.4")),
              expression(paste(alpha, " = 0.5")),
              expression(paste(alpha, " = 0.6"))))

for(a in seq(0.4, 0.6, by=0.1)) {
  y <- vapply(x, PSf, FUN.VALUE=1, alpha=a, K=100)
  lines(x, y, col="gray", lwd=5, lty=1)
  lines(x, dstable(x, alpha=a, beta=1, gamma=g(a), delta=0, pm=1), 
        col="black", lwd=2, lty=2)
}

Lô đầu ra

  1. Feller, W. (1971). Giới thiệu về Lý thuyết xác suất và các ứng dụng của nó , 2 , tái bản lần 2. New York: Wiley.
  2. Hougaard, P. (1986). Các mô hình sinh tồn cho các quần thể không đồng nhất có nguồn gốc từ các bản phân phối ổn định , Biometrika 73 , 387-394.
  3. Samorodnitsky, G., Taqqu, MS (1994). Các quy trình ngẫu nhiên phi Gauss ổn định , Chapman & Hall, New York, 1994.
  4. Ma quỷ, R. (2001). Các bản phân phối ổn định Levy được xem xét lại: chỉ số đuôi> 2 không loại trừ chế độ ổn định Levy , Tạp chí quốc tế về Vật lý hiện đại C, 2001, 12 (2), 209-223.

1
Hân hạnh. Chủ đề về các tham số ổn định tích cực đã gây ra rất nhiều vấn đề đau đầu cho tôi vào đầu năm nay (nó thực sự là một mớ hỗn độn), vì vậy tôi đang đăng những gì tôi nghĩ ra. Hình thức đặc biệt này rất hữu ích trong phân tích sinh tồn vì hình thức Laplacian cho phép mối quan hệ đơn giản giữa các tham số hồi quy có điều kiện và biên trong các mô hình mối nguy theo tỷ lệ khi có một thuật ngữ yếu sau phân phối ổn định dương (xem bài viết của Hougaard).
P Schnell

6

Những gì tôi nghĩ đang xảy ra là trong đầu ra deltacó thể báo cáo một giá trị vị trí nội bộ, trong khi ở đầu vào deltalà mô tả sự thay đổi. [Dường như có vấn đề tương tự với gammakhi nào pm=2.] Vì vậy, nếu bạn cố gắng tăng ca lên 2

> dstable(4, alpha=0.4, beta=1, gamma=0.4, delta=2, pm=1)
[1] 0.06569375
attr(,"control")
   dist alpha beta gamma    delta pm
 stable   0.4    1   0.4 2.290617  1

sau đó bạn thêm 2 vào giá trị vị trí.

Với beta=1pm=1bạn có một biến ngẫu nhiên dương với phân phối giới hạn dưới 0.

> min(rstable(100000, alpha=0.4, beta=1, gamma=0.4, delta=0, pm=1))
[1] 0.002666507

Thay đổi 2 và giới hạn dưới tăng cùng một lượng

> min(rstable(100000, alpha=0.4, beta=1, gamma=0.4, delta=2, pm=1))
[1] 2.003286

Nhưng nếu bạn muốn deltađầu vào là giá trị vị trí bên trong thay vì thay đổi hoặc giới hạn dưới, thì bạn cần sử dụng một thông số kỹ thuật khác cho các tham số. Ví dụ: nếu bạn thử các cách sau (với pm=3và thử delta=0delta=0.290617bạn đã tìm thấy trước đó), bạn dường như nhận được cùng một thông deltatin. Với pm=3delta=0.290617bạn có cùng mật độ 0,02700602 bạn đã tìm thấy trước đó và giới hạn dưới ở mức 0. Với pm=3delta=0bạn nhận được giới hạn dưới âm (thực tế là -0.290617).

> dstable(4, alpha=0.4, beta=1, gamma=0.4, delta=0, pm=3)
[1] 0.02464434
attr(,"control")
   dist alpha beta gamma delta pm
 stable   0.4    1   0.4     0  3
> dstable(4, alpha=0.4, beta=1, gamma=0.4, delta=0.290617, pm=3)
[1] 0.02700602
attr(,"control")
   dist alpha beta gamma    delta pm
 stable   0.4    1   0.4 0.290617  3
> min(rstable(100000, alpha=0.4, beta=1, gamma=0.4, delta=0, pm=3))
[1] -0.2876658
> min(rstable(100000, alpha=0.4, beta=1, gamma=0.4, delta=0.290617, pm=3))
[1] 0.004303485

Bạn có thể thấy dễ dàng hơn để bỏ qua deltađầu ra, và miễn là bạn giữ beta=1thì sử dụng pm=1phương tiện deltatrong đầu vào là phân phối giới hạn dưới, có vẻ như bạn muốn là 0.


4

Cũng cần lưu ý: Martin Maechler chỉ tái cấu trúc mã cho phân phối ổn định và thêm một số cải tiến.

Stablesist gói mới của anh ấy cũng sẽ được sử dụng bởi fBasics, vì vậy bạn có thể muốn cung cấp cho cái nhìn này là tốt.

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.