n <- 100
x <- rnorm(n)
y <- 2*x + rnorm(n)
out <- lm(y ~ x)
library(pander)
panderOptions("digits", 2)
pander(out)
Fitting linear model: y ~ x
  Estimate Std. Error t value Pr(>|t|)
x 2.1 0.099 21 1.9e-38
(Intercept) 0.046 0.11 0.43 0.66