MathLabs

Xác suất và Thống kê

Ước lượng

Dùng dữ liệu mẫu để suy luận các giá trị chưa biết của tổng thể, bằng ước lượng điểm và khoảng.

Trực giácTừ một nắm nhỏ suy ra cả tổng thể

Hãy tưởng tượng nếm một thìa nhỏ từ một nồi súp lớn để đoán độ mặn của cả nồi. Bạn không cần uống hết nồi súp -- chỉ một thìa đại diện (một mẫu) là đủ để đoán mức độ mặn của cả nồi (tổng thể). Ước lượng là toán học biến thìa súp đó thành một dự đoán đáng tin cậy, cùng với cảm nhận về việc dự đoán đó có thể sai lệch bao nhiêu.

Đường cong hình chuông minh họa phân phối mẫu thu hẹp quanh giá trị trung bình khi cỡ mẫu tăng.
Phân phối mẫu chuẩn hóa N(0,1)\mathcal{N}(0,1) với khoảng tin cậy trung tâm [−z,+z][-z, +z] tô màu xanh lá; kéo pp để mở rộng hoặc thu hẹp mức độ tin cậy.

Phổ thôngƯớc lượng điểm: một dự đoán tốt nhất duy nhất

Định nghĩa: Ước lượng điểm

Ước lượng điểm là một con số duy nhất, tính từ dữ liệu mẫu, dùng làm dự đoán tốt nhất cho một tham số chưa biết của tổng thể. Để ước lượng trung bình tổng thể μ\mu, ước lượng điểm tự nhiên là trung bình mẫu Xˉ\bar X: trung bình của các giá trị quan sát được X1,…,XnX_1,\dots,X_n.

Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i

Ở đây nn là cỡ mẫu, X1,…,XnX_1,\dots,X_n là các giá trị quan sát riêng lẻ của mẫu, và Xˉ\bar X là trung bình của chúng. Vì các mẫu khác nhau cho ra các giá trị Xˉ\bar X khác nhau, nên bản thân nó là một biến ngẫu nhiên có phân phối riêng -- gọi là phân phối mẫu của trung bình -- có độ phân tán thu hẹp dần khi nn tăng.

Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt{n}} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt{n}}

Đây là khoảng tin cậy 100(1−α)%100(1-\alpha)\% cho μ\mu khi độ lệch chuẩn tổng thể σ\sigma đã biết: zα/2z_{\alpha/2} là giá trị tới hạn chuẩn tắc để lại xác suất α/2\alpha/2 ở mỗi đuôi, và σn\frac{\sigma}{\sqrt n} là sai số chuẩn của Xˉ\bar X, khoảng cách điển hình giữa Xˉ\bar X và μ\mu.

Các mức độ tin cậy thường dùng và giá trị tới hạn zz tương ứng
Mức độ tin cậyzα/2z_{\alpha/2}
90%zα/2=1.645z_{\alpha/2}=1.645
95%zα/2=1.960z_{\alpha/2}=1.960
99%zα/2=2.576z_{\alpha/2}=2.576

Đại họcVì sao khoảng tin cậy hoạt động đúng

Nếu X1,…,XnX_1,\dots,X_n độc lập, cùng phân phối với E[Xi]=μE[X_i] = \mu với mọi ii, thì trung bình mẫu Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i thỏa mãn E[Xˉ]=μE[\bar X] = \mu.

Vì sao đúng?

Việc lấy trung bình loại bỏ sai số hệ thống: tính trung bình qua nhiều mẫu giả định, Xˉ\bar X rơi đúng vào trung bình tổng thể thực μ\mu, không lệch cao hay lệch thấp một cách có hệ thống.

Chứng minh

Theo định nghĩa, Xˉ=1n∑i=1nXi\bar X = \frac{1}{n}\sum_{i=1}^{n} X_i. Kỳ vọng là một toán tử tuyến tính nên phân phối được qua tổng và hệ số hằng 1n\frac{1}{n}: E[Xˉ]=1n∑i=1nE[Xi]E[\bar X] = \frac{1}{n}\sum_{i=1}^{n} E[X_i].

