% usefile c:\oupbuch\leukemia\leukemia_ch6.prg dataset d d.infile using c:\oupbuch\applications\leukemia\leuk.dat d.replace sex = 1*(sex=1) map m m.infile using c:\oupbuch\leukemia\nwengland_smoothed.bnd m.reorder % Empirical Bayes estimates remlreg r % Markov random field logopen using c:\oupbuch\leukemia\results\reml_mrf\remlmrf.log r.outfile = c:\oupbuch\leukemia\results\reml_mrf\remlmrf r.regress cens = time(baseline) + age(psplinerw2) + wbc(psplinerw2) + tpi(psplinerw2) + district(spatial,map=m) + sex, family=cox using d logclose % Bivariate spline logopen using c:\oupbuch\leukemia\results\reml_pspline\remlpspline.log r.outfile = c:\oupbuch\leukemia\results\reml_pspline\remlpspline r.regress cens = time(baseline) + age(psplinerw2) + wbc(psplinerw2) + tpi(psplinerw2) + xcoord*ycoord(pspline2dimrw1,nrknots=12) + sex, family=cox using d logclose % Kriging logopen using c:\oupbuch\leukemia\results\reml_kriging\remlkriging.log r.outfile = c:\oupbuch\leukemia\results\reml_kriging\remlkriging r.regress cens = time(baseline) + age(psplinerw2) + wbc(psplinerw2) + tpi(psplinerw2) + xcoord*ycoord(kriging,nrknots=100) + sex, family=cox using d logclose % Full Bayes estimates bayesreg b % Markov random field logopen using c:\oupbuch\leukemia\results\mcmc_mrf\mcmcmrf.log b.outfile = c:\oupbuch\leukemia\results\mcmc_mrf\mcmcmrf b.regress cens = time(baseline) + age(psplinerw2) + wbc(psplinerw2) + tpi(psplinerw2) + district(spatial,map=m) + sex, family=cox iterations=12000 burnin=2000 predict using d logclose % Bivariate spline logopen using c:\oupbuch\leukemia\results\mcmc_pspline\mcmcpspline.log b.outfile = c:\oupbuch\leukemia\results\mcmc_pspline\mcmcpspline b.regress cens = time(baseline) + age(psplinerw2) + wbc(psplinerw2) + tpi(psplinerw2) + xcoord*ycoord(pspline2dimrw1,nrknots=12) + sex, family=cox iterations=12000 burnin=2000 predict using d logclose drop d m r b