R中mgcv包画图——单位x对应y改变量图

我是看了最近发表在《Environmental Health》上的文章，采用的是R中的mgcv包，目的是探求血压和温度的非线性关系，我利用找到数据在R中学习操作有些问题。
ct=gam(Blood pressure~s(Temp,k=5)+s(BMI,k=4)+s(Age,k=6)+s(RH,k=4)+s(Pa,k=5)+Gender+Season,family=gaussian,data=mydata)                                                                                                                                                                                           （公式1）

1. 如何控制个体特征变量：性别（0,1）是分类变量，是否直接将其当成线性关系放进等式里，如公式1所示。年龄是连续性变量，用s平滑。是这样吗？我是把季节（1,2,3,4）当作分类变量，以线性关系放入公式控制的，这样作对吗？

2. 我在R中采用公式1，然后用plot（ct）画图，得到图1，这个纵坐标s（Tmean，3.87）是应变量血压吗？如果图1纵坐标为血压的话，我看文章中画出了降低单位温度血压改变值的图（图2）。想得到纵坐标为血压改变值的图是否就是求图1的导数，但是求导的程序代码一直尝试未果。这个猜想不知道对不对，希望各位大侠给出得到图2的具体代码，不甚感激。

1 个解决方案