質問 無作為抽出を実行できるプロシジャはありますか？ 回答 SAS 8.1以降では、SURVEYSELECTプロシジャを使用してさまざまな抽出を実行できます。   以下のサンプルプログラムは、706オブザベーションからなるサンプルデータから100件のオブザベーションを無作為抽出（非復元抽出、Random sampling without replacement）します。   /* サンプルデータの作成 */ data test; keep group type x; do group=1 to 10; do type=1 to 2; r=int(ranuni(12345)*50+10); do j=1 to r; x=rannor(12346); output; end; end; end; run; /* SURVEYSELECTプロシジャの実行 */ /* METHOD=SRSで単純無作為抽出（非復元抽出）を指定 */ proc surveyselect data=test method=srs rep=5 n=100 seed=12345 out=out; /* n=100の指定で100オブザベーションを抽出 */ run;   METHOD=オプションにおける指定をMETHOD=URSに変更すると、重複を許した無作為抽出（復元抽出）を行なうことも可能です。   また、層別抽出を行なうにはSTRATAステートメントを使用します。この場合、STRATAステートメントで指定する変数について、事前にソートしておく必要があります。   /* 同じサンプルデータを使用します */ /* SORTプロシジャ */ proc sort data=test out=sortout; by group type; run; proc surveyselect data=sortout method=srs rate=.1 seed=12345 out=out2; /* rate=.1で各層からそれぞれ10%のオブザベーションを抽出 */ strata group type; /* 層を表す変数を指定 */ run;   上記プログラムで指定しているステートメントおよびオプションは下記のとおりです。 STRATA 層を表す変数を指定する N= サンプリング数をオブザベーションの数で指定する オプション RATE= サンプリング数を割合で指定するオプション SEED= ランダムサンプリングを行なう際のシード値を 設定する METHOD= サンプリングの方法を指定する OUT= サンプリング後の出力データセットを指定する REP= 無作為抽出を行なう回数を指定する     各オプションの指定方法など、SURVEYSELECTプロシジャの詳細に関しては「SAS/STAT User's Guide」の「The SURVEYSELECT Procedure」の項をご参照ください。 ... View more

質問 SASデータセットの内容をCSVファイルにしています。その際、文字変数のデータに「"(ダブルクォーテーション)」をつけた形で出力したいのですが、何かよい方法はありますか。 回答 FILEステートメントにDSDオプションを指定できるようになったため、CSVファイルを簡単に作成できます。なお、文字変数のデータに「"(ダブルクォーテーション)」を追加するためには、$QUOTEフォーマットが便利です。次のサンプルプログラムを参考にしてください。   DATA test; INPUT key $ a1 a2; CARDS; 001 1 1 ; DATA _null_; SET test; FILE "c:\temp\test.csv" DSD; FORMAT key $QUOTE.; PUT key a1 a2; RUN;   ... View more

質問 SASデータセット中のある変数の最大値と最小値を簡単に求めるにはどうすればよいでしょうか。 回答 MEANSプロシジャを使用して簡単に求めることができます。 in_dataというSASデータセットの変数名hensuの最大値と最小値を求める例を下記に示します。   PROC MEANS DATA=in_data NOPRINT; VAR hensu; OUTPUT OUT=out_data(DROP=_TYPE_ _FREQ_) MAX=MAX MIN=MIN; RUN;   上記のプログラムを実行すると、out_dataというSASデータセット中の変数名MAXに最大値が、変数名MINに最小値がセットされます。 MEANSプロシジャでは、このほかにも平均値や標準偏差等の統計量を簡単に求めることができます。   参考: Base SAS® 9.4 Procedures Guide, Seventh Edition http://go.documentation.sas.com/?docsetId=proc&docsetTarget=p0f0fjpjeuco4gn1ri963f683mi4.htm&docsetVersion=9.4&locale=ja ... View more

