similar to: attributing names in predicted type="terms" gam object

Dear all, i generated several GLM objects, named myobject1 to myobject25. Now i'd like to update both of them. So i tried: for(ii in 1:25) { assign(paste("myobject.updated", ii, sep=""), update( myobject[ii] ,.~ + VAR2)) } Doesn't work I also tried to get all the names in a vector and update(names.myobject[ii],.~ + VAR2) Stiil doesn't work Any ideas

Hi, I'm using a loop to evaluate several models by taking adjacent variables from my dataframe. When i try to get predictions for new values, i get an error message about a missing variable in my new dataframe. Below is an example adapted from ?gam in mgcv package library(mgcv) set.seed(0) n<-400 sig<-2 x0 <- runif(n, 0, 1) x1 <- runif(n, 0, 1) x2 <- runif(n, 0, 1) x3 <-

I''m doing some simulations for which i need to use both S-plus and R. I generate in S+ some random normal distributions to define one dataset by iteration. I need to use the same dataset generated in S-plus in R. I was first thinking to generate in R the same dataset by using the same random number generator with a fixed seed. But It seems that S-plus and R don''t use the same

Dear List, I'm using GAMs in a multiple imputation project, and I want to be able to combine the parameter estimates and covariance matrices from each completed dataset's fitted model in the end. In order to do this, I need the knots to be uniform for each model with partially-imputed data. I want to specify these knots based on the quantiles of the unique values of the non-missing

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

Hi ! Is it plan to release a Browser History Management in Scriptaculous ? Or does anybody know a "prototype" library to manage browser history compatible with IE6, IE7 and FF ? We tried to use Really Simple History: http://code.google.com/p/reallysimplehistory/ But it seems to be incompatible with some Prototype event management. Regards, Jp -- Jean-Philippe Encausse - Veille /

Hello, I am trying to update the mgcv package on my Linux box and I keep getting an "Out of memory!" error. Does anyone know of a fix for this? Below is a snippet of the message that I keep getting: Thank you. Geoff ** R ** inst ** preparing package for lazy loading ** help *** installing help indices >>> Building/Updating help pages for package 'mgcv' Formats:

Bonjour, j'ai l'honneur de vous informer que je postule ? la liste de diffusion pour des fins utiles de cette formation "MOOC". Salutations ? toute l'?quipe d'encadrement. Cordialement. -- *DOUMBOUYA Moussa* *Epid?miologiste bio statisticien* ************************** Certifi? en FETP premi?re ligne en Guin?e *Chef de section pr?vention et lutte contre la maladie*

Hi, Is there migration for custom Effect between prototype/scriptaculous 1.4 and 1.5 ? Kevin Davis send me it''s cool Poof Effect but it seems to be tooooo fast in last version of scriptaculous. I think it comes from the code ? But I don''t know the best way to correct that ? Here is a copy/past of his mail for anybody interested in his Effect: > > Hey there, sorry for

Dear all, I have a serious problem to solve my model. I study over exploitation of fish in the bay of biscay (france). I know only the level of catch and the fishing effort (see data below) by year. My model is composed by the following equations: * the growth function Gt(St) = r*St*(1-St/sbar) with Gt the growth of each period t r intrinsec growth of the stock sbar carriyng capacity of the

Dear R list, I'm trying to calculate Mahalanobis distances for 'Species' of 'iris' data as obtained below: Squared Distance to Species From Species: Setosa Versicolor Virginica Setosa 0 89.86419 179.38471 Versicolor 89.86419 0 17.20107 Virginica 179.38471 17.20107 0 These distances were obtained with proc 'CANDISC'

Good day time! We've noticed the following issue with Samba 3.0.1 on Linux. When SMB client asks for the first time on FID for a locking or read request on an opened file, it causes an ACK to be sent from client side. Windows server also requests ACK BUT much more rarely. Sending ACK causes additional overhead when working in heavy locking and unlocking environment. Windows server

Good day time! We've noticed the following issue with Samba 3.0.1 on Linux. When SMB client asks for the first time on FID for a locking or read request on an opened file, it causes an ACK to be sent from client side. Windows server also requests ACK BUT much more rarely. Sending ACK causes additional overhead when working in heavy locking and unlocking environment. Windows server

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

Hi, I need to plot a spectogram for a speech signal. Is there a package that can do this? -- Rajesh.J [[alternative HTML version deleted]]

Hi, I need to plot spectrogram of a waveform.What package offers this? -- Rajesh.J [[alternative HTML version deleted]]

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