a=rnorm(300)
b=runif(300,-1,1)
x1=a[1:200]
x2=b[1:200]
y=1+2*x1+sin((x2^2)/3)+0.3*rnorm(200)
mp1.gam<-lm(y~x1+x2)
s=data.frame(x1,x2)
s1<-predict(mp1.gam,s)
天啊
这网速
刚才回复没看到
现在去信箱找你的信
居然没了
我网不好
也下不到相应的包
我就用OLS做了一下
就在上面.
你可以看看
结论:
你的数据框s=data.frame(Y,X1,X2)
貌似多了个Y
你把他搞掉在试试