/* ----------------------------------------------------------------
	FILE: logist.sas

	DATA: fabricated data on schizophrenia

	PURPOSE: Logistic regression and discrininamt function analysis
   ---------------------------------------------------------------- */

OPTIONS NODATA NOCENTER PAGENO=1 LINESIZE=80;

*- formats for printing;
proc format;
   value schiz 0='Normal'
               1='Schizophrenic';
run;

* - read in the data;
title 'Example of Logistic Regression and Discriminant Analysis';
data logist;
   infile '~carey/p7291dir/logist.dat';
   input schizo test1-test6;
   format schizo schiz.;
DATALINES;
 1 45 54 50 53 28 44
 1 33 38 55 58 31 37
 1 47 45 45 50 27 30
 1 55 33 43 49 47 45
 0 62 50 64 61 65 61
 0 53 49 58 52 71 56
 1 56 61 57 45 75 63
 1 48 43 52 57 54 50
 0 59 47 48 47 61 51
 1 26 46 31 37 28 42
 0 53 50 47 30 56 48
 0 52 74 71 62 53 57
 0 61 48 53 54 55 52
 0 54 57 60 73 42 60
 1 55 43 56 48 54 46
 0 58 57 40 52 54 69
 1 28 58 38 49 39 45
 0 58 57 58 50 61 70
 0 62 49 80 68 66 51
 1 53 44 45 37 40 53
 1 52 54 51 67 35 46
 0 62 39 50 40 53 51
 1 44 42 44 50 36 35
 0 64 54 48 46 66 63
 1 59 57 41 56 43 51
 1 48 65 47 62 45 55
 0 56 45 40 35 59 51
 0 58 69 74 71 69 65
 0 57 52 50 46 50 61
 1 41 57 38 46 41 40
 1 34 43 38 42 30 44
 1 38 44 56 46 57 41
 0 43 54 44 35 54 50
 1 59 66 61 53 62 60
 0 65 54 49 38 57 52
 1 46 40 43 60 41 39
 0 57 64 52 64 57 63
 0 51 40 63 58 42 48
 1 45 43 38 42 44 50
 1 45 50 62 61 45 50
 0 56 73 39 37 63 48
 0 64 70 47 54 67 59
 1 46 42 55 50 39 33
 0 61 53 55 58 69 65
 0 62 61 53 51 66 56
 1 35 44 58 53 50 44
 0 48 60 50 53 58 44
 0 52 32 34 37 40 52
 1 42 43 34 40 39 27
 1 36 49 58 53 38 35
 0 44 64 48 45 34 48
 1 42 52 41 41 50 49
 1 56 46 52 64 45 47
 1 41 39 40 44 38 35
 0 62 63 58 56 70 63
 0 49 73 45 39 55 41
 1 50 50 48 54 40 49
 1 42 49 49 53 41 51
 0 41 49 72 63 43 43
 0 52 46 49 41 53 36
 0 72 60 61 52 55 60
 1 26 46 44 42 45 43
 0 61 56 65 61 77 63
 0 51 47 60 59 48 51
 0 50 40 48 40 47 42
 0 64 60 65 60 62 68
 1 57 50 47 37 51 59
 1 38 47 44 33 48 52
 0 66 38 61 60 62 50
 1 43 37 37 44 43 50
 1 45 61 56 54 59 54
 1 50 37 49 42 52 48
 0 45 40 52 45 56 59
 0 68 66 82 78 67 64
 0 59 53 62 51 57 68
 1 52 54 38 38 50 52
 0 50 46 69 62 53 60
 1 41 50 34 37 50 44
 1 33 33 58 51 46 44
 0 52 67 66 64 65 59
 0 61 69 65 76 68 88
 1 33 46 24 22 38 31
 1 37 42 39 39 41 29
 1 36 26 60 55 67 37
 1 39 36 53 56 37 36
 1 51 32 52 49 48 32
 0 64 77 70 78 67 58
 0 54 49 61 59 52 57
 0 48 66 66 66 60 43
 1 40 67 43 52 45 65
 1 43 46 54 41 50 38
 1 53 53 52 56 52 63
 1 38 38 42 46 52 55
 1 34 58 43 48 52 46
 0 52 42 64 62 59 55
 0 40 37 41 37 42 30
 1 46 35 47 48 42 34
 0 50 52 51 47 54 54
 0 58 54 39 45 50 56
 0 50 63 49 53 54 40
run;


title2 'ANOVAs';
proc glm DATA=logist;
  classes schizo;
  model test1-test6=schizo;
  means schizo;
  manova h=schizo;
run;

title2 'Logistic Regression';
proc logist;
  model schizo=test1-test6;
run;

title2 'Classic Discriminant Analysis';
proc discrim pool=test pcorr manova;
   priors proportional;
   class schizo;
   var test1-test6;
run;

title2 'Canonical discriminant analysis';
proc candisc;
   class schizo;
   var test1-test6;
run;