SAS 9: 解析の新機能(reprise)
泉水（いずみ） 克之
SAS Institute Japan株式会社 技術本部

Copyright © 2004, SAS Institute Inc. All rights reserved.

Greetings, SAS 9 最新のリリースは SAS 9.1.3 です。

Copyright © 2004, SAS Institute Inc. All rights reserved.

2

Informations of SAS 9
 SAS 9.1 のドキュメントは下記のWEBサイトから。
http://support.sas.com/documentation/onlinedoc/sas9doc.html

Copyright © 2004, SAS Institute Inc. All rights reserved.

3

Informations of SAS 9 (continued)
 どなたでも、無料で登録なしにご覧になれます。
WEB版のオンラインドキュメント

Copyright © 2004, SAS Institute Inc. All rights reserved.

4

Informations of SAS 9 (continued)
 どなたでも、無料で登録なしにご覧になれます。
PDF版のオンラインドキュメント

Copyright © 2004, SAS Institute Inc. All rights reserved.

5

Informations of SAS 9 (continued)
 新機能については、「What’s New 」の項に記載されています。

Copyright © 2004, SAS Institute Inc. All rights reserved.

6

Shortcut to the SAS 9

Copyright © 2004, SAS Institute Inc. All rights reserved.

7

Power Analysis, Sample Size
Determination and PSS application
 SAS9.1から、POWER、及びGLMPOWERプロシジャがProduction となりました。

 これらの機能をGUI環境で利用するための、「Webブラウザアプリケーション」である「PSS Application」も開発されました。

Copyright © 2004, SAS Institute Inc. All rights reserved.

8

PROC POWER Supports:
 POWERプロシジャは、以下の検定における検出力の計算、サンプルサイズの決定を行うことができます。
• MULTREGステートメント− 回帰分析におけるTYPE III F検定• ONEWAYANOVAステートメント− 1元配置分散分析における対比に対する検定• ONECORRステートメント− Pearsonの相関係数に関する検定• ONESAMPLEFREQステートメント− 2項検定Copyright © 2004, SAS Institute Inc. All rights reserved.

9

PROC POWER Supports: (continued)• PAIREDFREQステートメント− McNemar検定• TWOSAMPLEFREQステートメント− 2×2表に対するPearsonのカイ2乗検定、尤度比検定、及びFisherの正確検定• ONESAMPLEMEANSステートメント− 1標本に対するt検定、同等性検定など• PAIREDMEANSステートメント− 対応のあるt検定など• TWOSAMPLEMEANSステートメント− 2標本に対するt検定、同等性検定など• TWOSAMPLESURVIVALステートメント− 生存時間解析におけるGehan検定（Wilcoxon検定）、LogRank検定、及びTarone-Ware検定

Copyright © 2004, SAS Institute Inc. All rights reserved.

10

PROC POWER: Syntax
 知りたい箇所を欠損値としてください。
/*1標本 t 検定における検出力を算出する例*/
proc power;
onesamplemeans
mean = 8 to 16 by 2
ntotal = 10 to 100 by 10
stddev = 20
power = . ;plot x=n max=100 min=5 ;run;
Copyright © 2004, SAS Institute Inc. All rights reserved.

11

PROC POWER: Graph Output

Copyright © 2004, SAS Institute Inc. All rights reserved.

12

PROC POWER: Syntax (2) 知りたい箇所を欠損値としてください。
/*1標本 t 検定におけるサンプルサイズを算出する例*/
proc power;
onesamplemeans
mean = 8 to 16 by 2
ntotal = .
stddev = 20
power = 0.6 to 0.9 by .05 ;plot ;run;
Copyright © 2004, SAS Institute Inc. All rights reserved.

13

PROC POWER: Output (2)Copyright © 2004, SAS Institute Inc. All rights reserved.

14

PROC GLMPOWER Supports:
 GLMPOWERプロシジャは、分散分析において、ユーザーが定義した対比にたいする検出力の計算、サンプルサイズの決定を行うことができます。

Copyright © 2004, SAS Institute Inc. All rights reserved.

