Tôi có một tập dữ liệu mẫu với 31 giá trị. Tôi đã chạy thử nghiệm t hai đầu bằng R để kiểm tra nếu giá trị trung bình thực bằng 10:
t.test(x=data, mu=10, conf.level=0.95)
Đầu ra:
t = 11.244, df = 30, p-value = 2.786e-12
alternative hypothesis: true mean is not equal to 10
95 percent confidence interval:
19.18980 23.26907
sample estimates:
mean of x
21.22944
Bây giờ tôi đang cố gắng làm điều tương tự bằng tay:
t.value = (mean(data) - 10) / (sd(data) / sqrt(length(data)))
p.value = dt(t.value, df=length(lengths-1))
Giá trị t được tính bằng phương pháp này giống như đầu ra của hàm t-test R. Giá trị p, tuy nhiên, đi ra là 3.025804e-12.
Có ai biết tôi đang làm gì sai không?
Cảm ơn!
BIÊN TẬP
Đây là mã R đầy đủ, bao gồm cả tập dữ liệu của tôi:
# Raw dataset -- 32 observations
data = c(21.75, 18.0875, 18.75, 23.5, 14.125, 16.75, 11.125, 11.125, 14.875, 15.5, 20.875,
17.125, 19.075, 25.125, 27.75, 29.825, 17.825, 28.375, 22.625, 28.75, 27, 12.825,
26, 32.825, 25.375, 24.825, 25.825, 15.625, 26.825, 24.625, 26.625, 19.625)
# Student t-Test
t.test(x=data, mu=10, conf.level=0.95)
# Manually calculate p-value
t.value = (mean(data) - 10) / (sd(data) / sqrt(length(data)))
p.value = dt(t.value, df=length(data) - 1)