質問 SQLプロシジャのUPDATE句を利用して、値の更新を行なっています。 条件によって代入する値を変更したいのですが、このような場合に簡潔に記述する方法はないでしょうか。 やはり、UPDATEステートメントを複数行記述するしかないのでしょうか。 回答 指定される選択条件が比較的単純であれば、次のサンプルプログラムのような記述で、代入する値を変更できます。 より複雑な条件の場合は、数回に分けてUPDATE句を実行してください。   UPDATE句とCASE式の構文   UPDATE テーブル名 SET 列名 = CASE WHEN 条件_1 THEN 値_1 ... WHEN 条件_n THEN 値_n ELSE 値_m END WHERE 抽出条件 ;   プログラム例   DATA sample; INPUT a b; DATALINES; 1 0 2 0 3 0 4 0 5 0 ; PROC SQL ; UPDATE sample SET b = CASE WHEN a = 1 THEN 1 WHEN a in (2, 3) THEN a * 2 ELSE 0 END WHERE a <= 4 ; QUIT;   ... View more

質問 SASデータセット内の変数名や属性などの情報をDATAステップ中で参照したいのですが方法はありますか？ 回答 SASデータセットの変数名や属性についての情報は、ビュー「SASHELP.VCOLUMN」に保存されています。 このビューの持つ情報は次の通りで、「SASHELP.VMEMBER」にないフォーマットなどの情報も含まれています。   LIBNAME ライブラリ参照名 MEMNAME データセット名 MEMTYPE ライブラリメンバのタイプ（DATAあるいはVIEW） NAME 変数名 TYPE 変数のタイプ（CHAR　文字変数　NUM　数値変数） LENGTH 変数の長さ NPOS データセット入力バッファ中の相対位置 VARNUM 変数につけられた通し番号 LABEL 変数ラベル FORMAT フォーマット INFORMAT インフォーマット IDXUSAGE インデックスの使用     DATAステップ中で任意のデータセットの変数名や属性情報を参照するには、SASHELP.VCOLUMNから、 参照するデータセットに対応したオブザベーションの値を用いるという方法が考えられます。   SASHELPライブラリには、上記以外にも、各種情報を格納したSQLビューが格納されています。     ... View more

質問 サブセット化IFステートメントとWHEREステートメントは非常によく似ています。 どう使い分ければよいか教えてください。 回答 サブセット化IFステートメントとWHEREステートメントは、どちらも条件にもとづいてオブザベーションを選択するのに使えます。また構文もよく似ています。しかしこれらのステートメントは、次の点で異なっています。   オブザベーション選択のタイミング 記述できる条件式 使用できるステップ ●オブザベーション選択のタイミング オブザベーションを読み込む場合、PDV(プログラムデータベクトル)に値が読み込まれます。 WHEREステートメントはPDVに読み込まれる前に作用します。一方サブセット化IFステートメントはPDVに読み込まれた後に作用します。 この違いにより、次のような相違点が生じます。   外部ファイルをINFILE/INPUTで読み込み場合、WHEREステートメントは使用できません。 グループ処理を行うときに使用するFIRST.by変数およびLAST.by変数の値の挙動が、WHEREステートメントとサブセット化IFステートメントでは異なることがあります。   ●記述できる条件式 サブセット化IFステートメントでは、SAS式のルールにもとづく任意の表現が使えます。 一方WHEREステートメントは表現にDATAステップ関数は使えません。 しかしWHEREステートメント特有の便利な各種演算子(CONTAINS、BETWEEN-AND、LIKEなど)が使えます。   変数NAMEに"John"が含まれていてAGEが20以下のオブザベーションを選択するWHEREステートメントは、次のように記述できます。 where name contains "John" and age <= 20; サブセット化IFステートメントでは、同じ条件式を表すのにINDEX関数を使用します。   ●使用できるステップ サブセット化IFステートメントはDATAステップでしか使えません。 一方WHEREステートメントは、DATAステップだけでなくPROCステップでも使えます。 PROCステップではWHEREステートメントが使えるので、1つのPROCステップで簡単にデータのサブセットを処理できます。 これは、CPU時間およびディスク容量の節約になります。   変数SEXの値が"女子"のオブザベーションについて散布図を作るには、WHEREステートメントを使えばPROCステップだけで処理できます。 proc plot data=sashelp.class; plot weight * height; where sex="女子"; quit;   サブセット化IFステートメントを使う場合、DATAステップとPROCステップの2つのステップが必要です。 data female; set sashelp.class; if sex="女子"; run; proc plot data=female; plot weight * height; quit;   ... View more

