data(air)
y=air$ozone # response
X=as.matrix(air[,3:4]) # single index term ;
Z=air[,2] # partially linear term ;
result <- gplsim(y,X,Z=Z,family = gaussian,k=10)
result$theta
result$coefficients
summary(result)
# Or you can try different spline basis
result <- gplsim(y,X,Z=Z,family = gaussian,bs="tp",k=10)
result$theta
result$coefficients
summary(result)
Run the code above in your browser using DataLab