無相関だけど独立でないデータの例

「講義で使える統計素材」シリーズ．今回は，相関の説明の際に，相関があるということとデータが独立であるということは違うという説明をする際に必要そうな図を作りました．

In [1]:

import numpy as np
import matplotlib.pyplot as plt

極座標で円を描き半径にノイズを加える¶

In [2]:

delta = 1/100
theta = np.arange(0, 2*np.pi, delta)
e = np.random.normal(0, 0.1, theta.shape[0])

r = 1 + e
x = r*np.cos(theta)
y = r*np.sin(theta)

#軸の設定
plt.axis([-2, 2, -2, 2])
plt.gca().set_aspect('equal', adjustable='box')
print(np.corrcoef(x,y))

#描画
plt.scatter(x, y, s=1)
plt.savefig('plot_out.svg')

[[ 1.         -0.00289982]
 [-0.00289982  1.        ]]