15

PROC GLMPOWER: Syntax
 知りたい箇所を欠損値としてください。
/*2-way ANOVAおける検出力を算出する例*/
proc glmpower data=sample;
class x1 x2;
model y = x1 | x2;
powerntotal = 50
stddev = 10
power =. ;run;
Copyright © 2004, SAS Institute Inc. All rights reserved.

16

Power and Sample Size (PSS) Application
 PSS Application は、Webアプリケーションです。利用するためには、Web Serverが必要です。

 Stand-aloneの環境でも、付属の「Tomcat」を利用して、どなたでも使用可能です。

Copyright © 2004, SAS Institute Inc. All rights reserved.

17

PSS Application: Interface
 こんな感じです。

Copyright © 2004, SAS Institute Inc. All rights reserved.

18

ODS Graphics (experimental)
 分析の結果を表現することに特化した新しいグラフ出力機能です。

 Base SAS, SAS/STAT, SAS/ETSの多くのプロシジャで対応しています。
• CORR• GLM, LIFETEST, LOGISTIC, MIXED, REG, …
• ARIMA, AUTOREG, VARMAX, X12, …

Copyright © 2004, SAS Institute Inc. All rights reserved.

19

ODS Graphics (experimental) (cont’d) 分析結果とグラフを、以下のファイルとして出力できます。
• HTML• PDF, PS, PCL• RTF• LaTeX（ソースファイル）

Copyright © 2004, SAS Institute Inc. All rights reserved.

20

ODS Graphics: How to Specify?
ODS HTML;ODS GRAPHICS ON;←挿入/*分析 PROC。 オプションの指定を行う必要があるときも。 */ODS GRAPHICS OFF; ←挿入
ODS HTML CLOSE;

Copyright © 2004, SAS Institute Inc. All rights reserved.

21

ODS Graphics Template Language
 ODS Graphicsで出力されるグラフをカスタマイズするためには、TEMPLATEプロシジャを使用します。

proc template;define statgraph StatGraph.PredActual;layout Gridded;
layout Overlay / yaxisopts=( label="Predicted and Actual Values" );Band ylimitlower=eval (PREDICTED-TVAL*SQRT((1+LEVERAGE)*MSE))
ylimitupper=eval (PREDICTED+TVAL*SQRT((1+LEVERAGE)*MSE)) x=ID /fill=true lines=false fillcolor=StatGraphConfidence:foregroundlegendlabel="95% Confidence Limits" name="Conf";
SeriesPlot y=PREDICTED x=ID / linecolor=
StatGraphPredictionLines:contrastcolor linepattern=
StatGraphPredictionLines:linestyle legendlabel="Predicted" name="Predicted";
Scatter y=ACTUAL x=ID / markersize=GraphDataDefault:markersize
markersymbol=GraphDataDefault:markersymbol markercolor=
GraphDataDefault:contrastcolor legendlabel="Actual" name="Actual";EndLayout;
DiscreteLegend "Actual" "Predicted" "Conf" / border=true;
Copyright © 2004, SAS Institute Inc. All rights reserved.
EndLayout;

22

ODS Graphics in future SAS9.1ではexperimentalですが、次の大きなアップデートでは正規版となる予定です。また、もっと多くのプロシジャで利用できるようになります。

 テンプレートをカスタマイズするのは大変そう???
• ＯＤＳ Graphicsによるグラフ作成をサポートするGUIアプリ
ケーションが次の大きなアップデートで用意されるかもしれません。（現状では不確定要素が多いですが、検討しているようです。)

 正規版になったときには、SAS/GRAPHのライセンスが必要です。
• 使用後のご意見、ご要望をお寄せください。
Copyright © 2004, SAS Institute Inc. All rights reserved.

23

Generalized Linear Mixed Models(experimental)
 一般化線形混合モデルをサポートするGLIMMIXプロシジャが評価版として提供されています。

 米国SAS Institute Inc.のWEBサイトからダウンロード（EXEファイルとPDFドキュメント）する必要があります。

