SAS（Statistical Analysis System）とは

Size: px

Start display at page:

Download "SAS（Statistical Analysis System）とは"

えりかさかわ
7 years ago
Views:

1 SAS 講習会資料 SAS(Statistical Analysis System)とは 1) 統計解析パッケージの1つ他に SPSS, BMDP など 2) 統計手法が豊富で先進的医薬系で多く使われている 3)データハンドリングが柔軟で容易 4)グラフ機能が優れている 5) 記述言語が簡潔で容易 Page 1/11 Last Edit: 吉永実行環境 : 記述方式とメニュー方式 1) 記述方式データの編集や統計計算の指示を記述する 2)メニュー方式メニュー/マウス操作を基本にデータ解析を行う(SAS Enterprise Guide) [ 記述方式 ] 3つの Window エディタ, ログ, アウトプット, 他に GRAPH 1)エディタデータの編集や統計計算の指示を記述する 2)ログ計算を実行した際の履歴 (エラーやデータ入力状況を確認する) 3)アウトプット計算結果 SAS プログラムの構成 1)2 種類のステップ DATA ステップ: データの入力加工等行う PROC フロシーシャ: 主に統計処理 2)ステートメント: セミコロン()までの文字列自由形式 (コマンドや変数名の行またがり不可 ) 大文字小文字の同一視 ( 文字データを除く) 3)データの型 : 文字型, 数値型, 他欠損値 : ピリオド(.) 4) 変数名 : 英字あるいはアンタースコア(_)で始まり 32 文字以内の英数アンタースコア * 印刷ページの設定 OPTIONS PAGENO=1 PAGESIZE=60 LINESIZE=100 /*データの読み込み BMIの計算 */ DATA taikaku INFILE 'd:\ex\ex1.data' LENGTH name $ 20 INPUT name $ sex $ height weight BMI=weight / height **2 PROC SORT DATA=taikaku BY sex PROC MEANS DATA=taikkau BY sex PROC UNIVARIATE DATA=taikkau NORMAL PLOT VAR weight 5)SAS データセット: 主に DATA ステップで作成される PROC ステップで参照される SAS 専用形式 1 度作成されると SAS を終了するまで有効 6)コメントの記述 (2 通り) a)コメントステートメント * b)コメントブロック /* */ 中にセミコロンが含まれてもよいステートメントの中にあってもよい

エディタ, ログ, アウトプット, 他に GRAPH 1)エディタデータの編集や統計計算の指示を記述する 2)ログ計算を実行した際の履歴 (エラーやデータ入力状況を確認する) 3)アウトプット計算結果 SAS プログラムの構成 1)2 種類のステップ DATA ステップ: データの入力加工等行う

2 Page 2/11 Ⅰ.DATAステップ 1)SASデータセットの作成 : データの入力, 新しい変数の作成,データ値の変換例. DATA ensyu1 ensyu1 というSASデータセットを作成する 2)データの入力 : CARDS, データファイル a)cards DATA ステップに続いてデータを記述するデータが少ないときに有効 b)データファイルデータファイルから読込 3) 入出力書式 : リスト入力,カラム入力, 他 a)リスト入力境界値 ( 標準はスペース)により順次割り当てる b)カラム入力桁位置指定により割り当てる演習 1: CARDS によるリスト入力 ( 境界値はスペース) 各データ項目は1 個以上のスペースで区切られている先頭のスペースは無視される INFILE CARDSは省略可能縦にデータを整列する必要はない DATA ensyu1 INFILE CARDS INPUT no height weight CARDS 演習 2: CARDS による CSV 形式 ( 境界値がカンマ(,))のリスト入力 (Excel の CSV 形式出力 ) 各データ項目は 1 個の境界値 (,)で区切られている縦にデータを整列する必要はない補足 : 一般に境界値は任意の文字で下記のように境界値を指定すればよい INFILE CARDS DELIMITER=',' 注意 :この場合項目間を1 個以上の区切りとして扱う演習 1では DELIMITER=' 'が省略されている DATA ensyu2 INFILE CARDS DSD INPUT no height weight CARDS 1, 160, 60 2, 165, 63 3, 168, 65 演習 3: CARDS によるカラム入力データの項目が縦に整列している場合に用いる桁範囲を指定して読込む( 右例 ) 桁位置 (@x) 桁数 (n.)による読込方法 height 3. weight no 5. 補足 : 入力桁位置を指定すれば入力変数の並びは任意 2 種類の指定を混合してもよい演習 4: データファイルからのリスト入力ファイルからデータを読込む(CSV 形式 ) 他は CARDS 入力と同様 DATA ensyu3 INFILE CARDS INPUT height 8-10 weight no 1-5 /* */ CARDS DATA ensyu4 INFILE 'H:\SAS 講習会 \ensyu.csv' DSD INPUT no height weight