質問 LAG関数をIFステートメントの中で使用している場合、結果が期待と異なる場合があります。なぜでしょうか。 回答 この現象は、LAG関数は、キューの操作、つまり値の保存を「実行されたときだけ行なう」という性質を持っているためです。LAG関数は、nオブザベーション前の値を返すのではなく、LAG関数を実行するたびに蓄積されたキューのn個前の値を返します。nオブザベーション前の値を取得するには、LAG関数を毎回実行して全ての値をキューに蓄積する必要があります。   プログラムを例に説明します。下記のプログラムは「前月との気温差をdif、偶数月のときにはその直前の偶数月(その月の2月前)との気温差をdif2」と出力することを目的として作ったプログラムです。   （例１） DATA lag ; INPUT chiku $ month kion ; dif = kion - LAG1(kion) ; IF MOD(month,2) = 0 THEN dif2 = kion - LAG2(kion) ; DATALINES; 北海道 1 -20 北海道 2 -18 北海道 3 0 北海道 4 10 北海道 5 20 北海道 6 15 北海道 7 20 北海道 8 28 北海道 9 22 北海道 10 17 北海道 11 10 北海道 12 2 ; RUN;   （例1で作成されたデータセットlagの内容） month kion dif dif2 (期待していたdif2の結果) 1 -20 . . 2 -18 2 . ← . = -18 - . 3 0 18 . 4 10 10 . ← 28 = 10 - (-18) 5 20 10 . 6 15 -5 33 ← 5 = 15 - 10 7 20 5 . 8 28 8 18 ← 13 = 28 - 15 9 22 -6 . 10 17 -5 2 ← -11 = 17 - 28 11 10 -7 . 12 2 -8 -26 ← -15 = 2 - 17     上記の結果からわかるように、dif2の結果は2月前のものではなく、「2回前にLAG関数により保存された値との計算結果」となっています。   month kion dif2 (実際のdif2の計算) 1 -20 . 2 -18 . = -18 - . 3 0 . 4 10 . = 10 - . (2回前のLAGは存在しない) 5 20 . 6 15 33 = 15 - (-18) (2回前のmonth=2で保存された値) 7 20 . 8 28 18 = 28 - 10 (2回前のmonth=4で保存された値) 9 22 . 10 17 2 = 17 - 15 (2回前のmonth=6で保存された値) 11 10 . 12 2 -26 = 2 - 28 (2回前のmonth=8で保存された値)   期待しているとおりの結果を出力させる場合は、次のようなプログラムを実行します。   （例2） DATA lag ; INPUT chiku $ month kion ; dif = kion - LAG1(kion) ; dif2 = kion - LAG2(kion) ; IF MOD(month,2) THEN dif2 = . ; /* 奇数月の場合は欠損値 */ DATALINES; 北海道 1 -20 北海道 2 -18 北海道 3 0 北海道 4 10 北海道 5 20 北海道 6 15 北海道 7 20 北海道 8 28 北海道 9 22 北海道 10 17 北海道 11 10 北海道 12 2 ; RUN;   （例2で作成されたデータセットlagの内容） month kion dif dif2 1 -20 . . 2 -18 2 . 3 0 18 . 4 10 10 28 5 20 10 . 6 15 -5 5 7 20 5 . 8 28 8 13 9 22 -6 . 10 17 -5 -11 11 10 -7 . 12 2 -8 -15   ... View more

