similar to: Update: Is there an implementation of loess with more than 3 parametric predictors or a trick to a similar effect?

I noticed a problem using R 2.7.1 on Windows XP SP2 with the precompiled Atlas Rblas.dll. Running the code below causes R to crash. I started R using Rgui --vanilla and am using the precompiled Atlas Rblas.dll from cran.fhcrc.org dated 17-Jul-2007 05:04 for Core2 Duo. The code that causes the crash: x <- rnorm(100) y <- rnorm(100) z <- rnorm(100) loess(z ~ x * y) loess(z ~ x) does

Dear R experts, I have a problem that is a related to the question raised in this earlier post ??? https://stat.ethz.ch/pipermail/r-help/2007-January/124064.html My situation is different in that I have only 2 predictors (coordinates x,y) for local regression but a number of global ("parametric") offsets that I need to consider. Essentially, I have a spatial distortion overlaid over a

Dear R experts, I have a problem that is a related to the question raised in this earlier post https://stat.ethz.ch/pipermail/r-help/2007-January/124064.html My situation is different in that I have only 2 predictors (coordinates x,y) for local regression but a number of global ("parametric") offsets that I need to consider. Essentially, I have a spatial distortion overlaid over a

Greetings All! After much trial and error, and an exhaustive search of the archives, I'm writing to solicit help with a problem. I want to dynamically modify variable names and values as function arguments inside a loop. I have a canned function from a package that takes ellipsis (...) arguments of the form varb=vval where varb is a name and vval is either numeric or character. I want to

Is there an R implementation of a scheme for automatic smoothing parameter selection with loess, e.g., by minimizing one of the AIC/GCV statistics discussed by Hurvich, Simonoff & Tsai (1998)? Below is a function that calculates the relevant values of AICC, AICC1 and GCV--- I think, because I to guess from the names of the components returned in a loess object. I guess I could use

I would like to replace a few varaibles within a data frame. For example, in the dataframe below (contrived) I would like to replace the current housesize value only if the Location is HSV. However, I would like to leave the other values intact. I tried "ifelse", but I don't really need the else condition. test_data2_df<-data.frame(Variables=c("SQR

Hi there, I would like to draw a scatter plot and fit a smooth line by loess. Below is the data. However, the curve line started from 0, which my "resid" list doesn't consist of 0 value. It returned some warnings which I don't know if this is the reason affecting such problem. Here I also attached the warning messages. Please let me know if there is a solution to fix this. Thank

Dear John, > I suggest that you look at the abilities of the mgcv package. > There are notes of mine at > > http://www.maths.anu.edu.au/%7Ejohnm/r-book/xtras/autosmooth.pdf > > that may help you get started. Thank?you very much for the suggestion and the link to your write-up, it was indeed very helpful! I have experimented with this library for a while now and am really happy

Hello list, I just started R today and tried something quite simple. I wanted to create a colored plot and eventually after hours of fiddling around got it working. However, my solution seems very suboptimal and I'd really appreciate your hints on how to improve. I believe that R already offers many functions I coded (e.g. distance between two vectors, vector length, vector normalization and

Hi, I am segfaulting when using predict.loess() (checked with r62092). I've traced the source with the help of valgrind (output pasted below) and it appears that this is due to int overflow when allocating an int work array in loess_workspace(): liv = 50 + ((int)pow((double)2, (double)D) + 4) * nvmax + 2 * N; where liv is an (global) int. For D=1 (one x variable), this overflows at

Good morning, I am investigating dll import in R under Windows XP. Using examples I found on the internet, I started with a very simple dll, e.g. including only the basic function: void { *x2 = x*x; }sqr(doublex, double*x2) I compiled it as a dll with Eclipse and Cygwin's gcc. It works when I call it with another simple .exe C program, compile with Eclipse and gcc as well. I can do what I

Hi Larry, Thanks so much! This seems to do the trick. I would have spun my wheels for > a long time before discovering all of this, wow. No worries. :) I'll try to keep this in mind and make sure I address it in future Kaleidoscope tutorial chapters - these issues tripped me up the first time I encountered them too. Do I even want to know what additional chickens need to be sacrificed

>> llvm::sys::DynamicLibrary::LoadLibraryPermanently(nullptr) This is one is a bit tricky and hard to find. I spent quiet some time digging into MC and ORC JIT execution engines trying to find what makes them work. The problem is that this trick (LoadLibraryPermanently) happens inside of EngineBuilder, despite that the functionality belongs to a JIT engine itself, not to the builder. I

Thanks so much! This seems to do the trick. I would have spun my wheels for a long time before discovering all of this, wow. Do I even want to know what additional chickens need to be sacrificed to get this to work on Windows? -- lg > On May 18, 2016, at 1:52 PM, Lang Hames <lhames at gmail.com> wrote: > > Hi Larry, > > You're basically there, but you're hitting

Following a discussion about numerical performance on comp.lang.functional recently I just tried running a simple C mandelbrot benchmark that uses C99's complex arithmetic using gcc and llvm-gcc on a 2.1GHz Opteron 2352 running Debian: gcc: 5.727s llvm-gcc: 1.393s There is still 20% room for improvement but LLVM is >4x faster than gcc here. Sweet. Here's the code: #include

I was having some problems persuading optim() to give me the answers I wanted, & simplified down to: sqr<-function(x){(x+1)^2} optim(1,sqr) I accept this is a hammer to crack a nut, but was still expecting the answer -1. I got: $par [1] -0.8 $value [1] 0.04 $counts function gradient 12 NA $convergence [1] 0 $message NULL so I've

When using ORC JIT, I'm having trouble with external function resolution (that is, of a function defined in the app, with C linkage). I add a declaration for the function to my IR, and when I use MCJIT, it finds it and all is well, But when I use ORC JIT (I *think* correctly, at least it closely matches what I see in the tutorial), I get an LLVM error, "Program used external function

Dear R-users, I have some trouble to perform a 2-samples KS test. Apparently my 2 samples are numerical (see below) but R complains that the "y" term is not.... ------------------------------------------------------------------------------------- > str(subset(mydata,Identified=="NO",select=KD)) 'data.frame': 2889 obs. of 1 variable: $ KD: num -0.272 -0.080

On gcc's side, this is a simple missed opt on the part of builtin lowering. As a result, the gcc code ends up with a call to muldc3 (complex = 2x2 multiply double) and the llvm code doesn't. GCC should be fixed in a second, and with that, there is no appreciable performance difference between the two. On Thu, Feb 26, 2009 at 4:07 PM, Jon Harrop <jon at ffconsultancy.com> wrote: >

Fit NLE - was: [R] computer algebra in R Original post: https://stat.ethz.ch/pipermail/r-help/2023-November/478619.html Dear Kornad, I think I have started to understand what you try to achieve. The problem is to fit a NLE and compute the parameters of the NL-Eq. I have included the R Help-list back in the loop, as I am not an expert in optimization. Goal: y ~ I0 + IHD * hd + ID * d; where: y