BIBD - การทำงานของไตในแมว

เพื่อให้แน่ใจว่าเข้าใจ BIBD อย่างถ่องแท้ มาสร้างชุดข้อมูลที่มีโครงสร้าง BIBD ขึ้นมาตั้งแต่ต้นกัน

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

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

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

ดูคอร์ส

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

ฟังก์ชัน lambda() ที่กำหนดเองถูกโหลดไว้ให้แล้ว คำนวณค่า lambda โดยใช้ t = 4, k = 3 และ r = 3 เพื่อตรวจสอบว่า BIBD เป็นไปได้
รันโค้ดเพื่อประกอบชุดข้อมูล จะเห็นลำดับการใช้ treatment ของอาหารในแต่ละบล็อก
สร้าง cat_model ด้วย aov() ตามคำอธิบายการทดลองข้างต้น แล้วตรวจสอบผลลัพธ์ด้วย summary() ประเภทของอาหารเปียกมีผลต่อระดับครีเอตินินหรือไม่

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

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

# Calculate lambda
___

# Build the data.frame
creatinine <- c(1.98, 1.97, 2.35, 2.09, 1.87, 1.95, 2.08, 2.01, 1.84, 2.06, 1.97, 2.22)
food <- as.factor(c("A", "C", "D", "A", "B", "C", "B", "C", "D", "A", "B", "D"))
color <- as.factor(rep(c("Black", "White", "Orange", "Spotted"), each = 3))
cat_experiment <- as.data.frame(cbind(creatinine, food, color))

# Create cat_model and examine with summary()
___ <- ___ (___ ~ ___ + color, data = ___)
___

แก้ไขและรันโค้ด