AssertionError with logcdf method for a Beta

Questions
6

I see, thanks. So with this code :

basic_model = pm.Model()

with basic_model:

    # Priors for unknown model parameters
    a = pm.HalfNormal('alpha', sd=10)
    b = 1.0
    sigma = pm.HalfNormal('sigma', sd=1)
    a.tag.test_value = 1

    m1 = pm.Beta.dist(alpha=a, beta=b).logcdf(0.1)


    # Likelihood (sampling distribution) of observations
    Y_obs = pm.Normal('Y_obs', mu=m1, sd=sigma, observed=Y)

with basic_model:
      trace = pm.sample(tune=500)
      pm.plot_posterior(trace)

I get :

Running on PyMC3 v3.6
Auto-assigning NUTS sampler…
Initializing NUTS using jitter+adapt_diag…
Traceback (most recent call last):
File “/Users/Colo/Google Drive/Projects/greghec/learning pymc2.py”, line 62, in
trace = pm.sample(tune=500)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/sampling.py”, line 395, in sample
progressbar=progressbar, **args)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/sampling.py”, line 1515, in init_nuts
step = pm.NUTS(potential=potential, model=model, **kwargs)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/step_methods/hmc/nuts.py”, line 154, in init
super(NUTS, self).init(vars, **kwargs)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/step_methods/hmc/base_hmc.py”, line 75, in init
vars, blocked=blocked, model=model, dtype=dtype, **theano_kwargs
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/step_methods/arraystep.py”, line 228, in init
vars, dtype=dtype, **theano_kwargs)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/model.py”, line 710, in logp_dlogp_function
return ValueGradFunction(self.logpt, grad_vars, extra_vars, **kwargs)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/pymc3/model.py”, line 443, in init
grad = tt.grad(self._cost_joined, self._vars_joined)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 605, in grad
grad_dict, wrt, cost_name)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1371, in _populate_grad_dict
rval = [access_grad_cache(elem) for elem in wrt]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1021, in access_term_cache
output_grads = [access_grad_cache(var) for var in node.outputs]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1326, in access_grad_cache
term = access_term_cache(node)[idx]
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/gradient.py”, line 1162, in access_term_cache
new_output_grads)
File “/Users/Colo/anaconda/lib/python2.7/site-packages/theano/scan_module/scan_op.py”, line 2222, in L_op
assert (x[::-1][:-1].tag.test_value.shape[0] == n)
AssertionError