วันอังคารที่ 10 พฤศจิกายน พ.ศ. 2558

t-test ด้วย R (How to do t-test in R)

       การวิเคราะห์ปริมาณสารที่สนใจด้วย HPLC เมื่อเติมสารที่ความเข้มข้น 2,999 mg/kg จำนวน 5 ซ้ำ เป็นดั้งนี้ 3019, 2982, 3007, 3068, 3033 mg/kg ต้องการทราบว่าค่าเฉลี่ยจากการวิเคราะห์มีค่าแตกต่างจากค่าที่เติมลงไปหรือไม่ เราสามารถใช้วิธีการทดลองแบบ t-test โดยเปรียบเทียบค่าที่วัดได้กับค่าอ้างอิงคือค่าที่เติมลงไป (2999 mg/kg) ซึ่งใน R มีวิธีการดังนี้


กรอกข้อมูลในชื่อ analyzed

> analyzed <- c(3019, 2982, 3007, 3068, 3033)

ใช้ฟังก์ชั่น t.test()

> t.test(analyzed, mu = 2999)

        One Sample t-test

data:  analyzed
t = 1.5981, df = 4, p-value = 0.1853
alternative hypothesis: true mean is not equal to 2999
95 percent confidence interval:
 2982.189 3061.411
sample estimates:
mean of x
   3021.8


      จากค่า p-value = 0.1853 แสดงให้เห็นว่าค่าเฉลี่ยของข้อมูลไม่ต่างกับค่าอ้างอิง 2999

ไม่มีความคิดเห็น:

แสดงความคิดเห็น