Pystan -Initilization failed-

Hello everyone.

I’m a begginer of pystan, so I can’t understand why this error has happened.
Please help me. I used Python in Google Colab.

# sample.stan
data {
  int<lower = 1> C1;
  int<lower = 1> C2;
  int<lower = 1> C3;
  int<lower = 1> C4;
  int<lower = 1> C5;
}

transformed data{
  int Count[5];
  Count[1] = C1;
  Count[2] = C2;
  Count[3] = C3;
  Count[4] = C4;
  Count[5] = C5;
}

parameters{
  real<lower=0,upper=1> p1;
  real<lower=0,upper=1> p2;
  real<lower=0,upper=1> p3;
  real<lower=0,upper=1> p4;
  real<lower=0,upper=1> p5;
}

transformed parameters{
  vector[5] p;
  p[1] = p1;
  p[2] = p2;
  p[3] = p3;
  p[4] = p4;
  p[5] = p5;
}

model{
  p1 ~ uniform(0,1);
  p2 ~ uniform(0,1);
  p3 ~ uniform(0,1);
  p4 ~ uniform(0,1);
  p5 ~ uniform(0,1);

  Count ~ multinomial(p);
}
# python
import pystan

stan_data = {'C1':36, 'C2':47, 'C3':37, 'C4':45, 'C5':40}

sm = pystan.StanModel(file='sample.stan')
fit = sm.sampling(data=stan_data,n_jobs=1)

---------------------------------------------------------------------------
RuntimeError                              Traceback (most recent call last)
<ipython-input-7-68dcd40d8b15> in <module>()
----> 1 fit = sm.sampling(data=stan_data,n_jobs=1)

1 frames
/usr/local/lib/python3.6/dist-packages/pystan/model.py in sampling(self, data, pars, chains, iter, warmup, thin, seed, init, sample_file, diagnostic_file, verbose, algorithm, control, n_jobs, **kwargs)
    811         call_sampler_args = izip(itertools.repeat(data), args_list, itertools.repeat(pars))
    812         call_sampler_star = self.module._call_sampler_star
--> 813         ret_and_samples = _map_parallel(call_sampler_star, call_sampler_args, n_jobs)
    814         samples = [smpl for _, smpl in ret_and_samples]
    815 

/usr/local/lib/python3.6/dist-packages/pystan/model.py in _map_parallel(function, args, n_jobs)
     88             pool.join()
     89     else:
---> 90         map_result = list(map(function, args))
     91     return map_result
     92 

stanfit4anon_model_531be191e7a7ea491d3706e6921b03ba_1149686610379806072.pyx in stanfit4anon_model_531be191e7a7ea491d3706e6921b03ba_1149686610379806072._call_sampler_star()

stanfit4anon_model_531be191e7a7ea491d3706e6921b03ba_1149686610379806072.pyx in stanfit4anon_model_531be191e7a7ea491d3706e6921b03ba_1149686610379806072._call_sampler()

RuntimeError: Initialization failed.
data {
  int<lower = 1> C1;
  int<lower = 1> C2;
  int<lower = 1> C3;
  int<lower = 1> C4;
  int<lower = 1> C5;
}

transformed data{
  int Count[5] = {C1, C2, C3, C4, C5};
}

parameters{
  simplex[5] p;
}
model{
  p ~ dirichlet([1,1,1,1,1]');
  Count ~ multinomial(p);
}
2 Likes

Thank you very much! It went well.