用R语言做限制性立方样条的交互作用,如何做呢?

我已经用R语言做了2个连续变量x1 x2分别和y的限制性立方样条函数,

还想把x1 x2结合在一起,想做x1 x2的交互作用,这个代码该怎么输入呀?

mgcv包可以用吗?

希望能得到老师的解答帮助!