BTYDplus icon indicating copy to clipboard operation
BTYDplus copied to clipboard

slice_sample_cpp loop does not finish

Open mplatzer opened this issue 6 years ago • 0 comments

structure(list(cust = c(103523979L, 106539747L, 109669649L, 114536229L, 
115893543L, 115959927L, 117271474L, 121815850L, 123091742L, 123545940L, 
126225453L, 126396424L, 126844622L, 126963045L, 132218991L, 133320753L, 
134979819L, 135327695L, 136081180L, 136651599L, 137087039L, 138027869L, 
138717722L, 139202404L, 141638965L, 142495136L, 143705614L, 145281013L, 
147186313L, 147408159L, 149208624L), x = c(0, 0, 0, 0, 1, 0, 
0, 0, 0, 0, 0, 1, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0, 4, 0, 0, 0, 
1, 0, 0, 0), t.x = c(0, 0, 0, 0, 20.1428571428571, 0, 0, 0, 0, 
0, 0, 3, 0, 0, 0.285714285714286, 0, 0, 16.8571428571429, 15.4285714285714, 
0, 0, 0, 0, 9.28571428571429, 0, 0, 0, 16.1428571428571, 0, 0, 
0), litt = c(0, 0, 0, 0, 3.00284974132285, 0, 0, 0, 0, 0, 0, 
1.09861228866811, 0, 0, -1.25276296849537, 0, 0, 2.82477447541035, 
2.73622107806891, 0, 0, 0, 0, 2.06327660482648, 0, 0, 0, 2.78147766965703, 
0, 0, 0), sales = c(154.98, 849.98, 269.94, 54.98, 499.43, 249.99, 
106.98, 944.97, 36.67, 349.97, 24.99, 2149.93, 257.71, 889.95, 
37.97, 29.99, 349.99, 1494.92, 594.92, 977.93, 479.99, 59.99, 
204.96, 7566.69, 54.98, 839.97, 14.99, 2914.91, 192.98, 208.94, 
1029.92), sales.x = c(0, 0, 0, 0, 162.98, 0, 0, 0, 0, 0, 0, 869.95, 
0, 0, 29.98, 0, 0, 274.94, 149.99, 0, 0, 0, 0, 1698.89, 0, 0, 
0, 139.99, 0, 0, 0), first = structure(c(12059, 12056, 12070, 
12079, 12075, 12056, 12077, 12066, 12058, 12070, 12059, 12078, 
12062, 12054, 12055, 12054, 12055, 12071, 12069, 12056, 12077, 
12064, 12062, 12066, 12082, 12080, 12082, 12069, 12062, 12069, 
12059), class = "Date"), T.cal = c(23.8571428571429, 24.2857142857143, 
22.2857142857143, 21, 21.5714285714286, 24.2857142857143, 21.2857142857143, 
22.8571428571429, 24, 22.2857142857143, 23.8571428571429, 21.1428571428571, 
23.4285714285714, 24.5714285714286, 24.4285714285714, 24.5714285714286, 
24.4285714285714, 22.1428571428571, 22.4285714285714, 24.2857142857143, 
21.2857142857143, 23.1428571428571, 23.4285714285714, 22.8571428571429, 
20.5714285714286, 20.8571428571429, 20.5714285714286, 22.4285714285714, 
23.4285714285714, 22.4285714285714, 23.8571428571429), T.star = c(74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857, 74.7142857142857, 74.7142857142857, 
74.7142857142857, 74.7142857142857), x.star = c(0L, 0L, 0L, 1L, 
0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 1L, 1L, 0L, 0L, 6L, 1L, 0L, 
0L, 1L, 0L, 2L, 0L, 2L, 0L, 2L, 0L, 0L, 1L), sales.star = c(0, 
0, 0, 54.9, 0, 155.97, 0, 0, 0, 299.94, 0, 379.98, 0, 210.76, 
149.99, 0, 0, 1597.82, 389.98, 0, 0, 393.92, 0, 399.97, 0, 1764.45, 
0, 1899.94, 0, 0, 49.99)), .Names = c("cust", "x", "t.x", "litt", 
"sales", "sales.x", "first", "T.cal", "T.star", "x.star", "sales.star"
), row.names = c(NA, -31L), class = "data.frame")

pggg.draws <- pggg.mcmc.DrawParameters(cbs, mcmc=2500, burnin=500, chains = 1, thin=50)

results in

Error in pggg_slice_sample("lambda", x = data$x, tx = data$t.x, Tcal = data$T.cal,  :
  slice_sample_cpp loop did not finish

mplatzer avatar Nov 15 '17 14:11 mplatzer