similar to: Poly-phase Filters in R

Hello, I've been using jpeg(), bmp() and pdf() to capture plots. I've used the parameter "at" in a plot, to define the tickmarks. It works fine on screen, but when I try to print it to a file, it gives a warning: "at" is not a graphical parameter and prints an empty figure. Can you help? Thanks, Eran. [[alternative HTML version deleted]]

Hello, I've switched to R studio from the StatET Eclipse plug-in. I have a question regarding navigating between plots. When I use x11() or windows() new devices are created and I know how to switch back and forth between them. However, when I plot on the device that stands for R-Studio's built-in plot browser, is there a way to switch back between plots? Each new "plot"

Hello, I have a time-series that has some missing samples. I was thinking on completing them using either zero-order hold or linear interpolation. I am looking for an efiicient way (other than a loop...) of identifiying the missing time slots and filling them. Can you think of any methods that might help here? (obviously which(diff(time)>min(diff(time))) will give the locations, but what

Hi, I am not sure if this is the right list to ask this question (though I did not find a more appropriate one). I've started using R a month ago, and small scripts work fine. However, when I start writing more complex code, it gets messy. 1. Is there any way to debug "normally", with breakpoints? 2. I am using the Eclipse plugin (StatET), and tried JGR(). Is there an IDE that

Hello, I wonder how I might create a package that only reveals some of the function in the package to the user. I've tried creating an R package using the following: f <- function(x,y) x+y g <- function(x,y) x-y h <- function(x,y) f(x,y)*g(x,y) package.skeleton(list=c("f","g","h"), name="mypkg") and would like only h() to be available when I

Hello all, I need to zoom in and out and "travel"(pan) inside a plot, like you can do on a Matlab plot. If possible, I would also like the option to use the mouse to set a marker on the graph and get the (x,y) data for it, again, like in Matlab. Is this possible in R with the regular packages, or do you maybe know a different package that will allow this? Eran. * * [[alternative

Hello, I have a function with a long list of parameters (of different types, numeric and string) myFunc <-function(p1, p2, p3, p4, p5...etc) { do.something(p1,p2,....) } I want to loop over this to provide a different set of parameters to the list every time. for (ii in 1:N) { myFunc(p1(ii), p2(ii),....etc) } I would like to simplify the notation and use some kind of structure, maybe

Hi all, Is there a way to create a "dictionary" in R, such that it has <key,value> pairs? Something to the effect of: x=dictionary(c("Hi","Why","water") , c(1,5,4)) x["Why"]=5 In truth I'm looking two categorial variables function. So that if x=dictionary(c("a","b"),c(5,2)) x val 1 a 5 2 b 2 I want to compute

To add to the wish-list for "convolve": For modeling processes that decay exponentially in time, e.g., fluorescence, it is desirable to have a function that convolves an arbitrary vector with an exponential using an iterative method. In the TIMP package (which won't be on CRAN till R 2.5.0 is official, but is for now at www.nat.vu.nl/~kate/TIMP) we implemented this special-purpose

I have been experimenting with convolve(). What I know about convolution I learned from engineering (they call this stuff linear systems theory), not statistics, so maybe this is all just a matter of different conventions. BUT I notice very weird things with convolve(). 1. First example, from the classic Bracewell The Fourier transform and its applications, chap 3 (p.32 in 2nd edition): {2 2 3 3

Dear R People: Here is something that I am sure is very simple. I'm just trying to re-create the C convolution example in the Extensions manual. Here is the subroutine: void convolve(double *a, int *na, double *b, int *nb, double *ab) { R_len_t i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j < *nb; j++) ab[i + j] += a[i] *

Dear R users, I want to call C code via the .C or .Call interface. This works fine with integer values but using doubles the array received in C will be set to zeros. I have tried the convolve examples (Writing R extensions, chapter 5.2) and still the resulting array consists of zeros. My code (shortened for my purposes. Original did not work either):

We obtained some disturbing results from convolve() (inaccuracies and negative probabilities). We'll try to make the context clear in as few lines as possible... Our function panjer() (code below) basically computes recursively the probability mass function of a compound Poisson distribution. When the Poisson parameter lambda is very large, the starting value of the recursive scheme ---

Hi, I compiled a C program file on Borland C++ 5.5 compiler to get one dll output (as instructed in the file readme.package). The C program file is just the example on page 31 of "writing R extensions": void convolve(double *a, int *na, double *b, int *nb, double *ab) { int i, j, nab = *na + *nb - 1; for(i = 0; i < nab; i++) ab[i] = 0.0; for(i = 0; i < *na; i++) for(j = 0; j

Hi, Man page for 'convolve' says: conj: logical; if 'TRUE', take the complex _conjugate_ before back-transforming (default, and used for usual convolution). The complex conjugate of 'x', of 'y', of both? In fact it seems that it takes the complex conjugate of 'y' only which is OK but might be worth mentioning because (1) conj=TRUE is the

Submitted to bugzilla as PR 15729 From: llvmdev-bounces at cs.uiuc.edu [mailto:llvmdev-bounces at cs.uiuc.edu] On Behalf Of Kaylor, Andrew Sent: Thursday, April 11, 2013 1:13 PM To: Weiss, Eran Cc: llvmdev at cs.uiuc.edu Subject: Re: [LLVMdev] Migration from JIT to MCJIT Thanks, Eran. I’m not sure how soon I’ll have a solution for you, but it’s on my to-do list now. I’ll also create a bugzilla

Hi, the function stats::convolve does not mention efficient usage of the underlying FFT algorithm, such as (a) if type="circular", then length(x)=length(y) should have many factors (e.g. length(x) = length(y) = 2^n) (b) if type="open" or "filter", then length(x)+length(y)-1 should have many factors (e.g. length(x)+length(y)-1 = 2^n) In particular the latter may

Hi, I'm new to R, and I'm a bit confused with the "convolve()" function. If I do: x<-c(1, 2, 3) convolve(x, rev(x), TRUE, "open") = 9 12 10 4 1 But I expected: 3 8 14 8 3 (like in Octave/MATLAB - conv(x, reverse(x)) ) 3 2 1 x 1 2 3 = 3 2 1 0 6 4 2 0 0 9 6 3 = 3 8 14 8 3 The thing is, that "convolve(x, x, TRUE, "open")" works. For me

Dear R-ophiles, I've found something very odd when I apply convolve to ever larger vectors. Here is an example below with vectors ranging from 2^11 to 2^17. There is a funny bump up at 2^12. Then it gets very slow at 2^16. > for( i in 11:20 )print( system.time(convolve(1:2^i,1:2^i,type="o"))) user system elapsed 0.002 0.000 0.002 user system elapsed 0.373

Thanks, Eran. I’m not sure how soon I’ll have a solution for you, but it’s on my to-do list now. I’ll also create a bugzilla record for this problem. -Andy From: Weiss, Eran [mailto:Eran.Weiss at emc.com] Sent: Thursday, April 11, 2013 12:40 AM To: Kaylor, Andrew Cc: llvmdev at cs.uiuc.edu; Jim Grosbach; Jiong Wang Subject: Re: [LLVMdev] Migration from JIT to MCJIT Andrew, I've attached