Tôi có một tập hợp các câu trả lời được thể hiện dưới dạng một khoảng như mẫu dưới đây.
> head(left)
[1] 860 516 430 1118 860 602
> head(right)
[1] 946 602 516 1204 946 688
trong đó bên trái là giới hạn dưới và bên phải là giới hạn trên của phản hồi. Tôi muốn ước tính các tham số theo phân phối logic.
Trong một thời gian khi tôi đang cố gắng tính toán khả năng trực tiếp, tôi đã phải vật lộn với thực tế là do hai giới hạn được phân phối dọc theo các bộ thông số khác nhau, tôi đã nhận được một số giá trị âm như dưới đây:
> Pr_high=plnorm(wta_high,meanlog_high,sdlog_high)
> Pr_low=plnorm(wta_low, meanlog_low,sdlog_low)
> Pr=Pr_high-Pr_low
>
> head(Pr)
[1] -0.0079951419 0.0001207749 0.0008002343 -0.0009705125 -0.0079951419 -0.0022395514
Tôi thực sự không thể tìm ra cách giải quyết nó và quyết định sử dụng điểm giữa của khoảng thay vào đó là một sự thỏa hiệp tốt cho đến khi tôi tìm thấy hàm mledist trích xuất khả năng loglik của một phản hồi khoảng, đây là tóm tắt tôi nhận được:
> mledist(int, distr="lnorm")
$estimate
meanlog sdlog
6.9092257 0.3120138
$convergence
[1] 0
$loglik
[1] -152.1236
$hessian
meanlog sdlog
meanlog 570.760358 7.183723
sdlog 7.183723 1112.098031
$optim.function
[1] "optim"
$fix.arg
NULL
Warning messages:
1: In plnorm(q = c(946L, 602L, 516L, 1204L, 946L, 688L, 1376L, 1376L, :
NaNs produced
2: In plnorm(q = c(860L, 516L, 430L, 1118L, 860L, 602L, 1290L, 1290L, :
NaNs produced
Các giá trị tham số dường như có ý nghĩa và loglikabilities lớn hơn bất kỳ phương pháp nào khác mà tôi đã sử dụng (phân phối trung điểm hoặc phân phối của một trong hai giới hạn).
Có một thông điệp cảnh báo mà tôi không hiểu vì vậy có ai có thể cho tôi biết nếu tôi đang làm đúng và thông điệp này có ý nghĩa gì không?
Đánh giá cao sự giúp đỡ!
fitdistrplus
.