วิธีสร้างพล็อต qq ใน sas

พล็อต QQ ย่อมาจาก “ควอนไทล์-ควอนไทล์” ใช้เพื่อประเมินว่าชุดข้อมูลอาจมาจากการแจกแจงทางทฤษฎีหรือไม่

ในกรณีส่วนใหญ่ การลงจุดประเภทนี้ใช้เพื่อพิจารณาว่าชุดข้อมูลเป็นไปตามการแจกแจงแบบปกติหรือไม่

หากข้อมูลมีการกระจายตามปกติ จุดบนพล็อต QQ จะอยู่บนเส้นทแยงมุม

ในทางกลับกัน ยิ่งจุดบนกราฟเบี่ยงเบนไปจากเส้นทแยงมุมที่มีนัยสำคัญมากเท่าใด ชุดข้อมูลก็จะเป็นไปตามการแจกแจงแบบปกติก็จะน้อยลงเท่านั้น

วิธีที่ง่ายที่สุดในการสร้างพล็อต QQ ใน SAS คือการใช้คำสั่ง PROC UNIVARIATE กับคำสั่ง QQPLOT :

 proc univariate data =my_data noprint ;
   qqplot my_variable;
run ;

ตัวอย่างต่อไปนี้แสดงวิธีใช้ไวยากรณ์นี้ในทางปฏิบัติ

หมายเหตุ : เราใช้คำสั่ง NOPRINT เพื่อลบสถิติสรุปและตารางอื่นๆ ทั้งหมดที่สร้างขึ้นโดยอัตโนมัติโดยคำสั่ง PROC UNIVARIATE

ตัวอย่างที่ 1: สร้างพล็อต QQ ใน SAS สำหรับข้อมูลปกติ

รหัสต่อไปนี้แสดงวิธีสร้างพล็อต QQ สำหรับชุดข้อมูลที่มี การสังเกต 1,000 ครั้งที่สร้างขึ้นจาก การแจกแจงแบบปกติ โดยมีค่าเฉลี่ย 10 และค่าเบี่ยงเบนมาตรฐาน 2:

 /*generate 1000 values that follow normal distribution with mean 10 and sd 2 */
data normal_data;
    do i = 1 to 1000;
	x = 10 + 2* rannor (1);
        output ;
    end ;
run ;

/*create some plot*/
proc univariate data =normal_data noprint ;
   qqplot x;
run ; 

เราจะเห็นได้ว่าจุดส่วนใหญ่อยู่ตามแนวเส้นทแยงมุมตรงโดยมีการเบี่ยงเบนเล็กน้อยตามหางแต่ละข้าง

จากกราฟนี้ เราสามารถสรุปได้อย่างปลอดภัยว่าชุดข้อมูลนี้มีการกระจายตามปกติ

ตัวอย่างที่ 2: พล็อต QQ สำหรับข้อมูลที่ไม่ปกติ

รหัสต่อไปนี้แสดงวิธีสร้างพล็อต QQ สำหรับชุดข้อมูลที่มีการสังเกต 1,000 รายการที่สร้างจาก การแจกแจงแบบเอ็กซ์โพเนนเชียล :

 /*generate 1000 values that follow an exponential distribution*/
data exp_data;
    do i = 1 to 1000;
	x = ranexp (1);
        output ;
    end ;
run ;

/*create some plot*/
proc univariate data =exp_data noprint ;
   qqplot x;
run ; 

ฉันเห็นว่าจุดเบี่ยงเบนไปจากเส้นทแยงมุมอย่างมาก สิ่งนี้แสดงให้เห็นอย่างชัดเจนว่าชุดข้อมูลไม่ได้กระจายตามปกติ

สิ่งนี้น่าจะสมเหตุสมผลเนื่องจากเราระบุว่าข้อมูลควรเป็นไปตามการแจกแจงแบบเอ็กซ์โพเนนเชียล

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

บทช่วยสอนต่อไปนี้จะอธิบายวิธีดำเนินการงานทั่วไปอื่นๆ ใน SAS:

วิธีใช้ Proc Univariate สำหรับการทดสอบความเป็นมาตรฐานใน SAS
วิธีดำเนินการแปลงบันทึกใน SAS
วิธีระบุค่าผิดปกติใน SAS