set.seed(129)
n <- 7; p <- 2
X <- matrix(rnorm(n * p), n, p) # no intercept!
y <- rnorm(n)
obj_1 <- .lm.fit (x = cbind(1, X), y = y)
obj_2 <- lm(y~X)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.