I have a data frame with 10 columns(features) COL0 to COL9 and a column RESP. How do I calculate a LinearRegression Model for each pair COL0 to COL9 ~ RESP?
I am expecting to get 10 graphs showing the Model and also a table with the coefficients of my model for each column.

What I tried do far:

x = df.Col0
y = df.RESP
data = DataFrame(X=x,Y=y)
model = lm(@formula(Y~X),data)

And I get what i want for this first pair COl0 ~Resp. But still i need to plot it.

Now i need to do the same steps more 9 times. I want to optimize this as i can have over 100 columns.

I am new to Julia and I really dont have a clue how to get this done. Any help?