Posterior sampling after model fitting

  In [12]:
    # posterior predictive distribution we can get samples from
    trace_pred = TracePredictive(wrapped_model,
                                 posterior,
                                 num_samples=1000)
    post_pred = trace_pred.run(x_data, None)

How do I generate samples with a different size? I tried changing “x_data -> x_data[:70,:]” in In [12] which seems to be reasonable to me but returns an error complaining about the shape of the input.

Also, I guess there is a typo in the ‘Linear Regression’ paragraph earlier in the example:
Our input X is a matrix of size N×2 and our output y is a vector of size 2×1.
Input is N x 2, then the output should be N x 1, because the weight matrix, as you point out later, is a 2 x 1 matrix (linear.weight [[-1.90511 -0.18619268]]). (N x 2) times (2 x 1) = (N x 1)