In stla/ocpuHello:

library(ggplot2)

dd <- setNames(dat, c("x","y"))
fit <- lm(y~x, data=dd)

Summary

summary(fit)

Confidence intervals

out <- confint(fit, level=conflevel)
rownames(out) <- c("Intercept", "slope")
pander::pandoc.table(out)

ggplot(dd, aes(x=x, y=y)) + geom_point() + geom_smooth(method="lm") +
  xlab(colnames(dat)[1]) + ylab(colnames(dat)[2]) +
  ggtitle(sprintf("Plot of %s vs %s.", colnames(dat)[2], colnames(dat)[1]))

stla/ocpuHello documentation built on May 30, 2019, 5:47 p.m.