Displaying 20 results from an estimated 3000 matches similar to: "GAM-binomial logit link"
2008 Mar 31
1
unexpected GAM result - at least for me!
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
2007 Jun 22
1
two basic question regarding model selection in GAM
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
2005 Sep 26
4
p-level in packages mgcv and gam
Hi,
I am fairly new to GAM and started using package mgcv. I like the
fact that optimal smoothing is automatically used (i.e. df are not
determined a priori but calculated by the gam procedure).
But the mgcv manual warns that p-level for the smooth can be
underestimated when df are estimated by the model. Most of the time
my p-levels are so small that even doubling them would not result
2002 Nov 13
2
Comparing GAM objects using ANOVA
Hi,
Is it possible to compare two GAM objects created with the gam() function from the mgcv package. I use a slightly modified version of anova.glm() named anova.gam(), modified from John Fox (2002). It often gives me some aberant responses, especially with "F" test. I use a quasibinomial model and scale (dispersion) is calculated and used in the calculation of the F value. Does someone
2004 Mar 12
1
GCV UBRE score in GAM models
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)
2003 Jun 04
2
gam()
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
2004 Jan 19
2
Relative risk using GAM
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 ~
2003 Jul 14
1
gam and step
hello,
I am looking for a step() function for GAM's.
In the book Statistical Computing by Crawley and a removal of predictors has
been done "by hand"
model <- gam(y ~s(x1) +s(x2) + s(x3))
summary(model)
model2 <- gam(y ~s(x2) + s(x3)) # removal of the unsignificant variable
#then comparing these two models if an significant increase occurs.
anova(model, model2,
2009 Mar 31
1
CV and GCV for finding smoothness parameter
I received an assignment that I have to do in R, but I'm absolutely not very
good at it.
The task is the following:
http://www.nabble.com/file/p22804957/question8.jpg
To do this, we also get the following pieces of code (not in correct order):
http://www.nabble.com/file/p22804957/hints.jpg
I'm terrible at this and I'm completely stuck. The model I chose can be
found in here:
2006 Dec 04
1
GAM model selection and dropping terms based on GCV
Hello,
I have a question regarding model selection and dropping of terms for GAMs fitted with package mgcv. I am following the approach suggested in Wood (2001), Wood and Augustin (2002).
I fitted a saturated model, and I find from the plots that for two of the covariates,
1. The confidence interval includes 0 almost everywhere
2. The degrees of freedom are NOT close to 1
3. The partial
2004 Oct 26
3
GLM model vs. GAM model
I have a question about how to compare a GLM with a GAM model using anova
function.
A GLM is performed for example:
model1 <-glm(formula = exitus ~ age+gender+diabetes, family = "binomial",
na.action = na.exclude)
A second nested model could be:
model2 <-glm(formula = exitus ~ age+gender, family = "binomial", na.action =
na.exclude)
To compare these two GLM
2005 Oct 04
1
time out error and email with attachments
I'm facing some problem with emails and attachments.
I'm not sure it's dealing with MTU (it is set to the maximum number
before a fragmentation occurs).
I have no errors from pop3, but the following message is the best I can
show you just to understand my problem.
The log shows this:
"imap(roberto.cipollini at allmediasolution.com): Oct 04 17:35:03 Error:
Timeout while
2008 Jan 08
3
GAM, GLM, Logit, infinite or missing values in 'x'
Hi,
I'm running gam (mgcv version 1.3-29) and glm (logit) (stats R 2.61) on
the same models/data, and I got error messages for the gam() model and
warnings for the glm() model.
R-help suggested that the glm() warning messages are due to the model
perfectly predicting binary output. Perhaps the model overfits the data? I
inspected my data and it was not immediately obvious to me (though I
2002 Sep 10
2
Hat values for generalized additive models
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
2004 Jun 03
3
Problem with mgcv PACKAGES file format?
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
.......... .......... .......... .......... ..........
.......... .......... .......... .......... ..........
.......... .......... ..........
2010 Aug 05
1
plot points using vis.gam
Hello,
I'm trying to illustrate the relationships between various trait and
environment data gathered from a number of sites. I've created a GAM to do
this: gam1=gam(trait~s(env1)+s(env2)+te(env1,env2)) and I know how to create
a 3D plot using vis.gam. I want to be able to show points on the 3D plot
indicating the sites that the data came from. I can do this on a 2D plot
when there is one
2003 May 06
2
Sum by categorical variable
Dear R-list,
I have two variables (numerical and categorical) and would like to have the
sum (and maybe some other statistics) of the numerical variable by the
categorical one.
Can you help me,
Thank you,
Fabrizio
---------------------------------------------------------------
Fabrizio De Amicis
IT Department
Generali Information Technologies - (GIT)
Centro Galleria 2,
Via Cantonale
CH -
2007 Apr 09
1
How to solve differential and integral equation using R?
Hello,
I want to know if there are some functions or packages to solve differential
and integral equation using R.
Thanks.
Shao chunxuan.
[[alternative HTML version deleted]]
2006 Nov 28
4
GAMS and Knots
Hi
I was wondering if anyone knew how to work out the number of knots that
should be applied to each variable when using gams in the mgcv library?
Any help or references would be much appreciated.
Thanks
Kathryn Baldwin
2003 May 06
3
FW: Sum by categorical variable
Another easy(/stupid) question:
with the following command
j2<-xf1[1:10,"V4"]
I have
> j2
[1] CHROMOLI LINEAGERMAI RINALDI GIUNTIMA
AUTOSTELLA
[6] CAIZZONE CENTRO B PEL E C SNC CONSORZI MAN NORD
PDM
1304 Levels: MACHIAVELLI Snc MENARINI MANUFACTURERS LOGISTICS SERVICES
... ZUCCHERIFICIO