I am estimating a bunch of linear models using `GLM.lm(X, y)`

, where `X`

is a matrix. I have column names for `X`

from a database, and I would like to have the result from `fit`

/`lm`

to include them. MWE:

```
using GLM
X = hcat(ones(100), randn(100, 2))
y = X * [1, 2, 3]
reg = lm(X, y)
colnames = ["intercept", "age", "distance"]
```

so that

```
julia> reg
LinearModel{GLM.LmResp{Array{Float64,1}},GLM.DensePredChol{Float64,LinearAlgebra.Cholesky{Float64,Array{Float64,2}}}}:
Coefficients:
─────────────────────────────────────────────────────────────────────
Estimate Std. Error t value Pr(>|t|) Lower 95% Upper 95%
─────────────────────────────────────────────────────────────────────
x1 1.0 9.41024e-17 1.06267e16 <1e-99 1.0 1.0
x2 2.0 1.09949e-16 1.81903e16 <1e-99 2.0 2.0
x3 3.0 9.99446e-17 3.00166e16 <1e-99 3.0 3.0
─────────────────────────────────────────────────────────────────────
```

and I would like to replace `x1`

, `x2`

, `x3`

with column names in `colnames`

.