はスペース)により順次割り当てる b)カラム入力桁位置指定により割り当てる演習 1: CARDS によるリスト入力 ( 境界値はスペース) 各データ項目は1 個以上のスペースで区切られている先頭のスペースは無視される INFILE CARDSは省略可能縦にデータを整列する必要はない DATA

3 補足 a) 読込行指定 FIRSTOBS: 読込開始行 OBS: 読込最終行例. INFILE FIRSTOBS=100 OBS= 行目から 200 行目まで読込む Page 3/11 b) 文字データの読込 INPUT 文で読込む際に変数名の後ろにドル記号 ($)を付ける以後ではドル記号は必要ない c)リスト入力の場合の行末の欠損値扱い入力項目に対するデータがない場合次の行を読みにいく INFILE CARDS MISSOVER とすると次の行へ行くのを禁止する後続を欠損値とする DATA ensyuaddb INFILE CARDS INPUT no sex$ CARDS 1 M 3 F d) 演算子 : 算術演算子 +( 加算 ),-( 減算 ),*( 乗算 ),/( 除算 ),**(べき乗 ) 文字結合 ( 結合 ) 比較演算子 =, <, >, <=, >=, ^=( 等しくない) 論理演算子 &(AND), (OR),^(NOT),~(NOT) (^と~は同じ) 4)プログラムの制御 : IF THEN ELSE a)if ( 条件 ) THEN ( 真の場合の処理 ) ELSE ( 偽の場合の処理 ) 注意 : 処理には1ステートメントしか書けない複数ステートメント書く場合には DO~END ブロックを用いる注意 :IF~THEN~ELSE DO~END などの対応を考えて分かりやすい構造に記述する b) 多岐選択 : SELECT WHEN END 注意 : 条件式を先頭から判定し該当した時点で SELECT ブロックを抜ける DATA ensyu4a INFILE CARDS INPUT no sex$ IF sex='m' THEN sexcode=1 ELSE sexcode=2 IF sex='m' THEN DO sexcode=1 sexcolor='blue' END ELSE DO sexcode=2 sexcolor='red' END CARDS 1 M 3 F DATA ensyu4b INFILE CARDS INPUT no height weight SELECT WHEN (150<=height<160) heigclas=5 WHEN (160<=height<170) heigclas=6 WHEN (170<=height<180) heigclas=7 WHEN (180<=height<190) heigclas=8 OTHERWISE END CARDS heigclas=. /*.は欠損値 */

