/* Lesson 02-03_10 */
 /*    File Name = les0203_10.sas   01/10/11   */
/* options linesize=72 pagesize=20; */
options linesize=72 pagesize=30;

data gakusei;
  infile 'all10ae.csv'
         firstobs=2
         dlm=',' 
         truncover 
         missover
         dsd
;

  input sex $ shintyou taijyuu kyoui 
        jitaku $ kodukai carryer $ tsuuwa;

if sex^='M' & sex^='F' then delete;
if shintyou=. | taijyuu=. then delete;
if taijyuu> 90            then delete;
if kodukai>=200000        then delete;

proc print data=gakusei(obs=10);
run;

proc means data=gakusei;
run;

proc chart data=gakusei;
  hbar shintyou taijyuu kyoui kodukai tsuuwa;
run;

proc chart data=gakusei;
  hbar shintyou taijyuu kyoui kodukai tsuuwa/group=sex;
run;


proc corr data=gakusei;
run;

proc plot data=gakusei;
  plot shintyou *(shintyou taijyuu kyoui kodukai tsuuwa);
  plot taijyuu  *(shintyou taijyuu kyoui kodukai tsuuwa);
  plot kyoui    *(shintyou taijyuu kyoui kodukai tsuuwa);
  plot kodukai  *(shintyou taijyuu kyoui kodukai tsuuwa);
  plot tsuuwa   *(shintyou taijyuu kyoui kodukai tsuuwa);
run;

proc reg data=gakusei;
  model taijyuu=shintyou;
  output out=outreg1 predicted=pred1 residual=resid1;
run;

proc print data=outreg1(obs=15);
run;

proc plot data=outreg1;
  plot taijyuu*shintyou;
  plot taijyuu*pred1;
  plot resid1*(pred1 shintyou taijyuu)/vref=0;
run;

proc univariate data=outreg1 plot normal;
  var resid1;
run;