similar to: Optimal parameters for Savitzky-Golay smoothing filter (loop)

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 When we use the sav_gol command in R , it shows an error which says: " error in as.matrix". We've downloaded the necessary packages. Kindly help us with this issue. If there is any other function to perform Savitzky-Golay smoothing in R, please let me know. With Regards Reynolds

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.

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

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 as the sample size increases?

Hi I have dataset that consists of two columns AB 0.102 AC -0.002 BA -0.102 BC 0.270 CA 0.002 CB -0.270 I wish to create a matrix so that I can eventually plot the data. A B C A 1 0.102 -0.002 B -0.102 1 0.27 C 0.002 -0.27 1 Any help or guidance would be greatly

My apologies, to elaborate I carried out a correlation analysis and I want to plot the data (similar to the graphs available in the corrplot package. However, the results are as follows (I have many more combinations): AB 0.102 AC -0.002 BA -0.102 BC 0.270 CA 0.002 CB -0.270 I now want to create a matrix that uses the names from the

Dear All, I have a pair of .cpp and .def file can be compiled using VC++ and works perfectly well in S-PLUS. I wanted to do the same for R; so I followed the guidline given in "Building from a source-code library under Windows" as much as possible and manage to compile them using VC++ and call it from R. But it gives different answer from the one called from S-Plus. I know that I did

Hello, I want an instance of CallGraph in my pass. By looking at -dot-callgraph source, I've tried something like this: CallGraphWrapperPass *CGWP = new CallGraphWrapperPass(); PM.add(CGWP); CallGraph *CG = &CGWP->getCallGraph(); PM.add(new MyPass(CG)); I get the following error: /home/riyad/installs/llvm-3.7.0/include/llvm/PassSupport.h:95:38: error: no matching constructor for

Hi As the title says, are all the debug nodes of LLVM IR avaialble in DWARF? If I add a new custom MDNode in LLVM IR, can I expect the debug data will be available in DWARF in the binary if I compile in debug mode? Thanks Riyad -------------- next part -------------- An HTML attachment was scrubbed... URL:

Hi All, I want to change debug information of an llvm instruction so that the modified debug info is subsequently passed to executable binary. So if I use "addr2line" utility on the binary, it will return my modified debug information. I've tried to change by using the following code snippet: MDNode *N = Inst->getMetadata("dbg"); DebugLoc Loc =

Hi, There is a "Root" private member in CallGraph <http://llvm.org/docs/doxygen/html/CallGraph_8h_source.html#l00192> class. But there is no way to get access to that member. How can I do that? I can get CallGraphNode of "main" method. But I am curious is there more cleaner way. Thanks, Riyad -------------- next part -------------- An HTML attachment was scrubbed...

I have some data coming from SQL sources that I wish to relate in various ways. For reasons only known to our IT people, this can't be done in SQL at present. I am looking for an R'ish technique for looking up new columns on a data frame. As a simple, hardwired example I have tried the following: # This gives me two columns, one the lookup value and the second one # the result column,

Hi, I want to add "mem2reg" pass similar to this: PassManager PM; PM.add(new LoopInfo()); PM.add(new Mem2Reg()); // What's the class name for this pass I couldn't find the name of corresponding class and header file of mem2reg pass except "mem2reg.cpp". So how can I add mem2reg pass into my pass manager? Thanks, Riyad -------------- next part -------------- An

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

Greetings. I am working on setting up an asterisk server (SIP only) and am running into a few issues getting RTP working correctly. Here is our setup: SIP Client (Public IP) <---> Asterisk Server (Public IP/Private IP) <--> Nortel CSG (Internal IP) <--> PSTN So far we have SIP to SIP working through Asterisk without any problems (using various sip clients). When I call from

http://bugzilla.mindrot.org/show_bug.cgi?id=1340 Summary: Support for Camellia block cipher to OpenSSH-portable. Product: Portable OpenSSH Version: -current Platform: All OS/Version: All Status: NEW Severity: enhancement Priority: P2 Component: Miscellaneous AssignedTo: bitbucket at mindrot.org

Have been running samba successfully authenticating to a windows 2003 domain since 3.0.1. Starting in 3.0.23 and 3.0.23b I can't do a "net ads join" on a HPUX itanium server running 11.23 ia64. I can kinit just fine w/ this userid. Samba was built with gcc 4.1.1. See below: root@serv00 # kinit jjurich_wa Password for jjurich_wa@DIVMS.UIOWA.EDU: root@serv00 #