title: Unprecise p-values in Stargazer
body: <p>I want the same stars for significancies in regression output in stargazer as in the "normal output".</p>
  
  <p>I produce data</p>
  
  <pre><code>library("stargazer"); library("lmtest"); library("sandwich")
  set.seed(1234)
  df &lt;- data.frame(y=1001:1100)
  df$x &lt;- c(1:70,-100:-71) + rnorm(100, 0, 74.8)
  model &lt;- lm(log(y) ~ x, data=df)
  </code></pre>
  
  <p>and get some model estimates where the coefficient on x has a p-value of 0.1023</p>
  
  <pre><code>coeftest(model, vcov = vcovHC(model, type="HC3"))
  </code></pre>
  
  <p>I want to have these results in LaTeX. Based on the same function I calculate heteroscedasticity consistent standard estimates and let stargazer use them.</p>
  
  <pre><code>stderr_HC3_model &lt;- sqrt(diag(vcovHC(model, type = "HC3")))
  stargazer(model, se=list(stderr_HC3_model))
  </code></pre>
  
  <p>The stargazer output has a star at the coefficient indicating significance when alpha=10%. I want stargazer to give the same as the coeftest. (Because of the comparability with Stata where reg L_y x, vce(hc3) gives exactly the coeftest results.)</p>
  
  <p>I played around with the stargazer options p.auto, t.auto which did not help. When I execute "stargazer" I cannot view the underlying code as it is possible in other cases. What to do?</p>
  
  <hr>
  
  <p>Richards answer helped me. I indicate the steps I used to give out more than one regression (let's say ols_a and ols_b).</p>
  
  <pre><code>ses &lt;- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,2],
          coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,2])
  pvals &lt;- list(coeftest(ols_a, vcov = vcovHC(ols_a, type="HC3"))[,4],
            coeftest(ols_b, vcov = vcovHC(ols_b, type="HC3"))[,4])
  stargazer(ols_a, ols_b, type="text", p=pvals, se=ses)
  </code></pre>