Displaying 20 results from an estimated 10000 matches similar to: "Newbie question--locally weighted regression"
2008 Feb 16
2
Possible overfitting of a GAM
The subject is a Generalized Additive Model. Experts caution us against
overfitting the data, which can cause inaccurate results. I am not a
statistician (my background is in Computer Science). Perhaps some kind soul
would take a look and vet the model for overfitting the data.
The study estimated the ebb and flow of traffic through a voting place. Just
one voting place was studied; the
2004 Sep 24
3
geographically weighted glm
Hi all,
I am interested in obtaining R code related to geographically weighted
regression.
In particular, I am interested in building geographically weighted
Poisson GLMs. The model will contain categorical and continuous x
independent variables, with interaction effects between categorical and
continuous variables.
Anybody have anything I can look at?
thanks,
Mark.
--
2010 Sep 03
2
density() with confidence intervals
Hello R users & R friends,
I just want to ask you if density() can produce a confidence interval, indicating how "certain" the density() line follows the true frequency distribution based on the sample you feed into density().
I've heard of loess.predict(loess(y ~ x), se=TRUE) which gives you a SE estimate of the smoothed scatterplot - but density() kernel smoothing is not the
2000 May 01
1
GAMs under R?
At 06:09 AM 5/1/00 +0100, Prof Brian D Ripley wrote:
>On Sun, 30 Apr 2000, Stephen R. Laniel wrote:
>
>> I was just now surprised to note that functions to go generalized additive
>> models don't appear to exist under R 1.000. In particular, the gam() and
>> loess() functions aren't there. Are they hidden somewhere and I just
>> haven't noticed?
>
2010 May 31
3
What does LOESS stand for?
Dear R-community,
maybe someone can help me with this:
I've been using the loess() smoother for quite a while now, and for
the matter of documentation I'd like to resolve the acronym LOESS.
Unfortunately there's no explanation in the help file, and I didn't
get anything convincing from google either.
I know that the predecessor LOWESS stands for "Locally Weighted
2003 Jun 03
1
S+ style implementation of GAM for R?
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,
2001 Dec 17
3
smoothing line and a pair of confidence intervals
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
2009 Jul 28
2
A hiccup when using anova on gam() fits.
I stumbled across a mild glitch when trying to compare the
result of gam() fitting with the result of lm() fitting.
The following code demonstrates the problem:
library(gam)
x <- rep(1:10,10)
set.seed(42)
y <- rnorm(100)
fit1 <- lm(y~x)
fit2 <- gam(y~lo(x))
fit3 <- lm(y~factor(x))
print(anova(fit1,fit2)) # No worries.
print(anova(fit1,fit3)) # Likewise.
print(anova(fit2,fit3)) #
2005 Oct 12
1
step.gam and number of tested smooth functions
Hi,
I'm working with step.gam in gam package. I'm interested both in spline and
lowess functions and when I define all the models that I'm interested in I get
something like that:
> gam.object.ALC<-gam(X143S~ALC,data=dane,family=binomial)
>
2004 Aug 06
2
gam --- a new contributed package
I have contributed a "gam" library to CRAN,
which implements "Generalized Additive Models".
This implementation follows closely the description in
the GAM chapter 7 of the "white" book "Statistical Models in S"
(Chambers & Hastie (eds), 1992, Wadsworth), as well as the philosophy
in "Generalized Additive Models" (Hastie & Tibshirani 1990,
2004 Aug 06
2
gam --- a new contributed package
I have contributed a "gam" library to CRAN,
which implements "Generalized Additive Models".
This implementation follows closely the description in
the GAM chapter 7 of the "white" book "Statistical Models in S"
(Chambers & Hastie (eds), 1992, Wadsworth), as well as the philosophy
in "Generalized Additive Models" (Hastie & Tibshirani 1990,
2007 Apr 03
3
Testing additive nonparametric model
I have estimated a multiple nonparametric regression using the loess
command in R. I have also estimated an additive version of the model using
the gam function. Is there a way of using the output of these two models to
test the restrictions imposed by the additive model?
2011 Jun 04
1
nonparametric logistic regression based on locally weighted scatterplot smoothing (lowess)
Dear UseRs:
Recently, I have read an article regarding the association between age and
lymph node metastases.
http://jco.ascopubs.org/content/27/18/2931.long
In statistical analysis, the authors stated "Because a nonlinear
relationship between age and lymph node involvement was expected based on
existing literature, lymph node involvement was also regressed on age using
nonparametric
2007 Apr 08
1
Relative GCV - poisson and negbin GAMs (mgcv)
I am using gam in mgcv (1.3-22) and trying to use gcv to help with model selection. However, I'm a little confused by the process of assessing GCV scores based on their magnitude (or on relative changes in magnitude).
Differences in GCV scores often seem "obvious" with my poisson gams but with negative binomial, the decision seems less clear.
My data represent a similar pattern as
2010 Jun 18
2
varIdent error using gam function in mgcv
Hello,
As I am relatively new to the R environment this question may be either
a) Really simple to answer
b) Or I am overlooking something relatively simple.
I am trying to add a VarIdent structure to my gam model which is fitting
smoothing functions to the time variables year and month for a particular
species. When I try to add the varIdent weights to variable Month I get this
error returned.
2002 Oct 31
3
Loess with glm ?
Hello,
I am wondering if there is an easy way to combine loess() with glm()
to produce a locally fitted generalised regression.
I have a data set of about 5,000 observations and 5 explanatory variables,
with a binary outcome. One of the explanatory variables (lets call it X)
is much more predictive than the others. A single glm() regression over
the entire data set produces rather poor results,
2017 Jun 12
2
plotting gamm results in lattice
Dear all,?
I hope that you can help me on this. I have been struggling to figure this out but I haven't found any solution.
I am running a generalised mixed effect model, gamm4, for an ecology project. Below is the code for the model:
model<-gamm4(LIFE.OE_spring~s(Q95, by=super.end.group)+Year+Hms_Rsctned+Hms_Poaching+X.broadleaved_woodland? ? ? ? ? ? ?+X.urban.suburban+X.CapWks,
2017 Jun 12
0
plotting gamm results in lattice
Hi Maria
If you have problems just start with a small model with predictions and then plot with xyplot
the same applies to xyplot
Try
library(gamm4)
spring <- dget(file = "G:/1/example.txt")
str(spring)
'data.frame': 11744 obs. of 11 variables:
$ WATERBODY_ID : Factor w/ 1994 levels "GB102021072830",..: 1 1 2 2 2 3 3 3 4 4 ...
$ SITE_ID
2009 Jun 08
3
Plotting two regression lines on one graph
Hi!
I have fitted two glms assuming a poisson distribution which are:
fit1 <- glm(Aids ~ Year, data=aids, family=poisson())
fit2 <- glm(Aids ~ Year+I(Year^2), data=aids, family=poisson())
I am trying to work out how to represent the fitted regression curves of fit1 and fit2 on the one graph. I have tried:
graphics.off()
plot(Aids ~ Year, data = aids)
line(glm(Aids ~ Year,
2006 Jan 02
1
An embarrassment of riches
I have a dataset which I am trying to smooth, using locally weighted
regression. The y values are count data, integers with Poisson
distribution, and it is important for the regression function to know
this, since assuming a Gaussian distribution will lead to substantial
errors. It is a time series; the x values have equal five minute
intervals.
Here is the problem: I have an embarrassment