假如你已将将Y x1 x2 x3的数据导入R了
lnY=ln(Y)
lnx1=ln(x1)
lnx2=ln(x2)
lnx3=ln(x3)
reg lnY c lnx1 lnx2 lnx3
lm.sol <- lm( lnY~1+lnx1+lnx2+lnx3)  ##lm()函数返回拟合结果的对象，可以用summary()函数查看其内容。
summary(lm.sol)

这样做的原因只要是 lny=lnx.  是源于y=x1e^x2  此时因为x在指数幂次和乘积上所以无法估计，因此等号两边同时取ln ，形成参数线性回归方程，进而对参数进行估计。另外，ln表示弹性的意思，经济学上有讲，可以查一下。