วิธีการใช้การแจกแจงแกมมาใน r (พร้อมตัวอย่าง)

ในสถิติ การแจกแจงแกมม่า มักใช้ในการสร้างแบบจำลองความน่าจะเป็นที่เกี่ยวข้องกับเวลาที่รอคอย

เราสามารถใช้ฟังก์ชันต่อไปนี้เพื่อทำงานกับการแจกแจงแกมมาใน R:

• dgamma(x, รูปร่าง, อัตรา) – ค้นหาค่าของฟังก์ชันความหนาแน่นของการแจกแจงแกมมาด้วยพารามิเตอร์รูปร่างและอัตราที่แน่นอน
• pgamma(q, รูปร่าง, อัตรา) – ค้นหาค่าของฟังก์ชันความหนาแน่นสะสมของการแจกแจงแกมมาด้วยพารามิเตอร์รูปร่างและอัตราที่แน่นอน
• qgamma(p, รูปร่าง, อัตรา) – ค้นหาค่าของฟังก์ชันความหนาแน่นสะสมผกผันของการแจกแจงแกมมาด้วยพารามิเตอร์รูปร่างและอัตราที่แน่นอน
• rgamma(n, รูปร่าง, อัตรา) – สร้างตัวแปรสุ่ม n ตัวที่ตามหลังการแจกแจงแกมมาด้วยพารามิเตอร์รูปร่างและอัตราที่แน่นอน

ตัวอย่างต่อไปนี้แสดงวิธีใช้แต่ละฟังก์ชันเหล่านี้ในทางปฏิบัติ

ตัวอย่างที่ 1: วิธีใช้ dgamma()

รหัสต่อไปนี้แสดงวิธีใช้ฟังก์ชัน dgamma() เพื่อสร้างพล็อตความหนาแน่นของความน่าจะเป็นของการแจกแจงแกมมาด้วยพารามิเตอร์บางอย่าง:

 #define x-values
x <- seq(0, 2, by=0.01)   
  
#calculate gamma density for each x-value
y <- dgamma(x, shape=5) 
  
#create density plot
plot(y)

ตัวอย่างที่ 2: วิธีใช้ pgamma()

รหัสต่อไปนี้แสดงวิธีใช้ฟังก์ชัน pgamma() เพื่อสร้างพล็อตความหนาแน่นสะสมของการแจกแจงแกมมาด้วยพารามิเตอร์บางอย่าง:

 #define x-values
x <- seq(0, 2, by=0.01)   
  
#calculate gamma density for each x-value
y <- pgamma(x, shape=5) 
  
#create cumulative density plot
plot(y) 

ตัวอย่างที่ 3: วิธีใช้ qgamma()

รหัสต่อไปนี้แสดงวิธีใช้ฟังก์ชัน qgamma() เพื่อสร้างพล็อตควอนไทล์ของการแจกแจงแกมมาด้วยพารามิเตอร์บางตัว:

 #define x-values
x <- seq(0, 1, by=0.01)   
  
#calculate gamma density for each x-value
y <- qgamma(x, shape=5) 
  
#create quantile plot
plot(y) 

ตัวอย่างที่ 4: วิธีใช้ rgamma()

รหัสต่อไปนี้แสดงวิธีใช้ฟังก์ชัน rgamma() เพื่อสร้างและแสดงภาพตัวแปรสุ่ม 1,000 ตัวที่เป็นไปตามการแจกแจงแกมมาด้วยพารามิเตอร์รูปร่างเป็น 5 และพารามิเตอร์อัตราเป็น 3:

 #make this example reproducible
set. seeds (0)

#generate 1,000 random values that follow gamma distribution
x <- rgamma(n=1000, shape=5, rate=3)

#create histogram to view distribution of values
hist(x)

แหล่งข้อมูลเพิ่มเติม

บทช่วยสอนต่อไปนี้จะอธิบายวิธีใช้การแจกแจงทางสถิติทั่วไปอื่นๆ ใน R:

วิธีใช้การแจกแจงแบบปกติใน R
วิธีใช้การแจกแจงแบบทวินามใน R
วิธีใช้การแจกแจงปัวซองใน R
วิธีใช้การกระจายทางเรขาคณิตใน R