質問 「2004/12/25」から「2005/01/25」、また「2004/12/02」から「2005/02/02」など、基準となる日付からnヵ月後の同日の日付を求めるにはどうしたらよいでしょうか。 回答 SAS日付を年月の単位で移動させる関数には、INTNX関数があります。 INTNX関数は、基準となる日付に対し、以下の構文で指定します。   1ヵ月後 INTNX('MONTH',基準日付,1); 2ヵ月後 INTNX('MONTH',基準日付,2);   ただしINTNX関数は、デフォルトではnヵ月後の月の初日を返します。 つまり、2004年11月26日を基準日付にした場合、以下のような値を返します。   1ヵ月後 2004年12月1日（に相当するSAS日付値） 2ヵ月後 2005年1月1日（に相当するSAS日付値）     したがって、nヵ月後の同日日付を求めるには少し工夫が必要になります。 以下にプログラム例と解説を記載しますので、参考にしてください。   /* サンプルデータ作成 */ DATA days; INPUT hizuke YYMMDD10.; CARDS; 2004/01/01 2000/02/29 2003/02/28 2004/12/02 2004/12/25 ; /* 1ヵ月後の日付計算 */ DATA temp; SET days; daydiff = DAY(hizuke)-1; monthafter = INTNX('MONTH', hizuke, 1) + daydiff; FORMAT hizuke YYMMDD10. monthafter YYMMDD10.; RUN;     上記のプログラムでは、与えられた基準日付 hizukeがその月の1日から何日後であるかを求め、daydiffに格納します。hizuke が2004年11月26日ならば、daydiffは25となります。   ここで、「1ヵ月後の日付 ＝ 1ヵ月後の月初 ＋ daydiff」と考えられますので、INTNX関数を使用して1ヵ月後の日付を求めています。 ただし上記のプログラムでは、うるう年でない年の1月29日の1ヶ月後の日付は、2月29日が存在しないために「3月1日」となりますのでご注意ください。   ... View more

質問 単一値代入法（single imputation）や多重代入法（multiple imputation）に基づいて、欠損値への代入（補完、埋め込み）を行なうには、どのようにしたらよいですか。 回答 多重代入法についてはMIプロシジャとMIANALYZEプロシジャが対応しています。これら2種のプロシジャは、SAS 8では評価版の扱いですが、SAS9では正規版となり、新たな機能も追加されています。  詳細については、SAS/STAT User's Guideをご参照ください。   一方、単一値代入法を包括的に扱うプロシジャはありませんが、代入方法に応じて以下のような対処策が考えられます。 ■ 全ての欠損値に同じ値（たとえば0）を代入する方法   DATAステップでMISSING関数を使用すれば可能です。下記のFAQをご参照ください。   • 「欠損値を任意の値に置き換える」   また、SAS/STATのSTDIZEプロシジャを利用することもできます。この場合、MISSING=オプションで指定した数値が欠損値に代入されます。   PROC STDIZE DATA=data1 REPONLY MISSING=0 OUT=out1; VAR _numeric_; RUN;   なお、STDIZEプロシジャでは中央値を代入することもできます。PROC STDIZEステートメントのオプションとしてMETHOD=MEDIANを指定してください。     ■ 各変数に指定した値を代入する方法 STDIZEプロシジャで対応可能です。代入する値を事前にデータセットとして用意し、オプションとして「METHOD=IN(  )」を指定してください。かっこの中には、用意したデータセットの名前を指定します。 以下の例は、変数xにおける欠損値に対して1を、変数yにおける欠損値に対して2を代入するものです。   /* 方法1 */ DATA impute; LENGTH _type_$ 8; INPUT _type_ x y; CARDS; LOCATION 1 2 SCALE 1 1 ; PROC STDIZE DATA=data3 METHOD=IN(impute) REPONLY OUT=out3; VAR x y; RUN; /* 方法2 */ DATA impute2; loc_x=1; loc_y=2; scal_x=1; scal_y=1; RUN; PROC STDIZE DATA=data3 METHOD=IN(impute2) REPONLY OUT=out4; VAR x y; LOCATION loc_x loc_y; SCALE scal_x scal_y; RUN;   ■ 線形回帰に基づく方法や、群ごとの平均を代入する方法 SAS/STATのREGプロシジャやGLMプロシジャでは、OUTPUTステートメントを利用すると、従属変数が欠損値である場合でもその予測値が算出されます。この機能を応用してください。     ■ LOCF（Last Observation Carried Forward）に基づく方法 LVCF（Last Value Carried Forward）とも呼ばれます。一般には、DATAステップを利用して対応することができます。下記のFAQのプログラムを参考にしてください。   「任意変数が欠損値のとき、前のオブザベーションの欠損値ではない値のセット方法」   また、米国で開かれているSASユーザー会（SUGI28）では、LOCFに関する下記の論文が発表されています。   「The DOW (not that DOW!!!) and the LOCF in Clinical Trials」（PDF）     ■ Hot Deck法、およびCold Deck法 一般には、DATAステップを用いて条件分けを行って対応することが考えられます。また、状況によっては何らかのプロシジャの併用が必要となる場合もあります。  ... View more

