similar to: Problem wiht mvbutils and timeDate in R 2.12

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

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

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))

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 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

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

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 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 =

Hello, Sorry if anyone gets this message twice, as my mailserver may not be working. Thanks for your response. Your idea makes a lot of sense to me, but I've been unable to get seconds to work. I ended up with this format finally: "2007-10-31_16:20:22" Problem is I am unable to get it recognized as a date using timeDate(): R>

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]]

For a model I am working on, I have samples organized by year and week of the year. For this model, the data (year and week) comes from the basic sample data, but I require a value representing the amount of time since the sample was taken (actually, for the purpose of the model, it is sufficient to use the number of weeks from the middle of the sample week to the present). What I have found so

Version 2.5.1 of 'mvbutils' is now on CRAN. This version offers improved support for easy package preparation and maintenance, plus minor changes required by the new version of the 'debug' package (see separate announcement). Package 'mvbutils' offers the following main features, as well as a number of miscellaneous goodies: - Hierarchical organization of projects (tasks)

Version 2.5.1 of 'mvbutils' is now on CRAN. This version offers improved support for easy package preparation and maintenance, plus minor changes required by the new version of the 'debug' package (see separate announcement). Package 'mvbutils' offers the following main features, as well as a number of miscellaneous goodies: - Hierarchical organization of projects (tasks)

Dear R users New versions of the 'mvbutils' and 'debug' packages are now available on CRAN, both in source form and as precompiled binaries. 'mvbutils' offers the following (as well as many miscellaneous utilities): ? hiearchical, searchable project organization, with workspaces switchable inside a single R session, and objects in "ancestor" projects always

Dear R users New versions of the 'mvbutils' and 'debug' packages are now available on CRAN, both in source form and as precompiled binaries. 'mvbutils' offers the following (as well as many miscellaneous utilities): ? hiearchical, searchable project organization, with workspaces switchable inside a single R session, and objects in "ancestor" projects always

There is a new version of the 'mvbutils' package (v1.1.1) available on CRAN. For existing users, the main new features are: (i) 'help' now works properly with with R2.0+ (!) (ii) 'mvbutils' is now NAMESPACEd, so you can avoid naming conflicts with other packages. (Many functions are no longer user-visible. (iii) Lazy-loading of individual objects is now available. Among

There is a new version of the 'mvbutils' package (v1.1.1) available on CRAN. For existing users, the main new features are: (i) 'help' now works properly with with R2.0+ (!) (ii) 'mvbutils' is now NAMESPACEd, so you can avoid naming conflicts with other packages. (Many functions are no longer user-visible. (iii) Lazy-loading of individual objects is now available. Among

I have a daily time series and have two questions to get some help with. Firs,t I have dates in simple numeric values. e.g. ymd [1] 20050104 20050105 20050106 20050107 20050110 20050111 20050113 20050114 [9] 20050118 20050120 20050121 20050124 20050125 20050126 20050127 20050128 [17] 20050201 20050202 20050203 20050204 Now, I'd like to compute statistics, e.g. acf, by business days. So, I