请教如何在广义可加模型（gam）中加入时间或者截面固定效应的估计

请教各位：
       我想用广义可加模型去拟合几个预测变量对一个反应变量的影响关系。由于数据结构是面板，也即同时包含了时间和截面两个维度的数据框。我想在模型中控制固定效应的影响。
问题1：是不是可以通过放入时间或者截面的因子进去实现：

modelfit<-gam(ly~s(real,bs="cr",k=37)+s(krate,bs="cr",k=37)+s(gdpreal, bs="cr",k=37)+regionfactor+yearfactor,family=gaussian(link=identity),

data= gugz,gamma=1.4)

其中，regionfactor和yearfactor是我定义的区域和时间因子，s(,,,)是模型的非参数部分

问题2：如果问题1可以实现，那么可不可以同时测算时间固定效应和截面固定效应，不会出现哪种低级的常识性错误吧？

请大家帮忙了，坐等回复，感激不尽！

这个问题难度好高啊

没有合适答案怎么处理啊