similar to: Time series

Dear R-help, How could a cross-correlation plot be optimized such that the relationship between seasonal time-series can be studied? We are working with strong seasonal time-series and derived a cross-correlation plot to study the relationship between time-series. The seasonal variation however strongly influences the cross-correlation plot and the plot seems to be ?rather? symmetrical (max

HI, May be this helps: ?ts1<- ts(1:20) ?ts2<- ts(1:25) ts1[-(1:3)]<- ts1[-(1:3)]+ts2[1:17] ?as.numeric(ts1) # [1]? 1? 2? 3? 5? 7? 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 A.K. Hey everyone, I`m an absolut beginner in R and need some help for an exercise: I want to do ordinary calculations with 2 time series. The issue with this, that I want to use different elements of time

Dear Helpers, I am looking for methods and tools to compare and then to concatenate two successive time series. They are both in the same frequency and they describe one phenomena. There is no time gap between them. The problem is that the method of measurements has changed between both time series and they are no statistically the same. I would like to merge them to receive one homogeneous

One thing that has given me trouble is the fact that the time series implementation in the class ts relies on the concept of a "start" to the time series. For example, if I have ts1 <- ts(c(1,2,3)) dts1 <- diff(ts1) then dts1 will be a vector c(1,1) but with the attribute start = 2. Similarly, when one takes lags the "start" is moved around and the underlying vector

Hi there, I'm rying to 'turn' an Schunk in an .Rnw file(Xemacs-21.4.13, ESS-5.2.8, R-2.1, miktex-2.4.1705). Has anyone got the isorot package to work with Sweave? JC example test.Rnw: \documentclass[a4paper]{article} \usepackage{Sweave} \usepackage{isorot} \rotdriver{dvips} \clockwise \title{Sweave Example 1} \author{apologies to Friedrich Leisch } \begin{document} \maketitle

I have a data file with a given time series of price data and I would like to split the time series into a test set and training set. I would then like to build an ARIMA model on the training set and apply this model on test set. Below is some code: [CODE] data= read.table("A.txt",sep=",") attach(data) training = data[1:120, 6] test = data[121:245, 6] ts1 = ts(training) ts2 =

HI: I need ideas on how to make this code shorter (maybe with a second loop?). The code as it is works, but in this case I only have 14 samples, but it will become insane with more, so I need a way to make it more automatic. The problem is that the output from ts1, ts2, and so on is a vector with more than one value, so I do not know how to solve this. Thanks Prenewbie The code is the

Hello R Gurus: I would like to take a character string and split at the $ sign. I thought that strsplit would do it, but here are the results: > vv [1] "whine$ts1" > vv [1] "whine$ts1" > strsplit(vv,"$") [[1]] [1] "whine$ts1" Does anyone have any suggestions, please? Thanks, Edna Bell

How can I get a single attribute value of an object ? I jhave the tiemSeries object > ts1 Open 2003-10-09 02:00:00 1.27 2003-10-10 02:00:00 1.25 2003-10-13 02:00:00 1.27 2003-10-14 02:00:00 1.29 When I unclass ts1 I get: > unclass(ts1) list() attr(,"Data") Open 2003-10-09 02:00:00 1.27 2003-10-10 02:00:00 1.25 2003-10-13 02:00:00 1.27

Dear list, I am trying to produce a pdf manual using Rtools (version 3.1.0.1936) on Windows 8 (64-bit). I am also aware of the inconsolata issue posted in the forum (https://stat.ethz.ch/pipermail/r-devel/2013-June/066850.html), and then I have installed the patched version of R-3.0.1. However I still have some errors with the following message: R CMD Rd2pdf mypack Hmm ... looks like a package

To the Lattice expert - I am new to the lattice package and I would like to add a curve to the xyplot(). I know that I need to use panel function to add it. But it doesn't work. Is there anyone can help me out on how to transfer the data into the function? The following is my code and I would like to add lines for xx,yy xx <- seq(0,1,length = 100); yy <-

Hi, I'm having problems trying to join my XPpro boxes to a Samba (2.2.7a) domain. I was able to get 1 XP box to join (the very first one) but the rest complain about a duplicate name on the network just after I enter the root/pw when joining the domain. None of the smb user names are in common with the machine names. The NT4 workstations have no problem joining the samba domain. All

Hi. I'm trying to delete a machine account using pdbedit -x -m and it fails. Any idea why that might happen? The machine account exists.. # pdbedit -Lv 'ts$' Unix username: TS$ NT username: Account Flags: [W ] User SID: S-1-5-21-2594621239-3917873753-3682423317-3286 Primary Group SID: S-1-5-21-2594621239-3917873753-3682423317-515 Full Name: Home Directory:

Dear list, We have encountered problems with the package mvbutils while creating (or reading from disk) timeSeries objects. Below is a very simple example of the problem. We are using R 2.12 in  a PC running Windows 7  64-bit Professional. ----------- library(timeSeries)   cv1 <- as.character(as.Date(1:10, origin = "1970-01-01")) ts1 = timeSeries(1:10, cv1, units = NULL)  

Hi all, please cc me, i'm not on the list. Second: All google findable information about problems setting up ntlm_auth for squid with winbind are read and checked more than three times. After breaking a running setup under debian squeeze, i go back to debian lenny to circumvent the actual MIT kerberos problem[1]. [1] http://bugs.debian.org/cgi-bin/bugreport.cgi?bug=566977#57 Now i

Hi, I have a problem with the appearance of legend in my plots. If I specify the legend positions in characters like "topright"..etc, it appears, if i specify it in terms of coordinates like "-1, 1" .. etc, it does not appear. Can anyone help me with this? script - x.date <- as.Date(paste(year, month, day, sep="-")) ts1.n.e3 <- ts(data.nemr.e3[,3]) z1.n.e3

Hi, All, This is my program ts1.sim <- arima.sim(list(order = c(1,1,0), ar = c(0.7)), n = 200) ts2.sim <- arima.sim(list(order = c(1,1,0), ar = c(0.5)), n = 200) tdata<-ts(c(ts1.sim[-1],ts2.sim[-1])) tre<-c(rep(0,200),rep(1,200)) gender<-rbinom(400,1,.5) x<-matrix(0,2,400) x[1,]<-tre x[2,]<-gender fit <- arima(tdata, c(1, 1, 0), method = "CSS",xreg=t(x))

I am interested in creating multiple (say 1000) time series, from a given stochastic process, of length 250. I want to refer to each realization with its own variable name, of the format say, tsn, where n is the n'th simulation. i.e. ts1, ts2, ts3, ts4, .... , ts1000 The way I am thinking of doing this is placing the following code within another loop, and the 'tsn' assignment should

i call to my tollfree number buy my CLI send the next error: Aug 14 08:15:22 WARNING[25931]: format_wav.c:169 check_header: Unexpected freqency 22050 Aug 14 08:15:22 WARNING[25931]: file.c:441 ast_filehelper: Unable to open file on /var/lib/asterisk/sounds/procall3.wav Aug 14 08:15:22 WARNING[25931]: file.c:828 ast_streamfile: Unable to open procall3 (format ulaw): No such file or directory Aug

Fix internal ppmtolss comment and man page describing ppmtolss syntax with color index Fix typo s/freqency/frequency/ in comment Not so sure it will not be mangled, so it's also as attachment. Signed-off-by: Gilles Espinasse <g.esp at free.fr> diff --git a/man/ppmtolss16.1 b/man/ppmtolss16.1 index 954e94b..5be3263 100644 --- a/man/ppmtolss16.1 +++ b/man/ppmtolss16.1 @@ -3,7 +3,7 @@