search for: gamobject

...le 2, but it shouldn't be. Once I figure out how to do this, I'll take the fitted coefficients from each model and average them, then take the vcv's from each model and average them, and add a correction to account for within and between imputation variability, then plug them into a gamObject$coeffient and gamObject$Vp matrix, plot/summarize, and have my result. Comments welcome on whether or not this would be somehow incorrect would be welcome as well. Still have a lot to learn! Thanks, Andrew [[alternative HTML version deleted]]

Hi, suppose i have a gam object gamobject<- gam( Y~ s(X1)+s(X2)+ X3) I would like to extract the predicted partial effect of X3 but selecting it by its name, as it's to be included in a function and i don't always know the exact position of X3. something like predict(gamobject,type="terms")["X3",]. But th...

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 model, without eliminating any variable. Anybody can help me? Thanks -- David...

...text html latex gam.selection text html latex example gam.side text html latex example gam.vcomp text html latex example gam2objective text html latex gamObject text html latex Out of memory! ERROR: building help failed for package ‘mgcv’ * Removing ‘/usr/lib/R/library/mgcv’ * Restoring previous ‘/usr/lib/R/library/mgcv’ The downloaded packages are in ‘/tmp/Rtmp9H7LN8/downloaded_packages’ Updating HTML index of packages i...

Dear users, I am trying to show the equation (including coefficients from the model estimates) for a gam model but do not understand how to. Slide 7 from one of the authors presentations (gam-theory.pdf URL: http://people.bath.ac.uk/sw283/mgcv/) shows a general equation log{E(yi )} = ?+ ?xi + f (zi ) . What I would like to do is put my model coefficients and present the equation used. I am an

Hello evyrone, I would be grateful if you could help me in (I hope) simple problem. I fit a gam model (from mgcv package) with several smooth functions . I don't know how to extract values of just one smooth function. Can you please help me in this? Kind regards, Daniel Rabczenko

Dear all, I am analysing growth data - response variable - using GAM and GAMM models, and 4 covariates: mean size, mean capture year, growth interval, having tumors vs. not The models work fine, and fit the data well, however when I try to compare models using AIC I cannot get an AIC value. This is the code for the gam model:

My data set consist of number of calls (lcin) across Day. I am looking for activity differences between three features (4 sites per feature). I am also looking for peaks of activity across time (Day). I am using a gamm since I believe these are nonlinear trends with nested data. gammdata<-gamm(lcin~Temp+s(Day)+fType+wind+fFeature+Forest+Water+Built, list=fSite,data=data, family=gaussian)

Hello, I'd like to understand 'what' is predicting the response for library(mgcv) gam? For example: library(mgcv) fit <- gam(y~s(x),data=as.data.frame(l_yx),family=binomial) xx <- seq(min(l_yx[,2]),max(l_yx[,2]),len=101) plot(xx,predict(fit,data.frame(x=xx),type="response"),type="l") I want to see the generalized function(s) used to predict the response