data rnafup ;
infile 'c:\temp\175rna_fup.dat' ;
input
@1 pidnum 6.
@7 specdate date7.
@15 assaydt date7.
@23 rna ;
timefrz = assaydt - specdate ;
format specdate assaydt mmddyy8. ;
run ;

proc print data = rnafup (obs=12 );
run ;

proc sort data = rnafup ;
by pidnum specdate ; 
run ;

data recent ;
set rnafup ;
by pidnum ;
if first.pidnum ;
run ;

proc print data = recent (obs=12) ;
run ;

proc means data = rnafup noprint nway ;
class pidnum ;
var rna timefrz ;
output out = rnasum mean = m_rna m_frz   stddev = sd_rna sd_frz;
run ;

proc print data = rnasum (obs = 20 ) ;
title 'rnasum dataset' ;
run ;

data newrnasum ;
set rnasum ;
count = _freq_ ;
drop _type_ _freq_ ;
run ;

proc print data = newrnasum (obs=12) ;
run ;

proc freq data = newrnasum ;
tables count ;
run ;




data scores ;
input id m1 m2 m3 m4 f ;
datalines ;
1 100 100 100 100 100
2 95   95  86  92  97
;
run ;

proc print data = scores ;
run ;



data reformscore ;
set scores ;
score = m1 ;
test = 'm1' ;
output ;
score = m2 ;
test = 'm2' ;
output ;
score = m3 ;
test = 'm3' ;
output ;
score = m4 ;
test = 'm4' ;
output ;
score = f ;
test = 'final' ;
output ;
drop m1 m2 m3 m4 f ;
run ;

proc print data = reformscore ;
run ;