我自己制作了60个数据如下，用广义加性模型之后有了以下的图，但不知如何解释，请教高手了：

应变量Y 1.5乘N X2的平方 X3方立 X4乘以5 随机变量*10

282.5874871 1.5 9 216 50 5.534079179

434.0977677 3 25 343 55 7.361607001

609.1677622 4.5 25 512 60 6.970692913

828.8728258 6 25 729 65 3.520750693

1105.338709 7.5 25 1000 70 2.58064455

1449.579293 9 25 1331 75 8.70844773

1848.852401 10.5 25 1728 80 4.865818735

2320.51964 12 25 2197 85 1.38149075

2879.538943 13.5 25 2744 90 6.399039327

3519.999912 15 25 3375 95 9.0908289

4244.50342 16.5 25 4096 100 6.366745152

5064.147181 18 25 4913 105 2.861073981

5987.366826 19.5 25 5832 110 0.788023418

7029.596295 21 25 6859 115 8.723904841

8149.154244 22.5 1 8000 120 5.140221646

9415.759973 24 4 9261 125 1.599975146

10816.92702 25.5 9 10648 130 4.024562024

12348.85329 27 16 12167 135 3.502990394

14021.79527 28.5 25 13824 140 3.904789597

15840.45908 30 36 15625 145 4.053712927

17807.16557 31.5 49 17576 150 0.605066127

19942.22324 33 64 19683 155 6.566585958

22230.59887 34.5 81 21952 160 2.81715778

24696.23078 36 100 24389 165 5.664345634

27331.83215 37.5 121 27000 170 3.029228778

30159.29745 39 144 29791 175 9.361317255

33166.78084 40.5 169 32768 180 8.437123434

430.2037178 42 196 1 185 5.639743492

476.9700769 43.5 225 8 190 9.518251765

532.3042691 45 256 27 195 8.458426428

604.0537992 46.5 289 64 200 4.139817453

507.734097 48 324 125 5 5.212815458

645.8374211 49.5 361 216 10 8.488564595

811.6825603 51 400 343 15 2.438691212

1031.834893 52.5 441 512 20 5.758993407

1297.463266 54 484 729 25 4.966605354

1624.21223 55.5 529 1000 30 8.829299788

2001.97129 57 576 1331 35 2.701172874

2452.423796 58.5 625 1728 40 0.83981419

2986.355787 60 676 2197 45 7.596170108

3588.732256 61.5 729 2744 50 3.847505156

4281.324553 63 784 3375 55 3.931412266

5062.029829 64.5 841 4096 60 0.481662416

5947.817369 66 900 4913 65 3.470335356

6931.72539 67.5 961 5832 70 1.113991176

8033.11115 69 1024 6859 75 5.55559082

9247.115554 70.5 1089 8000 80 6.923230804

10579.11515 72 1156 9261 85 4.650131975

12039.40292 73.5 1225 10648 90 2.63901496

13639.15996 75 1296 12167 95 5.599966634

15370.76132 76.5 1369 13824 100 1.146658615

17252.51216 78 1444 15625 105 0.465599063

19293.0136 79.5 1521 17576 110 5.921452802

21489.52891 81 1600 19683 115 9.57173934

23844.06026 82.5 1681 21952 120 7.782054418

26371.63309 84 1764 24389 125 8.757355583

29067.06032 85.5 1849 27000 130 2.327563141

31949.21552 87 1936 29791 135 0.195930038

35027.35387 88.5 2025 32768 140 5.321698888

38297.53222 90 2116 35937 145 8.665653364

代码如下：

rm(list=ls())

data<-read.table("d:/广义加性模型试验.csv",head=TRUE,sep=",")

head(data)

colnames(data)<-c("y","x1","x2","x3","x4","x5")

data<-as.data.frame(data)

data

head(data)

library(mgcv)

library(nlme)

whole.data<-gam(y~s(x1,bs="ts")+s(x2,bs="ts")+s(x3,bs="ts")

+s(x4,bs="ts")+s(x5,bs="ts"),

data=data,fimily="poisson")

par(mfrow=c(3,2),mar=c(4,4,0.5,0.5))

plot(whole.data,scale=0,pages=0,select=1,

xlab="x1",ylab="y",

residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

plot(whole.data,scale=0,pages=0,select=2,

xlab="x2",ylab="y",

residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

plot(whole.data,scale=0,pages=0,select=3,

xlab="x3",ylab="y",

residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

plot(whole.data,scale=0,pages=0,select=4,

xlab="x4",ylab="y",

residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

plot(whole.data,scale=0,pages=0,select=5,

xlab="x5",ylab="y",

residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

五个图如下：[s:12]

[s:12][s:12]，图怎么上不去，但仍要请高手解释，图的含义。

多谢了！