人常说,成名要趁早。那么我们将从数理统计的角度看,影响世界的科学家与年龄之间的相互关系。
假设此样本来自样本总体,回答以下问题:
1、求发现者当时的平均年龄
μ
\mu
μ的置信水平为95%的单侧置信下限。
2、求发现者当时的平均年龄
μ
\mu
μ的置信水平为95%的单侧置信上线。
3、求发现者当时的平均年龄
μ
\mu
μ置信水平为95%的双侧置信区间。
解答:
样本函数
t
=
X
ˉ
−
μ
S
/
n
∼
t
(
n
−
1
)
t=\frac{\bar X-\mu }{S/\sqrt{n} } \sim t\left ( n-1 \right )
t=S/nXˉ−μ∼t(n−1)。
1、
P
(
t
<
t
α
)
=
1
−
P
(
t
≥
t
α
)
=
1
−
α
P\left ( t< t_{\alpha } \right ) =1-P\left ( t\ge t_{\alpha } \right ) =1-\alpha
P(t<tα)=1−P(t≥tα)=1−α
即
P
(
X
ˉ
−
μ
S
/
n
<
t
α
)
=
1
−
α
P\left (\frac{\bar X-\mu }{S/\sqrt{n} }< t_{\alpha } \right ) =1-\alpha
P(S/nXˉ−μ<tα)=1−α
得
P
(
μ
>
X
ˉ
−
S
t
α
n
)
=
1
−
α
P\left (\mu > \bar X-\frac{St_{\alpha } }{\sqrt{n} } \right ) =1-\alpha
P(μ>Xˉ−nStα)=1−α
根据样本均值和样本方差有:
x
x
x=35.25,
s
2
s^2
s2=52.3864;
已知置信水平为1-
α
\alpha
α =0.95,
α
\alpha
α =0.05,查表得
t
α
(
n
−
1
)
=
t
0.05
(
11
)
=
1.796
t_{\alpha } \left ( n-1 \right )= t_{0.05 } \left (11 \right )= 1.796
tα(n−1)=t0.05(11)=1.796
解得
μ
l
=
35.25
−
7.2378
×
1.796
12
≈
31.497
\mu _{l} =35.25-\frac{7.2378\times 1.796}{\sqrt{12} } \approx 31.497
μl=35.25−127.2378×1.796≈31.497岁
发现者当时的平均年龄
μ
\mu
μ的置信水平为95%的单侧置信下限约为31岁。
2、
P
(
t
>
t
1
−
α
)
=
1
−
α
P\left ( t>t_{1-\alpha } \right ) =1-\alpha
P(t>t1−α)=1−α
即
P
(
X
ˉ
−
μ
S
/
n
>
t
1
−
α
)
=
1
−
α
P \left ( \frac{\bar X-\mu }{S/\sqrt{n}}>t_{1-\alpha } \right ) =1-\alpha
P(S/nXˉ−μ>t1−α)=1−α
得
P
(
μ
<
X
ˉ
−
S
t
α
n
)
=
1
−
α
P\left (\mu < \bar X-\frac{St_{\alpha } }{\sqrt{n} } \right ) =1-\alpha
P(μ<Xˉ−nStα)=1−α
t
1
−
α
=
−
t
α
(
11
)
=
−
1.796
t_{1-\alpha }=-t_{\alpha }\left ( 11 \right ) =-1.796
t1−α=−tα(11)=−1.796
t
α
(
n
−
1
)
=
t
0.05
(
11
)
=
1.796
t_{\alpha } \left ( n-1 \right )= t_{0.05 } \left (11 \right )= 1.796
tα(n−1)=t0.05(11)=1.796
解得
μ
l
=
35.25
+
7.2378
×
1.796
12
≈
39.0025
\mu _{l} =35.25+\frac{7.2378\times 1.796}{\sqrt{12} } \approx 39.0025
μl=35.25+127.2378×1.796≈39.0025岁
3、总体均值
μ
\mu
μ的置信水平为1-
α
\alpha
α的置信区间为
(
X
ˉ
−
S
t
α
2
n
,
X
ˉ
+
S
t
α
2
n
)
\left ( \bar X-\frac{St_{\frac{\alpha }{2} } }{\sqrt{n} } ,\bar X+\frac{St_{\frac{\alpha }{2} } }{\sqrt{n} } \right )
(Xˉ−nSt2α,Xˉ+nSt2α),带入可得到置信区间为(30.65,39.85)
发现者当时的平均年龄
μ
\mu
μ的置信水平为95%的单侧置信上限约为39岁。
总结与所想:发明者的重大发现所处的年龄服从X\sim N\left ( \mu ,\sigma^2 \right )
由于统计量
X
ˉ
\bar X
Xˉ与
S
2
S^2
S2是总体均值
μ
\mu
μ和总体方差
σ
2
\sigma^2
σ2的无偏、一致估计量。此时便可认为
X
∼
N
(
35.25
,
7.237
8
2
)
X\sim N\left ( 35.25 ,7.2378^2 \right )
X∼N(35.25,7.23782)
假设发明者的寿命为100岁,例用matlab软件绘制出重大发现年龄对应的概率密度函数。程序代码如下。
x=0:100;
age=[40,36,23,40,31,33,49,33,31,42,26,39];
mu=mean(age);
sigma=std(age);
y=normpdf(x,mu,sigma);
plot(x,y);
grid on;
运行后的图像为:
可以发现,科学家在30至40岁年龄之间,发现重大发明的几率最大。
置信限(Confidence Limit):你的年龄是否限制了你的创造力
最新推荐文章于 2022-11-22 22:13:26 发布