するデータがない場合次の行を読みにいく INFILE CARDS MISSOVER とすると次の行へ行くのを禁止する後続を欠損値とする DATA ensyuaddb INFILE CARDS INPUT no sex$ CARDS 1 M 3 F d) 演算子 : 算術演算子 +( 加算 ),-( 減算 ),*( 乗算 ),/(

4 c) 反復 DO ステートメント(ARRAY と一緒によく用いる) ARRAY: 変数リスト(q1 q2 q5)を配列 (q(5))の要素に順次割り当てる DO 変数 = 初期値 TO 終値 BY 増分値 DO 変数 = 値リスト DO 変数 WHILE ( 条件 ) DO 変数 UNTIL ( 条件 ) 変数リストの簡略化 q1 q2 q3 q4 q5 q1-q5 と記述できる INPUT q1-q5 ARRAY q(5) q1-q5 DATA ensyu4c ARRAY q(5) q1 q2 q3 q4 q5 INFILE CARDS INPUT q1 q2 q3 q4 q5 qs=0 DO i=1 TO 3, 5 qs=qs+q(i) END CARDS Page 4/11 5) 不必要なレコードの削除 : DELETE IF ( 条件 ) THEN DELETE 6) 既存 SAS データセットの読込 (SET ステートメント) BMI= 体重 (kg)/ 身長 (m) 2 DATA ensyu5 INFILE 'H:\SAS 講習会 \ensyu.csv' DSD INPUT no height weight IF no=2 THEN DELETE DATA ensyu6 set ensyu1 BMI=weight/(height/100)**2 7)SAS データセットの結合 SET ステートメント,MERGE( 順次,マッチ) a) 複数のSASデータセットを連結する DATA ensyu7a 変数名にて連結する ( 縦方向に連結する) SET ensyu1 ensyu2 b) 複数のSASデータセットをキーにより結合する注意 : 入力 SASデータセットはキーにより並べ替えられていなければならない PROC SORT DATA=ensyu1 OUT=ensyu1 BY no PROC SORT DATA=ensyuAddB OUT=ensyuAddB BY no DATA all MERGE ensyu1 ensyuaddb BY no 並べ替え(ソート)の記述 PROC SORT DATA= 入力 SASデータセット OUT=SASデータセット BY キー変数 8)インポート/エクスポートによる他のファイルとの読書き例 : Excel ファイルの読書きエクセルファイル(ensyu.xls)に2つのワークシート PROC IMPORT DBMS=EXCEL OUT= ensyu8a DATAFILE= "H:\SAS 講習会 \ensyu.xls" REPLACE SHEET="ensyu" PROC SORT BY no PROC IMPORT DBMS=EXCEL OUT= ensyu8b DATAFILE= "H:\SAS 講習会 \ensyu.xls" REPLACE SHEET="sex" PROC SORT BY no DATA all MERGE ensyu8a ensyu8b BY no 注意 : 第 1 行目が変数名となる PROC EXPORT DBMS=EXCEL DATA=all OUTFILE= "H:\SAS 講習会 \ensyu.xls" REPLACE SHEET="Merge"

$4/11 5) 不必要なレコードの削除 : DELETE IF ( 条件 ) THEN DELETE 6) 既存 SAS データセットの読込 (SET ステートメント) BMI= 体重 (kg)/ 身長 (m) 2 DATA ensyu5 INFILE 'H:\SAS 講習会 \ensyu.$

5 9) 書式設定 a)カテゴリラベル(format,value) 各カテゴリにラベルをつけるカテゴリ値の代わりにラベルが表示される PROC FORMAT VALUE sexfmt 1='MALE' 2='FEMALE' DATA ensyu10 INFILE CARDS INPUT no sex FORMAT sex sexfmt. LABEL sex=' 性別 ' CARDS PROC FREQ TABLES SEX Page 5/11 b) 変数ラベル(LABEL) 変数にラベルをつける統計量などの表示時にラベルが表示される FREQ プロシジャ性別累積累積 sex 度数パーセント度数パーセント MALE FEMALE Ⅱ.PROCステップ一般形式 PROC プロシジャ名オプション付随ステートメント一般に用いるオプション DATA=SASデータセット最後に作成されてものと同じであれば省略できるよく用いる付随ステートメント WHERE ( 条件 ) BY 変数名対象オブザベーションを制限するカテゴリ別に処理する事前に SORT されている必要がある簡単なものをいくつか紹介する 1)データのプリント例印刷対象の指定 : WHERE ステートメント性別 : BY ステートメント印刷変数の指定 : VAR ステートメント PROC PRINT DATA=ensyu1 WHERE (no>2) BY sex VAR no height 2)データの並べ替え昇順 / 降順, 複数キーの指定も出来る出力が入力と同じ場合は OUT を省略できる PROC SORT DATA=ensyu1 OUT=ensyu2 BY sex 記述統計と予備的解析として単変量統計量 (p.96 のサンプルデータの IMPORT による読込 ) /* p.96の身長データファイル(エクセル形式 )*/ PROC IMPORT DBMS=EXCEL OUT= height DATAFILE= "H:\SAS 講習会 \height.xls" REPLACE SHEET="height"

