我正在尝试使用“ mgcv”在R中创建变化系数GAM，其中针对不同水平的因子生成变化系数项。连续变量在空间上是显式的变系数项如下：
s(longitude, latitude, by=continuous_variable)
但是我找不到添加第二个“ by”变量以产生所需输出（例如）的方法：
s(longitude, latitude, by=continuous_variable, by=factor) ....据我所知
有没有一种方法可以针对不同水平的因子创建变化的系数项？
类似地，不是对不同因子水平进行全局参数估计，而是可能对因子水平具有空间上明确的变化系数？例如，如果我想估计诱饵类型对鱼类捕获的影响，则可以包括一个全局因素：
Catches ~ fBAIT_TYPE + ....
但是是否可以使fBAIT_TYPE在空间上具有明显的变化系数？如果只有两个级别（例如A鱼和B鱿鱼），是否可以查看B级别的系数相对于A在空间上是否发生变化？感兴趣的目标物种可能总是喜欢鱿鱼而不是鱼类（全局因子），或者对于鱿鱼而言，其相对偏好在空间上是明确的（系数变化）。
任何帮助将不胜感激！谢谢！