この例は例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に示されているとおり、これらの限界はデータが正規分散であることを前提にしています。この前提を確認してから、これらの信頼限界を使用する必要があります。PROC UNIVARIATEにおける正規性に関するShapiro-Wilk検定の詳細は、Shapiro-Wilk統計量のセクションを参照してください。正規性検定の使用例は、例4.19を参照してください。
Analysis of Female Heights |
Quantiles (Definition 5) | |||
---|---|---|---|
Quantile | Estimate | 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 |
PROC UNIVARIATEを使用して、正規分布を前提としない分位点の信頼限界を計算することもできます。次のステートメントはCIQUANTDFオプションを使用して、身長の母集団の分位点に対する分布に依存しない信頼限界を要求します。
title 'Analysis of Female Heights'; ods select Quantiles; proc univariate data=Heights ciquantdf(alpha=.1); var Height; run;
出力4.10.2に分布に依存しない信頼限界が示されています。
Analysis of Female Heights |
Quantiles (Definition 5) | ||||||
---|---|---|---|---|---|---|
Quantile | Estimate | 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サンプルライブラリに含まれています。