質問 レコード内に各レコード項目の識別子はあるが、レコードの長さが記載されたフィールドがない、つまり「改行コードのない」可変長のテキストファイルがあります。 このようなファイルをSASで読み込むには、どのようにすればよいでしょうか。 回答 改行コードが存在しない可変長データの代表的なケースには、メインフレームから転送されてきたデータがあります。   ホスト上にて作成されたデータであり、かつメインフレームからの転送時にコード変換がされておらず、制御コードが付加されたままであれば、INFILEステートメントでRECFM=オプションに S370V、またはS370VBを指定することで読み込むことが可能です。   しかし、一般的なファイル転送ソフトウェアなどで転送された場合、可変長レコードの長さを示すRDWやBDWと呼ばれるフィールドは省略されるため、RECFM=オプションは利用できません。   各レコードの長さが格納されたフィールドがない場合、データのレコードフォーマットをRECFM=N (unbufferd)として扱い、次に記載するような方法を利用することで読み込むことが可能です。   レコードの各フィールドを全てフォーマット指定して、フォーマット入力する レコードの各フィールドを、ポインタコントロールを利用してポインタ入力する （識別レコードの組み合わせごとに長さが固定であれば）レコードの各フィールドの読み込み位置に、レコードの長さを付与する   可変長レコードの読み込みについては、次のページで「ホストからのEBCDICデータを読み込む方法」として紹介しています。 TS-DOCS TS-642 Reading EBCDIC Files on ASCII Systems   なお、一例として簡単なサンプルを記載します。     前提条件 先頭１バイトのフラグで識別される長さの異なるデータがある フラグで識別された各オブザベーション（レコード）の 組み合わせの長さは固定 先頭1バイトが「1」の場合は200バイト 先頭1バイトが「2」であれば1000バイト  各オブザベーションには改行コードが存在しない データはパック十進データを除き、ASCII変換済み   プログラム概要 入力ファイル「VBREC1」をストリームレコード（1レコードのデータ）として読み込み、レコードレイアウトの異なる2つのデータを、それぞれ「REC1」「REC2」データセットとして分けて作成する。 次のデータレコードの読み込み位置を変数「strec」に保持し、該当のレコードの開始位置を特定させる。 /* 入力ファイル定義 */ FILENAME vbrec1 'C:¥test¥testvb.txt'; /* 出力データセット定義 */ DATA WORK.rec1 WORK.rec2; /* 先頭からの位置を記録した変数をRETAINステートメントで保存 */ RETAIN strec; /* １回目のみ0で初期化 */ IF _n_ = 1 THEN strec = 0 ; /* RECFM=n指定 */ INFILE vbrec1 RECFM=n ; /* １バイトの識別データの読み込み */ /* 行保持指示子"@"で行を保持 */ INPUT @1+strec w_flg 1. @; /* 識別変数の選択 */ SELECT(w_flg) ; /* 条件 =1 */ WHEN('1') DO; INPUT @2+strec w_num1 4.             @09+strec w_cdata1 $CHAR7. @16+strec w_cdata2 $CHAR10. @27+strec w_nopd1 S370FPD4. @32+strec w_nozd1 3. @35+strec w_alld1 $CHAR166. @; OUTPUT work.rec1; /* 条件のレコード長指定 */ endrec = 200; /* STRECへ次のレコード位置格納 */ strec = strec + endrec ; END; /* 条件 =2 */ WHEN('2') DO; INPUT @002+strec w_kdat1 $CHAR818. @820+strec w_cd1 4. w_cd2 3. @827+strec w_date YYMMDD10. @838+strec w_nopd1 S370FPD4. @843+strec w_nozd1 3. @846+strec w_alld2 $CHAR155. @; OUTPUT work.rec2 ; /* 条件のレコード長指定 */ endrec = 1000; /* STRECへ次のレコード位置格納 */ strec = strec + endrec ; END; OTHERWISE DELETE; END; /* DEBUG プログラム作成時のループ防止 */ IF _N_ > 5; STOP; RUN;         ... View more

