统计学(二):统计推断:速成包

公式表、每个专题的关键方法步骤和 10 道练习题,完整解答附在最后。

新题目返回课程

在打印窗口中选择“另存为 PDF”,即可在设备上保存一份副本。

1公式与方法

正态分布

z 分数
z = (x − μ)/σ
小于 x
P(X < x) = Φ(z)
大于 x
P(X > x) = 1 − Φ(z)
介于两者之间
Φ(z2) − Φ(z1)
百分位数
x = μ + z·σ

方法

  1. 把数值换算为 z 分数:z = (x − μ)/σ。
  2. 查 Φ(z),即 z 左侧的面积。
  3. 右侧面积用 1 − Φ(z);区间面积用两个左侧面积相减。
  4. 求百分位数时反向计算:先由面积求 z,再求 x = μ + z·σ。

置信区间

σ 已知
x̄ ± z·σ/√n
σ 未知
x̄ ± t·s/√n, df = n − 1
临界 z 值
90%: 1.645, 95%: 1.96, 99%: 2.576
样本量
n = (z·σ/E)², rounded up

方法

  1. 在 z(已知总体 σ)和 t(只有样本 s,自由度为 n − 1)之间作出选择。
  2. 求该置信水平下的临界值。
  3. 误差限:E = 临界值 × 标准差/√n。
  4. 区间:x̄ ± E。

假设检验与置信区间

均值的 z 检验
z = (x̄ − μ0)/(σ/√n)
均值的 t 检验
t = (x̄ − μ0)/(s/√n), df = n − 1
5% 双侧
若 |z| > 1.96 则拒绝 H0
置信区间
x̄ ± z·σ/√n

方法

  1. 写出 H₀ 和 H₁,并确定是单侧检验还是双侧检验。
  2. 计算检验统计量:(x̄ − μ₀)/(σ/√n),或者用 s 和 t。
  3. 与临界值比较,或求出 p 值。
  4. 若统计量超出临界值(p 小于 α),则拒绝 H₀,并用文字写出结论。

线性回归与相关

平方和
Sxx = Σ(x − x̄)² = Σx² − (Σx)²/n, Sxy = Σxy − ΣxΣy/n
斜率
b = Sxy/Sxx
截距
a = ȳ − b·x̄
相关系数
r = Sxy/√(SxxSyy)
残差
e = y − ŷ

卡方检验

检验统计量
χ² = Σ(O − E)²/E
拟合优度
E = n·p, df = k − 1 − (estimated parameters)
独立性
E = row total × column total / grand total, df = (r − 1)(c − 1)
临界值(α = 0.05)
df 1: 3.841, 2: 5.991, 3: 7.815, 4: 9.488
决策
若 χ² > χ²crit,则拒绝 H0

2练习题

1.正态分布中等
−3−2−10123zz 右侧的阴影面积 = 0.8599

拖动 z,使标准正态曲线下 z 右侧的阴影面积为 0.8599。求 z,保留两位小数。

A = 0.8599
z =
2.正态分布困难
235240245250255260265g顶部 15%

袋装重量 服从正态分布,均值为 250 g,标准差为 5 g。拖动截断值,使阴影的前 15% 位于其右侧。求截断值 x。

μ = 250 gσ = 5 g15%
x = g
3.置信区间中等
x̄E₁ = 4.590%, n = 64x̄E₂ = ?95%, n = 16

基于 64 个观测值且 σ 已知的均值 90% 置信区间,其误差范围为 4.5。 基于 16 个观测值的 95% 区间,其误差范围 E₂ 是多少?

90%n₁ = 64E₁ = 4.595%n₂ = 16
E₂ =
4.置信区间困难
p ≈ 0.8, n = ?p̂ − Ep̂ + Ep̂E = 0.0199%

一项预研究表明,约有 0.8 的产品有缺陷。 必须检验多少件产品 n,才能使比例的 99% 置信区间的误差范围不超过 0.01?

p ≈ 0.8E = 0.0199%
n =
5.假设检验与置信区间中等
t −2.01102.011

49 名学生的样本平均分为 54.11,样本标准差为 12。求 H₀: μ = 50 的检验统计量 t。

μ₀ = 50x̄ = 54.11s = 12n = 49
t =
6.假设检验与置信区间困难
−1.961.96z = ?

在一项试验中,接受治疗 A 的 300 名患者中有 60 人康复,接受治疗 B 的 150 名患者中有 39 人康复。用合并比例求 H₀: p₁ = p₂ 的检验统计量 z。

