/* Lesson 10-02b */
 /*    File Name = les1002b.sas  12/21/21   */

options nocenter linesize=78 pagesize=30;
options locale='en_US';
proc printto print = 'StatM21/les1002b-Results.txt' new;

data pref;
  infile 'StatM21/Prefectures20.csv'
    firstobs=8 dlm=',' dsd missover 
    encoding=sjis termstr=crlf;
/*
  input PName $ CName : $12. 
        population area density;
*/
  input PName : $20. CName : $20. 
        population area density;

lPName=length(PName);
lCName=length(CName);

title '*** proc print の出力 ***';
proc print data=pref(obs=100);
run;

title '*** proc means の出力 ***';
proc means data=pref;
run;

title '*** proc chart の出力 ***';
proc chart data=pref;
  hbar population area density;
  vbar population area density;
  hbar density / midpoints=0 to 6500 by 500;
  vbar density / midpoints=0 to 6500 by 500;
  hbar density / midpoints=0 to 6500 by 500 axis=0 to 25 by 5;
  vbar density / midpoints=0 to 6500 by 500 axis=0 to 25 by 5;
run;

/* キャラクタグラフ(文字でグラフを描画)が欲しい場合 */
ods listing;
ods graphics off;
title '*** proc univariate の出力(その1) ***';
proc univariate data=pref plot;
  var population area density;
run;
ods graphics on;
ods listing close;

/* グラフィック出力(キャラクタグラフではない)が欲しい場合 */
title '*** proc univariate の出力(その2) ***';
proc univariate data=pref plot;
  var population area density;
run;