search for: smoother

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 to do is plot the smoother for the x values that are only greater than or equal to 0. I don't want th...

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){...

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, va...

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 HTML ve...

Hi, I was reading the book on "Mixed Effects Models and Extensions in Ecology with R" by Zuur et al. In Section 6.2, an example is discussed where a gamm-model is fitted, with a smoother for time, which differs for each value of ID (4 different bird species). In earlier versions of R, the following code was used BM2<-gamm(Birds~Rain+ID+ s(Time,by=as.numeric(ID=="Stilt.Oahu"))+ s(Time,by=as.numeric(ID=="Stilt.Maui"))+ s(Time,by=as.numeri...

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

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

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 estimate the "background" intensity and then use this as the trend for a Strauss inhibition process. Thanks already in advance, Alejandro Veen

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 (oh323_write): H323:1637: Failed to fill smoother. WARNING[245776]: File file.c, Line 512 (ast_readaudio_callback): Failed to write frame -- Playing 'vm/109/unavail' ==...

...e 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 ast_readaudio_callback: Failed to write frame Asterisk is CVS Head, oh322+libs: [chan_oh323.so] => (OpenH323 Channel Driver) == Parsing '/etc/asterisk/rtp.conf': Found == Parsing '/etc/asterisk/oh323.conf': Found [1]WrapH323EndP...

...:19:05 DEBUG[524304]: res_features.c:422 ast_bridge_call: Read from SIP/xlite1-89a7 (4,4) Aug 13 10:19:05 WARNING[524304]: chan_oh323.c:1020 oh323_exception: OH323/R27469: Invalid format of RTP addresses. Aug 13 10:19:05 ERROR[524304]: chan_oh323.c:1933 oh323_write: OH323/R27469: Failed to create smoother. Aug 13 10:19:05 DEBUG[524304]: channel.c:2555 ast_channel_bridge: Bridge stops because we're zombie or need a soft hangup: c0=OH323/R27469, c1=SIP/xlite1-89a7, flags: No,Yes,No,No Aug 13 10:19:05 DEBUG[524304]: channel.c:2675 ast_channel_bridge: Bridge stops bridging channels OH323/R27469 an...

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

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, 1990 ? A sample code for fitting the model is: n <- 100 x <- runif...

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

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

...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: Model <- gam(var ~ s(var1, k=4) + s(var2, k=4), data = mydata) # Plotting the model: par(mfrow=c(1,1)) plot(Model ) #get a plot of the smoother plot(var ~ var1, data = mydata) #plotting the raw data ##predictions from M1 pred.gam <- predict.gam(M1, se.fit=TRUE, type="response") #Plotting the function and confidence intervals. I <- order(mydata$Ncell) lines(mydata$Ncell[I], pred.gam$fit[I], lwd = 2) lines(mydata$Ncell[I],...