You should have four files: this "installation-S-info.txt" file, and the files "combo.f", "usage-info.txt", and "Splus_code". The file "combo.f" contains Fortran programs which are designed to be called from within Splus. The file "Splus_code" contains Splus functions. To use these programs, carry out the following three steps. 1. Compile the program "combo.f" with the -c option (compile, but do not link) via the command: f77 -c combo.f or f90 -c combo.f This should create an object file named "combo.o" 2. Start Splus and then read in the required Splus code using: source("Splus_code") 3. Load the Fortran program into Splus using either dyn.open(), dyn.load(), dyn.load2(), dyn.load.shared(), or static loading with LOAD. The method used depends on what is available or preferred on your system. These are illustrated below. Our experience is with Unix-type systems. In the new version of S (version 6) you execute: dyn.open("combo.o") In older versions of Splus/S, you do one of the following. If dyn.load() or dyn.load2() are available, then after starting Splus you execute either: dyn.load("combo.o") or dyn.load2("combo.o") If dyn.load2() does not work for you, you may have to use the size option and do something like this: dyn.load2("combo.o", size=250000) If dyn.load.shared() is available, you issue the Unix command: Splus SHLIB -o combo.so combo.o Then you start Splus and execute the command: dyn.load.shared("./combo.so") If none of the above are available, you must use static loading as follows. Issue the Unix command: Splus LOAD combo.o This creates a local version of Splus named local.Sqpe in the current directory. This local version contains the Fortran subroutines in combo. If you execute Splus in this directory, the local version will be used. Usage ----- The following example should work in windows and also in Unix, if installed by root (else see below). library("npbayes") data(ovarian) ovarian.dat <- cbind(ovarian$futime, ovarian$fustat) ovarian.1000 <- gibbs1(data=ovarian.dat, prior=c(.1,.1,1000), nrep=10000) plot(ovarian.1000) summary(ovarian.1000) In Unix, if you installed in a local directory, instead of 'library("npbayes")' you do: library("npbayes", lib.loc="./lib") dyn.load("./lib/npbayes/libs/npbayes.so")