วิธีการพล็อตการแจกแจงแบบ weibull ใน r

ในการพล็อตฟังก์ชันความหนาแน่นของความน่าจะเป็นสำหรับการแจกแจงแบบไวบูลใน R เราสามารถใช้ฟังก์ชันต่อไปนี้:

• dweibull(x, รูปร่าง, สเกล = 1) เพื่อสร้างฟังก์ชันความหนาแน่นของความน่าจะเป็น
• curve(function, from = NULL, to = NULL) เพื่อพล็อตฟังก์ชันความหนาแน่นของความน่าจะเป็น

ในการพล็อตฟังก์ชันความหนาแน่นของความน่าจะเป็น เราจำเป็นต้องระบุค่าของพารามิเตอร์ รูปร่าง และ มาตราส่วน ในฟังก์ชัน dweibull รวมถึงค่า จาก และ ถึง ในฟังก์ชัน curve()

ตัวอย่างเช่น โค้ดต่อไปนี้แสดงวิธีการพล็อตฟังก์ชันความหนาแน่นของความน่าจะเป็นสำหรับการแจกแจงแบบ Weibull ด้วยพารามิเตอร์รูปร่าง = 2 และสเกล = 1 โดยที่แกน x ของพล็อตเปลี่ยนจาก 0 ถึง 4:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4)

เราสามารถเพิ่มชื่อเรื่อง เปลี่ยนป้ายกำกับแกน Y เพิ่มความกว้างของแถว และแม้แต่เปลี่ยนสีของแถวเพื่อทำให้โครงเรื่องดูสวยงามยิ่งขึ้น:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, 
    main = 'Weibull Distribution (shape = 2, scale = 1)', #add title
    ylab = 'Density', #change y-axis label
    lwd = 2, #increase line width to 2
    col = 'steelblue') #change line color to steelblue 

นอกจากนี้เรายังสามารถเพิ่มเส้นโค้งหลายๆ เส้นลงในกราฟเพื่อเปรียบเทียบการแจกแจงแบบ Weibull กับพารามิเตอร์รูปร่างและขนาดที่แตกต่างกัน:

 curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

เราสามารถเพิ่มคำอธิบายแผนภูมิให้กับโครงเรื่องได้โดยใช้ฟังก์ชัน legend() ซึ่งใช้ไวยากรณ์ต่อไปนี้:

ตำนาน (x, y=NULL, ตำนาน, เติม, col, bg, lty, cex)

ทอง:

• x, y: พิกัด x และ y ที่ใช้ในการวางตำแหน่งคำอธิบาย
• ตำนาน: ข้อความที่จะใส่ในตำนาน
• เติม: เติมสีในตำนาน
• col: รายการสีที่จะใช้สำหรับเส้นภายในคำอธิบายแผนภูมิ
• bg: สีพื้นหลังของตำนาน
• lty: สไตล์เส้น
• cex: ขนาดของข้อความในตำนาน

ในตัวอย่างของเรา เราจะใช้ไวยากรณ์ต่อไปนี้เพื่อสร้างคำอธิบาย:

 #create density plots
curve(dweibull(x, shape=2, scale = 1), from=0, to=4, col='red')
curve(dweibull(x, shape=1.5, scale = 1), from=0, to=4, col='blue', add=TRUE)

#add legend
legend(2, .7, legend=c("shape=2, scale=1", "shape=1.5, scale=1"),
       col=c("red", "blue"), lty=1, cex=1.2)