x₁ = 60n₁ = 300x₂ = 39n₂ = 150
z =
7.线性回归与相关中等
xy128226322419517617xy2462025

表中给出 6 对观测值 x = 车龄(年) 和 y = 价格 (k€)。 求皮尔逊相关系数 r。

n = 6
r =
8.线性回归与相关困难
xy351463564664772878xy46850607080x = 4.5

表中给出 6 对观测值 x = 学习小时数 和 y = 考试成绩。 使用最小二乘直线 ŷ = a + b·x。 新得到一个观测值 (4.5, 57)。求其残差 e = y − ŷ。

n = 6(4.5, 57)
e =
9.卡方检验中等
好转未好转合计新药233760安慰剂334780合计5684140

列联表给出了观测频数。计算行变量与列变量独立性检验的卡方统计量 χ²。不使用连续性校正。

N = 1402 × 2
χ² =
10.卡方检验困难
dfα = 0.05α = 0.0113.8416.63525.9919.21037.81511.34549.48813.277511.07015.086612.59216.812714.06718.475815.50720.090916.91921.6661018.30723.209

对一个 3 行 4 列(不含合计)的列联表进行卡方独立性检验。 利用表格,求显著性水平 α = 0.01 下的临界值 χ²crit。

α = 0.01
χ²crit =

3溶液

1.正态分布

答案: z = −1.08

由面积反求 z:Φ(z) 是 z 左侧的面积。对于左尾面积 A,在表中找出使 Φ(z) = A 的 z。右尾面积 A 意味着 Φ(z) = 1 − A。−z 与 z 之间的中间面积 A 使每个尾部剩下 (1 − A)/2,所以 Φ(z) = (1 + A)/2。对于数值 x,用 x = μ + zσ。

z 左侧的面积

Φ(z) = 1 − 0.8599 = 0.1401

反向查 z 表

在表中找到 0.1401,读出 z。

z = −1.08

答案

z = −1.08
2.正态分布

答案: x = 255.2 g

由面积反求 z:Φ(z) 是 z 左侧的面积。对于左尾面积 A,在表中找出使 Φ(z) = A 的 z。右尾面积 A 意味着 Φ(z) = 1 − A。−z 与 z 之间的中间面积 A 使每个尾部剩下 (1 − A)/2,所以 Φ(z) = (1 + A)/2。对于数值 x,用 x = μ + zσ。

左侧的面积

Φ(z) = 0.85

由表得到的 z

z ≈ 1.04

回到原来的单位

x = 250 + 1.04·5 = 255.2 g

答案

x = 255.2 g
3.置信区间

答案: E₂ = 10.72

Confidence interval for a mean: x̄ ± margin of error. With a known population σ the margin is z·σ/√n; with only the sample standard deviation s, use t with n − 1 degrees of freedom: t·s/√n. Critical values: z = 1.645 (90%), 1.96 (95%), 2.576 (99%). To get a margin E or smaller, you need n ≥ (z·σ/E)², rounded up.

误差范围之比

E₂/E₁ = (1.96/1.645)·√(64/16) = 2.383

新的误差范围

E₂ = 4.5·2.383 = 10.72

答案

E₂ = 10.72
4.置信区间

答案: n = 10618

Confidence interval for a proportion: p̂ ± z·√(p̂(1 − p̂)/n), where p̂ = x/n. To reach a margin E you need n ≥ z²·p(1 − p)/E², rounded up; without a prior estimate use p = 0.5, which gives the largest n.

解出 n

n = 2.576²·0.8·0.2/0.01² = 10617

向上取整

n = 10618

答案

n = 10618
5.假设检验与置信区间

答案: t = 2.397

均值检验:检验统计量衡量样本均值离 H₀ 所声称的值有多少个标准误:σ 已知时 z = (x̄ − μ₀)/(σ/√n),由样本估计时 t = (x̄ − μ₀)/(s/√n)。置信区间为 x̄ ±(临界值)·(标准误)。

标准误

SE = 12/√49 = 12/7 = 1.714

检验统计量

t = (54.11 − 50)/1.714 = 2.397

答案

t = 2.397
t0.025; df = 48 = 2.011

|t| > 2.011,所以在 5% 水平下拒绝 H₀(双侧)。

6.假设检验与置信区间

答案: z = −1.448

