similar to: glm and gam confidence intervals

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 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

Hello Im a beginner on R and I would like how to develop a ROC statistic to evaluate a GAM model with a binomial distribution (Im using mgcv package) 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)

Hi everybody Im triying to develop a GAM using GraspeR, but i have found some problem in grasp.model I always obtain same error: RESPONSE NAME: Aves Error in "[.data.frame"(XXX, , c(1, sX)) : invalid subscript type The data are a vector for response and a vector for predictors. They must be data.frame? -- David Nogu?s Bravo Functional Ecology and Biodiversity

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

This is in regards to the suggested use of type="lpmatrix" in the documentation for mgcv::predict.gam. Could one not get the same result more simply by using type="terms" and interpolating each term directly? What is the advantage of the lpmatrix approach for prediction outside R? Thanks. -- View this message in context:

I have downloaded the source code for Linux. the ./configure procedure can find gcc but cannot find f77 or f2c I have looked for the above compiler with Test and installed all what pertains to Fortran. Still not even the man pages for f77 anf fc2 work >. Where is f77 and f2c with linux 10.3 ?? Thank you do much -- Maura E.M [[alternative HTML version deleted]]

Dear R-Users, Dr. Wood replied to a similar topic before where confidence intervals were for a ratio of two treatments ( https://stat.ethz.ch/pipermail/r-help/2011-June/282190.html). But my question is more complicated than that one. In my case, log(E(y)) = s(x) where y is a smooth function of x. What I want is the confidence interval of a ratio of log[(E(y2))/E(y1)] given two fixed x values of

Hello, Does anybody know how to make the following tutorial work on Typo? http://theory.isthereason.com/?p=499 Thanks! Jonas -- Posted via http://www.ruby-forum.com/.

Please use dput() to post your example matrix. Rambler1 wrote > > I have run into a problem in my code. What I want to accomplish is this: > I have a user input stock symbols into a list and from there I run the > quantmod package to get historical data. I compute the correlation matrix > and then turn that matrix into a simple matrix with 1's or 0's depending > on

Hi, I would be very grateful for advice on getting confidence intervals for the ordinary (non smoothed) parameter estimates from a gam. Motivation I am studying hospital outcomes in a large data set. The outcomes of interest to me are all binary variables. The one in the example here, Dead30d, is death within 30 days of admission. Sexf is gender (M or F), Age is age in years at the start

Here is my problem, On my network there is a subversion server to which everybody connects to and checks in/out code and documents. The workstations on the network are all based on Fedora core 3. And everybody's home directory is on a nfs share. This nfs share is mounted via the fstab. So no matter to which box a user logs in, his/her home dir is gonna be the same. Again, /home is a nfs

Hello, How do I plot a gam fit object on probability (Y axis) vs raw values (X axis) axis and include the confidence plot lines? Details... I'm using the gam function like this: l_yx[,2] = log(l_yx[,2] + .0004) fit <- gam(y~s(x),data=as.data.frame(l_yx),family=binomial) And I want to plot it so that probability is on the Y axis and values are on the X axis (i.e. I don't want log

Hi all, I am using GAM to model time trends in a logistic regression. Yet I would like to extract the the fitted spline from it to add it to another model, that cannot be fitted in GAM or GAMM. Thus I have 2 questions: 1) How can I fit a smoother over time so that I force one knot to be at a particular location while letting the model to find the other knots? 2) how can I extract the matrix

R-help, Sorry for posting again the same question (dated 26-03-2007) but all my mails have been sent to the recycle bin without possibility of recovering and thus I don't know if anyone has answer my query. Here is the original message: I'm applying a gam model (package mgcv) to predict relative abundances of a fish species. The covariates are year, month, vessel and statistical

Dear All Why NAT rules stays valid even if I flush nat anf table chains?? I have: iptables -P FORWARD DROP iptables -A FORWARD -m state --state ESTABLISHED,RELATED -j ACCEPT iptables -A FORWARD -s SOME_IP -d SOME_BCP_5_IP --dport 1234 -j ACCEPT iptables -i nat -A PREROUTING -s SOME_IP -d MY_INTERNET_IP \\ --dport 1234 -j DNAT --to-destination SOME_BCP_5_IP The conection is

Hai all, Iam a newbie to Asterisk. I want to configure my Asterisk thru Command Line Interface to connect two internal extensions and two external numbers and calls should occur between any of the two numbers. Can anybody kindly send me the configyration details for extensions.conf anf sip.conf file.. and if anything else needed to serve my purpose. Iam in a great need of this.... A reply word

Hi, I would like to have some advise about how to install Dopt on R for Windows. I can not manage to properly compile Dopt.f anf to get a suitable dll for R. Can someone explain me step by step what should I do ? Thanks in advance Isabelle Zabalza-Mezghani

Hello, I am a beginner in R and statistics, so my question may be trivial. Sorry in advance. I performed a Cox proportion hazard regression with 2 categorical variables with cph{design}. Then an anova on the results. the output is > anova(cph(surv(survival, censor) ~ plant + leaf.age + plant*leaf.age, > Mpnymph) Wald Statistics Response: Surv(survival, censored)

Dear R-helpers, I am trying to construct a confidence interval on a prediction of a gam fit. I have the Wood (2006) book, and section 5.2.7 seems relevant but I am not able to apply that to this, different, problem. Any help is appreciated! Basically I have a function Y = f(X) for two different treatments A and B. I am interested in the treatment ratios : Y(treatment = B) / Y(treatment = A) as