predict.spm {SemiPar} | R Documentation |
Takes a fitted spm
object produced by
spm()
and obtains predictions at new data values.
predict.spm(object,newdata,se,...)
object |
a fitted spm object as produced by spm() . |
newdata |
a data frame containing the values of the predictors at which predictions are required. The columns should have the same name as the predictors. |
se |
when this is TRUE standard error estimates are returned for each prediction. The default is FALSE. |
... |
other arguments. |
Takes a fitted spm
object produced by
spm()
and obtains predictions at new data values
as specified by the `newdata' argument. If `se=TRUE' then
standard error estimates are also obtained.
If se=FALSE then a vector of predictions at `newdata' is returned. If se=TRUE then a list with components named `fit' and `se' is returned. The `fit' component contains the predictions. The `se' component contains standard error estimates.
M.P. Wand wand@maths.unsw.edu.au (other contributors listed in SemiPar Users' Manual).
Ruppert, D., Wand, M.P. and Carroll, R.J. (2003)
Semiparametric Regression Cambridge University Press.
http://stat.tamu.edu/~carroll/semiregbook/
Ganguli, B. and Wand, M.P. (2005)
SemiPar 1.0 Users' Manual.
http://www.maths.unsw.edu.au/~wand/papers.html
spm
lines.spm
plot.spm
summary.spm
library(SemiPar) data(fossil) attach(fossil) fit <- spm(strontium.ratio~f(age)) newdata.age <- data.frame(age=c(90,100,110,120,130)) preds <- predict(fit,newdata=newdata.age,se=TRUE) print(preds) plot(fit,xlim=c(90,130)) points(unlist(newdata.age),preds$fit,col="red") points(unlist(newdata.age),preds$fit+2*preds$se,col="blue") points(unlist(newdata.age),preds$fit-2*preds$se,col="green")