---------------------------------------------------------- MALE 1 50.00 1 50.00 FEMALE 1 50.00 2 100.00 Ⅱ.

6 Page 6/11 1)グラフによる視覚的観察 (ヒストグラム:GCHART)(p.97) PROC CHART HVAR sex/type=freq PROC GCHART HBAR sex/type=freq キャラクター図グラフィカル図 2) 単純集計表 (FREQ)(p.103) 性別の単純集計 PROC FREQ TABLES sex 3) 要約統計量 a) MEANS 平均値, 標準偏差, 最小値, 最大値など基本統計量 (p.107~) sex 累積累積 sex 度数パーセント度数パーセント F M PROC MEANS VAR student MEANS プロシジャ分析変数 : student N 平均値標準偏差最小値最大値 b) UNIVARIATE MEANS に加え分布の検定などを行う(p.109~) PROC UNIVARIATE NORMAL の NORMAL は正規分布検定を行うオプション HISTGRAM student / NORMAL の NORMAL は正規分布グラフを重ねるオプション PROC UNIVARIATE NORMAL VAR student HISTGRAM student / NORMAL

107~) sex 累積累積 sex 度数パーセント度数パーセント ------------------------------------------------------- F 19 47.50 19 47.50 M 21 52.50 40 100.

7 Page 7/11 UNIVARIATE プロシジャ変数 : student (student) モーメント標本数 40 重み変数の合計 40 平均値合計 6615 標準偏差分散歪度尖度無修正平方和修正平方和変動係数平均の標準誤差基本統計量位置ばらつき平均値標準偏差中央値分散最頻値範囲四分位範囲正規性の検定検定 ---- 統計量 p 値 Shapiro-Wilk W Pr < W Kolmogorov-Smirnov D Pr > D Cramer-von Mises W-Sq Pr > W-Sq Anderson-Darling A-Sq Pr > A-Sq :: :: :: :: 2 変量間の関連 1) 相関係数 (CORR)(p.125~) Pearson の相関係数, N = 21 帰無仮説 Rho=0 に対する Prob > r student father mother student student father father mother mother PROC CORR PEARSON SPEARMAN WHERE sex='m' Spearman の相関係数, N = 21 帰無仮説 Rho=0 に対する Prob > r student father mother student student father father mother mother

00000 四分位範囲 13.50000 正規性の検定検定 ---- 統計量 ---- -------p 値 ------- Shapiro-Wilk W 0.940563 Pr < W 0.0361 Kolmogorov-Smirnov D 0.156253 Pr > D 0.0151 Cramer-von Mises W-Sq 0.

8 Page 8/11 統計的検定 ( 第 6 章 ) 2 群間の平均値の差の検定 (p.142~) 1) 対応なし:TTEST 等分散性のチェック:F-test PROC TTEST CLASS sex VAR student The TTEST Procedure Statistics Lower CL Upper CL Lower CL Upper CL Variable sex N Mean Mean Mean Std Dev Std Dev Std Dev Std Err Minimum Maximum student F student M student Diff (1-2) T-Tests Variable Method Variances DF t Value Pr > t student Pooled Equal <.0001 student Satterthwaite Unequal <.0001 Equality of Variances Variable Method Num DF Den DF F Value Pr > F student Folded F 等分散の場合は上段, 不等分散の場合は Satterthwaite の補正検定 2) 対応あり a)paired Ttest PROC TTEST WHERE sex='m' PAIRED student*father The TTEST Procedure Statistics Lower CL Upper CL Lower CL Upper CL Difference N Mean Mean Mean Std Dev Std Dev Std Dev Std Err Minimum Maximum student - father T-Tests Difference DF t Value Pr > t student - father <.0001 b) : 各差に対して平均値 =0 の検定 :UNIVARIAVTE ( 省略 )

