การสำรวจข้อมูลเบื้องต้น (EDA) ของ NHANES
มาสำรวจชุดข้อมูลที่สร้างขึ้นใหม่โดยใช้แนวทาง EDA กัน เช่นเดียวกับบทที่แล้ว การดูทั้งตัวเลขสรุปและกราฟจะช่วยให้เข้าใจข้อมูลได้ดีขึ้น และยังเป็นวิธีที่ดีในการตรวจหาขั้นตอนการทำความสะอาดข้อมูลที่อาจพลาดไป ชุดข้อมูล nhanes_combined ถูกโหลดไว้ให้แล้ว
สมมติว่าเรามีสิทธิ์เข้าถึงข้อมูลผู้ป่วย NHANES และต้องการศึกษาผลของการที่แพทย์แนะนำให้ผู้ป่วยลดแคลอรี/ไขมันในอาหารต่อน้ำหนักตัว นี่คือ treatment ของเรา โดยสมมติว่าแทนที่จะเป็นคำถามที่ถามผู้ป่วย เราจะให้แพทย์สุ่มให้คำปรึกษาด้านโภชนาการแก่ผู้ป่วยบางราย อย่างไรก็ตาม เราสงสัยว่าอาจมีความแตกต่างของน้ำหนักตัวตามเพศของผู้ป่วย ซึ่งถือเป็น blocking factor!
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การออกแบบการทดลองใน R
คำแนะนำการฝึกหัด
- เติมโค้ด
dplyrแล้วรันเพื่อหาค่าเฉลี่ยน้ำหนัก (bmxwt) หน่วยกิโลกรัม จำแนกตาม treatment (mcq365d) มีอะไรน่าสังเกตเกี่ยวกับผู้ป่วยกลุ่มNAหรือไม่? - เติมโค้ด
ggplot2เพื่อสร้าง boxplot แสดง IQR ของน้ำหนักผู้ป่วย จำแนกตามตัวแปร treatment
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Fill in the dplyr code
___ %>%
group_by(___) %>%
summarize(mean = mean(___, na.rm = TRUE))
# Fill in the ggplot2 code
___ %>%
ggplot(aes(as.factor(___), ___)) +
geom_boxplot() +
labs(x = "Treatment",
y = "Weight")