Hypothesis test: the test statistic = (estimate − value under H₀)/(standard error under H₀). The p-value is the probability, assuming H₀ is true, of a statistic at least as extreme as the one observed: the area in one tail for a one-sided H₁, twice the tail area for a two-sided H₁. Reject H₀ when p ≤ α.

样本比例

p̂₁ = 60/300 = 0.2
p̂₂ = 39/150 = 0.26
p̂ = 99/450 = 0.22

H₀ 下的标准误

SE = √(0.22·0.78·(1/300 + 1/150)) = 0.04142

检验统计量

z = (0.2 − 0.26)/0.04142 = −1.448

答案

z = −1.448

|z| < 1.96,所以在 5% 水平下不拒绝 H₀(双侧)。

7.线性回归与相关

答案: r = −0.9709

最小二乘回归:设 x̄ 和 ȳ 为均值,Sxx = Σ(x − x̄)²,Syy = Σ(y − ȳ)²,Sxy = Σ(x − x̄)(y − ȳ),则直线 ŷ = a + b·x 的斜率 b = Sxy/Sxx,截距 a = ȳ − b·x̄,因此它总是通过 (x̄, ȳ)。皮尔逊相关系数 r = Sxy/√(Sxx·Syy) 介于 −1 与 1 之间;r² 是 y 的变异中能被该直线解释的比例。

均值

x̄ = 21/6 = 3.5
ȳ = 129/6 = 21.5

平方和与乘积和

Sxx = Σ(x − x̄)² = 17.5
Sxy = Σ(x − x̄)(y − ȳ) = −42.5
Syy = Σ(y − ȳ)² = 109.5

相关系数

r = −42.5/√(17.5·109.5) = −0.9709

答案

r = −0.9709
8.线性回归与相关

答案: e = −3.705

最小二乘回归:设 x̄ 和 ȳ 为均值,Sxx = Σ(x − x̄)²,Syy = Σ(y − ȳ)²,Sxy = Σ(x − x̄)(y − ȳ),则直线 ŷ = a + b·x 的斜率 b = Sxy/Sxx,截距 a = ȳ − b·x̄,因此它总是通过 (x̄, ȳ)。皮尔逊相关系数 r = Sxy/√(Sxx·Syy) 介于 −1 与 1 之间;r² 是 y 的变异中能被该直线解释的比例。

均值

x̄ = 33/6 = 5.5
ȳ = 392/6 = 65.33

平方和与乘积和

Sxx = Σ(x − x̄)² = 17.5
Sxy = Σ(x − x̄)(y − ȳ) = 81

回归直线

b = 81/17.5 = 4.629
a = 65.33 − 4.629·5.5 = 39.88
ŷ = 39.88 + 4.629·x

预测值

ŷ = 39.88 + 4.629·4.5 = 60.7

残差

e = 57 − 60.7 = −3.705

答案

e = −3.705

该点位于直线下方。

9.卡方检验

答案: χ² = 0.1215

卡方检验:用 χ² = Σ(O − E)²/E 比较观测频数 O 与 H₀ 下的期望频数 E。拟合优度检验:每个类别 E = n·p,df = k − 1(每由数据估计一个参数再减一)。r × c 列联表的独立性检验:E = 行合计 × 列合计 / 总计,df = (r − 1)(c − 1)。当 χ² 超过所选 α 下的临界值时拒绝 H₀。

期望频数

新药: 60·56/140 = 24, 60·84/140 = 36
安慰剂: 80·56/140 = 32, 80·84/140 = 48

卡方统计量

χ² = (23 − 24)²/24 + (37 − 36)²/36 + (33 − 32)²/32 + (47 − 48)²/48
χ² = 0.1215

答案

χ² = 0.1215

df = 1。在 α = 0.05 时临界值为 3.841,所以 H₀(独立)不被拒绝。

10.卡方检验

答案: χ²crit = 16.81

卡方检验:用 χ² = Σ(O − E)²/E 比较观测频数 O 与 H₀ 下的期望频数 E。拟合优度检验:每个类别 E = n·p,df = k − 1(每由数据估计一个参数再减一)。r × c 列联表的独立性检验:E = 行合计 × 列合计 / 总计,df = (r − 1)(c − 1)。当 χ² 超过所选 α 下的临界值时拒绝 H₀。

自由度

df = (3 − 1)·(4 − 1) = 6

表

χ²crit(df = 6, α = 0.01) = 16.812

答案

χ²crit = 16.812

若计算出的 χ² 大于该值,则拒绝 H₀。

ExamProblems · examproblems.com/courses/statistics-2-inference