UNIVARIATEプロシジャ

 

例4.10 分位点とパーセント点に対する信頼限界の計算

この例は例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を参照してください。

出力4.10.1 正規分布に基づく分位点の信頼限界
Analysis of Female Heights

The UNIVARIATE Procedure
Variable: Height (Height (in))

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に分布に依存しない信頼限界が示されています。

出力4.10.2 分布に依存しない分位点の信頼限界
Analysis of Female Heights

The UNIVARIATE Procedure
Variable: Height (Height (in))

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サンプルライブラリに含まれています。