ANOVA
จนถึงตอนนี้ คุณได้สร้างออบเจกต์เอาต์พุตเพื่อบันทึกและดึงองค์ประกอบต่าง ๆ สำหรับการนำเสนอผลลัพธ์ที่ปรับแต่งได้และละเอียดยิ่งขึ้น อย่างไรก็ตาม ออบเจกต์เอาต์พุตบางประเภท เช่น ออบเจกต์จากฟังก์ชัน aov() ที่ใช้รันการวิเคราะห์ความแปรปรวน (ANOVA) สามารถนำไปใช้เป็น input ของฟังก์ชันอื่น ๆ ได้ เช่น summary() และ plot() เพื่อรายงานผลและสร้างการแสดงผลแบบกราฟที่ปรับแต่งได้ นอกจากนี้ เอาต์พุตของ aov() ยังสามารถใช้เป็น input ให้กับ TukeyHSD() เพื่อทำการทดสอบเชิงคู่แบบ post hoc ได้อีกด้วย
ชุดข้อมูล abaloneKeep และแพ็กเกจ dplyr ถูกโหลดให้แล้ว
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
R สำหรับผู้ใช้ SAS
คำแนะนำการฝึกหัด
- คำนวณค่า
n, ค่าเฉลี่ย, ส่วนเบี่ยงเบนมาตรฐาน และความแปรปรวนของ age จำแนกตาม sex - รัน ANOVA ของ
ageตามกลุ่มsexแล้วบันทึกเอาต์พุตเป็นabaov - รัน
summary()ของabaov - รันการทดสอบเชิงคู่แบบ post hoc ด้วย
TukeyHSD()บนabaov
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Compute n(), mean(), sd() and var() of age by sex
abaloneKeep %>% group_by(___) %>% select(___, ___) %>%
summarise(across(everything(), list(mean = ___,
___ = ~sd(.x),
var = ___)),
N = ___)
# Run aov() of age by sex, save as abaov
abaov <- aov(___ ~ ___, data = ___)
# Run summary() of abaov
___
# Perform TukeyHSD posthoc pairwise tests on abaov
TukeyHSD(___)