[img]file:///C:\Documents and Settings\fang\Application Data\Tencent\Users\153421304\QQ\WinTemp\RichOle\7P{V[_UG4DUXF88KT2%4]PM.jpg[/img] 这个是我要做的回归方程，其中β我想用可变系数来估计，要用到mgcv包。    下面是回归用到的数据。        执行ct1<-gam(s~s(b-1)+Yvar+Gvar,data=mydata)；         plot(ct1)
之后出来的图像横坐标并不是时间，而是b，请问这个代码应该怎么写。谢谢啦。




                   s       b            Gvar          Yvar
1980 -0.02400 0.22577   3.06121   2.023190
1981 -0.07578 0.26739   1.05229   1.399608
1982 -0.05751 0.29374  -0.54425   0.792793
1983 -0.06503 0.33665  -2.14653   0.246519
1984 -0.07638 0.40149  -0.84130   0.051629
1985 -0.10312 0.46724   0.87540   0.084501
1986 -0.09240 0.47246   0.50350  -0.467000
1987 -0.08628 0.52528  -2.84640  -1.385860
1988 -0.10394 0.57195   0.23400  -1.729390
1989 -0.12981 0.59952   2.83660  -1.811040
1990 -0.14506 0.73316   1.33500  -0.060970
1991 -0.10370 0.74847   3.01720  -0.894230
1992 -0.11504 0.80144   1.92420   0.038730
1993 -0.12397 1.00510  -2.01360   0.779420
1994 -0.09016 0.98512  -2.29960  -0.760040
1995 -0.06986 0.99209  -2.84640   0.191660
1996 -0.06791 1.01603  -3.50640  -0.010550
1997 -0.06024 0.98736  -3.06770   0.610330
1998 -0.03905 0.96581  -3.25720  -0.902580
1999 -0.03137 1.02509  -3.63620  -2.676290
2000 -0.03692 1.03442  -2.80130   1.164250
2001 -0.04328 1.03716  -2.27500  -1.560180
2002 -0.04805 1.01451  -2.62750  -2.836410
2003 -0.05668 0.97270   0.99960  -1.930700
2004 -0.07398 0.98844   3.02200  -0.336920
2005 -0.05298 1.00291   2.47140  -4.516660
2006 -0.06121 1.06107   7.47520   0.029900
2007 -0.06686 1.05412  12.18530   5.639700
2008 -0.09804 1.10721  12.15270  12.980400
2009 -0.15513 1.27100   6.75950  15.766300
2010 -0.10420 1.42757  -1.53070   2.270800
2011 -0.07982 1.65559 -10.21970  -8.338200
2012 -0.06874 1.89149 -13.44530 -13.852700