การใช้งานผลลัพธ์จาก R (2)
เมื่อพิจารณาคำถามวิจัยทางวิทยาศาสตร์ หาก IQ เกิดจากพันธุกรรมเพียงอย่างเดียว เราจะคาดว่าความชันของเส้นระหว่างแฝดทั้งสองกลุ่มจะเท่ากับ 1 การทดสอบสมมติฐานที่กำหนดค่าความชันเป็น 1 ทำได้โดยสร้างค่าสถิติทดสอบใหม่ที่วัดว่าความชันที่สังเกตได้ห่างจากค่าสมมติฐาน 1 มากน้อยเพียงใด
$$new_t = \frac{slope - 1}{SE}$$
หากสมมติฐานที่ว่าความชันเท่ากับ 1 เป็นจริง ค่าสถิติทดสอบใหม่นี้จะมีการแจกแจงแบบ t ซึ่งสามารถนำไปใช้คำนวณ p-value ได้
ข้อมูลจากโมเดลในส่วนของชีววิทยาพร้อมใช้งานในชื่อ biological_term
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การอนุมานสำหรับการถดถอยเชิงเส้นใน R
คำแนะนำการฝึกหัด
- คำนวณดีกรีแห่งความอิสระของชุดข้อมูล
twins - คำนวณ p-value แบบสองด้านสำหรับสมมติฐานทางเลือกที่ว่าความชันที่แท้จริงแตกต่างจาก 1 โดยสร้างการคำนวณเป็นขั้นตอน ดังนี้
- คำนวณค่าสถิติทดสอบโดยนำ
estimateของความชันลบด้วย1แล้วหารด้วยค่าความคลาดเคลื่อนมาตรฐาน - คำนวณ p-value แบบด้านเดียวของค่าสถิติทดสอบโดยใช้ฟังก์ชันการแจกแจงสะสมของการแจกแจงแบบ t ด้วย
pt()ที่test_statisticพร้อมระบุดีกรีแห่งความอิสระที่คำนวณไว้ - คำนวณ p-value แบบสองด้านโดยนำ p-value แบบด้านเดียวคูณด้วยสอง
- คำนวณค่าสถิติทดสอบโดยนำ
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Calculate the degrees of freedom of twins
degrees_of_freedom <- nrow(twins) - 2
biological_term %>%
mutate(
# Calculate the test statistic
test_statistic = ___,
# Calculate its one-sided p-value
one_sided_p_value_of_test_statistic = ___,
# ... and its two-sided p-value
two_sided_p_value_of_test_statistic = ___
)