Phát hiện ngoại lệ đa biến
100 người sống trong cùng một khu vực đã nộp yêu cầu bồi thường vì nhà của họ bị hư hại do mưa đá vào đêm Chủ nhật. Bộ dữ liệu hailinsurance có 100 quan sát và 2 biến. Cột thứ nhất là số tiền chi trả mà công ty bảo hiểm đã thanh toán cho mỗi khách hàng, còn cột thứ hai là giá nhà gần nhất.
Trong bài tập này, trước tiên bạn sẽ dùng các ước lượng cổ điển trên bộ dữ liệu. Sau đó, bạn sẽ so sánh kết quả với các ước lượng vững (robust).
Bài tập này là một phần của khóa học
Phát hiện gian lận với R
Bài tập tương tác thực hành trực tiếp
Hãy thử làm bài tập này bằng cách hoàn thành đoạn mã mẫu này.
# Create a scatterplot
plot(hailinsurance, xlab = "price house", ylab = "claim")
# Compute the sample mean and sample covariance matrix
clcenter <- colMeans(___)
clcov <- cov(___)
# Add 97.5% tolerance ellipsoid
rad <- sqrt(qchisq(___, ___))
ellipse(center = clcenter, shape = clcov, radius = rad,col = "blue", lty = 2)