options ls=78;
title "Partial Correlations - Wechsler Data";

 /* The first two lines define the name of the data set with the name 'wechsler'
  * and specify the path where the contents of the data set are read from.
  * Since we have a header row, the first observation begins on the 2nd row,
  * and the delimiter option is needed because columns are separated by commas.
  * The input statement is where we provide names for the variables in order 
  * of the columns in the data set. If any were categorical (not the case here), 
  * we would need to put a '$' character after its name.
  */

data wechsler;
  infile "D:\Statistics\STAT 505\data\wechsler.csv" firstobs=2 delimiter=',';
  input id info sim arith pict;
  run;

 /* glm stands for 'general linear model'
  * the model statement specifies info and sim
  * as responses and arith and pict as predictors
  * the 'nouni' option suppresses univariate stats
  * the 'manova' statement requests multivariate
  * statistics for info and sim jointly
  * the 'printe' option provides the sum of squares
  * and cross products matrix for error
  */

proc glm; 
  model info sim = arith pict / nouni;
  manova / printe;
  run;