similar to: concatenating two successive time series

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

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

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

Hi again. Perhaps a simple question this time.... I am analysing data with a dependent variable of insect counts, a fixed effect of site and two random effects, day, which is the same set of 10 days for each site, and then transect, which is nested within site (5 each). I am trying to fit the cross classified model using GLMM in lme4. I have, for potential use, created a second coding

Muy buenas, necesito crear un vector a partir de los valores de una tabla como la de abajo. Debe ser algo muy fácil pero no lo encuentro en la web. Gracias, Manuel 246, 345, 401, 131,125, 69 a partir de: TS1 TS2 TS3 TS4 TS5 TS6 246 345 401 131 125 69 [[alternative HTML version deleted]]

Sorry to keep posting but I want to do this right and I'm hoping for some pointers I now have two time series objects which I need to subtract. Unfortunatly the two series dont have the same sample rates. When I try to subtract them avgSub<-avg1-avg2 The time series object is clever enough to object. So I guess I need to write a function for subtraction of the time series objects which

Hi, I have got four correlation matrix. They are the same set of variables under different conditions. Is there a way to test whether the correlation matrix are significently different among each other? Could anyone give me some advice? -- View this message in context: http://www.nabble.com/correlation-matrix-difference-tf4073868.html#a11578046 Sent from the R help mailing list archive at

I would like to add a rule allowing only the address 192.168.150.20 and the range of addresses from 192.169.150.100 to 192.168.150.150 in zone dmz0 to connect to two terminal servers in the local zone. Is there a syntax that can specify a range of addresses in the rules file? Do I have to enter each one separately? -- Stephen Carville Unix and Network Adminstrator DPSI 6033 W.Century Blvd.

I am trying to extract only the winters (defined to be 01-Dec through 28-Feb) of daily data from 1948-2002. There are 90 days in each winter season. I wrote the following code to gather the winter dates into a single vector: DJF <- NULL for(year in 1949:1999) { temp.begin <- strptime(paste("01/12", year-1, sep="/"), "%d/%m/%Y") temp.end <-

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)  

Hallo everyone I have two univariate time series (class ts) describing the same variable. They have the same resolution, but span different periods in time with a big gap in between. I need to append one to the other such that they are one object, with the gap filled with NAs. The method ts.union produces a multivariate time series where the time axis is correct, but the individual time series

Hello, I would like to apply the EF function (a goodness of fit test between predicted (ts2) and observed (Tw) values in the qualV package) to each of the Sites in my data frame. However, when I try the following script, R gives an error message that lists unused arguments, which include the Tw and ts2 that are explicitly in the function: nsc = by(tsdem16, tsdem16['Site'],

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

Hi people, does someone know how to only log the name of the machine and not the name of the machine AND the ip? Let me explain this: -rw-r--r-- 1 root root 0 Apr 28 02:47 0.0.0.0.log -rw-r--r-- 1 root root 0 Apr 28 02:49 10.3.0.10.log -rw-r--r-- 1 root root 0 Apr 28 02:47 10.3.0.21.log -rw-r--r-- 1 root root 0 Apr 28 02:46 10.3.0.22.log -rw-r--r-- 1 root root 0 Apr 28 02:48

Hi. I am sure there is a better way in R to do this then using a loop but I am new to it and not sure what to do. I think it might be something about using a function as an argument but not sure. I have a 1 x 2000 vector TS2 which has entries from the set {x: x is in Z and 0<x<8} (where Z is the set of Integers). Then I also have a 5050 x 7 matrix called 'perm' whose entries are

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 =

I have been trying to use the ets function in the forecast package on a daily time series (ts2 is a ts object with frequency =7). However when I run the following code I get an error related to etsmodel. I have looked at ets and I can see that there is a call to the function etsmodel, but I cant seem to find info on the ets function anywhere. Does anyone know anything about the etsmodel function?

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