自定义函数线性回归

Question

I am trying to implement the following "R" code in python:

fit = lm(log(y) ~ log(x1) + log(x2) + 
x3, data=data);

I know in sklearn, you can make a linear regression with multiple variables. However, I specifically want to make the formula above.

Any guidance would be appreciated.

Answer 1

Apply a log transformation to x1 and x2 and then run the linear regression:

import numpy as np
from sklearn.linear_model import LinearRegression

log_x1 = np.log(x1)
log_x2 = np.log(x2)

log_y = np.log( y)

log_model = LinearRegression().fit( np.c_[log_x1, log_x2, x3], log_y)