各位大佬,我想问一下,我做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