Course 960-542 Home Page
Instructor
John E. Kolassa,
kolassa@stat.rutgers.edu
Announcements
- Office hours for 11 Apr moved to 3:30-4:20.
- The data set for
hw 4.
- Hint for HW3: For question 1 part c, you need to test whether
a factor with five levels has an effect, stratifying on stage.
We've done something like this in class, when the factor only has two
levels. In this case, you use
proc lifetest with
strata and an indicator
for one of the factors used in the test statement. What I've asked you
to do is a bit harder. In this case, create dummy variables representing
four of the five levels of the factor of interest (in this case, cell type).
Now use all four variables after the test statement. That is, if
ct is cell type, coded as zero through 4, create
ct1 taking the values 1 for cell type 1 and 0 otherwise, and
similarly construct ct2, ct3, and ct4
for types 2 through 4 respectively, and use test ct1 ct2 ct3
ct4; in conjunction with strata.
- Also, I didn't specify alpha for problem 2. Use alpha=.025.
- Office hours for Mon 28 Feb canceled.
- The format of the
data file for HW3 appears to have been changed.
Refer to the
documentation for the correct position of the variables.
- The documentation ordered for the course and listed on the course
information page is unfortunately no longer available.
You can get similar information online.
- A student called to my attention that SAS versions 8 and 9 give
different results for some calculations I offered in class. Specifically,
Version 9 allows SAS to calculate internally transformed confidence
intervals for survival curves. In class I did this using PROC PHREG.
Unfortunately, the default that Version 9 uses is different from what you get
from Version 8. Specifically, Version 9 uses as the default the complementary
log-log transform, while Version 8 used no transformation. Hence if the
data set mydata contains two variables, t a time variable and s a status
variable, taking the values 0 for censored and 1 for event, then the code
proc lifetest data=mydata outsurv=a; time t*s(0); run;
puts into a confidence intervals on the original scale in version 8,
but on the log log scale in version 9. To get confidence intervals on the
original scale in version 8,do
proc lifetest data=mydata; survival out=a conftype=linear; time t*s(0); run;
- Solutions for exam 1 are available.
- Office hours are 4:30-5:20 Monday and Thursday, and 10:30-11:20
Wednesday, starting immediately.
- A final exam from last year is
available to help students know what format to expect on the final.
Course information
Course Information Student
Survey
Lectures
Lecture 1 (small format)
Lecture 2 (small format)
(SAS input)
(SAS macro)
(SAS output)
(graphs)
(R code)
Lecture 3 (small format)
(SAS input)
(SAS macro)
(SAS output)
(graphs)
(R code)
(R function to do actuarial estimate)
(R code to simulate simultaneous coverage)
Lecture 4 (small format)
(SAS input)
(SAS macro)
(SAS output)
(R code)
(R function to do Renyi and Cramer-von Mises testing)
(graphs)
Lecture 5 (small format)
(Critical values for the Renyi statistic)
Lecture 6 (small format)
(examples) (SAS output)
(R code)
Lecture 7 (small format)
(examples)
(SAS output)
(R code)
(graphs)
Lecture 8 (small format)
(examples)
(graphs)
(R code for some of the examples)
Lecture 9 (small format)
(examples) (SAS macros)
(R code) (R macros)
(graphic output)
Lecture 10 (small format)
(examples) (SAS output)
(graphic output)
Lecture 11 (small format)
(examples) (SAS output)
(graphs) (macros)
Lecture 12 (small format)
(examples) (SAS output)
Lecture 13 (examples)
(graphs)
Review (examples)
(SAS output)
(graphs)
Homeworks
- Solutions
HW 1
HW 2 , plus the data set
HW 3
HW 4
HW 5
HW 6, plus the data set ; note that
- you'll need to add descending to the call to proc genmod, and
- use freq count in proc genmod.
- Already assigned
- Yet to be assigned
Prerequisites
Level V statistics is presumed. This prerequisite may be met by
563 or 586 or 593.
Supplimentary Materials
Other resources
- Software:
- The R program.
You have to separately download and install the survival analysis library.
-
If you have a hard time reading the pdf files I put up, and
are using Windows, you might try downloading
gs5.50,
a postscript reader that will also read pdf, or Acrobat Reader
- Data:
-
The home page for the text,
Survival Analysis: Techniques for Censored and Truncated Data, has links to
data sets
that it references, including
the the breast feeding data ,
the the kidney
transplant data for HW 2,
and
the larynx cancer data of lecture 2, and the
rat data
of lecture 4.
If you use these data sets,
you'll have to clip off the text at the top and set the file extension
yourself.
- The census data from lecture 3.
-
Another source of data sets
that will be used in homework problems., including the raw version of the
Emerson data
for lecture 8.
- Still another
data source, the text Data by Andrews and Herzberg. I have
reformatted their
VA bladder tumor
recurrence data .
- The
Multi-Center AIDS cohort study data, also for lecture 8.
- The ALR data, provided by Prof. Westfall at Texas Tech, and collected by the
Texas Tech Industrial Engineering department.
- The
Stanford Heart Transplant data set. There are two versions of this
data set; I use the first one. To make the data file to read into sas,
- down load the above file,
- cut off the descriptive text at the
beginning,
- cut off the part of the file after the line of all "=="
You'll find that some of the entries are empty. This is ok, since
they all occur at the end of the line and SAS treats them as missing
(which they are).
- The arsenic data for lecture 12,
an ASCII version of
the Word file at
the Statlib archive.
- The ulcer data for lecture 12,
from Collett's book.
- A cleaned up version of the
gamma frailty macro
found in Klein and Klein Moeschberger's text.
-->
Return to:
SAS "feature"