[R] can I get standard error from predict.gam()?
Andrew Robinson
A.Robinson at ms.unimelb.edu.au
Wed Nov 15 22:54:09 CET 2006
Hi Kevin,
sadly, you have (again) not provided enough information to enable us
to make a constructive response - for example, what package are you
using? I assume that it's not mgcv, as mgcv does have a predict.gam
with an se.fit argument:
predict.gam package:mgcv R Documentation
Prediction from fitted GAM model
Description:
...
Usage:
## S3 method for class 'gam':
predict(object,newdata,type="link",se.fit=FALSE,terms=NULL,
block.size=1000,newdata.guaranteed=FALSE,na.action=na.pass,...)
But -- maybe I'm wrong.
Furthermore, we can't tell from your note in what way se=T or se.fit=T
is "no use". What is the error message? Please try to read the
posting guide http://www.R-project.org/posting-guide.html and provide
commented, minimal, self-contained, reproducible code.
Cheers,
Andrew
On Wed, Nov 15, 2006 at 04:35:21PM -0500, seeTigers wrote:
> Hi everybody,
> I am using predict.gam() now. I but it seems there is no such option to get
> standard errors of the predicted values. I tried to set se=T or se.fit=T but
> no use.
>
> If you know anything about that please let me know. Thanks very much.
> Kevin.
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
--
Andrew Robinson
Department of Mathematics and Statistics Tel: +61-3-8344-9763
University of Melbourne, VIC 3010 Australia Fax: +61-3-8344-4599
http://www.ms.unimelb.edu.au/~andrewpr
http://blogs.mbs.edu/fishing-in-the-bay/
More information about the R-help
mailing list