similar to: mgcv 1.0

mgcv version 1.2 is on CRAN now. mgcv provides generalized additive models and generalized additive mixed models with automatic estimation of the smoothness of model components. Changes in this version are: * A new gam fitting method is implemented for the generalized case. It provides more reliable convergence than the previous default, but can be a little slower. See ?gam.method,

mgcv version 1.2 is on CRAN now. mgcv provides generalized additive models and generalized additive mixed models with automatic estimation of the smoothness of model components. Changes in this version are: * A new gam fitting method is implemented for the generalized case. It provides more reliable convergence than the previous default, but can be a little slower. See ?gam.method,

Hi, I think there is a problem with the LAPACK version of qr() in version 1.7.0. (version below). 1. The documentation states that LAPACK=TRUE is the default, but the code has LAPACK=FALSE. 2. With LAPACK=TRUE qr() is never pivoting, even in cases where it very clearly should be. e.g. set.seed(0) X<-matrix(rnorm(40),10,4);X[,1]<-X[,2] qrx<-qr(X,LAPACK=TRUE) qrx$pivot # note, no

The rpanel package builds on the tcltk package to provide simple interactive controls for R functions, in particular to provide simple forms of dynamic graphics. The intention is to make this form of control particularly easy for R users to implement, with full documentation. The necessary tcltk variables are managed behind the scenes so that users need not be concerned with any

Hola! I have some problems with gam in mgcv. Firts a detail: it would be nice igf gam would accept an na.action argument, but that not the main point. I want to have a smooth term for time over a year, the same pattern repeating in succesive years. It would be natural then to impose the condition s(0)=s(12). Is this possible within mgcv? I tried to obtain this with trigonometric terms, aca:

Hi list, I'm working with gamm models of this sort, using Simon Wood's mgcv library: gm<- gamm(Z~te(x,y),data=DATA,random=list(Group=~1)) gm1<-gamm(Z~te(x,y,by=Factor)+Factor,data=DATA,random=list(Group=~1)) with a dataset of about 70000 rows and 110 levels for Group in order to test whether tensor product smooths vary across factor levels. I was wondering if comparing those two

Dear All, I am trying to produce log-log plots in R and I was wondering if any of you have a 'template' for generating these with 'nice' labels and log-log grids? I know I can set up axes individually and use the intervals I want, however, I will be producing a large number of these plots and would not like to do this manually for each of them + I am very new to R and at the

In R version 2.7.0 (2008-04-22) as.numeric(".") returns zero. > as.numeric(".") [1] 0 This must be a bug. Splus and previous versions of R (<= 2.6.0) return NA, as you might expect. I'm running R version 2.7.0 (2008-04-22) on Windows XP. Paul _____________________________________________________ Paul Johnson Robertson Centre for Biostatistics University of

I fitted a GAM model with Poisson distribution using the function gam() in the mgcv package. My model is of the form: mod<-gam(y~s(x0)+s(x1)+s(x2),family=poisson). To extract estimates at a specified set of covariate values I used the gam `predict' method. But I want to get estimate and standard error of the difference of two fitted values. Can someone explain what should I do? Thank

Hi all! I?ve got some problems with the function gam (library mgcv). For some models I get the error message : Error: no valid set of coefficients has been found:please supply starting values In addition: Warning message: NaNs produced in: log(x) This is a shortened code I used: gam(y ~ M1 + M3 + M4 + M5 + M6 + sex + M1*M3 + s(age), family=Gamma(link ="identity"), weights=days) If

hi, i'm working with mgcv packages and specially gam. My exemple is: >test<-gam(B~s(pred1)+s(pred2)) >plot(test,pages=1) when ploting test, you can view pred1 vs s(pred1, edf[1] ) & pred2 vs s(pred2, edf[2] ) I would like to know if there is a way to access to those terms (s(pred1) & s(pred2)). Does someone know how? the purpose is to access to equation of smooths terms

here is a distillation of a problem encountered in transfering some working code from R-0.63 to R-1.1.0 a1 <- 1:10 b1 <- tapply(a1,a1, sum) c1 <- ts(b1) c1 Error in if (NCOL(x) == 1) { : missing value where logical needed note that the error is returned as the value of calling ts() and is not automatically displayed problem seems to be that is.array(b1) returns TRUE yet dim(b1)

Hi! I'm learning mgcv, and reading Simon Wood's book on GAMs, as recommended to me earlier by some folks on this list. I've run into a question to which I can't find the answer in his book, so I'm hoping somebody here knows. My outcome variable is binary, so I'm doing a binomial fit with gam(). I have five independent variables, all continuous, all uniformly

Hi I wonder if anyone can throw some light on the * console message. This only occurs when I register a phone on the end of a BT ADSL line, with a Draytec router/modem. The phone registers okay but cannot dial out. Console message: Notice[1125329600]: File chan_sip.c, Line 1759 (copy_via_headers): No field 'Via' present to copy Thanks Steven *****************************************

the following occurred inadvertently and brought R-1.1.0 down rep(1:3, c(4,2,-6)) Segmentation Fault (core dumped) platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 1 minor 1.0 year 2000 month June day 15 language R

I have come across three problems in the past few days, in spell of heavy R ( version R0.50-a4/Sun Solaris2.5.1 ) use. 1. I was using lwd=2 to get thicker lines on plots for printing, but although the 'lwd' parameter works fine with x11(), the thickened lines do not print with print.plot, or by using postscript() directly. 2. Try the following, plot(1:10, -(1:10))

Setup: Sparcstation 4 with Solaris 2.5.1 and Sun compilers After downloading R-0.16.1 and runnning configure, make fails as follows cd main; make cc -o R.binary arithmetic.o array.o attrib.o bind.o builtin.o character.o coerce.o colors.o context.o cov.o cum.o dataframe.o debug.o deparse.o dotcode.o dstruct.o duplicate.o edit.o envir.o errors.o eval.o format.o gram.o

The following works as expected plot(1:5, panel.first=grid(2,2)) and if my.panel<-function() grid(2,2) then plot(1:5, panel.first=my.panel() ) is also OK but, do.call("plot", list(x=1:5, panel.first=grid(2,2)) do.call("plot", list(x=1:5, panel.first=my.panel)) do.plot("plot", list(x=1:5, panel.first=my.panel() ))

Hi, I've got the R library "mgcv" for GAM written by Simon Wood which works well in many instances. However, over the years I got attached to the S+ implementation of GAM which allows loess smoothing in more than 1 dimension as well as spline smoothing. Has anyone ported the S+ GAM library to R? Regards, Doug Beare. Fisheries Research Services, Marine Laboratory, Victoria Road,

Hi I have installed Festival on the same box as asterisk and followed the instructions to integrate it with asterisk. Festival seems to work fine on its own performing text to speech from the command line or via a file. Asterisk answers the call but there is no speech. I can see no errors in the Festival log file The asterisk console shows --Executing Answer("SIP/81801-c091",