質問 長いプログラムを実行する際、「ログ WINDOW FULL：ウィンドウのテキスト域が一杯のため削除します」というウィンドウが表示されて、処理が中断します。 これを回避するにはどうすればよいでしょうか。 回答 上記メッセージは、ログの出力行数が、ウィンドウの最大表示行数である99,999行を超えたことを示しています。 この現象を回避するには、下記の3種類の対応法があります。   1. PRINTTOプロシジャでログの出力先を外部ファイルに切り替える PRINTTOプロシジャを使用して、ログの出力先をログウィンドウから外部ファイルに切り替えることができます。 この場合は、行数の制限はありません。     例 /* ログの出力先を外部ファイルに変更 */ PROC PRINTTO LOG="C:\SASLog\print.log" NEW; RUN; DATA class; SET sashelp.class; RUN; /* ログの出力先をウィンドウに戻す */ PROC PRINTTO; RUN;     PRINTTOプロシジャの LOG=オプションに、出力先ファイルのパスまたはファイル参照名を記述します。以降、ログはすべてファイルに記録されます。   ログの出力先をウィンドウに戻すには、オプションなしでPRINTTOプロシジャを実行します。これにより、ログの出力先がデフォルト（ウィンドウ）に戻ります。   すでに同名のファイルが存在する場合、上書きするにはNEWオプションを指定します（NEWオプションを指定しない場合は、追加出力になります）。     2. NONOTES, NOSOURCEシステムオプションを指定して出力行数を減らす プログラムソースやNOTEメッセージの出力が不要な場合は、NONOTES, NOSOURCEシステムオプションを指定して、出力しないようにできます。   NOSOURCE ソースプログラムを出力しない （デフォルトはSOURCE） NONOTES 「NOTE:」で始まるメッセージを出力しない （デフォルトはNOTES）     例 OPTIONS NOSOURCE NONOTES;     デフォルトに戻すときは、「OPTIONS SOURCE NOTES;」を指定します。   3. DMSLOGSIZEシステムオプションを指定して、出力行数の制限を拡張する   SAS 9.1以降では、SASコンフィグレーションファイルに、 DMSLOGSIZEシステムオプションを指定することで、出力行数の制限を拡張できます。 以下の例では999,999行まで拡張しております。   例   -DMSLOGSIZE 999999     ... View more

