similar to: GAM: Remedial measures

I am analyzing particulate matter data (PM10) on a small data set (147 observations). I fitted a semi-parametric model and am worried about overfitting. How can one check for model fit in GAM? Jean G. Orelien

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20070404/0f61f54a/attachment.pl

Dear list members, while I appreciate the possibility to deal with overdispersion for count data either by specifying the family argument to be quasipoisson() or negative.binomial(), it estimates just one overdispersion parameter for the entire data set. In my applications I often would like the estimate for overdispersion to depend on the covariates in the same manner as the mean. For example,

Hi, there are a lot of messages dealing with overdispersion, but I couldn't find anything about how to test for overdispersion. I applied a GAM with binomial distribution on my presence/absence data, and would like to check for overdispersion. Does anyone know the command? Many thanks, Anna -- View this message in context:

I am new to R. I'm using the function GAM and wanted to get standard errors and p-values for the parametric terms (I fitted a semi-parametric models). Using the function anova() on the object from GAM, I only get p-values for the nonparametric terms. Does anyone know if and how to get standard errors for the parametric terms? Thanks. Jean G. Orelien

Dear R-users: Im trying (using gam package) to develop a stepwise analysis. My gam object contains five pedictor variables (a,b,c,d,e,f). I define the step.gam: step.gam(gamobject, scope=list("a"= ~s(a,4), "b"= ~s(b,4), "c"= ~s(c,4), "d"= ~s(d,4), "e"= ~s(e,4), "f"= ~s(f,4))) However, the result shows a formula containing the whole

How can I obtain the values of confidence intervals from gam anf glm objects? Thanks in advance -- David Nogu?s Bravo Functional Ecology and Biodiversity Department Pyrenean Institute of Ecology Spanish Research Council Av. Monta?ana 1005 Zaragoza - CP 50059 976716030 - 976716019 (fax)

Get some question about the function "gam". Suppose I have a semiparametric model, Y~x1+x2+s(z1). Using "gam", how could I get the estimates for the parametric part and nonparametric part respectively? And another question: we could find the coefficients for both parametric term and nonparametric term, what do these coefficients for the nonparametric term stand for, the

hello, I am looking for a step() function for GAM's. In the book Statistical Computing by Crawley and a removal of predictors has been done "by hand" model <- gam(y ~s(x1) +s(x2) + s(x3)) summary(model) model2 <- gam(y ~s(x2) + s(x3)) # removal of the unsignificant variable #then comparing these two models if an significant increase occurs. anova(model, model2,

Hi, I've got the R library "mgcv" for GAM written by Simon Wood which works well in many instances. However, over the years I got attached to the S+ implementation of GAM which allows loess smoothing in more than 1 dimension as well as spline smoothing. Has anyone ported the S+ GAM library to R? Regards, Doug Beare. Fisheries Research Services, Marine Laboratory, Victoria Road,

Hi, Is there any way to change the y axis range of values in a plot.gam()? I need that two different GAM plots to be of the same scale. Also, it is possible to change the labels? I tried with "ylab" and "ylim" and did not work Thanks in advance Ricardo Lopes Ricardo Lopes ............................................. Instituto do Mar Departamento de Zoologia

All, Sorry for bombarding you with GAM related questions, but... I know a partial residual option in plot.gam() is on Simon Wood's todo list, but since I'm in the midst of a project and not yet having acquired sufficient R knowledge to code something usable myself I'll have to put my trust in you. Anybody got some code lying around for doing this? Or if someone can supply me with

hi, i'm working with mgcv packages and specially gam. My exemple is: >test<-gam(B~s(pred1)+s(pred2)) >plot(test,pages=1) when ploting test, you can view pred1 vs s(pred1, edf[1] ) & pred2 vs s(pred2, edf[2] ) I would like to know if there is a way to access to those terms (s(pred1) & s(pred2)). Does someone know how? the purpose is to access to equation of smooths terms

I fitted a GAM model with Poisson distribution using the function gam() in the mgcv package. My model is of the form: mod<-gam(y~s(x0)+s(x1)+s(x2),family=poisson). To extract estimates at a specified set of covariate values I used the gam `predict' method. But I want to get estimate and standard error of the difference of two fitted values. Can someone explain what should I do? Thank

Dear all, I'm a fairly new R user having two questions regarding gam: 1. The prediction example on p. 38 in the mgcv manual. In order to get predictions based on the original data set, by leaving out the 'newdata' argument ("newd" in the example), I get an error message "Warning message: the condition has length > 1 and only the first element will be used in: if

Hi all! I?ve got some problems with the function gam (library mgcv). For some models I get the error message : Error: no valid set of coefficients has been found:please supply starting values In addition: Warning message: NaNs produced in: log(x) This is a shortened code I used: gam(y ~ M1 + M3 + M4 + M5 + M6 + sex + M1*M3 + s(age), family=Gamma(link ="identity"), weights=days) If

Hi.. I'm using gam() to fit a spline model for a data set that has two predictor variables (say A and B). The results indicate that the higher order interaction terms are significant. The R^2 jumps from .5 to .9 when I change the maximum order for the interaction from 10 to 15 (i.e. (AB)^10 to (AB)^15). Is there any way of finding out which of the terms in the model are really

hello to everybody: I would to know with ranges of GCV or UBRE values can be considered as adequate to consider a GAM as correct Thanks in advance -- David Nogu?s Bravo Functional Ecology and Biodiversity Department Pyrenean Institute of Ecology Spanish Research Council Av. Monta?ana 1005 Zaragoza - CP 50059 976716030 - 976716019 (fax)

I have a question about how to compare a GLM with a GAM model using anova function. A GLM is performed for example: model1 <-glm(formula = exitus ~ age+gender+diabetes, family = "binomial", na.action = na.exclude) A second nested model could be: model2 <-glm(formula = exitus ~ age+gender, family = "binomial", na.action = na.exclude) To compare these two GLM

Hola! I have some problems with gam in mgcv. Firts a detail: it would be nice igf gam would accept an na.action argument, but that not the main point. I want to have a smooth term for time over a year, the same pattern repeating in succesive years. It would be natural then to impose the condition s(0)=s(12). Is this possible within mgcv? I tried to obtain this with trigonometric terms, aca: