title: Bayes predict, subscript out of bounds
body: <p>I'm having some problems with the predict function when using bayesglm. I've read some posts that say this problem may arise when the out of sample data has more levels than the in sample data, but I'm using the same data for the fit and predict functions. Predict works fine with regular glm, but not with bayesglm. Example:</p>
  
  <pre><code>control &lt;- y ~ x1 + x2
  
  # this works fine:
  glmObject &lt;- glm(control, myData, family = binomial())
  predicted1 &lt;- predict.glm(glmObject , myData, type = "response")
  
  # this gives an error: 
  bayesglmObject &lt;- bayesglm(control, myData, family = binomial())
  predicted2 &lt;- predict.bayesglm(bayesglmObject , myData, type = "response") 
  Error in X[, piv, drop = FALSE] : subscript out of bounds
  
  # Edit... I just discovered this works. 
  # Should I be concerned about using these results?
  # Not sure why is fails when I specify the dataset
  predicted3 &lt;- predict(bayesglmObject, type = "response")
  </code></pre>
  
  <p>Can't figure out how to predict with a bayesglm object. Any ideas? Thanks!</p>