この例は例4.9の続きで、分位点とパーセント点に対する信頼限界の計算方法を示しています。あるもう1人の調査員は、平均や標準偏差ではなく分位点を使用して身長を集計し、 分位点に対する90%信頼区間も計算してみようと考えています。次のステートメントは母集団の分位点の推定された分位点と信頼限界を作成します。
title 'Analysis of Female Heights'; ods select Quantiles; proc univariate data=Heights ciquantnormal(alpha=.1); var Height; run;
ODS SELECTステートメントは出力を"Quantiles"テーブルに制限します。ODSテーブル名のセクションを参照してください。CIQUANTNORMALオプションは分位点に対する信頼限界を作成します。出力4.10.1に示されているとおり、これらの限界はデータが正規分散であることを前提にしています。この前提を確認してから、これらの信頼限界を使用する必要があります。UNIVARIATEプロシジャにおける正規性のShapiro-Wilk検定の詳細は、Shapiro-Wilk統計量のセクションを参照してください。正規性検定の使用例は、例4.19を参照してください。
出力4.10.1: 正規分布に基づく分位点の信頼限界
Analysis of Female Heights |
Quantiles (Definition 5) | |||
---|---|---|---|
Level | Quantile | 90% Confidence Limits Assuming Normality |
|
100% Max | 70.0 | ||
99% | 70.0 | 68.94553 | 70.58228 |
95% | 68.6 | 67.59184 | 68.89311 |
90% | 67.5 | 66.85981 | 68.00273 |
75% Q3 | 66.0 | 65.60757 | 66.54262 |
50% Median | 64.4 | 64.14564 | 64.98770 |
25% Q1 | 63.1 | 62.59071 | 63.52576 |
10% | 61.6 | 61.13060 | 62.27352 |
5% | 60.6 | 60.24022 | 61.54149 |
1% | 60.0 | 58.55106 | 60.18781 |
0% Min | 60.0 |
UNIVARIATEプロシジャを使用して、正規分布を前提としない分位点の信頼限界を計算することもできます。次のステートメントはCIQUANTDFオプションを使用して、身長の母集団の分位点に対する分布に依存しない信頼限界を要求します。
title 'Analysis of Female Heights'; ods select Quantiles; proc univariate data=Heights ciquantdf(alpha=.1); var Height; run;
分布に依存しない信頼限界を出力4.10.2に示します。
出力4.10.2: 分布に依存しない分位点の信頼限界
Analysis of Female Heights |
Quantiles (Definition 5) | ||||||
---|---|---|---|---|---|---|
Level | Quantile | Order Statistics | ||||
90% Confidence Limits Distribution Free |
LCL Rank | UCL Rank | Coverage | |||
100% Max | 70.0 | |||||
99% | 70.0 | 68.6 | 70.0 | 73 | 75 | 48.97 |
95% | 68.6 | 67.5 | 70.0 | 68 | 75 | 94.50 |
90% | 67.5 | 66.6 | 68.6 | 63 | 72 | 91.53 |
75% Q3 | 66.0 | 65.7 | 66.6 | 50 | 63 | 91.77 |
50% Median | 64.4 | 64.1 | 65.1 | 31 | 46 | 91.54 |
25% Q1 | 63.1 | 62.7 | 63.7 | 13 | 26 | 91.77 |
10% | 61.6 | 60.6 | 62.7 | 4 | 13 | 91.53 |
5% | 60.6 | 60.0 | 61.6 | 1 | 8 | 94.50 |
1% | 60.0 | 60.0 | 60.5 | 1 | 3 | 48.97 |
0% Min | 60.0 |
出力4.10.2のテーブルには、計算された信頼限界の順位も示されています。信頼限界の計算方法の詳細は、パーセント点の信頼限界のセクションを参照してください。分位点の信頼限界は、WEIGHTステートメントが使用されている場合は作成されません。
この例のサンプルプログラムuniex07.sasは、Base SASソフトウェアのSASサンプルライブラリに含まれています。