-
-
Notifications
You must be signed in to change notification settings - Fork 77
Closed
Description
Summary:
I have a model, and it compiles. Yay! But when I try to sample, I get error code 70 on each chain, with no explanation about what that means.
Description:
Describe the issue as clearly as possible. Code is below. The stack trace is ~useless.
import cmdstanpy
model_code = """
data {
int <lower=0> num_points;
int <lower=0> num_features;
vector[num_features] locations[num_points];
real <lower=0> durations[num_points];
int <lower=0> observations[num_points];
}
parameters {
real <lower=0> amplitude;
real <lower=0> length_scale;
vector[num_points] log_intensities;
}
transformed parameters {
vector[num_points] zero;
matrix[num_points, num_points] L_K;
real mean_log_intensity;
{
matrix[num_points, num_points] K = gp_matern32_cov(locations, amplitude, length_scale);
L_K = cholesky_decompose(K);
for (i in 1:num_points) mean_log_intensity += (log(observations[i]) - log(durations[i]));
mean_log_intensity /= num_points;
}
}
model {
log_intensities ~ multi_normal_cholesky(zero, L_K);
for (i in 1:num_points) observations[i] ~ poisson_log(mean_log_intensity + log_intensities[i] + log(durations[i]));
}
"""
path='/tmp/foo.stan'
with open(path, 'w') as f:
f.write(model_code)
m = cmdstanpy.Model(path)
m.compile()
fit = m.sample(chains=4, cores=2, data=dict(num_points=3, num_features=2, locations=[[1., 2.], [2., 3], [3., 4]], durations=[2., 4., 3.], observations=[100., 175., 165.]))
fit.summary()
Error log here:
INFO:cmdstanpy:compiling c++
INFO:cmdstanpy:compiled model file: /tmp/foo
INFO:cmdstanpy:start chain 1
INFO:cmdstanpy:start chain 2
INFO:cmdstanpy:finish chain 1
INFO:cmdstanpy:start chain 3
INFO:cmdstanpy:finish chain 2
INFO:cmdstanpy:start chain 4
INFO:cmdstanpy:finish chain 3
INFO:cmdstanpy:finish chain 4
---------------------------------------------------------------------------
RuntimeError Traceback (most recent call last)
<ipython-input-59-09ad9f23ba87> in <module>()
39 m.compile()
40 # print(m.code())
---> 41 fit = m.sample(chains=4, cores=2, data=dict(num_points=3, num_features=2, locations=[[1., 2.], [2., 3], [3., 4]], durations=[2., 4., 3.], observations=[100., 175., 165.]))
42 fit.summary()
/usr/local/lib/python3.6/dist-packages/cmdstanpy/model.py in sample(self, data, chains, cores, seed, chain_ids, inits, warmup_iters, sampling_iters, save_warmup, thin, max_treedepth, metric, step_size, adapt_engaged, adapt_delta, fixed_param, csv_basename, show_progress)
610 msg, i, runset._retcode(i)
611 )
--> 612 raise RuntimeError(msg)
613 stanfit = StanFit(runset, fixed_param)
614 stanfit._validate_csv_files()
RuntimeError: Error during sampling, chain 0 returned error code 70, chain 1 returned error code 70, chain 2 returned error code 70, chain 3 returned error code 70
Additional Information:
Provide any additional information here.
Current Version:
'2.19.0.0'
davidar
Metadata
Metadata
Assignees
Labels
No labels