Introduction

Let's simulate some data

x <- rnorm(100)
y <- x + rnorm(100, sd = 0.5)

Here is a scatter plot of the above data

par(mar = c(5, 4, 1, 1), las = 1)
plot(x, y, main = "My Simulated Data")

Making Tables with xtable

library(datasets)
data(airquality)
fit <- lm(Ozone ~ Wind + Temp + Solar.R, data = airquality)

Here is table of regression coefficients

library(xtable)
xt <- xtable(summary(fit))
print(xt, type = "html")