similar to: residuals in plot.gam (mgcv)

Hi Allstaters, Does anybody know if the R package can fit Generalised Additive Models?. Thanks indeed, Aurelio. atobias at ole.com ___________________________________________________________________ Consigue tu e-mail gratuito TERRA.ES Haz click en http://www.terra.es/correo/ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hi helpteam, may you send me the exact syntax for getting the estimated degrees of freedom for a gam? I tried to do it by name <- edf(gam_name) but this didn't work. Thanks a lot, greetings from Munich, Thomas Nittner -- Thomas Nittner Tel: ++49 89 2180 3197 Akademiestra?e 1 Fax: ++49 89 2180 5042 D - 80799 M?nchen http://www.stat.uni-muenchen.de

Hi, Is there any difference in the gam() function in R and Splus? I just tried to do an Ace Plot (plotting the result of the gam() ) on a dataset. R gave me four (the data set has one response and four explanatory variables) ace plot with perfectly straight lines. I kind of felt that it is highly unlikely, so I tried it in Splus. Using the exact commands Splus gave me four completely

Would anyone be able to provide insight for the following question, please? Setting: estimation of prediction intervals for age-period-cohort models using GAMs (rate ~ s(age,period)) Method: bootstrap (Davison and Hinkley, 1997) Issue: standardisation of the residuals for resampling requires an adjustment using the diagonals of the hat matrix. Is there a simple way to get the hat values out of a

I hope this is the right place for my question, my apologise if it is not. I am a long time Mac/UNIX user and my local machine is a Mac. At present I run an orphan version of R which a worried out of (*****). On reflection I guess I'd better not say, but it changed my life. It is howerver a little out of date ( say 0.6ish) So the obvious question, do I move to MacLinux or is there likely

Hello Is there a simple way of fitting a 'smooth curve' to a time series with the additional constraint that the first derivative be not negative? I know that if I choose a "large enough" window for a moving average sort of filter on my data I end up with a non-decreasing curve anyway, but I'm wondering if there is a pre-built method that can include the constraint from

To R users: I am running an R program with .Fortran() calls, and discovered that the memory, as reported by gc(), continues to grow with new loop iterations. I hadn't been able to find information on how to 'destroy' objects in earlier loops to keep the program size from continuous growth; calling rm() to remove some what I thought to be large arrays didn't seem to help. Any

Hello: I sent this question the other day with the wrong subject heading and couple typos, with no response. So, here I go again, having made those corrections. I would like to estimate, for lack of a better description, a partially additive non-parametric model with the following structure: z~ f(x,y):w1 + g(x,y):w2 + e In other words, I'd like to estimate the marginals with respect to

Dear R users, Using the library(mgcv) and running R under MacOSX, I have fitted a generalised additive model with binomial errors in order to check the linearity of two continuous variables ap2mm and diffdaysm in a glm: > mymodel.gam <- gam(diedhos~ s(ap2mm) + Dweekm + s(diffdaysm) + Dweekm:diffdaysm + ap2mm:Dweekm, binomial) I would like postscript gam plots for the two smoothed

Dear R-users, I've a question regardiing multiple linear regression. Is it possible to regress a function of the type y=b0+b1*x1+b2*x2 constraining the partial derivatives dy/dx1 and dy/dx2 to be greater than 0 ??? Thank you very much for any possible suggestion. Best regards Marco -- Marco Pagani - PhD Student Dipartimento di Scienze della Terra - Sezione Geofisica Universita' degli

Hi, I need help regarding estimating a linear model where restrictions are imposed on the coefficients. An example is as follows: Y_{t+2}=a1Y_{t+1} + a2 Y_t + b x_t + e_t restriction a1+ a2 =1 Is there a function or a package that can estimate the coefficient of a model like this? I want to estimate the coefficients rather than test them. Thank you for your help Ahmad Abu Hammour --------------

Hello everybody, I think my question has been asked before but I am posing it once again since I need it. Is there any way to find the rank of a matrix in R or Splus? ______________________ Indrajit SenGupta Department Of Statistics St. Xavier's College Calcutta University mailto:indra_calisto at yahoo.com mailto:indrajitsg at vsnl.net

I'm trying to figure out a good way to smooth binary data. The ideal approach appears to be the "sm.logit" function in library "sm", but I haven't had success with it. Below, is some code illustrating what I've come up with so far, but I'm hoping there is a better approach. I'm using R 1.2 (development) under Windows 98. -Bill library(MASS) data(birthwt)

Hi All, Thought I would let those that are interested know that I had success in running 24G on an Asus P6T with 24G kit of Kingston DDR3. While I was putting this together I saw lots of forum posts asking if anyone had tried it. Well we did here at our work and all looks great including running "memtest86" overnight. I have a fluid dynamics simulation running on it with 90% memory

Has anyone got a function for smooth monotonic interpolation of a univariate function? I'm after something like the NAG function PCHIM which does monotonic Hermite interpolation. Alternatively, montononic cubic spline interpolation. Please reply directly. Rob Hyndman ___________________________________________________ Rob J Hyndman Associate Professor & Director of Consulting

I noticed that there is no generalised additive model functions in R (1.1.1) ... is there a package that implements them? Thanks Prasad ***************************************************************** Mr. Anantha Prasad, Ecologist/GIS Specialist USDA Forest Service, 359 Main Rd. Delaware OHIO 43015 USA Ph: 740-368-0103 Email: aprasad at fs.fed.us Web:

Dear List I am using R with mcgv package to model spatial variation in density estimates of dorcas gazelle in Sinai. I have 59 points of data and 4 explanatory variables(distance from mountain edge, camel presence, Latitude & Longitude). I want to test the model fir via bootstraping. I have used the jacknife bootstraping but it have the limitation of allowing only 58 trials. I tried to use the

I was having some problems persuading optim() to give me the answers I wanted, & simplified down to: sqr<-function(x){(x+1)^2} optim(1,sqr) I accept this is a hammer to crack a nut, but was still expecting the answer -1. I got: $par [1] -0.8 $value [1] 0.04 $counts function gradient 12 NA $convergence [1] 0 $message NULL so I've

Hi R Users, I am very new to R and would like to do something quick if possible, please help! Suppose I have a data set of y versus x, how can I generate a smoothing line of y versus x (for example, using loess) and at the same time, generate a pair of confidence intervals for the smoothing or mean plus/minus standard deviation? Yi Zhu Golder Associates Inc. USA

Hi there, I'm using R1.0.1 Windows 98. This file contains some inputs and an aov function code. Can someone check it for me? Somehow I got completely different answer when typing them in R and in Splus. Splus gives me this: > summary( Turnip.aov ) Error: Blocks Df Sum of Sq Mean Sq F Value Pr(F) Residuals 3 163.7367 54.57891 Error: Plots %in% Blocks