Làm cách nào tôi có thể ước tính khoảng tin cậy 95% bằng cách sử dụng cấu hình cho các tham số được ước tính bằng cách tối đa hóa hàm khả năng ghi nhật ký bằng cách sử dụng tối ưu trong R?

Làm cách nào tôi có thể ước tính khoảng tin cậy 95% bằng cách sử dụng cấu hình cho các tham số được ước tính bằng cách tối đa hóa hàm khả năng ghi nhật ký bằng cách sử dụng tối ưu trong R?

Tôi biết rằng tôi có thể ước lượng bất thường ma trận hiệp phương sai bằng cách đảo ngược hessian , nhưng tôi lo ngại rằng dữ liệu của tôi không đáp ứng các giả định cần thiết cho phương pháp này là hợp lệ. Tôi muốn ước tính khoảng tin cậy bằng một số phương pháp khác.

Là hồ sơ khả năng phương pháp thích hợp, như đã thảo luận trong Stryhn và Christensen , và trong cuốn sách MASS Venables và Ripley, §8.4, pp. 220-221?

Nếu vậy, có gói nào có thể giúp tôi làm điều này trong R không? Nếu không, mã giả cho một phương thức như vậy sẽ như thế nào?

Các mlechức năng từ gói stats4 là một wrapper của optim, mà làm cho nó khá dễ dàng để tạo ra tính toán hồ sơ khả năng. Xem help("profile,mle-method", package = "stats4")để biết thêm thông tin.

Có gói ProfileLikabilities nếu bạn sử dụng nlme. Cá nhân, tôi đã không quản lý để sử dụng nó.

Sử dụng gói lme4a hoặc lmeEigen có một hàm profile () nhằm mục đích chính xác là làm những gì bạn muốn. Hãy thử một cái gì đó như thế để cài đặt các gói này:

install.packages("lme4a",repos="http://lme4.r-forge.r-project.org/repos") 

hoặc vào trang web để lấy kho lưu trữ zip. Tương tự và thật không may, tôi đã không quản lý để sử dụng nó :) Có lẽ chúng ta nên chờ đợi một bản cập nhật của lme4.

Phương pháp này được trình bày chi tiết trong bản thảo cuốn sách của Douglas Bates

EDIT: Tuyệt! Hàm profile () cho các mô hình lmer hiện có sẵn trong phiên bản mới nhất của lme4, được cài đặt bằng cách gõ:

install.packages("lme4",repos="http://r-forge.r-project.org")