similar to: Hat values for generalized additive models

Updating my 1.8.0 R installation (>update.packages() ) I obtain the following (SORRY FOR THE LENGTH OF THE LOG BUT IT HELPS!!!): ................ downloaded 135Kb KernSmooth : Version 2.22-11 in /usr/lib/R/library Version 2.22-12 on CRAN Update (y/N)? y mgcv : Version 0.9-3.1 in /usr/lib/R/library Version 0.9-6 on CRAN Update (y/N)? y trying URL

Trying following example from Pinheiro and Bates in order to fit an ARMA(1,1) model: library(nlme) fm1Ovary.lme<-lme(follicles~sin(2*pi*Time)+cos(*pi*Time),data=Ovary,random=p dDiag(~sin(2*pi*Time))) fm5Ovary.lme<-update(fm1Ovary.lme,corr=corARMA(p=1,q=1)) I get follwing error message: Error in "coef<-.corARMA"(`*tmp*`, value = c(62.3428455941166, 62.3428517930051 :

Hello, I'm using R 1.40 on MacOS X X.1.2 (installed via the fink package manager). To upgrade my installed packages, I tried to use update.packages() today. All went well for most packages, with the exception of KernSmooth and cluster. In both cases, libraries were not found although I think they are present. Here's what happened: ---------------------------------- >

Qusetion #1 ********* Model selection in GAM can be done by using: 1. step.gam {gam} : A directional stepwise search 2. gam {mgcv} : Smoothness estimation using GCV or UBRE/AIC criterion Suppose my model starts with a additive model (linear part + spline part). Using gam() {mgcv} i got estimated degrees of freedom(edf) for the smoothing splines. Now I want to use the functional form of my model

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)

Dear Simon, Thank you for your quick reply! I used to perform the GAMsetup in the following manner: GAMsetup sintax: x.summer: vector used for construct the spline knots<-14 N<-length(x.summer) x<-array(x.summer,dim=c(1,N)) G<-list(m=1,n=N,nsdf=0,df=knots+1,dim=1,s.type=0,by=0,by.exists=FALSE,p.order=0,x=x,n.knots=knots,fit.method="mgcv") H<-GAMsetup(G) with the

I'm looking for the best way to do the following: run a set of GAM models, and then make predictions with new data. My problem is the size of the gam model object, I would like to strip it down to the bare minimum of information needed to apply the model to new data. For example, if this were a linear model, I would just keep the betas. If this were an ordinary spline fit, I think I

Dear all yesterday I got around to install R-1.5.1 on my machine at home which, admittedly, is running a very old version of RedHat Linux. I noticed two things: 1) My pdflatex version is too old to produce the pdf version of the manuals. `./configure' checks whether a recent version of makeinfo is installed and tells you whether you can produce the info version or not. (I was told I

I am a new user of R. I am trying to fit gam model with our air pollution data. I used Foreign package to call data from SPSS and used MGCV package to fit gam. The following are the steps I used: > dust<- read.spss("a:dust9600jan.sav") > c<-gam(MRESPALL~s(DUSTM)+s(TEMP)+s(RH),family=poisson,data=dust) > summary(c) Family: poisson Link function: log Formula: MRESPALL ~

Hi I am afraid i am not understanding something very fundamental.... and does not matter how much i am looking into the book "Generalized Additive Models" of S. Wood i still don't understand my result. I am trying to model presence / absence (presence = 1, absence = 0) of a species using some lidar metrics (i have 4 of these). I am using different models and such .... and when i

Hi I am trying to plot the following gam with termplot but keep getting the error message: Error in order(xx) : unimplemented type 'list' in 'orderVector1' Is there anyway I can rectify this to get my parametric coefficients plotted? Thanks Will Family: binomial Link function: logit Formula: fgha$pa ~ s(fgha$wspd) + fgha$depth + fgha$slha +

I have asked this question on Stackoverflow and was told it does not relate to the sites' mission as it is statistical question, thus I brought it here. I am fitting a gam mode in the mgcv package to study associations of environmental pollutants and mortality. The aim is to choose a model with lowest mgcv and also to reduce the PACF to less than < |0.1|. library(gamair) library(mgcv)

Hi, I'm adjusting a binomial gam and I use the inverse variance of each proportion as weights in the gam function (it's a model for the sex ratio of hake). My doubt is what are these weights doing ? I get a negative UBRE and if I use GCV the scale parameter is close to 0 ... Regards EJ

Dear all, I'm using a binomial distribution with a logit link function to fit a GAM model. I have 2 questions about it. First i am not sure if i've chosen the most adequate distribution. I don't have presence/absence data (0/1) but I do have a rate which values vary between 0 and 1. This means the response variable is continuous even if within a limited interval. Should i use

Dear all, I've now spent a couple of days trying to learn R and, in particular, the gam() function, and I now have a few questions and reflections regarding the latter. Maybe these things are implemented in some way that I'm not yet aware of or have perhaps been decided by the R community to not be what's wanted. Of course, my lack of complete theoretical understanding of what

Hello All, I'm getting this error (Version: 1.9.0-1 on a debian system) > update.packages("mgcv") trying URL `ftp://mirror.aarnet.edu.au/pub/cran/src/contrib/PACKAGES' ftp data connection made, file length 169516 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ..........

I always followed http://cran.r-project.org/bin/linux/ubuntu/ to install R on Ubuntu 8.1. I had no errors before! > install.packages("XML") Warning in install.packages("XML") : argument 'lib' is missing: using '/usr/local/lib/R/site-library' --- Please select a CRAN mirror for use in this session --- Loading Tcl/Tk interface ... done trying URL

From: Tom Rix <trix at redhat.com> This is a upcoming change to clean up a new warning treewide. I am wondering if the change could be one mega patch (see below) or normal patch per file about 100 patches or somewhere half way by collecting early acks. clang has a number of useful, new warnings see

From: Tom Rix <trix at redhat.com> This is a upcoming change to clean up a new warning treewide. I am wondering if the change could be one mega patch (see below) or normal patch per file about 100 patches or somewhere half way by collecting early acks. clang has a number of useful, new warnings see

Hello. I found that anonymous structures does not work on gcc-2.95. If you guys want to support a bit older platforms I suggest fixing it. You can check out patch I created to fix this issue. I just added 2 extra structures to remove anonymous structs inside connection_status_t and node_status_t. Patch is here: ftp://borg.uu3.net/pub/unix/tinc/tinc.patch Attaching it as well. Regards, Borg