Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

slice_sample_cpp loop does not finish #61

Open
mplatzer opened this issue Nov 15, 2017 · 0 comments
Open

slice_sample_cpp loop does not finish #61

mplatzer opened this issue Nov 15, 2017 · 0 comments

Comments

@mplatzer
Copy link
Owner

mplatzer commented Nov 15, 2017

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
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant