On Sep 11, 2012, at 12:10 PM, exu wrote: > Thanks for the response. But I'm relatively new to R, could you elaborate a > bit? I tried > > P <-predict.gam(model,type=terms") > > But I'm having a hard time to trying to make sense of the resulting data.
And I'm having even greater difficulty because I cannot see any of it. (My prediction: It gave you linear predictors at your original data points, since you did not give it a newdata argument and you only have one term. I suspect that the chosen smoothing basis might have constructed a matrix output.) The help page ?mgcv::predict.gam has worked examples. In fact the first example illustrated what I thought you would be attempting. Working through all the examples on a help page is a great way to learn a function's capabilities. -- David. > > Regards, > > Eugene Xu > LIBREMAX CAPITAL, LLC. > 830 Third Avenue > Suite 801 > NewYork, NY 10022 > (212) 612-1598 (O) > (908) 656-0306 (M) > (646) 380-2689 (F) > > From: David Winsemius [via R] > [mailto:ml-node+s789695n4642822...@n4.nabble.com] > Sent: Tuesday, September 11, 2012 3:00 PM > To: Eugene Xu > Subject: Re: How to get curve information from a gam model > > > On Sep 11, 2012, at 8:08 AM, exu wrote: > >> I'm using library(mgcv) >> >> I did a gam regression like this >> >> model <- gam(formula = y ~ s(x), family="binomial") >> >> Now I'd like to get 2 vectors x and z (with z=s(x)), so that I can implement >> the model outside R. How can I do this? > > ?predict.gam > > -- > > David Winsemius, MD > Alameda, CA, USA > David Winsemius, MD Alameda, CA, USA ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
