title: Extracting z values from hurdle model output, 'pscl' package in R
body: <p>I am using the <code>hurdle()</code> function from the <code>pscl</code> package to run negative binomial hurdle models in R.</p>
  <p>I would like to extract only the list of z values from the model output and assign those values to a vector. But I cannot figure out how to do this?</p>
  <p>Here is the basic outline of my model:</p>
  <pre><code>real.mod &lt;- summary(hurdle(y ~ x1 + x2 + x3 + x4  |  x1 + x2 + x3 + x4, data=mydata, dist = &quot;negbin&quot;, zero.dist = &quot;negbin&quot;, link = &quot;logit&quot;))
  
  real.mod 
  
  </code></pre>
  <p>For example, I have tried:</p>
  <p><code>coef(summary(real.mod))[,3]</code></p>
  <p>But this doesn't work.</p>
  <p>Any guidance on how to do this would be greatly appreciated.</p>