similar to: question on creating zoo object from indexed time series: re-arranges dates

Hello, My R skills are somewhere between novice and intermediary, and I am hoping that some of you very helpful forum members, whom I've seen work your magic on other peoples' problems/questions, can help me here. I have a matrix with the following format: (i) individual plants comprising many different genotype groups (i.e., a plant is genotype 1 or genotype 2 or genotype 3, etc). The

day<-structure(c(7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7, 7.69, 7.68, 7.68, 7.67, 7.67, 7.67, 7.66, 7.65, 7.65, 7.65, 7.64, 7.64, 7.63, 7.63, 7.63, 7.62, 7.62, 7.62, 7.62, 7.63, 7.63, 7.63, 7.63, 7.63, 7.64, 7.64, 7.65, 7.65, 7.65, 7.66, 7.66, 7.67, 7.67, 7.67, 7.68, 7.68, 7.69, 7.69, 7.69, 7.69, 7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7,

NOTE: I will provide data if necessary, but I didn't want clutter everyones mailbox All: I have a time series with level and temperature data for 11 sites for each of three bases. I will have to do this more than once is what I am saying here. OK, The time series are zoo objects with index values in chron format. The problem is that the date and times should be at even 15 min intervals,

Dear list, I have to look examine hourly time - series and would like to plot variable section of them using plot.zoo. Hourly time series data which looks like this: YYYY MM DD HH P-uk P-kor P-SME EPOT EREA RO R1 R2 RGES S-SNO SI SSM SUZ SLZ 2003 1 1 1 0.385 0.456 0.021 0.000 0.000 0.000 0.013 0.223 0.235 0.01 0.38

New user here. My goal is pull daily averages from a long dataset. I've been working with some code I got from this list from https://stat.ethz.ch/pipermail/r-help/2009-March/191302.html The code how I have been using it is as follows: library(zoo) library(chron) DB<-read.table("/Users/me/Desktop/R/data.csv", sep=",", header=TRUE, as.is =TRUE) z<-zoo(LTER6$temp,

Is there a easy way to create the time index for a zoo/xts object for every 100 milliseconds. eg. time Index would be: 10:00:00:100 10:00:00:200 10:00:00:300 10:00:00:400 I am looking to build an empty zoo/xts object with time index from 10am to 3pm, index jumps by 100ms each row. Thanks, Chris -- View this message in context:

Goal: get time series data interpolated on to desired time stamps. I have two or more data sets that have time stamps that vary from 5 mins to 3-5 hours. I want to get all the data put on common time stamps e.g. "00:05:00" intervals. I asked Gabor and got some very good code ( zoo aggregate, na.spline, na.approx) but I'm having trouble getting the csv file read in and converted to a

Hi guys, I have a list of 250 numbers as a result of using the ?by function! List of 246  $ 0   : num [1:28] 22 11 31...  $ 1   : num [1:15] 12 14 9 ... .. .. ..  - attr(*, "dim")= int 250  - attr(*, "dimnames")=List of 1 The problem is that each list of 250 has different length! I would like to get the values of each list in a vector like vector_0 = (22,11,31,..), is this

I'm trying to export a large data frame to a text file for permanent storage. The only thing I could find was the treeglia Package but that didn't work. Any suggestions? Thanks! Christina D Smith PhD Student, GRA Statistics Department Kansas State University

Hi, I am running R 2.0.1.1. on Windows. It is a Dell Dimension with a 3.2 Ghz Processor and 4Gb RAM. When using the ReadAffy() function to read in 97 arrays, I get the below error messages: Error: cannot allocate vector of size 393529 Reached total allocation of 1024Mb: see help(memory.size) When I use the comman "memory.limit(size=4000)" to increase the memory size to the

Full_Name: Christina Merz Version: R version 2.5.0 (2007-04-23) OS: mingw32 Submission from: (NULL) (213.70.209.132) R> version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

Hi everyone I'm using Tukey HSD as post-hoc test following a lme analysis. I'm measuring hemicelluloses in different species treated with three different CO2 concentrations (l=low, m=medium, h=high). The whole experiment is a split-plot design and the Tukey-function from the package multcomp is suitable for lme-analysis with random factors. The analysis works fine but I get a non

Hi, I'd like to make a time series at an annual frequency. > a<-xts(x=c(2,4,5), order.by=c("1991","1992","1993")) Error in xts(x = c(2, 4, 5), order.by = c("1991", "1992", "1993")) : order.by requires an appropriate time-based object > a<-xts(x=c(2,4,5), order.by=1991:1993) Error in xts(x = c(2, 4, 5), order.by =

Dear list members, How can I shade the background area of a zoo time series object between specific dates? eg. library(tseries) library(zoo) SP500<-get.hist.quote("^GSPC", start = "1990-01-01", quote = "Close") plot(SP500) How can I produce the same plot but with a (say) red background between 2007-04-12 and 2008-05-14 ?

Hi! I am looking for a tool (application or webinterface) which shows me the current status of an Asterisk server, e.g.: - Status of the SIP peers (registered/offline) - current incoming and outgoing calls - start-time, numbers, some history - history (calls stopped in the last 15 minutes, who hang up?) - should be possible to link those calls to the relevant SIP peers -

Hi, I just started using speex recently and I've been having problem with Speex ver. 1.1.10 decoder. I used Speex to encode a wav file to spx then decode it back to wav file to see the voice quality that results from different speex options. Somehow, I couldn't get the original wav file back. I used both Speex ver 1.1.6 and ver 1.1.10 to both encode and decode, despite the same output

Hi guys, I still did not solve my problem properly! I have to compare the values of two lists of 250 numbers as a result of using the ?by function! List1 of 250  $ 0   : num [1:28] 22 11 31...  $ 1   : num [1:15] 12 14 9 ... .. .. ..  - attr(*, "dim")= int 250  - attr(*, "dimnames")=List of 1 List2 of 250  $ 0   : num [1:24] 20 12 22...  $ 1   : num [1:17] 11 12 19 ... .. ..

OK, that was a stupid error in the le_int prototype. It's fixed in svn now. Jean-Marc Le mercredi 28 septembre 2005 ? 11:53 +1000, Jean-Marc Valin a ?crit : > OK, it seems like I screwed up the wav header code between 1.1.8 and > 1.1.9. Everything seems fine (playing to soundcard also works), except > that the length reported for the wav is too short, which is why it stops >

Hi, I'm installing an Asterisk box with a TDM2400P in Brazil. I can make analog phones work while lines are not working. Since I do not know anything about brazilian lines, is there anybody who can tell me what is wrong/missing in my conf files (below)? TIA Giorgio _zaptel.conf:_ fxoks=9-16 fxsks=17-24 defaultzone=br loadzone=br* * _zapata.conf:_ context = inbound_zap echocancel = 128

Salve, sar? assente dal 18 Agosto al 9 Settembre compreso. Per importanti escalation utilizzi la coda Clarify EUA_CC NTWKST. Se desidera parlare col mio manager contatti Alexander Weigl Cordialmente, Sandro ------------------------------------------------------------------------ Dear sender, I'm absent from August 18th till September 9th included. During this time I wont be able to