Std Err Minimum Maximum student F 19 157.16 158.53 159.89 2.1427 2.8357 4.1934 0.6505 153 164 student M 21 169.31 171.57 173.84 3.8067 4.9757 7.1852 1.0858 162 181 student Diff (1-2) -15.68-13.05-10.

9 独立性の検定, 一致性の検定 (p.149~) 1)χ 2 検定 ( 最小期待値が5 以上 ) Fisher の直接確率 ( 最小期待値が 5 未満 ) FREQ プロシジャ表 : treat * resp treat resp 度数パーセント行のパーセント列のパーセント + - 合計 control treat 合計 Page 9/11 DATA frq INPUT treat$ resp$ w CARDS control - 13 control + 7 treat - 8 treat + 12 PROC FREQ WEIGHT w TABLES treat * resp / CHISQ treat と resp の統計量統計量自由度値 p 値 χ 2 乗値尤度比 χ 2 乗値連続性補正 χ 2 乗値 Mantel-Haenszel のχ 2 乗値 φ 係数一致係数 Cramer の V 統計量 Fisher の正確検定セル (1,1) 度数 (F) 7 左側 Pr <= F 右側 Pr >= F 表の確率 (P) 両側 Pr <= P サンプルサイズ = 40

20 17.50 32.50 50.00 35.00 65.00 36.84 61.90 --------------+--------+--------+ treat 12 8 20 30.00 20.00 50.00 60.00 40.00 63.16 38.10 --------------+--------+--------+ 合計 19 21 40 47.50 52.50 100.

10 SAS プロダクトの構成 (SAS ヘルプ) SAS ヘルプメイン SAS/ BASE SAS/ STAT( 統計 ) Page 10/11

11 Page 11/11 メニュー/マウス操作によるデータ解析 [メニュー方式 ] SAS Enterprise Guide(SAS EG と記す) Statview や SPSS のようにメニュー形式で指示し統計計算を行うインストール(SAS とは別にインストールする) 1)インストールする前にコンピュータが最低限のシステム必要条件を満たしているか確認 System Requirements Wizard for SAS Enterprise Guide(CD)により事前設定 2)SAS Enterprise Guideのインストール SAS Enterprise Guide(CD)によりインストール体験してみよう 1)Enterprise Guide の起動スタートすべてのプログラム SAS Enterprise Guide xx 2) 新規プロジェクトの作成 3)エクセルデータを読込む(さきの height ファイルを読込む) ファイルデータのインポートローカルコンピュータ 4) 要約統計量の計算男女別に student の統計量を計算メニューから要約統計量データの新規作成 (SAS EG でデータを新規作成し保存する場合 ) SAS EG 上でデータの新規作成を行うとデータは一時ファイル扱いとなり保存されない保存する場合は SAS 形式にエクスポートして保存する 1) 新規作成 :ファイル新規作成データ SAS の作業フォルダ(WORK)に作成する変数名, 変数ラベル,データ形式等を設定する 2) データのエクスポート:ファイルデータのエクスポート 3) エクスポートしたデータと置換 :フロー上のデータを右クリックプロパティファイルを変更

SAS Enterprise Guide xx 2) 新規プロジェクトの作成 3)エクセルデータを読込む(さきの height ファイルを読込む) ファイルデータのインポートローカルコンピュータ 4) 要約統計量の計算男女別に student の統計量を計算メニューから要約統計量データの新規作成 (SAS EG

Microsoft Word - 新Excel&SAS資料.doc

Microsoft Word - 新Excel&SAS資料.doc 1 SAS と EXCEL による統計処理 1. はじめにアンケート調査などで集めたデータを加工解析するためには統計ソフトの利用が必須である本演習では EXCEL と SAS を使ったデータ解析の方法について説明する 2. 前準備 EXCEL EXCEL での統計処理は関数を用いる方法と分