確率プロットとQ-Qプロット

データセットが与えられた分布に従うかどうか検定するために確率プロットを使用します。これは、X軸上に観測した累積パーセント、Y軸上に期待累積パーセントを持つグラフを表示します。すべてのデータポイントが参照線に近い場合、データセットは与えられた分布に従うものと結論付けできます。

Q-Q(Quantile-Quantile) プロットは、データセットが与えられた分布に従うかどうかを検定する手法です。これは、確率プロットとは異なり、XとY軸上でパーセントの変わりに実測値と期待値を示します。すべてのデータポイントが参照線に近い場合、データセットは与えられた分布に従うものと結論付けできます。

Originでは、4つの分布(正規, 正規対数, 指数, ワイブル)と5つのパーセンタイル近似のプロット手法(Blom, Benard, Hazen, Van der Waerden, Kaplan-Meier)をサポートしています。

確率プロット/Q-Qプロットを作成するには、

確率プロット、または、Q-Qプロットを作成するには、

  1. 1つのY列を選択します。
  2. 確率/Q-Qプロットダイアログを開くには、
    確率プロット:メインメニューで、作図をクリックして、統計を選択して、確率プロット:P-P図をクリックします。 あるいは、 2Dグラフツールバーからの確率プロットボタン(P-P図)Button Probability Plot.pngをクリックします。
    Q-Qプロット:メインメニューで、作図をクリックして、統計を選択して、Q-Q図をクリックします。 あるいは、 2DグラフツールバーからのQ-Q図ボタンButton Q Q Plot.png をクリックします。
  3. plot_prob Xファンクションダイアログで、分布と方法を指定します。
  4. OKをクリックして、確率プロット、または、Q-Qプロットを作成します。
    Probability Plot and Q Q Plot 2.png

plot_probXファンクションのダイアログ

Probability Plot and Q Q Plot 1.png

入力データ

入力データを指定します。

分布

データに対する分布の種類を選択します。分布についてのより詳細な情報は、分布セクションを参照してください。

分布
4つの分布が利用できます。
  • 正規分布
  • 正規対数
  • 指数分布
  • Weibull
  • ガンマ
データから推定
入力データから分布パラメータを推定するかどうかを指定します。チェックを付けない場合、パラメータを手動で指定することができます。
u
正規分布の平均
シグマ
正規分布の標準偏差。
形状
指定した分布の形状。正規対数、Weibull、ガンマの分布で利用可能です。
スケール
指定した分布のスケール。正規対数、指数分布、Weibull、ガンマの分布で利用可能です。

スコア法

パーセンタイル近似をプロットする方法を選択します。手法についてのより詳細な情報は、スコア法セクションを参照してください。

  • Blom
  • Benard
  • Hazen
  • Van der Waerden
  • Kaplan-Meier

信頼帯

確率プロットの信頼帯を出力するかどうか指定します。計算の詳細は、アルゴリズムを確認してください。

信頼度 (%)

信頼帯が選択されている場合のみ利用できます。選択した分布に対するパーセントで、信頼水準を指定します。

XY軸の交換

X軸とY軸を変更するかどうか指定します。

X最小
X最大

自動の値は、X最小 = 1 および X最大 = 99.5です。自動のチェックを外すと、出力の参照線列の最小および最大値を使用します。

X最大が自動の値より大きい場合、X最大のパーセント値p1を計算して、 パーセンタイル列には、デフォルトリストにp1およびp1より大きい値のみ含める必要があります。X最大自動の値より小さい場合、X最大のパーセント値p2を計算して、 パーセンタイル列には、デフォルトリストにp2およびp2より小さい値のみ含める必要があります。

X最小が自動の値よりも小さい場合、X最小値のパーセンタイル値p1を計算します。p1<1e-5, p1=1e-5の場合、p1よりも大きい最小値10^(-m)を見つけ、パーセンタイル列には p1, 10^(-m), 10^(-m+1), ,,,,1, 2,...が含まれます。

X最大が自動の値よりも大きい場合、X最大値のパーセンタイル値p2を計算します。p2>99.99, p2=99.99の場合、リスト (99.9, 99.99) から2より小さい最大値をみつけ、パーセンタイル列には99, 99.5, 99.9,..p2.が含まれます。

出力範囲

グラフのために計算されたデータを保存する場所を設定します。

分布

Originには、確率プロットとQ-Qプロットに対して、4つの分布があります。 以下の表に、密度関数を示します。

分布 密度関数p(x) 範囲 パラメータ

正規分布

\frac 1{\sigma \sqrt{2\pi }}\exp \left( -\frac{\left( x-\mu \right) ^2}{2\sigma ^2}\right)

すべてx

\mu平均は、位置パラメータ
で、\sigma(>0)標準偏差は、スケールパラメータです。

正規対数

\frac 1{\sigma x\sqrt{2\pi }}\exp \left( -\frac{\left( \ln \left( x\right) -\mu \right) ^2}{2\sigma ^2}\right)

x>0

\muは形状パラメータ、
\sigma(>0)はスケールパラメータです。

指数分布

\frac 1\sigma \exp \left( -\frac x\sigma \right)

x>0 \sigma(>0)は、スケールパラメータです。

Weibull

\frac c\sigma \left( \frac x\sigma \right) ^{c-1}\exp \left( -\left( \frac x\sigma \right) ^c\right)

x>0

\sigma(>0)はスケールパラメータ、
c(>0) は形状パラメータです。

ガンマ

\frac{1}{\Gamma(c)\sigma^c}x^{c -1} exp(-x/\sigma),

x>0

\sigma(>0)はスケールパラメータ、
c(>0) は形状パラメータです。

確率プロット作成の詳細

確率プロットを作成して、小さいデータセットから大きいものの順にソートするには、

x[1]\le x[2]\le x[3]\le \cdots \le x[n-1]\le x[n], nは観測データセットの合計数です。

ソートした観測値は、X座標がx[i]\ 、Y座標が スコア法で計算されたポイントで代表されています。

Probability Plot and Q Q Plot 7.png

確率プロットのスケールタイプは、分布により異なります。

分布 Xスケールタイプ Yスケールタイプ

正規分布

線形

確率

正規対数

Ln

確率

指数分布

Ln

2重対数逆数

Weibull

Log10

2重対数逆数

ガンマ

Log10

確率

Q-Qプロット作成の詳細

Q-Qプロットを作成して、小さいデータセットから大きいものの順にソートするには、

x[1]\le x[2]\le x[3]\le \cdots \le x[n-1]\le x[n], ここで、nは、 観測データセットの合計数です。

Y値は、スコア法の逆累積分布関数です。

Probability Plot and Q Q Plot 6.png

スコア法

以下に示した方法のひとつを使って、入力値を小さい値から大きい値の順に並べ、ソートデータのシリアル番号がスコア付けされます。 この表では、iがシリアル番号で、nは、非欠損入力データです。

手法 プロット位置method(i,n)

Blom

(i-0.375)/(n+0.25)

Benard

(i-0.3)/(n+0.4)

Hazen

(i-0.5)/n

Van der Waerden

i/(n+1)

Kaplan-Meier

i/n

参考文献

  • Samuel Kotz , Campbell B. Read , N. Balakrishnan, Brani Vidakovic, 2005.Encyclopedia of statistical sciences., NewYork: John Wiley & Sons, Inc.
  • Thode, Henry C. 2002, Testing for Normality, CRC Press