Hãy xem xét các mã và đầu ra sau đây:
par(mfrow=c(3,2))
# generate random data from weibull distribution
x = rweibull(20, 8, 2)
# Quantile-Quantile Plot for different distributions
qqPlot(x, "log-normal")
qqPlot(x, "normal")
qqPlot(x, "exponential", DB = TRUE)
qqPlot(x, "cauchy")
qqPlot(x, "weibull")
qqPlot(x, "logistic")
Có vẻ như âm mưu QQ cho log-normal gần giống như âm mưu QQ cho weibull. Làm thế nào chúng ta có thể phân biệt chúng? Ngoài ra nếu các điểm nằm trong vùng được xác định bởi hai đường màu đen bên ngoài, điều đó có cho biết rằng chúng tuân theo phân phối đã chỉ định không?
library(car)
trong mã của mình để mọi người dễ theo dõi hơn. Nói chung, bạn cũng có thể muốn đặt hạt giống (ví dụset.seed(1)
:) để làm cho ví dụ có thể lặp lại, để bất kỳ ai cũng có thể nhận được chính xác các điểm dữ liệu bạn đã nhận được, mặc dù điều này có thể không quan trọng ở đây.