Vì mỗi XiX_i đều đến từ cùng một tổng thể, nên E[Xi]=μE[X_i] = \mu với mỗi trong số nn số hạng, nên tổng rút gọn thành 1n∑i=1nμ=1n(nμ)=μ\frac{1}{n}\sum_{i=1}^{n} \mu = \frac{1}{n}(n\mu) = \mu.

Do đó E[Xˉ]=μE[\bar X] = \mu với mọi cỡ mẫu nn: trung bình mẫu không chệch dù mẫu nhỏ hay lớn. Đây chỉ là phát biểu về tâm của phân phối mẫu -- độ phân tán của phân phối đó, Var⁡(Xˉ)=σ2/n\operatorname{Var}(\bar X) = \sigma^2/n, vẫn thu hẹp dần khi nn tăng, và đó chính là điều khoảng tin cậy bên dưới sử dụng.

Nếu X1,…,XnX_1,\dots,X_n độc lập, cùng phân phối với trung bình μ\mu và phương sai hữu hạn σ2\sigma^2, thì với nn lớn, khoảng tin cậy Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt n} chứa μ\mu với xác suất xấp xỉ 1−α1-\alpha.

Vì sao đúng?

Định lý giới hạn trung tâm nói rằng trung bình mẫu đã chuẩn hóa xử sự như một biến chuẩn tắc khi nn lớn, nên ta có thể dùng các phân vị chuẩn zα/2z_{\alpha/2} để kẹp μ\mu với một tỉ lệ thành công đã biết về lâu dài, bất kể hình dạng của tổng thể gốc.

Chứng minh

Theo định lý giới hạn trung tâm, Xˉ−μσ/n\frac{\bar X - \mu}{\sigma/\sqrt n} hội tụ theo phân phối về N(0,1)N(0,1) khi nn→∞\to\infty. Vậy với nn lớn, −zα/2≤Xˉ−μσ/n≤zα/2-z_{\alpha/2} \le \frac{\bar X - \mu}{\sigma/\sqrt n} \le z_{\alpha/2} đúng với xác suất xấp xỉ 1−α1-\alpha, trong đó zα/2z_{\alpha/2} là giá trị để lại xác suất α/2\alpha/2 ở mỗi đuôi của N(0,1)N(0,1).

Nhân cả ba vế của bất đẳng thức với σn\frac{\sigma}{\sqrt n}>0>0 (giữ nguyên chiều bất đẳng thức): −zα/2σn≤Xˉ−μ≤zα/2σn-z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \bar X - \mu \le z_{\alpha/2}\frac{\sigma}{\sqrt n}.

Trừ Xˉ\bar X ở mọi vế, rồi nhân với −1-1 (đảo chiều bất đẳng thức): Xˉ−zα/2σn≤μ≤Xˉ+zα/2σn\bar X - z_{\alpha/2}\frac{\sigma}{\sqrt n} \le \mu \le \bar X + z_{\alpha/2}\frac{\sigma}{\sqrt n}. Đây chính xác là công thức khoảng tin cậy: vì phát biểu xác suất ban đầu đúng với xác suất xấp xỉ 1−α1-\alpha, nên khoảng đã sắp xếp lại này cũng chứa μ\mu với xác suất đó.

Sai số biên E=zα/2σnE = z_{\alpha/2}\frac{\sigma}{\sqrt n} chỉ thu hẹp theo 1/n1/\sqrt n, không phải theo 1/n1/n: để giảm sai số biên xuống một nửa, bạn phải tăng gấp bốn cỡ mẫu nn, và để giảm còn một phần ba thì cần gấp chín lần số quan sát. Độ chính xác rất tốn kém: mỗi chữ số chính xác thêm vào tốn nhiều dữ liệu hơn hẳn chữ số trước đó.

Đại họcỨng dụng thực tiễn và Ví dụ minh họa

