search for: yts

...Search("multiple series plot") i found this post by Gabor Grothendieck: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42281.html Exactly what i need except for one detail. I want one series to be made of points and the other by a line. When I simply try: xts <- ts(x$b,start=x$a[1]) yts <- ts(y$b,start=y$a[1]) ts.plot(xts,yts,type=c("p","l")) or if i use the 'type' argument inside gpars=list() i get an error of "invalid plot type". Then if i try ts.plot(xts,yts,type1="p",type2="l") i get warnings about NAs introduced by...

Hi: I need to plot two data sets and then I need to change labels of xaxis. I can manage to achieve something like it with this: x<-data.frame(a=c(1,2,3,4,5),b=c(2,4,6,8,10)) y<-data.frame(a=c(3,4,5),b=c(1.5,2,2.5)) xts <- ts(x$b,start=x$a[1]) yts <- ts(y$b,start=y$a[1]) ts.plot(xts,yts,col=c("red","blue")) But ts.plot does not allow to change x labels. Besides, i need dots. That's the cause I would lke to know if there is some way to trick plot to do this. Thanks in advance. -- Felipe [[alternative HTML vers...

...tings, I would like to estimate a spectral coherence between two timeseries. The stats : spectrum() returns a coh matrix which estimates coherence (squared). A basic test which from which i expect near-zero coherence: x = rnorm(500) y = rnorm(500) xts = ts(x, frequency = 10) yts = ts(y, frequency = 10) gxy = spectrum( cbind( xts, yts ) ) plot( gxy $ freq, gxy $ coh ) yields a white spectrum of 1. Clearly i'm not using this correctly... or i mis-interpret the coh as a cross-spectral density estimate of coherence |Gxy|^2/(Gxx Gyy) Thanks in advance!

.... The data files was enclosed in the email. I'm sorry the errors message appeared in chinese, but it says "plot.xy(xy.coords(x, y), type = type, ...) : errors in argument has more than 3 matches" > y=scan('F:/Math355/AIRLINE.txt',skip=2,nlines=18) Read 144 items > yts=ts(log(y)) > plot(yts,main="LOG AIRLINE TOTALS",type='p',col=2) > n=length(y) > time=seq(1:n) > > month=c(rep(seq(1:12),12)) > fmonth=as.factor(month) > ymod1=lm(yts~time+fmonth) > summary(ymod1) > yfit1=ymod1$fitted > lines(yfit1,c=3) ??? plot.x...

Hi, I am fairly new to R so this may be easy for some. I have one column that is the Type="A, B' T" and I have another column that lists Price='894, 895, 896' I wanted to plot three different series based upon the Type...If B than xts..If T yts... How exactly do I do this? Thanks in advance -- View this message in context: http://www.nabble.com/Plotting-3-Time-Series-tp17405509p17405509.html Sent from the R help mailing list archive at Nabble.com.

...unning CentOS 6, all updates applied, maybe a week or two ago. ### HUME ### Hume has a backup saved in /home/spfs.450, the directory that I'm trying to export. [root at hume ~]# ll /home/spfs.450/ drwxr-xr-x 3 apache apache 3 Oct 8 2009 y.spfs drwxr-xr-x 3 apache apache 3 Feb 1 2010 yts.spfs --SNIP-- Hume is exporting with /etc/exports /home/spfs.450 192.168.254.0/255.255.255.0(rw,async,no_subtree_check,mp,no_acl,insecure,no_root_squash) Hume has appropriate /etc/passwd entries: [root at hume ~]# grep -i apache /etc/passwd apache:x:48:48:Apache:/var/www:/sbin/nologin To be sur...

I think this is another one of the same kind of bugs in ts: Version 0.64.1 (May 8, 1999) ... > z <- ts(matrix(1:20,10,2), start=c(1969,1), frequency=12) > (z > 5) | (z < 2) Error: invalid time series parameters specified > Paul -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read

say we have: > x<-data.frame(a=c(1,2,3,4,5),b=c(1,1,1.5,2,2)) > y<-data.frame(a=c(1,2,3,4,5),b=c(1,2,2,3,3)) How would I plot this so that, with the shared $a as the x-axis values, I have both $b columns plotted together? (a comparison of the two?) thanks. DISCLAIMER: This e-mail message and any attachments are inte...{{dropped}}