PRML Old Faithful間欠泉データ

PRMLでOld Faithful間欠泉データのクラスタリングをしたいけど、データがリンク切れでなかなか落ちていなかった。faithful.txtで検索すると出てきた。
2つの変数は「噴出継続時間」と「次の噴出までの時間」らしい。
一応プロットしてみる。

import numpy as np
import matplotlib.pyplot as plt

data = np.loadtxt("faithful.txt")

plt.plot(data[:,0], data[:,1], "o", color=(0.9, 0.0, 0.1))
plt.show()

colorでプロットの色もしていできるので、濃度の変更もよしと。