Khoảng tin cậy xuất hiện ở bất cứ đâu cần ra quyết định từ một mẫu không đầy đủ: kiểm soát chất lượng trên dây chuyền nhà máy, thăm dò ý kiến trước bầu cử, thử nghiệm liều lượng trong y học, hay hiệu chỉnh kính thiên văn từ các phép đo lặp lại. Trong mọi trường hợp, cùng một công thức biến trung bình mẫu thành một khoảng thành thật về độ bất định của chính nó.

Ví dụ: Khoảng tin cậy cho đường kính trung bình của bu lông

Một kỹ sư kiểm soát chất lượng lấy mẫu n=64n=64 bu lông và tìm được đường kính trung bình mẫu xˉ=12.02\bar x = 12.02 mm, với độ lệch chuẩn tổng thể đã biết σ=0.16\sigma = 0.16 mm. Hãy xây dựng khoảng tin cậy 95% cho đường kính trung bình thực μ\mu, dùng z0.025=1.96z_{0.025}=1.96.

Lời giải

Sai số chuẩn là 0.1664=0.02\frac{0.16}{\sqrt{64}} = 0.02 mm.

Sai số biên là 1.96×0.02=0.03921.96 \times 0.02 = 0.0392 mm.

Khoảng tin cậy 95% là (12.02−0.0392, 12.02+0.0392)=(11.981, 12.059)(12.02 - 0.0392,\ 12.02 + 0.0392) = (11.981,\ 12.059) mm: ta tin cậy 95% rằng đường kính bu lông trung bình thực nằm trong khoảng này.

Ví dụ: Ước lượng tỉ lệ cử tri ủng hộ một chính sách

Một nhà thăm dò khảo sát n=400n=400 cử tri được chọn ngẫu nhiên và thấy tỉ lệ mẫu p^=0.53\hat p = 0.53 ủng hộ một chính sách. Coi tỉ lệ mẫu như trung bình mẫu của các phản hồi 0/1, hãy xây dựng khoảng tin cậy 95% cho tỉ lệ tổng thể thực μ\mu ủng hộ, dùng z0.025=1.96z_{0.025}=1.96.

Lời giải

Độ lệch chuẩn ước lượng của một phản hồi 0/1 đơn lẻ là p^(1−p^)=0.53×0.47≈0.499\sqrt{\hat p(1-\hat p)} = \sqrt{0.53\times 0.47} \approx 0.499, nên sai số chuẩn của tỉ lệ mẫu là 0.499400≈0.025\frac{0.499}{\sqrt{400}} \approx 0.025.

Sai số biên là 1.96×0.025≈0.0491.96 \times 0.025 \approx 0.049.

Khoảng tin cậy 95% là (0.53−0.049, 0.53+0.049)=(0.481, 0.579)(0.53 - 0.049,\ 0.53 + 0.049) = (0.481,\ 0.579): cuộc thăm dò tin cậy 95% rằng tỉ lệ ủng hộ thực nằm trong khoảng 48,1% đến 57,9%, quá rộng để khẳng định chắc chắn chính sách được đa số ủng hộ -- cần mẫu lớn hơn để có kết luận chặt hơn.

Dùng công thức sai số biên E=zα/2σnE = z_{\alpha/2}\frac{\sigma}{\sqrt n}, sai số biên là bao nhiêu khi zα/2=1.96z_{\alpha/2}=1.96, σ=8\sigma=8, và n=64n=64?

Một khoảng tin cậy 95% cho nhiệt độ trung bình ngày được tính từ một mẫu. Cách hiểu đúng của "tin cậy 95%" là gì?

Tính chất E[Xˉ]=μE[\bar X] = \mu với mọi cỡ mẫu nn mô tả đặc điểm nào của ước lượng Xˉ\bar X?

Nếu cỡ mẫu nn tăng từ 100 lên 400, với σ\sigma và zα/2z_{\alpha/2} không đổi, sai số biên được nhân với ...

Tài liệu tham khảo

  1. NIST/SEMATECH (2013). Confidence Limits for the Mean
  2. Diez, D.; Cetinkaya-Rundel, M.; Barr, C. (2019). OpenIntro Statistics