M Sethi

Considering that it's been over a day and it's apparently still not starting the simulation, I probably need to do some rethinking, yes! This entire project has been driving me insane because the model and the dataset in combination are so computationally demanding. Again, I appreciate your responsiveness immensely.

Matt!! This is so helpful, and I appreciate it more than you know. I'm trying a test with just one simulation that drops 10% of the data, and I was wondering how you can tell when it's reached the simulation stage. This is the code I'm running now: datafit <- run.jags.study(simulations=1, model=model, datafunction=datafun_drop10pc, targets=list(y=y), monitor=params, data=datalist_study, inits=list(chain1, chain2, chain3)) It's been about an hour and I haven't yet seen the "Testing the model and data...

(Hmm, it looks like the attachment is still too big, so please let me know if it would be helpful to have it and I will send it to you another way. )

Coming back to say that since the model takes an extremely long time to run (though I suppose you could reduce the number of iterations), I'm also attaching a zip file containing results that you can use to test the drop.k() function. Thank you again!

Hi Matt, you're a rock star; thanks for being willing to take a look. I'm attaching the script and a data file (loaded in the script as jags_data. I used file.choose() here to load it for simplicity's sake. The data file is called pheno_tidy_disp.csv.)

Hi Matt, I'm still struggling to fix this problem, and I know it's not your job to help everyone who uses your package (which is awesome and I'm so grateful for it) but if you could even give me a hint of how to approach it I would be in such a better position. The complete model is attached:

Thanks for the quick response. Sorry for my ignorance, but I am not sure how I can address the issue. If you have any idea how I can fix this, please let me know.

Hmm. Ok, I need to think on that. In any case, this is eventually ended with an error: Error: The following error occured when compiling and adapting the model using rjags: Error in rjags::jags.model(model, data = dataenv, inits = inits, n.chains = length(runjags.object$end.state), : RUNTIME ERROR: Compilation error on line 86. Cannot evaluate range of counter j "j" is my datapoints. I'm not sure what is going on here, since the model itself compiled and ran successfully to begin with. Any insight?...