Làm thế nào để tính toán các dải dự đoán cho hồi quy phi tuyến tính?


15

Các trang trợ giúp cho Prism cung cấp cho các giải thích sau đây để biết nó như thế nào tính toán các ban nhạc dự đoán cho hồi quy phi tuyến tính. Xin thứ lỗi, nhưng tôi không theo đoạn thứ hai (điều đó giải thích cách được định nghĩa và được tính toán). Mọi sự trợ giúp sẽ rất được trân trọng.d Y / d PG|xdY/dP

Việc tính toán các dải tin cậy và dự đoán là khá chuẩn. Đọc để biết chi tiết về cách Prism tính toán dự đoán và các dải tin cậy của hồi quy phi tuyến.

Trước tiên, hãy xác định G | x, là độ dốc của các tham số tại một giá trị cụ thể của X và sử dụng tất cả các giá trị phù hợp nhất của các tham số. Kết quả là một vectơ, với một phần tử cho mỗi tham số. Đối với mỗi tham số, nó được định nghĩa là dY / dP, trong đó Y là giá trị Y của đường cong cho giá trị cụ thể của X và tất cả các giá trị tham số phù hợp nhất và P là một trong các tham số.)

G '| x là vectơ gradient được hoán vị, vì vậy nó là một cột chứ không phải là một hàng các giá trị.

Cov là ma trận hiệp phương sai (đảo ngược Hessian từ lần lặp cuối cùng). Nó là một ma trận vuông với số lượng hàng và cột bằng số lượng tham số. Mỗi mục trong ma trận là hiệp phương sai giữa hai tham số.

Bây giờ hãy tính c = G '| x * Cov * G | x. Kết quả là một số duy nhất cho bất kỳ giá trị nào của X.

Các dải tin cậy và dự đoán được tập trung vào đường cong phù hợp nhất và mở rộng trên và dưới đường cong một lượng bằng nhau.

Các dải tin cậy mở rộng trên và dưới đường cong bởi: = sqrt (c) * sqrt (SS / DF) * CriticalT (Confidence%, DF)

Các dải dự đoán mở rộng khoảng cách xa hơn trên và dưới đường cong, bằng: = sqrt (c + 1) * sqrt (SS / DF) * CriticalT (Confidence%, DF)


Hy vọng điều này sẽ giúp: stats.stackexchange.com/questions/74334/ từ
Bipi

Hy vọng điều này sẽ giúp: stats.stackexchange.com/questions/74334/ từ
Bipi

Đây thực sự được gọi là phương thức delta và sử dụng xấp xỉ Taylor bậc một. Mặc dù tốt hơn là sử dụng xấp xỉ Taylor bậc 2 cho việc này - hàm dự đoánNLS trong gói tuyên truyền sẽ làm điều đó nếu bạn quan tâm!
Tom Wenseleers

Câu trả lời:


18

Đây được gọi là Phương pháp Delta.

Giả sử rằng bạn có một số chức năng ; lưu ý rằng G ( ) là một hàm của các tham số mà bạn ước tính, β và các giá trị của các yếu tố dự đoán của bạn, x . Trước tiên, tìm đạo hàm của chức năng này đối với vector của bạn các thông số, với β : G ' ( β , x )y=G(β,x)+ϵG()βxβG(β,x). Điều này nói rằng, nếu bạn thay đổi một tham số một chút, chức năng của bạn thay đổi bao nhiêu? Lưu ý rằng đạo hàm này có thể là một hàm của chính các tham số của bạn cũng như các yếu tố dự đoán. Ví dụ, nếu , sau đó đạo hàm là x exp ( β x ) , mà phụ thuộc vào giá trị của β và giá trị của x . Để đánh giá này, bạn cắm trong dự toán của β rằng thủ tục của bạn mang đến, β , và giá trị của các yếu tố dự báo xG(β,x)=exp(βx)xexp(βx)βxββ^x nơi bạn muốn dự đoán.

Đồng bằng sông Cửu Phương pháp, xuất phát từ thủ tục maximum likelihood, tiểu bang rằng phương sai của sẽ là G ' ( β , x ) T Var ( β ) G ' ( β , x ) , nơi var ( β )G(β^,x)

G(β^,x)TVar(β^)G(β^,x),
Var(β^)là ma trận phương sai hiệp phương sai của các ước tính của bạn (giá trị này bằng với nghịch đảo của Hessian --- các đạo hàm thứ hai của hàm khả năng theo ước tính của bạn). Hàm mà các gói thống kê của bạn sử dụng sẽ tính toán giá trị này cho từng giá trị khác nhau của bộ dự đoán . Đây chỉ là một số, không phải là một vectơ, cho mỗi giá trị của x .xx

Điều này đưa ra phương sai của giá trị của hàm tại mỗi điểm và điều này được sử dụng giống như bất kỳ phương sai nào khác trong việc tính khoảng tin cậy: lấy căn bậc hai của giá trị này, nhân với giá trị tới hạn cho phân phối t thông thường hoặc áp dụng có liên quan cho a mức độ tin cậy cụ thể, và cộng và trừ giá trị này vào ước tính của tại điểm.G()

Đối với những khoảng thời gian dự đoán, chúng ta cần phải thực hiện sự thay đổi của kết quả đưa ra dự đoán , Var ( y | x ) σ 2 , vào tài khoản. Do đó, chúng ta phải đẩy mạnh sai của chúng tôi từ phương pháp Delta theo ước tính của chúng ta về phương sai của ε , σ 2 , để có được sự thay đổi của y , chứ không phải là sự thay đổi của giá trị kỳ vọng của y được sử dụng cho khoảng tin cậy. Lưu ý rằng σ 2 là tổng của bình phương lỗi ( trong tập tin ký hiệu giúp đỡ) chia cho bậc tự do ( ).xVar(yx)σ2ϵσ^2yyσ^2SSDF

cσ2σ2σc*SS/DF

c(xx)1Var(β^)=σ2(xx)1


Bạn có thể giải thích tính toán ci? Không giống như điểm quan trọng của t * sqrt (var)
B_Miner

Tôi nghĩ rằng tôi hiểu tính toán của họ; Tôi cập nhật phản hồi của tôi.
Charlie

Charlie, cảm ơn rất nhiều vì đã trả lời chi tiết. Tôi dự định viết mã để có thể tính toán dải dự đoán 95%. Tôi sẽ cho bạn biết làm thế nào đi.
Joe Listerr

@Charlie - rất rất tốt đẹp!
B_Miner

2
@Charlie. Cảm ơn. Tôi đã thêm một câu vào Câu hỏi thường gặp về GraphPad Prism giải thích rằng chúng tôi sử dụng cov để có nghĩa là ma trận hiệp phương sai chuẩn hóa (mỗi giá trị nằm trong khoảng từ -1 đến 1). Tôi cũng đã thêm một liên kết đến trang này, rất tốt cho bất cứ ai tìm kiếm chi tiết toán học.
Harvey Motulsky
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.