http://support.sas.com/rnd/app/da/glimmix.html

Copyright © 2004, SAS Institute Inc. All rights reserved.

24

Multiple Imputation
 正規版となりました。
 それにつれて、指定方法も変更（改善）されました。
 欠測パターンがmonotoneのとき• REGPMM• ロジスティックモデルに基づく補完（CLASSステートメントの利用、exeperimental）• 判別分析に基づく補完（CLASSステートメントの利用、exeperimental）

Copyright © 2004, SAS Institute Inc. All rights reserved.

25

Survival Analysis
 LIFETESTプロシジャ• 生存関数に対する検定手法が追加されました。（STRATAステートメント）− Tarone-Ware− Peto− Modified Peto
− Fleming-Harrington− トレンド検定との組み合わせ（全ての検定手法に対して）

 信頼帯(confidence band)の計算
Copyright © 2004, SAS Institute Inc. All rights reserved.

26

Survival Analysis (continued)
 TPHREGプロシジャ (experimental)• カテゴリ変数をモデルに• カテゴリ変数に関する対比の計算

 PHREGプロシジャ• WEIGHTステートメント• Proportional rate/mean モデル（再発事象データに対して）• 累積残差(cumulative residual)に基づいたモデルアセスメント(experimental)• etc
Copyright © 2004, SAS Institute Inc. All rights reserved.

27

Logistic Regression
 条件付きロジット回帰（STRATAステートメント）
 スコアリング（SCOREステートメント）
 EXACT な解析の処理速度を改善
 細々とした拡張• OFFSET=• OUTDESIGN=• カテゴリ変数の展開手法の拡張

Copyright © 2004, SAS Institute Inc. All rights reserved.

28

Other Enhancements
 DISTANCEプロシジャ• 様々な距離（類似度）を計算します。

 MIXEDプロシジャ• 診断統計量 (experimental)• 共分散構造の追加

 BOXPLOTプロシジャ• “What’s New”には記載されていない拡張点があります。

 ROBUSTREGプロシジャ• ロバストな手法を適用した「線形回帰」

Copyright © 2004, SAS Institute Inc. All rights reserved.

29

Other Enhancements (continued)
 FREQプロシジャ• ZEROSオプションで、度数が0のセルも分析に含めることができるようになりました。（2項検定やカッパ統計量の計算時に効果を発揮） 調査データの解析• SURVEYLOGISTIC, SURVEYFREQプロシジャが追加されました。
• 旧来のSURVEYREG, SURVEYMEANS,SURVEYSELECTプロシジャも強化されました。

 計算処理の並列化• REG, GLMプロシジャなどでは、CPUが複数存在するときに、一部の計算を並列処理できるようになりました。
Copyright © 2004, SAS Institute Inc. All rights reserved.

30

Random Number Generator------ Mersenne Twister 松本・西村(1998)による”Mersenne Twister”に基づく乱数列を生成する関数RANDが、SAS9.1ではproductionとなりました。

 旧来の乱数生成関数（RANUNI, RANNORなど）と比べて、様々な面で（特に大量の乱数を生成するとき）アドバンテージがあります。
• 生成速度 （のはず）• 周期が長い（2の19937乗 -1）• 生成される乱数列の「性質」が良い(623次元超立方体の中に均一に分布)

Copyright © 2004, SAS Institute Inc. All rights reserved.

31

Function compiler----- FCMP Procedure Base SASのプロシジャとして関数をコンパイルする
FCMPプロシジャが新たに用意されました。

 VBで関数を作成するときと同じような構文を使用します。
 コンパイルされた関数は、SAS/STATやSAS/ETS、SAS/ORの幾つかのプロシジャで利用可能です。（自動微分には対応していないようです。また、DATAステップ
では現在のところ使用不可。）

Copyright © 2004, SAS Institute Inc. All rights reserved.

32これで、おしまいです。
 本日はご参加いただきありがとうございました。

 ご質問？

Copyright © 2004, SAS Institute Inc. All rights reserved.

33

Copyright © 2004, SAS Institute Inc. All rights reserved.

34