เริ่มต้นใช้งานเริ่มต้นใช้งานได้ฟรี

RCBD กับข้อมูล NHANES

ทบทวนว่าการทดลองแบบบล็อกของเราใช้การรักษาที่แพทย์แนะนำให้ผู้ป่วยลดไขมันหรือแคลอรีในอาหาร และเราทดสอบผลที่มีต่อน้ำหนัก (bmxwt) โดยวางแผนจะบล็อกตามเพศ ซึ่งในชุดข้อมูล NHANES เก็บไว้ในตัวแปร riagendr การบล็อกมีจุดประสงค์เพื่อสร้างกลุ่มทดลองที่มีความคล้ายคลึงกันมากที่สุด การบล็อกการทดลองนี้ตามเพศหมายความว่า หากพบผลของการรักษาต่อ bmxwt ผลนั้นมีแนวโน้มเกิดจากการรักษาจริง ไม่ใช่จากเพศของผู้เข้าร่วม

ในโค้ด R การระบุการทดลองแบบบล็อกทำโดยใช้สูตรในรูปแบบ outcome ~ treatment + blocking_factor ในฟังก์ชันโมเดลที่เหมาะสม

nhanes_final พร้อมใช้งานแล้ว

แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร

การออกแบบการทดลองใน R

ดูคอร์ส

คำแนะนำการฝึกหัด

  • ใช้ aov() เพื่อสร้าง nhanes_rcbd โดยตัวแปรการรักษาคือ mcq365d ผลลัพธ์ที่ทดสอบคือ bmxwt และปัจจัยบล็อกคือ riagendr
  • ตรวจสอบผลลัพธ์ของ nhanes_rcbd ด้วย summary()
  • ใช้ฟังก์ชันจาก dplyr เพื่อดูค่าเฉลี่ยน้ำหนักแยกตาม mcq365d และ riagendr

แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ

ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์

# Use aov() to create nhanes_rcbd
___

# Check results of nhanes_rcbd with summary()
___

# Print mean weights by mcq365d and riagendr
___ %>% 
	group_by(___, ___) %>% 
	summarize(mean_wt = ___(___))
แก้ไขและรันโค้ด