UNIVARIATEプロシジャ

例4.3 極値のオブザベーションと極値の識別

この例では、例4.1で説明したデータセットBPressureを使用し、データセット内の極値のオブザベーションのテーブルと極値のテーブルを作成する方法を示します。次のステートメントは、SystolicとDiastolicに対する"Extreme Observations"テーブルを生成します。このテーブルで、各変数の極値のオブザベーションを識別できます。

title 'Extreme Blood Pressure Observations';
ods select ExtremeObs;
proc univariate data=BPressure;
   var Systolic Diastolic;
   id PatientID;
run;

ODS SELECTステートメントは出力を"ExtremeObs"テーブルに制限します。ODSテーブル名のセクションを参照してください。IDステートメントは、極値のオブザベーションがオブザベーション番号に加えてPatientIDの値で識別されるよう要求します。デフォルトでは、上位5件および下位5件のオブザベーションが表示されます。NEXTROBS=オプションを使用すると、異なる数の極値のオブザベーションを要求できます。

出力4.3.1はSystolicおよびDiastolicともに、'CP' (オブザベーション7)と識別される患者の値が最も高いことを示しています。極値のオブザベーションを視覚的に示すため、ヒストグラムを作成できます。例4.14を参照してください。

次のステートメントは、SystolicとDiastolicに対する"Extreme Values"テーブルを生成します。これは分布の裾の表です。

title 'Extreme Blood Pressure Values';
ods select ExtremeValues;
proc univariate data=BPressure nextrval=5;
   var Systolic Diastolic;
run;

ODS SELECTステートメントは出力を"ExtremeValues"テーブルに制限します。ODSテーブル名のセクションを参照してください。NEXTRVAL=オプションは、分布の両端の極値の数を出力4.3.2のテーブルに示されているように指定します。

出力4.3.2は、値78および80がDiastolicで2回出現することと、Diastolicの最大値が110であることを示しています。なお、出力4.3.1ではDiastolicに値80が2回表示されています。これはこの値を持つオブザベーションが2つあるためです。出力4.3.2では、値80は1回だけ表示されています。

この例のサンプルプログラムuniex01.sasは、Base SASソフトウェアのSASサンプルライブラリに含まれています。