誤差について（メモ）

（[ブログ本編](https://datadriven-rnd.com/about-error/)にも投稿してます）

測定には誤差が必ず伴うため、測定によって得られた値には不確かさがあることに留意する。
誤差が測定ごとにランダムであると仮定すれば、データを統計処理することで不確かさを小さくすることができる。

## 誤差（測定誤差）

<strong>誤差</strong>：測定値と真値との差（[誤差]ー[真値] で求められる。絶対誤差とも呼ぶこともある）
<strong>相対誤差</strong>：誤差と真値との比。（実際には真の値は分からないので、実際は測定値から真の値を推定し相対誤差の大きさを推定する）

ある量$x$を測定した時、$x$ の測定値$x$(measured)は<u>必ず</u>誤差を含む。
この時、$x$(measured)は、最良推定値（$x_{best}$)、誤差（不確かさ｜$\delta x$）を用いて、

$$ x(measured) = x_{best} ± \delta x$$

と表される。また、相対誤差（relative uncertainty）はは以下の式で表される。

$$ {relative \, uncertainty} = \frac{\delta x}{x_{best}} $$

<h2>誤差の種類、正確さと精密さについて</h2>

<h4>誤差の種類</h4>

- <strong>系統誤差（systematic errors）</strong>：再現性がある誤差。一方向。繰り返し測定を行っても原因は明らかにならないことが多い。誤差の原因がわかれば測定や解析の改良（測定値補正など）で誤差を小さくすることができる。<br> 例）理論（モデル）の誤り、測定器固有の特性など。
- <strong>ランダム誤差（random errors）
</strong>：制御できない種類の不規則な誤差のこと。偶発誤差ともいう。繰り返し測定やそれと統計的手法を組み合わせることで誤差の大きさを評価したり、誤差を小さくすることが可能。<br>例）ランダムな測定ムラ、測定器の精度限界など

<h4>正確さと精密さについて</h4>
<ul>
<li><strong>accuracy</strong>：正確さ。真の値にいかに近いかの尺度<br>
（accuracyが高い=系統誤差が小さい）</li>
<li><strong>precision</strong>：精密さ。一連の実験結果がどの程度一致しているかの尺度。真の値との関係は示していない。<br>
（precisionが高い=ランダム誤差が小さい）</li>
</ul>

<u>真の値が分かっているか否かで、accuracy, precision, 系統誤差を明らかにできるかどうかが決まる（下表参照）</u>

![スクリーンショット 2021-06-23 23.03.39.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/348376/b9963c1e-1845-0a39-f6bc-1f08ee27b6c6.png)

## ランダム誤差を小さくするには
ランダム誤差のみ考えた場合、誤差を評価し、測定精度を高めるためには<u>同じ測定を何回も実施し、その値から統計処理によって真値やランダム誤差を見積もる</u>ことが必要となる。

### （母分散が未知の場合の）母平均の信頼区間の求め方
母集団の平均を $\mu$、$x$ 標本平均を $\bar{x}$、不偏分散を $s^2$、抽出したサンプルサイズを $n$、信頼係数を $100(1 - \alpha)$%とすると、下式から母平均 $\mu$ の $100(1-\alpha))$% 信頼区間を求めることができる。
信頼係数：母数を区間推定するときに母数が信頼区間内に含まれる確率のこと。（[参考](https://bellcurve.jp/statistics/glossary/2009.html)）

<p>$$ \bar{x}-t_{\alpha/2}(n-1)\times\sqrt{\frac{s^2}{n}}\leq\mu\leq\bar{x} +t_{\alpha/2}(nー1)\times\sqrt{\frac{s^2}{n}} $$</p>


### 解析例）
ある天秤で、ある物質のその重さを10回測定したところ、下のデータが得られた。母分散（≒電子天秤の誤差,分散）は分かっておらず、天秤の誤差は正規分布に従うものとした場合の、このデータの母平均の95%信頼区間を求める。<br>
（上記の式を使い、母平均の信頼区間を求めてみる）

![スクリーンショット 2021-06-23 23.15.29.png](https://qiita-image-store.s3.ap-northeast-1.amazonaws.com/0/348376/6ff15b24-34e7-f5ce-fabe-9629e5ecf858.png)

1. 標本平均 $ \bar{x} $ と不偏分散 $s^2$ を求める。
  $$  \bar{x} = (100.5 + 99.5 + \cdots + 100.4) / 10 = 100.0 \\\
  s^2 = \frac{1}{10-1} \\{ (100.5-100.0)^2 + (99.5-100.0)^2 + \cdots + (100.4-100.0)^2 \\} = 0.545 $$  
  

2. t分布の上側 $\alpha/2$% 点を求める。
この場合、$\alpha/2 = (100-95)/2=2.5 $% なので、自由度9のt分布の上側2.5%点を求める。
$$ t_{\alpha/2} = 2.262 $$
  

3. 式より95%信頼区間を求める
<p>$$ 100-2.262 \times \sqrt{\frac{0.545}{10}}\leq\mu\leq100-2.262 \times \sqrt{\frac{0.545}{10}} $$</p>
よって、
<p>$$ 99.5 \leq\mu\leq  100.5 $$</p>

#### pythonで計算してみた

```python
import numpy as np
from scipy import stats

data = np.array([100.5, 99.5, 98.9, 101.5,100.1, 99.8, 99.9, 99.2, 100.3, 100.4])

#平均と不偏分散を求める
mean = round(data.mean(),2)
var = round(data.var(ddof=1),2)
print('平均：',mean)
print('分散：',var)

#t分布の上側 α/2 % 点を求める
upper_0025 = round(stats.t.ppf(q=0.975, df=len(data)-1),3)
print('上側2.5%点：',upper_0025,3)

#95%信頼区間を求める
lower = round(mean - upper_0025 * np.sqrt(var/len(data)),3)
upper = round(mean + upper_0025 * np.sqrt(var/len(data)),3)

print('95%信頼区間：',lower,'≦',mean,'≦',upper)
```

<h2>参考</h2>

- [母平均の信頼区間の求め方（母分散未知）|統計WEB](https://bellcurve.jp/statistics/course/8972.html)
- [誤差解析 Error Anlysis](http://www.fm.ehcc.kyoto-u.ac.jp/web_material/data_analysis.pdf)
- [はじめての誤差論](https://dora.bk.tsukuba.ac.jp/~takeuchi/?%E3%81%AF%E3%81%98%E3%82%81%E3%81%A6%E3%81%AE%E8%AA%A4%E5%B7%AE%E8%AB%96)