พิจารณารหัสและผลลัพธ์ต่อไปนี้:
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")
ดูเหมือนว่าพล็อต QQ สำหรับบันทึกปกติเกือบจะเหมือนกับพล็อต QQ สำหรับ weibull เราจะแยกแยะพวกมันได้อย่างไร นอกจากนี้หากคะแนนอยู่ในพื้นที่ที่กำหนดโดยเส้นสีดำด้านนอกสองเส้นนั่นแสดงว่ามันเป็นไปตามการแจกแจงที่ระบุหรือไม่?
library(car)
ในรหัสของคุณเพื่อให้ผู้คนติดตามได้ง่ายขึ้น โดยทั่วไปคุณอาจต้องการตั้งค่าเมล็ดพันธุ์ (เช่นset.seed(1)
) เพื่อทำซ้ำตัวอย่างเพื่อให้ทุกคนสามารถได้รับจุดข้อมูลเดียวกันกับที่คุณได้รับแม้ว่ามันอาจจะไม่สำคัญเท่านี้