search for: narzo

This code gives an error: a <- ts(1:10, start=0, freq=10) b <- ts(1:10, start=1, freq=10) ts.intersect(a,b) This one works normally (and correctly): a <- ts(1:10, start=0) b <- ts(1:10, start=1) ts.intersect(a,b) Antonio, Fabio Di Narzo. P.S. How to switch off italian error messages to post on r-help? > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 1.0 year 2005 month 04 day 18 language R

(Moved from r-devel to r-bugs) On 3/24/2006 5:03 AM, Antonio, Fabio Di Narzo wrote: > Hi all. > There's a bug in plot.acf, when plotting acf for multivariate time series. > Here a reproducible example: > > X <- rnorm(1000) > Y <- -X + rnorm(1000, sd=0.6) > Z <- cbind(X,Y) > > In > acf(Z) > cross-correlation plot y-axis is limi...

...ng deriv or numericDeriv (or something else). My difficulties here are that in deriv and numericDeric the function is passed as an expression, and one have to pass the list of variables involved as a char vector... So, it's a pure R programming question. Have a nice sunday, Antonio, Fabio Di Narzo.

...'] == d1[d1$time==2,'b'] #This is TRUE, as expected ##Try to compare d and d1 now: they look right Any hint on what's wrong here? By now, my workarond is changing variable names before reshaping, than re-assign old variable names back after reshape. Best regards, Antonio, Fabio Di Narzo. > R.version _ platform i686-pc-linux-gnu arch i686 os linux-gnu system i686, linux-gnu status major 2 minor 6.0 year 2007 month 10 day 03 svn rev 43063 language R version.string R ve...

Full_Name: Antonio, Fabio Di Narzo Version: 2.6.1 OS: linux Submission from: (NULL) (213.140.16.187) In 'outer' man page, there is no mention of the Vectorize function. Moreover, I think it isn't underlined enough that the FUN argument to 'outer' must be a vectorized function (doc speaks about a function which...

...ind(coords, data, day), realisations = 4) : realisations and coords have incompatible dimensions The problem disappear if I remove the NA manually from the dataset before passing to as.geodata. I.e.: > as.geodata(cbind(coords, data, day)[2:10,], realisations=4) works. Antonio, Fabio Di Narzo. [[alternative HTML version deleted]]

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060525/ef972e9f/attachment.pl

...xing data.frame), but this really isn't the point. The point is that I have somewhere a ready sequence of IDs from which obtaining corresponding data. Note that I don't know if some ID combination exists in the database. Some suggestions on a good way to do this? Tnx all, Antonio, Fabio Di Narzo. [[alternative HTML version deleted]]

...<- rep(times, each=2) vals <- rep(1:2,5) tapply(vals, id, mean) Under R-2.2.1 this gives: 1970-01-02 1970-01-03 1970-01-04 1970-01-05 1970-01-06 1.5 1.5 1.5 1.5 1.5 But under R-2.3.0 the output is: 1 2 3 4 5 1.5 1.5 1.5 1.5 1.5 Antonio, Fabio Di Narzo. [[alternative HTML version deleted]]

...<- rep(times, each=2) vals <- rep(1:2,5) tapply(vals, id, mean) Under R-2.2.1 this gives: 1970-01-02 1970-01-03 1970-01-04 1970-01-05 1970-01-06 1.5 1.5 1.5 1.5 1.5 But under R-2.3.0 the output is: 1 2 3 4 5 1.5 1.5 1.5 1.5 1.5 Antonio, Fabio Di Narzo. [[alternative HTML version deleted]]

...nt speed (w.r.t. JAGS) on some example models, and I'm writing a proof-of-concept, reproducible Sweave file about it, to be published online shortly. What do you think about it in general? What do you think about developing an R package of it as a GSoC project? Best regards, Antonio, Fabio Di Narzo. -- Antonio, Fabio Di Narzo Ph.D. student at Department of Statistical Sciences University of Bologna, Italy

...1) + X(2) + ...+ X(n) from an observed sample, where X(i) are *correlated* and not necessarly identically distributed. Someone can suggest a simple strategy (I hope by exploiting some already present R package) for obtaining such estimate from an observed vector X[1:n]? Tnx all, Antonio, Fabio Di Narzo. [[alternative HTML version deleted]]

Dear all, When using gdb to debug my C code, I use R_PV to show the content of SEXP variables: SEXP sexp; // it is a data.frame SEXP colNames = getAttrib(sexp, R_NameSymbol); A strange thing is that after halting the program: (gdb) p R_PV(colNames) does not show the content of colNames. I am positive my code is right because if I insert "PrintValue(colNames);" in the c code,

...i686 os linux-gnu system i686, linux-gnu status Patched major 2 minor 8.0 year 2008 month 12 day 04 svn rev 47063 language R version.string R version 2.8.0 Patched (2008-12-04 r47063) -- Antonio, Fabio Di Narzo Ph.D. student at Department of Statistical Sciences University of Bologna, Italy

...working. However, I was thinking if there is a smarter/cleaner/more robust way to do the job. For instance, with the above function I get dataframe rows re-ordering as a side-effect (anyway this is of no concern in my current analysis)... Any suggestion? All the bests, Fabio. -- Antonio, Fabio Di Narzo Ph.D. student at Department of Statistical Sciences University of Bologna, Italy

Dear all, I have the following problem. Given an expression object 'expr' containing a certain set of symbols (say 'a', 'b', 'c'), I would like to translate the expression object in an R function of, say, 'a', programmatically. Here an example of what I mean. Given: > expr <- expression(a+b+c) a call like: > asFunctionOf(expr, 'a',

Hello everybody, I am wondering if it is possible to create a gui to plot a time series that is very big, it's an EEG signal of 20mins. What I would like to do is plot the first 5mins, then have a button on the gui that plots the next 5mins when pushed. Is it possible? Thanks in advance ! Gael.

...function `_start': ../sysdeps/i386/elf/start.S:98: undefined reference to `main' collect2: ld returned 1 exit status make: *** [MCMCpack.so] Error 1 ERROR: compilation failed for package 'MCMCpack' I don't know why it searches a reference to 'main'... Antonio, Fabio Di Narzo. > version _ platform i386-pc-linux-gnu arch i386 os linux-gnu system i386, linux-gnu status major 2 minor 1.1 year 2005 month 06 day 20 language R

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060604/a8854663/attachment.pl

...("temp.dat") In other words, my real problem should be restated as: how to safely dump (and then load) an object which has a (perhaps valid) "tsp" attribute? More generally: can someone suggest me a safer way to dump/restoring R objects? Tnx all, Antonio. -- Antonio, Fabio Di Narzo Ph.D. student at Department of Statistical Sciences University of Bologna, Italy