各位大佬，我想问一下，我做GAM模型，里面包含了三个虚拟变量，最终出图的时候只能出来平滑项的图，虚拟变量的响应图怎么能出来呢？这是结果

Family: gaussian
Link function: identity

Formula:
sl ~ s(x, y) + s(Asp, bs = "cr") + s(Slo, bs = "cr") +
    s(abs, bs = "cr") + s(sl.1, bs = "cr") + s(sl.1.1,
    bs = "cr") + s(abs.1, bs = "cr") + T + s(abs.1.1,
    bs = "cr") + s(Mix, bs = "cr") + s(Dec, bs = "cr") +
    s(Eve, bs = "cr") + s(Bro, bs = "cr") + s(Riv,
    bs = "cr") + s(Cro, bs = "cr") + s(Res, bs = "cr") +
    ID + age + T

Parametric coefficients:
             Estimate Std. Error t value Pr(>|t|)   
(Intercept)  1.251715   0.008787 142.445  < 2e-16 ***
Tnight      -0.018519   0.010399  -1.781 0.074975 .  
ID3         -0.050561   0.019854  -2.547 0.010898 *  
ID4         -0.087514   0.025059  -3.492 0.000482 ***
ID18         0.459881   0.016350  28.128  < 2e-16 ***
ID19         0.309096   0.021335  14.488  < 2e-16 ***
agesub       0.482738   0.014407  33.506  < 2e-16 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Approximate significance of smooth terms:
             edf Ref.df       F  p-value   
s(x,y)     4.452  6.010   0.955 0.456077   
s(Asp)     6.883  7.899   2.529 0.008041 **
s(Slo)     1.000  1.000   8.562 0.003443 **
s(abs)     5.683  6.808   8.404  < 2e-16 ***
s(sl.1)    6.706  7.732 163.949  < 2e-16 ***
s(sl.1.1)  6.722  7.747 191.112  < 2e-16 ***
s(abs.1)   5.745  6.874  14.881  < 2e-16 ***
s(abs.1.1) 6.479  7.595   5.007 5.84e-06 ***
s(Mix)     8.799  8.981   2.954 0.001799 **
s(Dec)     8.253  8.807   3.159 0.000908 ***
s(Eve)     8.031  8.688   1.532 0.086646 .  
s(Bro)     7.840  8.460   1.671 0.070804 .  
s(Riv)     8.891  8.993   2.705 0.004099 **
s(Cro)     1.000  1.000   0.045 0.831671   
s(Res)     7.506  8.196   1.743 0.078944 .  
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

Rank: 161/162
R-sq.(adj) =  0.481   Deviance explained = 48.8%
GCV = 0.17994  Scale est. = 0.17746   n = 7240