Simulated data for the two-dimensional Lexis prior
Tommi Härkänen and Anna But
The Lexis prior for two-dimensional hazard surface and this example has been presented in (Härkänen, But and Haukka, 2017). Before rendering this file with Rmarkdown, please run the MCMC simulation using the Linux shell command $ bite < l1.opt > out.txt . Further instructions can be found in the manual.
Plot the figures
Figures of Bite output can be generated running the splus.S file in R e.g. using the source("splus.S") function.
Another option is to drop the pdf(), dev.off() and q() functions from the file and including the source in a code block of an R markdown document.
cmd <- scan("splus.S", what="", sep="\n")
## Drop the function calls for the pdf output:
cmd <- cmd[!grepl("pdf\\(|dev.off()|q()", cmd)]
## Run the plotting commands:
eval(parse(text=paste(cmd, collapse="\n")))## [1] 2 2
## [1] 2 2
## [1] 2 2
## [1] 2 2
Scalar parameters
Regression coefficients
Read the MCMC sample of the \(\phi\) coefficient from file phi1.tmp. Note that the results differ from those in the article (Härkänen, But and Haukka, 2017) as the seed of the random number generator is not the same.
## Read the output data files and construct a matrix 'betas':
betas <- c()
for (i in dir(pattern="phi1"))
betas <- cbind(betas, scan(i))
dimnames(betas) <- list(1:nrow(betas), sub(".tmp", "", dir(pattern="phi1")))
## Convert the matrix into MCMC object of the coda package:
library(coda)
betas <- mcmc(betas)
## Summary of betas:
summary(betas)##
## Iterations = 1:450
## Thinning interval = 1
## Number of chains = 1
## Sample size per chain = 450
##
## 1. Empirical mean and standard deviation for each variable,
## plus standard error of the mean:
##
## Mean SD Naive SE Time-series SE
## 3.12197 1.13469 0.05349 0.06295
##
## 2. Quantiles for each variable:
##
## 2.5% 25% 50% 75% 97.5%
## 1.384 2.317 2.934 3.726 5.737Hazard function values at predefined points
The expectation command can be applied to calculate posterior expectations of different functionals of parameters e.g. hazard functions. The save expectation command saves the values at the MCMC iterations. For example, expectation elambda150 : value(lambda0(),t) range t (5, 70, 13); calculates the posterior expectations of \(\lambda_0(t)\) at 13 time points 5, 10, , 65, and the output can be found in file elambda150.0.tmp. The MCMC iterations are saved using the command save expectation elambda150;, and the values are store in file elambda150.0.sav.
## Read the output data files and construct a matrix 'betas':
haz <- NULL
for (i in paste0("elambda150.0.sav")) {
tmp <- read.table(i)
colnames(tmp) <- paste0(sub(".0.sav", "", i), "_", 1:ncol(tmp))
if (!is.null(haz))
haz <- cbind(haz, tmp) else haz <- tmp
}
## Convert the matrix into MCMC object of the coda package:
haz <- mcmc(haz)
## Summary of hazard functions:
summary(haz)##
## Iterations = 1:450
## Thinning interval = 1
## Number of chains = 1
## Sample size per chain = 450
##
## 1. Empirical mean and standard deviation for each variable,
## plus standard error of the mean:
##
## Mean SD Naive SE Time-series SE
## elambda150_1 0.004060 0.0032975 1.554e-04 2.213e-04
## elambda150_2 0.004026 0.0032320 1.524e-04 2.211e-04
## elambda150_3 0.003953 0.0031013 1.462e-04 2.118e-04
## elambda150_4 0.003896 0.0030370 1.432e-04 2.165e-04
## elambda150_5 0.003852 0.0030194 1.423e-04 2.093e-04
## elambda150_6 0.003822 0.0028939 1.364e-04 2.035e-04
## elambda150_7 0.003795 0.0028791 1.357e-04 1.979e-04
## elambda150_8 0.003781 0.0028709 1.353e-04 1.946e-04
## elambda150_9 0.003643 0.0024830 1.170e-04 1.767e-04
## elambda150_10 0.003505 0.0019939 9.399e-05 1.613e-04
## elambda150_11 0.003397 0.0016921 7.977e-05 1.470e-04
## elambda150_12 0.003327 0.0016333 7.699e-05 1.510e-04
## elambda150_13 0.003263 0.0015422 7.270e-05 1.452e-04
## elambda150_14 0.003238 0.0015449 7.283e-05 1.555e-04
## elambda150_15 0.003185 0.0015136 7.135e-05 1.298e-04
## elambda150_16 0.003174 0.0015061 7.100e-05 1.254e-04
## elambda150_17 0.003180 0.0015030 7.085e-05 1.254e-04
## elambda150_18 0.003179 0.0015043 7.092e-05 1.260e-04
## elambda150_19 0.003157 0.0014984 7.064e-05 1.358e-04
## elambda150_20 0.003113 0.0014040 6.618e-05 1.229e-04
## elambda150_21 0.003101 0.0013993 6.596e-05 1.248e-04
## elambda150_22 0.003110 0.0014183 6.686e-05 1.253e-04
## elambda150_23 0.003120 0.0015299 7.212e-05 1.310e-04
## elambda150_24 0.003162 0.0016033 7.558e-05 1.501e-04
## elambda150_25 0.003181 0.0016060 7.571e-05 1.453e-04
## elambda150_26 0.003171 0.0016123 7.600e-05 1.469e-04
## elambda150_27 0.003141 0.0016131 7.604e-05 1.449e-04
## elambda150_28 0.003122 0.0016248 7.659e-05 1.376e-04
## elambda150_29 0.003031 0.0015230 7.180e-05 1.302e-04
## elambda150_30 0.002909 0.0014099 6.646e-05 9.777e-05
## elambda150_31 0.002884 0.0013895 6.550e-05 9.007e-05
## elambda150_32 0.002776 0.0012673 5.974e-05 7.873e-05
## elambda150_33 0.002735 0.0012632 5.955e-05 7.955e-05
## elambda150_34 0.002699 0.0012072 5.691e-05 6.941e-05
## elambda150_35 0.002684 0.0011979 5.647e-05 6.939e-05
## elambda150_36 0.002654 0.0011689 5.510e-05 6.811e-05
## elambda150_37 0.002611 0.0011176 5.268e-05 6.580e-05
## elambda150_38 0.002574 0.0010795 5.089e-05 6.845e-05
## elambda150_39 0.002564 0.0010601 4.997e-05 6.808e-05
## elambda150_40 0.002545 0.0010358 4.883e-05 6.811e-05
## elambda150_41 0.002512 0.0009950 4.691e-05 6.176e-05
## elambda150_42 0.002494 0.0009762 4.602e-05 6.081e-05
## elambda150_43 0.002476 0.0009547 4.501e-05 6.045e-05
## elambda150_44 0.002468 0.0009713 4.579e-05 6.315e-05
## elambda150_45 0.002453 0.0009333 4.400e-05 6.283e-05
## elambda150_46 0.002439 0.0008948 4.218e-05 5.772e-05
## elambda150_47 0.002435 0.0008982 4.234e-05 6.579e-05
## elambda150_48 0.002439 0.0008955 4.222e-05 6.496e-05
## elambda150_49 0.002435 0.0008995 4.240e-05 5.929e-05
## elambda150_50 0.002426 0.0008931 4.210e-05 6.847e-05
## elambda150_51 0.002425 0.0009035 4.259e-05 6.173e-05
## elambda150_52 0.002430 0.0008975 4.231e-05 6.040e-05
## elambda150_53 0.002421 0.0008963 4.225e-05 5.832e-05
## elambda150_54 0.002419 0.0008779 4.139e-05 5.683e-05
## elambda150_55 0.002422 0.0008788 4.142e-05 6.084e-05
## elambda150_56 0.002426 0.0008859 4.176e-05 5.475e-05
## elambda150_57 0.002419 0.0008713 4.107e-05 5.922e-05
## elambda150_58 0.002431 0.0008812 4.154e-05 6.033e-05
## elambda150_59 0.002435 0.0008970 4.228e-05 5.945e-05
## elambda150_60 0.002448 0.0009084 4.282e-05 5.997e-05
## elambda150_61 0.002453 0.0009141 4.309e-05 6.065e-05
## elambda150_62 0.002461 0.0009101 4.290e-05 5.978e-05
## elambda150_63 0.002477 0.0009135 4.306e-05 5.519e-05
## elambda150_64 0.002488 0.0009309 4.389e-05 6.017e-05
## elambda150_65 0.002501 0.0009332 4.399e-05 5.639e-05
## elambda150_66 0.002521 0.0009397 4.430e-05 5.460e-05
## elambda150_67 0.002535 0.0009499 4.478e-05 5.462e-05
## elambda150_68 0.002548 0.0009620 4.535e-05 7.493e-05
## elambda150_69 0.002557 0.0009524 4.490e-05 7.323e-05
## elambda150_70 0.002567 0.0009464 4.461e-05 5.436e-05
## elambda150_71 0.002570 0.0009411 4.436e-05 5.432e-05
## elambda150_72 0.002573 0.0009299 4.384e-05 5.522e-05
## elambda150_73 0.002573 0.0009264 4.367e-05 5.466e-05
## elambda150_74 0.002578 0.0009410 4.436e-05 7.442e-05
## elambda150_75 0.002572 0.0009401 4.432e-05 7.488e-05
## elambda150_76 0.002575 0.0009407 4.435e-05 7.530e-05
## elambda150_77 0.002584 0.0009411 4.436e-05 7.410e-05
## elambda150_78 0.002587 0.0009482 4.470e-05 7.650e-05
## elambda150_79 0.002589 0.0009469 4.464e-05 7.636e-05
## elambda150_80 0.002595 0.0009439 4.449e-05 6.614e-05
## elambda150_81 0.002597 0.0009407 4.434e-05 6.622e-05
## elambda150_82 0.002614 0.0009381 4.422e-05 6.645e-05
## elambda150_83 0.002613 0.0009389 4.426e-05 6.628e-05
## elambda150_84 0.002621 0.0009420 4.441e-05 6.736e-05
## elambda150_85 0.002647 0.0009615 4.532e-05 5.519e-05
## elambda150_86 0.002641 0.0009605 4.528e-05 5.875e-05
## elambda150_87 0.002641 0.0009621 4.535e-05 5.886e-05
## elambda150_88 0.002647 0.0009592 4.522e-05 5.857e-05
## elambda150_89 0.002637 0.0009517 4.486e-05 6.001e-05
## elambda150_90 0.002642 0.0009476 4.467e-05 5.979e-05
## elambda150_91 0.002639 0.0009474 4.466e-05 5.903e-05
## elambda150_92 0.002641 0.0009522 4.489e-05 6.006e-05
## elambda150_93 0.002647 0.0009513 4.484e-05 5.415e-05
## elambda150_94 0.002652 0.0009505 4.481e-05 6.138e-05
## elambda150_95 0.002653 0.0009630 4.540e-05 6.757e-05
## elambda150_96 0.002651 0.0009628 4.539e-05 6.752e-05
## elambda150_97 0.002668 0.0009794 4.617e-05 6.669e-05
## elambda150_98 0.002672 0.0009802 4.621e-05 6.600e-05
## elambda150_99 0.002676 0.0009895 4.665e-05 6.620e-05
## elambda150_100 0.002682 0.0010073 4.749e-05 6.628e-05
## elambda150_101 0.002690 0.0010205 4.810e-05 6.695e-05
## elambda150_102 0.002684 0.0010192 4.804e-05 6.712e-05
## elambda150_103 0.002678 0.0010101 4.762e-05 6.771e-05
## elambda150_104 0.002679 0.0010110 4.766e-05 6.799e-05
## elambda150_105 0.002673 0.0010625 5.009e-05 6.569e-05
## elambda150_106 0.002660 0.0010631 5.011e-05 5.929e-05
## elambda150_107 0.002651 0.0010709 5.048e-05 6.061e-05
## elambda150_108 0.002647 0.0010618 5.006e-05 5.927e-05
## elambda150_109 0.002639 0.0010759 5.072e-05 5.937e-05
## elambda150_110 0.002649 0.0010727 5.057e-05 5.970e-05
## elambda150_111 0.002624 0.0010222 4.819e-05 6.042e-05
## elambda150_112 0.002639 0.0010402 4.903e-05 6.241e-05
## elambda150_113 0.002631 0.0010475 4.938e-05 6.460e-05
## elambda150_114 0.002632 0.0010385 4.896e-05 6.433e-05
## elambda150_115 0.002626 0.0010377 4.892e-05 6.445e-05
## elambda150_116 0.002609 0.0010115 4.768e-05 6.225e-05
## elambda150_117 0.002612 0.0010112 4.767e-05 6.280e-05
## elambda150_118 0.002629 0.0010299 4.855e-05 6.614e-05
## elambda150_119 0.002621 0.0010339 4.874e-05 6.690e-05
## elambda150_120 0.002614 0.0010566 4.981e-05 5.980e-05
## elambda150_121 0.002598 0.0010330 4.869e-05 6.055e-05
## elambda150_122 0.002575 0.0010194 4.805e-05 5.789e-05
## elambda150_123 0.002568 0.0010182 4.800e-05 5.702e-05
## elambda150_124 0.002559 0.0010162 4.790e-05 5.647e-05
## elambda150_125 0.002556 0.0010197 4.807e-05 5.635e-05
## elambda150_126 0.002557 0.0010214 4.815e-05 5.614e-05
## elambda150_127 0.002568 0.0010234 4.824e-05 5.654e-05
## elambda150_128 0.002573 0.0010318 4.864e-05 5.727e-05
## elambda150_129 0.002585 0.0010434 4.919e-05 5.701e-05
## elambda150_130 0.002608 0.0010807 5.095e-05 6.235e-05
## elambda150_131 0.002605 0.0010624 5.008e-05 6.217e-05
## elambda150_132 0.002605 0.0010608 5.001e-05 6.214e-05
## elambda150_133 0.002624 0.0010952 5.163e-05 6.631e-05
## elambda150_134 0.002637 0.0011441 5.394e-05 6.771e-05
## elambda150_135 0.002589 0.0011632 5.484e-05 7.040e-05
## elambda150_136 0.002560 0.0011972 5.644e-05 7.530e-05
## elambda150_137 0.002556 0.0012096 5.702e-05 7.639e-05
## elambda150_138 0.002530 0.0012274 5.786e-05 7.954e-05
## elambda150_139 0.002491 0.0012352 5.823e-05 8.184e-05
## elambda150_140 0.002456 0.0012020 5.667e-05 7.790e-05
## elambda150_141 0.002400 0.0012066 5.688e-05 7.827e-05
## elambda150_142 0.002360 0.0012290 5.793e-05 7.948e-05
## elambda150_143 0.002335 0.0012439 5.864e-05 8.365e-05
## elambda150_144 0.002302 0.0012366 5.829e-05 8.811e-05
## elambda150_145 0.002287 0.0012553 5.918e-05 8.684e-05
## elambda150_146 0.002290 0.0012746 6.009e-05 8.600e-05
## elambda150_147 0.002258 0.0012674 5.975e-05 8.398e-05
## elambda150_148 0.002230 0.0012705 5.989e-05 9.051e-05
## elambda150_149 0.002228 0.0012771 6.021e-05 8.982e-05
## elambda150_150 0.002228 0.0012843 6.054e-05 9.041e-05
##
## 2. Quantiles for each variable:
##
## 2.5% 25% 50% 75% 97.5%
## elambda150_1 0.0013224 0.002383 0.003105 0.004253 0.011584
## elambda150_2 0.0013532 0.002383 0.003098 0.004213 0.011077
## elambda150_3 0.0013532 0.002383 0.003088 0.004175 0.010761
## elambda150_4 0.0013532 0.002383 0.003085 0.004164 0.010113
## elambda150_5 0.0013532 0.002361 0.003076 0.004117 0.010113
## elambda150_6 0.0013572 0.002361 0.003076 0.004117 0.009830
## elambda150_7 0.0013532 0.002361 0.003076 0.004081 0.009830
## elambda150_8 0.0013532 0.002361 0.003071 0.004081 0.009830
## elambda150_9 0.0013224 0.002339 0.003047 0.003987 0.008749
## elambda150_10 0.0013224 0.002324 0.003017 0.003968 0.007862
## elambda150_11 0.0013082 0.002324 0.002967 0.003926 0.007439
## elambda150_12 0.0012894 0.002272 0.002927 0.003903 0.007404
## elambda150_13 0.0012894 0.002228 0.002908 0.003869 0.007274
## elambda150_14 0.0012630 0.002213 0.002896 0.003846 0.007274
## elambda150_15 0.0012630 0.002159 0.002881 0.003773 0.007274
## elambda150_16 0.0012894 0.002158 0.002875 0.003751 0.007274
## elambda150_17 0.0012894 0.002167 0.002878 0.003751 0.007274
## elambda150_18 0.0012894 0.002167 0.002878 0.003751 0.007274
## elambda150_19 0.0012630 0.002158 0.002875 0.003750 0.007274
## elambda150_20 0.0012630 0.002158 0.002875 0.003725 0.006621
## elambda150_21 0.0012338 0.002152 0.002850 0.003743 0.006587
## elambda150_22 0.0012338 0.002152 0.002841 0.003750 0.006621
## elambda150_23 0.0012058 0.002151 0.002838 0.003743 0.007099
## elambda150_24 0.0012058 0.002151 0.002838 0.003751 0.007284
## elambda150_25 0.0012058 0.002152 0.002866 0.003773 0.007284
## elambda150_26 0.0012058 0.002152 0.002841 0.003773 0.007284
## elambda150_27 0.0011322 0.002138 0.002837 0.003743 0.007284
## elambda150_28 0.0011126 0.002131 0.002816 0.003719 0.007284
## elambda150_29 0.0011126 0.002045 0.002751 0.003646 0.007257
## elambda150_30 0.0010961 0.001978 0.002653 0.003521 0.006466
## elambda150_31 0.0010961 0.001973 0.002596 0.003510 0.006436
## elambda150_32 0.0010154 0.001948 0.002530 0.003382 0.005961
## elambda150_33 0.0009355 0.001878 0.002495 0.003280 0.005951
## elambda150_34 0.0009438 0.001854 0.002473 0.003279 0.005951
## elambda150_35 0.0009438 0.001843 0.002445 0.003226 0.005951
## elambda150_36 0.0009438 0.001831 0.002435 0.003207 0.005861
## elambda150_37 0.0009753 0.001831 0.002415 0.003166 0.005317
## elambda150_38 0.0009849 0.001815 0.002401 0.003144 0.005186
## elambda150_39 0.0009876 0.001822 0.002401 0.003105 0.005087
## elambda150_40 0.0009876 0.001822 0.002392 0.003085 0.004926
## elambda150_41 0.0009849 0.001796 0.002385 0.003063 0.004589
## elambda150_42 0.0009753 0.001796 0.002383 0.003048 0.004280
## elambda150_43 0.0009876 0.001790 0.002383 0.003015 0.004203
## elambda150_44 0.0009753 0.001790 0.002383 0.002993 0.004203
## elambda150_45 0.0009753 0.001790 0.002383 0.002961 0.004192
## elambda150_46 0.0009753 0.001790 0.002383 0.002943 0.004175
## elambda150_47 0.0009753 0.001790 0.002374 0.002943 0.004175
## elambda150_48 0.0009753 0.001790 0.002383 0.002943 0.004175
## elambda150_49 0.0009753 0.001789 0.002383 0.002943 0.004175
## elambda150_50 0.0009753 0.001789 0.002361 0.002928 0.004175
## elambda150_51 0.0009172 0.001765 0.002374 0.002938 0.004192
## elambda150_52 0.0009297 0.001790 0.002383 0.002928 0.004192
## elambda150_53 0.0008286 0.001790 0.002383 0.002927 0.004175
## elambda150_54 0.0009297 0.001790 0.002383 0.002927 0.004164
## elambda150_55 0.0009297 0.001790 0.002383 0.002927 0.004175
## elambda150_56 0.0009297 0.001789 0.002383 0.002936 0.004192
## elambda150_57 0.0009297 0.001785 0.002383 0.002928 0.004175
## elambda150_58 0.0009753 0.001790 0.002363 0.002928 0.004192
## elambda150_59 0.0009753 0.001790 0.002356 0.002936 0.004203
## elambda150_60 0.0009753 0.001796 0.002363 0.002943 0.004243
## elambda150_61 0.0009753 0.001815 0.002363 0.002948 0.004280
## elambda150_62 0.0009849 0.001825 0.002374 0.002948 0.004280
## elambda150_63 0.0009849 0.001831 0.002383 0.002964 0.004298
## elambda150_64 0.0009924 0.001831 0.002383 0.002965 0.004403
## elambda150_65 0.0010198 0.001831 0.002389 0.002968 0.004448
## elambda150_66 0.0010363 0.001833 0.002402 0.003024 0.004565
## elambda150_67 0.0010363 0.001843 0.002408 0.003039 0.004643
## elambda150_68 0.0010363 0.001854 0.002416 0.003048 0.004703
## elambda150_69 0.0010960 0.001864 0.002419 0.003052 0.004703
## elambda150_70 0.0012058 0.001869 0.002438 0.003053 0.004703
## elambda150_71 0.0012864 0.001869 0.002438 0.003053 0.004703
## elambda150_72 0.0012864 0.001877 0.002439 0.003066 0.004643
## elambda150_73 0.0012864 0.001893 0.002445 0.003063 0.004643
## elambda150_74 0.0012864 0.001893 0.002445 0.003053 0.004703
## elambda150_75 0.0012864 0.001877 0.002439 0.003052 0.004703
## elambda150_76 0.0012864 0.001877 0.002445 0.003053 0.004703
## elambda150_77 0.0012864 0.001893 0.002469 0.003066 0.004703
## elambda150_78 0.0012797 0.001893 0.002469 0.003066 0.004840
## elambda150_79 0.0012864 0.001895 0.002469 0.003066 0.004840
## elambda150_80 0.0012864 0.001905 0.002482 0.003079 0.004840
## elambda150_81 0.0012864 0.001927 0.002487 0.003079 0.004840
## elambda150_82 0.0013000 0.001972 0.002492 0.003085 0.004840
## elambda150_83 0.0012421 0.001973 0.002495 0.003085 0.004869
## elambda150_84 0.0012421 0.001976 0.002499 0.003089 0.004869
## elambda150_85 0.0012421 0.001989 0.002558 0.003104 0.004934
## elambda150_86 0.0011840 0.001989 0.002544 0.003104 0.004869
## elambda150_87 0.0012421 0.001981 0.002502 0.003131 0.004869
## elambda150_88 0.0012421 0.001989 0.002544 0.003121 0.004869
## elambda150_89 0.0011840 0.001989 0.002502 0.003104 0.004817
## elambda150_90 0.0012421 0.001989 0.002502 0.003104 0.004817
## elambda150_91 0.0012421 0.001989 0.002502 0.003099 0.004817
## elambda150_92 0.0013000 0.001989 0.002502 0.003089 0.004934
## elambda150_93 0.0012948 0.002003 0.002502 0.003097 0.004934
## elambda150_94 0.0012948 0.002010 0.002509 0.003097 0.005005
## elambda150_95 0.0012948 0.002008 0.002509 0.003097 0.005005
## elambda150_96 0.0012948 0.002008 0.002509 0.003089 0.005005
## elambda150_97 0.0013000 0.002008 0.002522 0.003099 0.005260
## elambda150_98 0.0013000 0.002008 0.002522 0.003104 0.005260
## elambda150_99 0.0013000 0.002010 0.002509 0.003104 0.005340
## elambda150_100 0.0012948 0.002008 0.002522 0.003112 0.005445
## elambda150_101 0.0012821 0.002008 0.002503 0.003123 0.005528
## elambda150_102 0.0012424 0.002008 0.002509 0.003123 0.005528
## elambda150_103 0.0012424 0.002008 0.002503 0.003112 0.005528
## elambda150_104 0.0012424 0.002008 0.002503 0.003112 0.005528
## elambda150_105 0.0011721 0.001992 0.002499 0.003104 0.005554
## elambda150_106 0.0011258 0.001989 0.002492 0.003097 0.005554
## elambda150_107 0.0010988 0.001983 0.002490 0.003089 0.005554
## elambda150_108 0.0011121 0.001981 0.002490 0.003085 0.005528
## elambda150_109 0.0010839 0.001974 0.002482 0.003089 0.005554
## elambda150_110 0.0011121 0.001976 0.002490 0.003099 0.005554
## elambda150_111 0.0010907 0.001974 0.002482 0.003097 0.005445
## elambda150_112 0.0010839 0.001976 0.002492 0.003112 0.005445
## elambda150_113 0.0010611 0.001974 0.002487 0.003123 0.005445
## elambda150_114 0.0010611 0.001974 0.002492 0.003127 0.005340
## elambda150_115 0.0010611 0.001972 0.002487 0.003123 0.005340
## elambda150_116 0.0010143 0.001972 0.002482 0.003123 0.005184
## elambda150_117 0.0010143 0.001972 0.002482 0.003127 0.005184
## elambda150_118 0.0010008 0.001976 0.002492 0.003136 0.005249
## elambda150_119 0.0009789 0.001972 0.002482 0.003127 0.005249
## elambda150_120 0.0008995 0.001963 0.002460 0.003123 0.005249
## elambda150_121 0.0008995 0.001957 0.002428 0.003105 0.005184
## elambda150_122 0.0008419 0.001905 0.002392 0.003101 0.005155
## elambda150_123 0.0008419 0.001905 0.002389 0.003085 0.005155
## elambda150_124 0.0008419 0.001899 0.002383 0.003080 0.005155
## elambda150_125 0.0008163 0.001896 0.002374 0.003080 0.005155
## elambda150_126 0.0008163 0.001896 0.002383 0.003085 0.005155
## elambda150_127 0.0008163 0.001896 0.002392 0.003101 0.005155
## elambda150_128 0.0007812 0.001896 0.002399 0.003105 0.005155
## elambda150_129 0.0007812 0.001899 0.002428 0.003127 0.005184
## elambda150_130 0.0007755 0.001905 0.002438 0.003139 0.005276
## elambda150_131 0.0007755 0.001942 0.002444 0.003139 0.005249
## elambda150_132 0.0007892 0.001942 0.002438 0.003139 0.005249
## elambda150_133 0.0007892 0.001950 0.002438 0.003144 0.005276
## elambda150_134 0.0007609 0.001946 0.002444 0.003167 0.005440
## elambda150_135 0.0006641 0.001891 0.002393 0.003139 0.005440
## elambda150_136 0.0005662 0.001860 0.002374 0.003127 0.005528
## elambda150_137 0.0005645 0.001855 0.002364 0.003127 0.005528
## elambda150_138 0.0002829 0.001815 0.002357 0.003127 0.005528
## elambda150_139 0.0002200 0.001780 0.002337 0.003089 0.005528
## elambda150_140 0.0002200 0.001777 0.002325 0.003066 0.005313
## elambda150_141 0.0001577 0.001684 0.002291 0.003037 0.005276
## elambda150_142 0.0001433 0.001603 0.002270 0.003035 0.005276
## elambda150_143 0.0001395 0.001549 0.002259 0.003028 0.005276
## elambda150_144 0.0001239 0.001528 0.002235 0.003008 0.005249
## elambda150_145 0.0001239 0.001501 0.002224 0.003008 0.005276
## elambda150_146 0.0001239 0.001493 0.002224 0.003008 0.005313
## elambda150_147 0.0001239 0.001463 0.002210 0.002983 0.005276
## elambda150_148 0.0001018 0.001439 0.002199 0.002940 0.005249
## elambda150_149 0.0001018 0.001430 0.002189 0.002983 0.005249
## elambda150_150 0.0001030 0.001427 0.002183 0.002983 0.005249The above output read the iterations saved for the first strip. The next strip can be found in file elambda150.1.sav etc. Respectively, file elambda150.14.sav corresponds to the 15th strip.
References
Härkänen T, But A and Haukka J (2017). Non-parametric Bayesian Intensity Model: Exploring Time-to-Event Data on Two Time Scales. To appear in the Scandinavian Journal of Statistics.