similar to: Regarding Savitzky-Golay Smoothing Filter

I got a question from a fellow PhD student that work with spectrum analysis in Excel and now he has lots of spectrums that needs to be smoothed, which would be nice to be able to do in batch. Is there an R package that can do: Savitzky-Golay smoothing for reflectance spectral data or a function that does something similar. _______________________________________ Henrik Andersson

As the request for the Savitzky-Golay Algorithm in R has come up several times, I here include my implementation based on code written for Matlab. Savitzky-Golay uses the pseudo-inverse pinv() of a matrix. There is an 'generalized inverse' ginv() in the MASS package, but I use a simpler form because I didn't want to 'require' MASS any time I apply Savitzky-Golay.

If I'm not mistaken, that's sort of local polynomial with even degree and fixed bandwidth (based on my own interpretation of description in Numerical Recipes). You can do that with functions in the KernSmooth package. HTH, Andy > -----Original Message----- > From: wolf at micro-biolytics.com [mailto:wolf at micro-biolytics.com] > Sent: Friday, March 21, 2003 1:43 PM > To:

Hi, Does anyone know how to use weights and generate error bounds for Savitzky-Golay? I have a (smallish) set of points y equally spaced each with a known error and would like to smooth them using S-G but so as to take into account the error already have and construct new error bounds around them that take into account the errors they had at the beginning and the erros they get as a result

Dear all, Sorry if this is slightly off the track as far as R is concerned, but I have been using the Savitzky-Golay filter to estimate some derivatives of interest. I am wondering however, if anyone has seen anything in the literature (or has any ideas) of how these estimates perform asymptotically. Does anyone know what the rate of convergence is for these? Thanks, matt.

Hi I am running a Savitzky-Golay smoothing filter (http://tolstoy.newcastle.edu.au/R/help/04/02/0385.html) for variables in my dataset, dim (272:90). I managed to run the code for individual variables in the dataset and then combine the results into a single dataset. My novice attempt at this task is shown below csg<-NULL for (i in 1:ncol(data.all)) {

Dear all, Sorry if this is slightly off the track as far as R is concerned, but I have been using the Savitzky-Golay filter to estimate some derivatives of interest. I am wondering however, if anyone has seen anything in the literature (or has any ideas) of how these estimates perform asymptotically. Does anyone know what the rate of convergence is for these as the sample size increases?

Hi I'm a B.E student pursuing my Project in the CEERI unit of the Council of Scientific and Industrial Research, Chennai, TN, India. I'm working on R -Language for my project. I'm in need of the functionality of SIMCA for the project. Is there any in-built function for SIMCA in R? Or are there people working on it? It would be a great help if you let us know as soon as possible.

Hi, I wonder if someone could help. I needed to transfer (copy) a workspace file that had been generated in linux (R 2.11) to windows running the same version of R 2.11 (but of course windows binary). Usually, there is no problem in doing this and all objects work as expected. I am often doing this to be able to produce wmf or emf graphic files that I need. This time I had some spectra that I

Hi the list. I do not understand the philosophy behind numeric and integer. - 1 is numeric (which I find surprising) - 2 is numeric. - 1:2 is integer. Why is that ? Christophe

Hi R-users, I'm trying to produce decompositions of a multiple time-series, grouped by a factor (called "area"). I'm modifying the code in the STLperArea function of package ndvits, as this function only plots produces stl plots, it does not return the underlying data. I want to extract the trend component of each decomposition ("x$time.series[,trend]), assign a name

Is the project on creating R GUIs using QT interfaces still on? Any plans of using PyQT Regards Ajay Ohri Websites- http://decisionstats.com http://dudeofdata.com Linkedin- www.linkedin.com/in/ajayohri On Tue, Nov 9, 2010 at 8:33 PM, Kari Ruohonen <kari.ruohonen at utu.fi> wrote: > Hi, > I wonder if someone could help. I needed to transfer (copy) a workspace > file that

Suppose some environment variable (say MY_R_INC) has a number of paths. I want to source some file relative to some path in $MY_R_INC (just as #include in C++ does, which looks for header file in a number of directories). I RSiteSearch'ed, but I don't find any function that satisfies my need. Could somebody let me know if I overlooked something?

Hi, I?m trying to use the STLperArea function in the ndvits package. I can run the sample data (?SLPSAs_full?) without any problem. However, when I come to run the function on my own data, I get the following message. Waiting to confirm page change... Error in .checkPath(path) : no TISEAN executables found in that directory. Please set a proper TISEAN executables path using

I want to run deconvolution of a time series by an impulse or point-spread function through Wiener filter, regularized filter, Lucy-Richardson method, or any other approaches. I searched the CRAN website and the mailing list archive, but could not find any package for such a deconvolution analysis. Does anybody know an existing R function for deconvolution? TIA, Gang

Dear All, I am a novice when it comes to time-series analysis and at the moment I am actually interested in calculating the Hurst exponent of a time series. This question has already been asked quite some time ago http://bit.ly/98dZsi and I trust some progress has been made ever since. I was able to find some functions in the packages http://cran.r-project.org/web/packages/Rwave/index.html

Hi there, How may I smooth spline two vectors with the smoothing parameter selected by generalized maximum likelihood (GML) .? Thanks a lot. -- View this message in context: http://r.789695.n4.nabble.com/Smoothing-spline-with-smoothing-parameters-selected-by-generalized-maximum-likelihood-tp4210679p4210679.html Sent from the R help mailing list archive at Nabble.com.

Dear all, I am new to R, and I am a biotechnologist, I want to fit a smoothing spline with smoothing parameter selected by generalized maximum likelihood. I was wondering what function implement this, and, if possible how I can find the fitted results for a certain point (or predict from the fitted spline if this is the correct language) -- View this message in context:

Hello all, I am looking for some help in understanding how to change the way R optimizes the smoothing constant selection process for the HoltWinters function. I'm a SAS veteran but very new to R and still learning my way around. Here is some sample data and the current HoltWinters code I'm using: rawdata <- c(294, 316, 427, 487, 441, 395, 473, 423, 389, 422, 458, 411, 433, 454,

Hello everyone, I was wondering if anyone was aware of a way in which I could use ggplot's stat_smooth() function for add an exponential moving average. I was thinking that I could maybe use something like: >myggplot + stat_smooth (method = 'HoltWinters( data , .9 , 0, 0)') but my efforts were futile. Perhaps there is a way to write my own custom method to throw into