質問 SASデータセットの内容を、既存のExcelファイルに新しいシートとして出力できますか。 また、既存のシートにデータを追加できますか。 回答 SAS/ACCESS Interface to PC Filesは、Excelエンジンをサポートしています。 このエンジンを使用して、既存のExcelファイル(xlsファイル)にシートを追加できます。   出力方法は、Excelファイルを新規に出力する場合と同様です。次のサンプルプログラムでは、既存のExcelファイル「c:¥DATA¥book1.xls」ファイルに対して、「excalss2」というシートを追加出力しています。   /* ライブラリ参照名の割り当て */ LIBNAME ex EXCEL 'c:¥DATA¥book1.xls'; /* データの出力 */ DATA ex.exclass2 ; SET sashelp.class; RUN; /* ライブラリ参照名のクリア */ LIBNAME ex CLEAR;     また、APPENDプロシジャを使用して、既存のExcelシートにデータを追加することも可能です。既存のExcelシートにデータを追加する際には、LIBNAMEステートメントで「SCANTEXT=NO」を指定する必要があります。   次のサンプルプログラムでは、「c:¥data¥book1.xls」ファイルに含まれるシート「exclass2」にデータを追加しています。   /* ライブラリの割り当て */ LIBNAME ex EXCEL 'c:¥data¥book1.xls' SCANTEXT=NO; /* APPENDプロシジャ */ PROC APPEND BASE=ex.exclass2 DATA=llibref.data FORCE; RUN; /* ライブラリ参照名のクリア */ LIBNAME ex CLEAR;     上記のサンプルプログラムを実行すると、「変数の長さが一致しない」という下記のWARNINGメッセージがログに出力されますが、データは追加されます。   WARNING: 変数 Name は BASE と DATA ファイル上で長さが一致しません。 (BASE 255 DATA 8).     Windows版SAS 8.2をご利用の場合は、弊社ホームページの下記FAQをご参照ください。   「Microsoft Excelへのデータの複数シート出力について」 http://www.sas.com/japan/service/technical/faq/list/body/pc054.html   ... View more

質問 複数のSASデータセットを、1つのMicrosoft Excelファイル（ブック）内で複数のシートに分けて出力する方法を教えてください。 回答 EXPORTプロシジャでは、同一ブック内の複数シートへの出力はできません。次の2通りの方法で出力する方法があります。 DDE の機能（Base SAS）を使用してExcelシートに書き込む方法 SAS/ACCESS to ODBCを用いてExcelシートに書き込む方法 具体的な説明とプログラム例をご紹介します。    ■ DDE の機能を使用してExcelシートに書き込む方法      DDEは、Dynamic Data Exchange の略で、Microsoft社が提供しているWindowsアプリケーション間でのデータ交換に広く    使用されているプロトコルです（DDEの詳細については、Microsoft社提供、市販されている書籍などでご確認ください）。    SASでは、このDDEプロトコルを利用して、MS-Excelにデータを出力することが可能です（※ MS-ExcelがSASと同じPCに    インストールされていることが必要です）。    以下のプログラム例では、2つのシートに出力していますが、この部分を適宜変更することにより、さらに複数のシートに    対して書き込むことができます。   OPTIONS NOXWAIT NOXSYNC; X "start excel"; /* MS Excel を起動 */ DATA _null_; x=SLEEP(10); RUN; /* Excelの起動を10秒待機 */ /* 待機時間は環境に応じて変えてください */ FILENAME ddedata DDE 'excel|sheet1!r1c1:r100c3'; /* DDEを利用したエクセル出力 */ /* シート１枚目 */ /* 乱数を使用して、サンプルを作成します */ DATA _null_; FILE ddedata; DO i = 1 TO 100; x = RANUNI(i); y = 10 + x; z = x - 10; PUT x y z; END; RUN; /* Excelへの書き込みを待機 */ DATA _null_; x=SLEEP(10); RUN; FILENAME ddedata2 DDE 'excel|sheet2!r1c1:r100c3'; /* シート2枚目 */ /* 乱数を使用して、サンプルを作成します */ DATA _null_; FILE ddedata2; DO i = 1 TO 100; x = RANUNI(i); y = 20 + x; z = x - 30; PUT x y z; END; RUN;   ■ SAS/ACCESS to ODBCを用いてExcelシートに書き込む方法   データソースとしてシステムまたはユーザーDSNを事前に定義します。 データソースの定義は、[ODBC データ ソース アドミニストレーター]を開き、データソースを定義します。なお、詳細の定義方法につきましては、Microsoft提供のヘルプ等をご参照ください。 [ODBC データ ソース アドミニストレーター]は、SASのビット数に合わせて、定義する必要があります。 一般的に、以下の実行ファイルより、[ODBC データ ソース アドミニストレーター]は、開くことができます。 32ビットの[ODBC データ ソース アドミニストレーター] C:\Windows\SysWOW64\odbcad32.exe 64ビットの[ODBC データ ソース アドミニストレーター] C:\WINDOWS\system32\odbcad32.exe   ライブラリを割り当てます。 libname mylib odbc dsn=Excel; /*データソース名をExcel とします。*/ /*上記データソースの設定で定義した名前を指定してください。*/   複数のシート別に出力します。下記の例では、sashelp.classのデータを利用しています。 DATA mylib.sheet1a; SET sashelp.class; RUN; /* 新しくsheet1a を作成 */ DATA mylib.sheet2a; SET sashelp.class; RUN; /* 新しくsheet2a を作成 */   SASを終了せずに作成したExcelファイルをExcelから利用するには、暗黙的に適用されたロックを解除するために、ライブラリの割り当てを開放する必要があります。 libname mylib clear;   なお、Excelへの複数シート出力に関するサンプルのマクロを、下記のURLで公開しています。あわせてご参照ください。 URL：http://support.sas.com/kb/6/695.html   ... View more

