search for: smoothers

Hi, I am using ggplot2 to with the following code: gmathk2 <- qplot(time,math,colour=Kids,data=kids.ach.lm.k5,geom="smooth",method="lm",formula=y~ns(x,1)) + opts(title="Smoother Plot: Math K-5") + xlab("Time") + ylab("Math") + scale_colour_brewer(pal="Set1"); gmathk2 This plots all the smoother for all the x values. What I'd like

Dear R help list, I'm fitting a 'variable coefficient model' in the MGCV package and I want to plot the different smoothers I get for each factor level in one graph. So, I do something like this to fit the gam: Mtest <- gam(outcome ~ s(age, by=as.numeric(gender==0)) + s(age,by=as.numeric(gender==1))+factor(Gender)) Then I can plot the smoother for gender=0: plot(Mtest,select=1) And for gender=1: plot(Mtest,select=...

Can anyone tell me the trick for obtaining the smoother matrix from smooth.spline when there are non-unique values for x. I have the following code but, of course, it only works when all values of x are unique. ## get the smoother matrix (x having unique values smooth.matrix = function(x, df){ n = length(x); A = matrix(0, n, n); for(i in 1:n){ y = rep(0, n); y[i]=1; yi =

Greetings, I would like to use a cubic spline or smoother to model the fixed effects within nlme. So far the only smoother I have been able to get to run successfully in nlme is smooth(). I tried smooth.spline: fixed=list(lKa~1,lCL~smooth.spline(BSA, df=3)) the error I got was the following. Error in model.frame(formula, rownames, variables, varnames, extras, extranames, : invalid

Hi there, I'm plotting some glass RI values just by plotting plot(x) then I put on my lowess smoother lines(lowess(x)) now I want to put on some 95% Confidence Interval bands of the lowess smoother, but don't know how?? Thanks -- Gareth Campbell PhD Candidate The University of Auckland P +649 815 3670 M +6421 256 3511 E gareth.campbell@esr.cri.nz gcam032@gmail.com [[alternative

...e |ID ), weights=varIdent(form=~1|ID)) However, in the current version of R, this does not work anymore, and should be changed into BM2<-gamm(Birds~Rain+ID+ s(Time,by=ID), correlation=corAR1(form=~Time |ID ), weights=varIdent(form=~1|ID)) It turns out that 2 of the 4 smoothers have estimated degrees of freedom of 1, so a linear effect would be sufficient. Now my question is how I need to change the code in order to have a time smoother for ID=Coot.Oahu and ID=Coot.Maui, and a linear time effect for ID=Stilt.Oahu and ID=Stilt.Maui. With the "old" R-code, this s...

Hello. I have data (plant weights over time) that are non-linear and in which the variance increases over time. I have to estimate the first derivatives of plant weight given time (i.e. growth rate) and their se, using a regression smoother, and I have been considering cubic spline smoothers. However, I do not know if this can be done given that the error variance would increase over time. Does anyone know what the effect of a non-constant error variance has on the estimates of the 1st derivative and its se? Bill Shipley Departement de biologie Universite de Sherbrooke Sherbrooke...

Dear all, Anybody knows about the difference amongst various spline smoothers, specifically in R, 'bs' (by default a cubic spline), 'ns', smooth.spline with roughness penalty along with many other smoothers? I've consulted serveral books like 'S-plus Guide to Statistics' by Mathsoft, Ripley and Venables' Splus book and Eubank's Spline Smoo...

Hello. I want to fit a smoother spline (or an equivalent local regression method) to a series of data in which the initial value of the 1st derivative (slope) is constrained to a specific value. Is it possible to do this? If so, how? Bill Shipley [[alternative HTML version deleted]]

Dear r-help, I would like to use the 'ppm' function of the 'spatstat' package to fit a Strauss inhibition model. I understand that I can specify a parametric model for the "background" trend, but how would I specify a trend which is estimated using a Kernel density smoother? In particular, I would like to use the 'kde' function of the 'ks' package to

Hi All I am having this problem when setting up a H323 call. Can anybody tell me what is going on? Thanks Serge ------------------ NOTICE[245776]: File chan_oh323.c, Line 1293 (oh323_write): H323:1637: Format changed from 4 to 8. WARNING[245776]: File frame.c, Line 76 (ast_smoother_feed): Smoother was working on 4 format frames, now trying to feed 8? ERROR[245776]: File chan_oh323.c, Line 1380

Hello, I tried to get the oh323 drivers running. The driver loads, but as soon as a H323 voice communication should be started, following error occurs: -- Executing Playback("OH323/R1", "invalid") in new stack Jun 3 01:26:20 ERROR[294931]: chan_oh323.c:1933 oh323_write: OH323/R1: Failed to create smoother. Jun 3 01:26:20 WARNING[294931]: file.c:539

hello, Im tryin to make Calls from MS Netmeeting(h323) to Xlite(SIP) it rings, but as soon as i answered it dissconnects!!!! This is what i get from the Asterisk console: -- Executing Dial("OH323/R27469", "SIP/xlite1|10") in new stack Aug 13 10:19:03 DEBUG[524304]: chan_sip.c:1265 create_addr: Setting NAT on RTP to 0 Aug 13 10:19:03 DEBUG[524304]: chan_sip.c:1500 sip_call:

I get the following error when I run qplot() qplot(grade, read,data = hhm.long.m, geom = c("point", "smooth")) Error in smooth.construct.cr.smooth.spec(object, data, knots) : x has insufficient unique values to support 10 knots: reduce k. I am not sure how to tackle this problem. When I take a subsample (< 1000) than I am able to run that function but with my sample

Hello. I want to estimate the predicted values and standard errors of Y=f(t) and its first derivative at each unique value of t using the smooth.spline function. However, the data (plant growth as a function of time) show substantial heterogeneity of variance since the variance of plant mass increases over time. What is the consequence of such heterogeneity of variance in terms of bias in the

Good morning all! I am attempting to superimpose a moving-average smoother onto a graph of daily plots. These plots (in table[,2] below) span about 350 days and looks very noisy. I'd like for this smoother to plot the average of each group of 7 consecutive days (weekly) and show a line which joins these series of averages. Given the definition of MA, the first and last points will

Dear R-list members, Once a gam (package gam) model has been fitted with family=poisson, is there some R function that could extract the diagonal elements of the smoother matrix S, to be used in a cross-validation for the selection of the best smoothing parameter, following equation 3.19 on page 48 of the book T.J. Hastie and R.J. Tibshirani, Generalized Additive Models, Chapman and Hall/CRC,

Hallo, I want to use gam from the mgcv package with a mrf smoother. This is my data set (`x`) ?????? y????????? id ??? 1? 0.6684496? 1 ??? 2? 0.6684496? 2 ??? 3? 0.6684496? 3 ??? 4? 0.6684496? 4 ??? 5? 0.6684496? 5 ??? 6? 0.6684496? 6 ??? 7? 0.6684496? 7 ??? 8? 0.5879492? 8 ??? 9? 0.5879492? 9 ??? 10 0.5879492 10 ??? 11 0.5879492 11 ??? 12 0.5879492 12 ??? 13 0.5879492 13 ??? 14

Good evening,I am Marta Colombo, student of "Politecnico di Milano". I'm studying Local Regression Techniques such as loess, smoothing splines and kernel smoothers. Choosing "symmetric" for the argument "family" in loess function it is possible to produce a robust estimate , in function smooth.spline and ksmooth I didn't find this possibility. Well, is there a way to produce a robust estimate using smoothing splines or kernel smoothers...

Hi, I have used the mgcv library to generate a simple additive model. I want to know how to plot the function on the raw data with confidence intervals whan I have TWO variables in the model. I get it to work with one variable but not with two. I am on the limit for what I understand in R, so be gentle. I have read the help file on predict.gam, but did not get any help out of it. #My model: