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

公式表、每个专题的关键方法步骤和 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.正态分布中等
340 kg328 kgT = X₁ + … + X₄?

带行李的乘客的质量相互独立,服从均值 85 kg、标准差 12 kg 的正态分布。求其中 4 个的总质量低于 328 kg 的概率 P。z 取两位小数。

μ = 85 kgσ = 12 kgn = 4< 328 kg
P =
2.正态分布困难
35455565758595分底部 80%

考试成绩 服从正态分布,均值为 65 分,标准差为 10 分。拖动截断值,使阴影的后 80% 位于其左侧。求截断值 x。

μ = 65 分σ = 10 分80%
x = 分
3.置信区间中等
n = 20, s = 8LUx̄ = 107.5E = ?95%

一个包含 20 个个体的样本,均值为 107.5,标准差 s = 8。总体近似服从正态分布。 求均值的 95% 置信区间的误差范围 E。

n = 20x̄ = 107.5s = 895%t(19) = 2.093
E =
4.置信区间困难
p = 0.5, n = ?p̂ − Ep̂ + Ep̂E = 0.0595%

事先对缺陷产品的比例一无所知。 必须检验多少件产品 n,才能使比例的 95% 置信区间的误差范围不超过 0.05?

p = 0.5E = 0.0595%
n =
5.假设检验与置信区间中等
1.96p = ?

据称某总体的均值为 50;已知总体标准差为 36。一个含 144 个观测值的随机样本均值为 55.88。检验 H₀: μ = 50 对 H₁: μ > 50。求 p 值 p。

μ₀ = 50σ = 36n = 144x̄ = 55.88H₁: μ > μ₀
p =
6.假设检验与置信区间困难
σ = 10, n = ?x̄ − Ex̄ + Ex̄E = 0.595%

样本至少要多大,才能使均值的 95% 置信区间的误差范围不超过 0.5?标准差为 10。求最小样本量 n。

σ = 10E = 0.595%
n =
7.线性回归与相关中等
xy1373505617749851193xy510406080

表中给出 6 对观测值 x = 学习小时数 和 y = 考试成绩。 求皮尔逊相关系数 r。

n = 6
r =
8.线性回归与相关困难
xy2941261782410301234xy510102030ŷ = 10

表中给出 6 对观测值 x = 荷载 (kN) 和 y = 伸长量 (mm)。 使用最小二乘直线 ŷ = a + b·x。 x 取何值 x 时,直线预测 ŷ = 10?

n = 6ŷ = 10
x =
9.卡方检验中等
好转未好转合计新药5446100安慰剂2674100合计80120200

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

N = 2002 × 2
χ² =
10.卡方检验困难
合格不合格合计工厂 A295180工厂 B322860工厂 C275380合计88132220dfα = 0.05α = 0.0113.8416.63525.9919.21037.81511.34549.48813.277511.07015.086612.59216.812714.06718.475815.50720.090916.91921.6661018.30723.209

表中给出了观测频数。在 α = 0.05 下检验两个变量是否独立(不使用连续性校正)。计算检验统计量 χ²,再与表中的临界值比较。

N = 2203 × 2α = 0.05
χ² =

3溶液

1.正态分布

答案: P = 0.3085

正态变量的和与差:若 X 和 Y 相互独立且服从正态分布,则 X + Y 与 X − Y 也服从正态分布。均值相加或相减:E(X ± Y) = μX ± μY。方差总是相加,差也是如此:Var(X ± Y) = σX² + σY²,所以 σ = √(σX² + σY²)。标准差从不直接相加。n 个独立同分布项的总和,均值为 n·μ,标准差为 σ·√n。然后照常标准化:z = (x − 均值)/σ,再查 Φ(z)。

总和的分布

μT = 340 kg, σT = 12·√4 = 24 kg

标准化

z = (328 − 340)/24 = −0.5

概率

P = Φ(−0.5) = 0.3085

答案

P = 0.3085
2.正态分布

答案: x = 73.4 分

由面积反求 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.8

由表得到的 z

z ≈ 0.84

回到原来的单位

x = 65 + 0.84·10 = 73.4 分

答案

x = 73.4 分
3.置信区间

答案: E = 3.744

均值的置信区间:x̄ ± 误差范围。总体 σ 已知时,误差范围为 z·σ/√n;只有样本标准差 s 时,用自由度为 n − 1 的 t:t·s/√n。临界值:z = 1.645(90%)、1.96(95%)、2.576(99%)。要使误差范围不超过 E,需要 n ≥ (z·σ/E)²,并向上取整。

临界值

df = n − 1 = 19, t = 2.093

误差范围

E = 2.093·8/√20 = 3.744

答案

E = 3.744
4.置信区间

答案: n = 385

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 = 1.96²·0.5·0.5/0.05² = 384.2

向上取整

n = 385

答案

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

答案: p = 0.025

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 ≤ α.

检验统计量

SE = σ/√n = 36/√144 = 3
z = (55.88 − 50)/3 = 1.96

尾部面积

Φ(1.96) = 0.9750
1 − Φ(1.96) = 0.0250

p 值

p = 0.025

答案

p = 0.025

p ≤ 0.05,所以在 5% 水平下拒绝 H₀。

6.假设检验与置信区间

答案: n = 1537

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

解出 n

n = (1.96·10/0.5)² = 1537

向上取整

n = 1537

答案

n = 1537
7.线性回归与相关

答案: r = 0.9977

最小二乘回归:设 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̄ = 36/6 = 6
ȳ = 400/6 = 66.67

平方和与乘积和

Sxx = Σ(x − x̄)² = 70
Sxy = Σ(x − x̄)(y − ȳ) = 398
Syy = Σ(y − ȳ)² = 2273

相关系数

r = 398/√(70·2273) = 0.9977

答案

r = 0.9977
8.线性回归与相关

答案: x = 2.86

最小二乘回归:设 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̄ = 42/6 = 7
ȳ = 126/6 = 21

平方和与乘积和

Sxx = Σ(x − x̄)² = 70
Sxy = Σ(x − x̄)(y − ȳ) = 186

回归直线

b = 186/70 = 2.657
a = 21 − 2.657·7 = 2.4
ŷ = 2.4 + 2.657·x

解出 x

x = (10 − 2.4)/2.657 = 2.86

答案

x = 2.86
9.卡方检验

答案: χ² = 16.33

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

期望频数

新药: 100·80/200 = 40, 100·120/200 = 60
安慰剂: 100·80/200 = 40, 100·120/200 = 60

卡方统计量

χ² = (54 − 40)²/40 + (46 − 60)²/60 + (26 − 40)²/40 + (74 − 60)²/60
χ² = 16.33

答案

χ² = 16.33

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

10.卡方检验

答案: χ² = 6.215

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

期望频数

工厂 A: 32, 48
工厂 B: 24, 36
工厂 C: 32, 48

卡方统计量

χ² = (29 − 32)²/32 + (51 − 48)²/48 + (32 − 24)²/24 + (28 − 36)²/36 + (27 − 32)²/32 + (53 − 48)²/48
χ² = 6.215

决策

df = 2, χ²crit = 5.991
6.215 > 5.991

答案

χ² = 6.215

χ² 高于临界值,所以拒绝 H₀:两个变量看起来是相关的。

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