* Record output to a Word file;

ODS RTF FILE = 'C:\Documents and Settings\Rich\Desktop\Rich\LYNAM\Nov13WorkshopOutput.rtf';


* Read data into SAS;

PROC IMPORT DATAFILE = 'C:\Users\Rich\Desktop\Rich\LYNAM\PROCTRAJData.xls'
     OUT = One   DBMS = XLS REPLACE;
* Supply path and name of Excel file with data, supply new name for data in SAS, change XLS to XLSX if using later version of Excel;
     SHEET = Sheet1;   
* Identify tab of Excel file on which data reside;
     GETNAMES = YES;
RUN;


* List all variables;

PROC CONTENTS DATA = One;
RUN;


* Print values of selected variables;

PROC PRINT DATA = One;
VAR ID Y1 Y2 Y3 Y4 Y5 X;
RUN;


* Print descriptive statistics for selected variables;

PROC UNIVARIATE DATA = One;
VAR Y1 Y2 Y3 Y4 Y5 X;
RUN;
 

* Add time indices to data set;

DATA One;
SET  One;
T1 = -3;
T2 = -2;
T3 = -1;
T4 =  0;
T5 =  1;
RUN;


* First version of trajectory analysis;

PROC TRAJ DATA=One OUTPLOT=OP OUTSTAT=OS OUT=OF OUTEST=OE;  
ID ID;						    		
VAR Y1 Y2 Y3 Y4 Y5; INDEP T1 T2 T3 T4 T5; MODEL ZIP; NGROUPS 3; ORDER 1 1 1;
RISK X;
RUN;

%TRAJPLOT(OP,OS,'Marijuana vs. Time','ZIP Model #1','Marijuana','Time')


* Examine numerical output;

PROC PRINT DATA=OF;
RUN;

PROC PRINT DATA=OE;
RUN;

PROC PRINT DATA=OP;
RUN;

PROC PRINT DATA=OS;
RUN;


* Second version of trajectory analysis;

PROC TRAJ DATA=One OUTPLOT=OP OUTSTAT=OS OUT=OF OUTEST=OE;  
ID ID;						    		
VAR Y1 Y2 Y3 Y4 Y5; INDEP T1 T2 T3 T4 T5; MODEL ZIP; NGROUPS 3; ORDER 2 2 2;
RISK X;
RUN;

%TRAJPLOT(OP,OS,'Marijuana vs. Time','ZIP Model #2','Marijuana','Time')



* Examine numerical output;

PROC PRINT DATA=OF;
RUN;

PROC PRINT DATA=OE;
RUN;

PROC PRINT DATA=OP;
RUN;

PROC PRINT DATA=OS;
RUN;


* Third version of trajectory analysis;

PROC TRAJ DATA=One OUTPLOT=OP OUTSTAT=OS OUT=OF OUTEST=OE;  
ID ID;						    		
VAR Y1 Y2 Y3 Y4 Y5; INDEP T1 T2 T3 T4 T5; MODEL ZIP; NGROUPS 3; ORDER 2 2 2; IORDER 0 0 0;
RISK X;
RUN;

%TRAJPLOT(OP,OS,'Marijuana vs. Time','ZIP Model #3','Marijuana','Time')


* Examine numerical output;

PROC PRINT DATA=OF;
RUN;

PROC PRINT DATA=OE;
RUN;

PROC PRINT DATA=OP;
RUN;

PROC PRINT DATA=OS;
RUN;


* Cease recording output to a Word file;

ODS RTF CLOSE;
RUN;