similar to: GAM selection error msgs (mgcv & gam packages)

Hi We are provisioning some default features to our customers like "automated answer when outside labor time", "rerouting when Subscriber Absent", and so. These are macro calls embedded in key points inside the dialplan. Since not all customers need/want all features and they (the features) are customized, I named them like [macro-feature1-ClientA], [macro-feature4-ClientF],

Hi Could anyone please tell me what to do to resolve this error message? I tried to run a gam with the mgcv package and got the following error: "Error in qr.qty(qrc, sm$S[[1]]): NA/NaN/Inf in foreign function call (arg 5)" (I have 116 covariates, I'm using the "cr" basis to speed things up, the binomial family and, where necessary, have set the required k to lower than

Dear Trevor, I've noticed a discrepancy in the degrees of freedom reported by gam() from the gam package in R vs. gam() in S-PLUS. The nonparametric df differ by 1; otherwise (except for things that depend upon the df), the output is the same: --------- snip ------------ *** From R (gam version 0.93): > mod.gam <- gam(prestige ~ lo(income, span=.6), data=Prestige) >

I am trying to compare and contrast the smoothing in the {mgcv} version of gam vs. the {gam} version of gam but I get a strange side effects when I try to alternate calls to these routines, even though I detach and unload namespaces. Specifically when I start up R the following code runs successfully until the last line i.e. plot(g4,se=TRUE) when I get "Error in dim(data) <- dim :

Hi all, Has anyone tried to plot math symbols with vis.gam() - a function that uses (I think) plot.gam() internally? This is from the package mgcv - I am new to gam(). For example, with the normal plot(), the expression() works fine for Math symbols: plot(xx, ylab=expression(paste(Delta,"H",sep=""))) # works fine But, with vis.gam(), the same

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

Sometimes, for specific models, I get this error from predict.gam in library mgcv: Error in complete.cases(object) : negative length vectors are not allowed Here's an example: model.calibrate <- gam(meansalesw ~ s(tscore,bs="cs",k=4), data=toplot, weights=weight, gam.method="perf.magic") > test <- predict(model.calibrate,newdata) Error in

Hallo I need some help with the output provided by gam.check after a gam fit (using the package mgcv). To give a brief description of my data, I have claims: a vector of values, which include NA's and one large negative value - otherwise all positive (55 values in total that are not NA). origin: a factor with 10 levels j : taking the values 1, 2, ...., 10 I have fitted a gam, with: >

Dear R-list members, Once a gam (package gam) model has been fitted with family=poisson, is there some R function that could extract the diagonal elements of the smoother matrix S, to be used in a cross-validation for the selection of the best smoothing parameter, following equation 3.19 on page 48 of the book T.J. Hastie and R.J. Tibshirani, Generalized Additive Models, Chapman and Hall/CRC,

Hello, I am using the function step.gam() from the 'gam' package (header info from library(help=gam) included below) and have come across some behavior that I cannot understand. In short, I have written a function that 1) creates a dataframe, 2) calls gam() to create a gam object, then 3) calls step.gam() to run stepwise selection on the output from gam(). When I do this, gam()

R-users E-mail: r-help@r-project.org I have a quenstion on "gam()" in "gam" package. The help of gam() says: 'gam' uses the _backfitting algorithm_ to combine different smoothing or fitting methods. On the other hand, lm.wfit(), which is a routine of gam.fit() contains: z <- .Fortran("dqrls", qr = x * wts, n = n, p = p, y = y *

Dear List, I'm just teaching myself semi-parametric techniques. Apologies in advance for the long post. I've got observational data and a longitudinal, semi-parametric model that I want to fit in GAM (or potentially something equivalent), and I'm not sure how to do it. I'm posting this to ask whether it is possible to do what I want to do using "canned" commands

This is in regards to the suggested use of type="lpmatrix" in the documentation for mgcv::predict.gam. Could one not get the same result more simply by using type="terms" and interpolating each term directly? What is the advantage of the lpmatrix approach for prediction outside R? Thanks. -- View this message in context:

As I understand it,  predict.lm(l ,newdata=nd ,interval="confidence") yields confidence bands for the predicted mean of new observations and lm.predict(l ,newdata=nd ,interval="prediction") yields confidence bands for new observations themselves, given an lm object l.   However with regard to {mgcv} although  predict.gam (g ,se.fit=TRUE ,interval= "prediction")

I am struggling to get GAM model predictions from the top models calculated using model.avg in the package "MuMIn". My model looks something like the following: gamp <- gam(log10(y)~s(x1,bs="tp",k=3)+s(x2,bs="tp",k=3)+ s(x3,bs="tp",k=3)+s(x4,bs="tp",k=3)+s(x5,bs="tp",k=3)+ s(x6,bs="tp",k=3)+x7,data=dat,

Hi, Got my issues sorted. Error message solved: I heard from the guy who developed MuMIn and his suggestion worked. "As for the error you get, it seems you are running an old version of MuMIn. Please update the package first." I did (I was only 1 version behind in both R and in MuMIn) and error disappeared! Running quasipoisson GAM in MuMIn: As for my questions on GAM and " to

Hi, I used GAM function with this code: library(mgcv) modgam=gam(Z~X+Y,data=table) pred=predict(modgam,type="response") vis.gam(modgam,view=c("X","Y"),plot.type="contour",zlim=c(1,16)) My problem is on my figure, the predicted values begin by 0 to 16 with an interval equal 2: 1- I want that the predicted values (Z) begin by 1 until 16 with an interval

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

The documentation for predict.gam in library mgcv gives an example of using an "lpmatrix" to do approximate prediction via interpolation. However, the code is specific to the example wrt the number of smooth terms, df's for each,etc. (which is entirely appropriate for an example) Has anyone generalized this to directly generate code from a gam object (eg SAS or C code)? I wanted to

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: