我想在R中用一个因变量进行多项式回归 y 和两个独立的变量 x1 和 x2。在我看来，模型应如下所示，
y=b0 + b1x1+ b2x2+ b3x1^2+ b4x2^2+ b5x1x2
我试过了 lm(y~x1+x2+poly(x1,2,raw=TRUE)+poly(x2,2,raw=TRUE)) 并且 lm(y~x1+x2+I(x1^2)+I(x2^2))。但这只给出了正方形，而不是两个变量的乘积。
我能做到的 lm(y~x1+ x2+ x1^2+ x2^2+ x1x2)。但我想知道是否有一个比写出整个方程更容易的方法。我也想做3和4型号的力量，这是更长的。