質問 計算結果やプロシジャからの出力値の桁数が異なる際、表示時に 同一桁数にならない場合があります。 その際、有効桁数を揃えて表示するには，どうしたらよいのでしょうか。 回答 ROUND関数で有効数字の桁を指定することにより、有効桁数を揃えることができます。 指定する有効桁数は、ABS関数、LOG10関数、およびINT関数を使用して以下のように求めることが可能です。   　　1. INT 関数により値が整数であるか小数であるかで処理を分ける。 　　2. LOG10関数により何桁の数値かを算出する 　　3. 2の値を用いて設定した有効桁となるように、累乗の値を調整する。 　　4. 3の値を用いてROUND関数により指定した有効桁での丸めを行なう。   以下の例は、前述の手順で有効桁数が3になるように桁数を揃えるプログラムです。    例：有効桁数を3桁で揃えるプログラム例    DATA one; INPUT x; IF int(x) ne 0 THEN DO; /* 整数の場合 */ _3sigdigit=ROUND(x,10**(INT(LOG10(ABS(x)))-2)); END; ELSE DO; /* 小数の場合 */ _3sigdigit=ROUND(x,10**(-1*(ABS(INT(LOG10(ABS(x))))+3))); END; DATALINES; 1234 12345 0.01234 ; RUN; PROC PRINT; FORMAT x _3sigdigit 12.5; RUN;   同様のサンプルプログラムが、以下のURLにも記載されています。 http://support.sas.com/kb/24/728.html   ... View more

質問 MERGEステートメントを使用して2つのデータセットを結合する際、両方のデータセットにキーがあるときだけデータを出力するようにできますか？   回答 両方のデータセットにキーが存在するときだけデータを出力するには、MERGEステートメントでIN=オプションを指定します。   「IN=変数名」はブーリアンフラグで、BYグループが替わるたびにリセットされます。 両方のデータセットにデータが存在するデータだけを出力する場合は、BYグループをチェックします。 下記の例を参考にしてください。   プログラム例：    DATA file1; INPUT var name $; CARDS; 100 john 200 joe 400 bill 600 bob ; RUN; PROC SORT data=file1; BY var; RUN; DATA file2; INFILE cards dsd truncover; INPUT var address $ 13.; CARDS; 100,34 Smith Road 200,67 Burt Ave 300,12 You St 400,45 Younge St 500,79 Wellington 600,23 Done Road ; RUN; PROC SORT data= file2; BY var; RUN; DATA three; MERGE file1 (in=a) file2 (in=b); BY var; IF a and b; RUN; PROC PRINT; RUN;     結果： Obs var name address 1 100 john 34 Smith Road 2 200 joe 67 Burt Ave 3 400 bill 45 Younge St 4 600 bob 23 Done Road     ... View more