CENSRE2 - Mixed effects models for censored data. See censre2.s for arguments and 
          calling structure

Reference: Hughes JP: Mixed effects models with censored data with application to HIV RNA levels.  Biometrics, 55:625-629, 1999.

2-Feb-2007

To install this software, do the following:

1) Compile the fortran code (need to do only once):

e.g.

% f77 -c censre2.f
  or
% Splus SHLIB -o censre2.so censre2.f
  or
 ...

2) In Splus, source the function censre2() (need to do only once):

> source("censre2.s")

3) In Splus, load the compiled fortran code (need to do everytime
   Splus is started):

e.g.

> dyn.load("censre2.o")
  or
> dyn.load2("censre2.o")
  or
> dyn.load.shared("censre2.so")
  or
 ...

4) Run the test dataset

source("censre.test")

5) Compare the output with censre.test.out


Report problems to Jim Hughes (jphughes@u.washington.edu)

KNOWN BUGS:

1) Values in the extreme tails of the normal distribution can cause
   numerical problems which lead to a NaN in the fortran routines.
   Such values sometimes arise during the gibbs sampling procedure.
   The program cannot recover and all parameter estimates and the
   loglikelihood become NaN's (printed as "NA"). If this occurs,
   try rerunning the program with a different random seed or different
   initial values.
