(2)SASは簡単プログラム、超高速計算

 簡単なSASプログラムの例を下に示してみましょう。

国語、数学、理科、社会、英語、体育です。試験得点について、289人分のデータがあって、その平均点、標準偏差、最高点、最低点、各テスト間の相関係数を求めるプログラムの例です。例えば、プログラムファイル名は TEST.SAS とつけたとします。

data test;
infile "test.dat";            ←データをtest.datという
ファイルから読み込む
input kokugo suugaku rika shakai eigo taiiku;
              ↑国、算、理、社、英、体の順にデータがなっている。
proc corr;  ←平均、相関などを求めるプログラムの指示
run;

たった、この5行のプログラムでこれだけのことが、簡単に、しかも高速にできてしまいます。

 データはこの場合でいえば、TEST.DATというファイル名で、次のように1行に科目数分のデータを、スペースを挟んで並べればよいのです。

データの例

45 55 80 90 60 50
48 50 80 60 80 20
42 50 60 80 55 33
** ** ** ** ** **

45 60 80 95 50 60

以上のようなプログラムファイルと、データファイルを作って、SASワークウインドウの中のアイコンをクリックすれば、計算できてしまいます。 次の図は、上のプログラムを実行した例です。

                                        SAS システム                   2001年 *月24日 土曜日 23:37   2

                                                CORR プロシジャ

                            5  変数 :    kokugo   suugaku  rika     shakai   eigo

                                                   要約統計量

           変数               N        平均値      標準偏差          合計        最小値        最大値

           kokugo            44      58.22727      23.46186          2562      21.00000      97.00000

           suugaku           44      35.95455      25.03482          1582             0      80.00000

           rika              44      66.81818      19.04374          2940      31.00000      96.00000

           shakai            44      55.09091      28.31946          2424       2.00000     100.00000

           eigo              44      40.95455      28.91524          1802             0      99.00000

                                          Pearson の 相関係数, N = 44

                                      帰無仮説 Rho=0 に対する Prob > |r|

                                kokugo       suugaku          rika        shakai          eigo

                 kokugo        1.00000      -0.12280      -0.08449       0.22247      -0.07773

                                              0.4271        0.5856        0.1466        0.6160

                 suugaku      -0.12280       1.00000       0.13515      -0.16850       0.19452

                                0.4271                      0.3817        0.2742        0.2058

                 rika         -0.08449       0.13515       1.00000       0.05143      -0.27297

                                0.5856        0.3817                      0.7402        0.0730

                 shakai        0.22247      -0.16850       0.05143       1.00000      -0.08247

                                0.1466        0.2742        0.7402                      0.5946

                 eigo         -0.07773       0.19452      -0.27297      -0.08247       1.00000

                                0.6160        0.2058        0.0730        0.5946

 どうですか?このように簡単に平均、標準偏差、相関行列などができてしまいます。しかも計算時間は、数秒です。教育に、ゼミに、研究にと大いに活用してください。もちろん、主な統計解析のほとんどがこのパッケージには入っています。

(3)複雑な因子分析も簡単、高速、いろいろ試せる

 こんどは、因子分析という多変量解析のプログラム例です。ご存じの方は多いと思いますが、これをFORTRANなどのプログラム言語でプログラムしますと、単純なものでも数百行以上のプログラムになります。ところが下のように、SASではデータ部分を除いた計算の命令の本体は数行ですんでしまいます。

DATA PSYCH9 (TYPE=CORR);

 DATA X1-X9 

  CARDS;

   1.00   .     .     .     .     .     .     .     .  ← データの始まり

    .75  1.00   .     .     .     .     .     .     .  (相関係数)

    .78   .72  1.00   .     .     .     .     .     .

    .44   .52   .47  1.00   .     .     .     .     .

    .45   .53   .48   .82  1.00   .     .     .     .

    .51   .58   .54   .82   .74  1.00   .     .     .

    .21   .23   .28   .33   .37   .35  1.00   .     .

    .30   .32   .37   .33   .36   .38   .45  1.00   .

    .31   .30   .37   .31   .36   .38   .52   .67  1.00

;                         ← データの終わり

PROC FACTOR METHOD=ML NFACT=3 ALL MSA SCREE PREPLOT ←因子分析の

     OUTSTAT=FACT;                                         第1回目

RUN;

PROC FACTOR ROTATE=PROMAX NORM=WEIGHT REORDER ROUND FUZZ=.1 PLOT

     OUTSTAT=FACT;                 ←↑因子分析の

RUN;                                                        第2回目

次のページへ