Function to calculate CI of linear regression coefficients

Amateur code

```

 jad<-function(x, y) {

  model<-lm(y~x)

  std.err<-coef(summary(model))[, 2]

  coef.model1<-coef(summary(model))[, 1]

  upper.ci<-coef.model1+1.96*std.err

  lower.ci<-coef.model1-1.96*std.err

  print(upper.ci)

  print(lower.ci)

}

```

Pro code

jad <- function(x, y) {

  `colnames<-`(

    coef(summary(lm(y ~ x)))[,1:2] %*% matrix(c(1, 1.96, 1, -1.96), nrow = 2), c("upper.ci", "lower.ci")

  )

}