similar to: plot.ts bad default labelling and lag.plot incomplete doc

Hi all. Currently, if you try: > lag.plot(1:10) you get superposed labels '1' and '10'. Things go worse in more extreme cases: x <- ts(1:10) x1 <- lag(x, 4) plot(x1, x) This is due to a mistake in plot.ts. My suggestion is the following really minimal patch to plot.ts: @@ -530,7 +530,7 @@ plot.ts <- text(xy, labels =

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

Hi all. I'm looking for robust ways of building lagged variables in a dataset with multiple individuals. Consider a dataset with variables like the following: ## set.seed(123) d <- data.frame(id = rep(1:2, each=3), time=rep(1:3, 2), value=rnorm(6)) ## >d id time value 1 1 1 -0.56047565 2 1 2 -0.23017749 3 1 3 1.55870831 4 2 1 0.07050839 5 2 2 0.12928774 6

This in an enhancement requestion, to add the capability to pass arguments through barplot() to legend(). I've created a version that does this. # Run this code; notice that the legend # is superposed over one bar: barplot(height = cbind(x = c(465,91) / 465 * 100, y = c(840,200) / 840 * 100, z = c(37,17) / 37 * 100), beside = FALSE, width =

Dear R users, I am trying to put labels on the data points on a ts.plot. I have tried: ts.plot(df.ts, gpars=list(xy.labels=colnames(df.ts), xlab="group", ylab="level")) which plots fine but gives no labels xy.labels=TRUE also plots without labels I also tried using text() but got an error Error in xy.coords(x, y, recycle = TRUE) : Argument "x" is

Dear R developers, these days I'm working on some R code for fitting completely generic Bayesian Hierarchical Models in R, a la OpenBUGS and JAGS. A key feature of OpenBUGS and JAGS is that they automatically build an appropriate MCMC sampler from a generic model, specified as a directed acyclic graph (DAG). The spirit of my (would-be) implementation is instead more focused on experimentation

Dear R-helpers! Using logarithmic axes sometimes results in an odd labelling. In some cases, only half of the axes do have labels. Hmmmm. I've played around with par(lab=c(5,5,7)) but that didn't help. An example: ---- begin example ---- x <- c(0.12, 0.3, 0.53, 1.1, 1.8, 2.6, 4.5) y <- c(0.012, 0.021, 0.032, 0.055, 0.1, 0.2, 0.35) plot(x, y, log="xy",

(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 limited to 0-1. But: >

Example: Will ts objects be obsolete or modified? > a [,1] 1983 Q1 2.747365190 1983 Q2 2.791594762 1983 Q3 -0.009953715 1983 Q4 -0.015059485 1984 Q1 -1.190061246 1984 Q2 -0.553031799 1984 Q3 0.686874720 1984 Q4 0.953911035> lag(a,4) [,1] 1983 Q1 NA 1983 Q2 NA 1983 Q3 NA 1983 Q4 NA 1984 Q1 2.747365190 1984 Q2

On Tue, 2016-02-09 at 16:56 -0600, Spencer Graves wrote: > Hello: > > > I'm having trouble plotting an object of class "ts" that is in a > data.frame. I can do it with(data.frame, plot(...)) but not with > plot(..., data.frame); see the example below. The plot function is generic so the actual function call depends on what arguments you give it:

This is my reproducible example tv.ms<-structure(list(inq = structure(4:17, .Label = c("D4", "D5", "D6a", "D6b", "D6c", "D7", "D8", "F4", "F5a", "F5b", "F6a", "F6b", "F6c", "F6d", "F7a", "F7b", "F8"), class =

Hi All, I am busy using R to do some regression modelling and have been using plot(x,y,"") to visualise my variables. I would now like to label my points using data stored in the data.frame used for the regression analysis. For example each of my data points is made up of a field measured forest volume value and a remotely sensed vegetation estimate (NDVI). Each point is an enumeration

I've enconuntered this problem with the last cran version of geoR: > library(geoR) > day <- rep(1:2, each=5) > coords <- matrix(rep(runif(10),2), 10, 2) > data <- rnorm(10) > data[1] <- NA > as.geodata(cbind(coords, data, day), realisations=4) as.geodata: 1 points removed due to NA in the data Errore in as.geodata(cbind(coords, data, day), realisations = 4) :

Dear all, I think I have found a buglet in plot.ts plot.ts(x=1,type="n") # correct: does not show the plot plot.ts(x=1,y=1,type="n") # not correct: does show the plot I did not investigate the problem in depth but it could be related to the switch xy.labels, in fact plot.ts(x=1,y=1,type="n",xy.labels=TRUE) # does show the plot

Dear list, I've noticed that, when writing examples in an Rd file, you need to escape the '%' character, even if it is valid R code. I can see maybe this is the intended behaviour, but I found it a bit surprising, and leads to cryptic error messages from 'R CMD check'. The relevant section in the "writing R extensions" manual isn't clear on that point either. To

I've rolled up R-2.6.1.tar.gz a short while ago. This is a maintenance release and fixes a number of mostly minor bugs and platform issues, as well as a couple of major problems. Probably, the most important issues were - save-to-postscript on Windows - tcltk event loop with X11 - build issue on recent MacOSX See the full list of changes below. You can get it (in a short while) from

I've rolled up R-2.6.1.tar.gz a short while ago. This is a maintenance release and fixes a number of mostly minor bugs and platform issues, as well as a couple of major problems. Probably, the most important issues were - save-to-postscript on Windows - tcltk event loop with X11 - build issue on recent MacOSX See the full list of changes below. You can get it (in a short while) from

Hello: I'm having trouble plotting an object of class "ts" that is in a data.frame. I can do it with(data.frame, plot(...)) but not with plot(..., data.frame); see the example below. This work around gets me past this problem. However, I thought the R Core team might want to know about this if they don't already. Thanks for all your work in making

Hi all. A statistical question. I have to estimate the variance of the sum: X(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.

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 has to 'operate elementwise'). A cross-reference from