similar to: Help with vectors and rollapply

Hi everyone, I have a rollapply statement that applies a function, in steps, over a data matrix as follows: #Code start testm<-rollapply(mat, 100, by=100, min, na.rm=F) #Code end This moves down matrix 'mat' and calculates the minimum value over a 100 row range, every 100 rows (i.e. no overlaps). NAs are not removed. I want to modify this statement somehow so that the rollapply()

zoo's rollapply() function appears to be extremely useful for plugging in a function on-the-fly to run over a window. With inline, there is a lot more coding and room for error, and the code is less portable because the user has to have R compiling set up or it won't work. However, rollapply() seems to be really slow. Several orders of magnitude slower than inline, in fact. I don't

Dear Colleagues, I have faced with the problem that function rollaply with rolling window for calculation of volatility doesn't give the all results of calculations. I have run the rolling window for calculation in Excel and obtained that the number of outputs for Excel is 36 and for R is 18. The total number of observations is 37. In the attachment you can find pdf of the Excel and Excel

Hi, Sorry ahead of time for not including data with this question. Using rollapply to calculate mean values for 5 day blocks, I'd use this: Roll5mean <- rollapply(data, 5, mean, by=5, align = c("left")) My question is, can someone tell me how to fill in the days between each of these means with the previously calculated mean? If this doesn't make sense, I will clarify and

Hi, I faced with problem when start using function - rollapply(returns, 3 , mean) Error in UseMethod("rollapply") : No suitable Method for "rollaply" How can I fix the problem? Thank you for help. -- Best regards, Andy Fetsun [[alternative HTML version deleted]]

Hi all- Thank you for reading my post. Please bear in mind that I'm very much a newbie with R! My question is this: I'm trying to use rollapply() on an irregular time series so I can't simply use the width parameter (I don't think). Rather than last 5 entries, I'd like to rollapply on last 6 months (for example). What would be the proper course of action for this? Thanks!

Hi I would like to understand how to extend the function (FUN) I am using in rollapply below. ###################################### With the following simplified data, test1 yields parameters for a rolling regression data = data.frame(Xvar=c(70.67,70.54,69.87,69.51,70.69,72.66,72.65,73.36), Yvar =c(78.01,77.07,77.35,76.72,77.49,78.70,77.78,79.58)) data.z = zoo(d) test1 =

Hi. I'm comparing output from rollapply.zoo, as produced by two versions of R and package zoo. I'm illustrating with an example from a R-help posting 'Zoo - bug ???' dated 2010-07-13. My question is not about the first version, or the questions raised in that posting, because the behaviour is as documented. I'm puzzled as to why na.rm no longer is passed to mean, i.e. why

# Hi there, # I am trying to apply a function over a moving-window for a large number of multivariate time-series that are grouped in a nested set of factors. I have spent a few days searching for solutions with no luck, so any suggestions are much appreciated. # The data I have are for the abundance dynamics of multiple species observed in multiple fixed plots at multiple sites. (I total I

Dear R users. I have zoo object "size_june" containing market-capital values: > dim(size_june) # market-cap data of 625 firms for 20 years [1] 20 625 > class(size_june) [1] "zoo" > size_june # colnames = "size.firmcode" size.34020 size.4710 size.11050 size.10660 size.9540 size.8060 size.16160 size.8080 size.9280 1988-06-30 NA

Here is a relatively simple script (with comments as to the logic interspersed): # Some of these libraries are probably not needed here, but leaving them in place harms nothing: library(tseries) library(xts) library(quantmod) library(fGarch) library(fTrading) library(ggplot2) # Set the working directory, where the data file is located, and read the raw data

Hi Howard I was reading your post on Syslinux of the PXE: EB No PXE err Did you ever find a resolution to what this error was? I'm looking, but not having much luck. I used http://www.rom-o-matic.net/ to make the boot disk. Cheers in advance ________________________________ David Pexton | Systems Administrator Complinet Limited, Vintners Place, 68 Upper Thames Street London EC4V

Dear List members I want to do the sliding window analysis of some specific values. Here is my code: require(zoo) dat <- read.table("chr1.txt", header = TRUE, sep="\t") dat2 <- cbind(dat[1,3]) #The first column is also important. It represents the position of the site on the chromosome. TS <- zoo(c(dat2)) a <- rollapply(TS, width=1000000, by=200000, FUN=mean,

#using "rollapply" to calculate a moving sum or running sum? #I am tryign to use rollapply to calcualte a moving sum? #I tried rollapply and get the error message #"Error in seq.default(start.at, NROW(data), by = by) : # wrong sign in 'by' argument" #example: mymatrix <- ( matrix(data=1:100, nrow=5, ncol=20) ) mymatrix_cumsum <- ( matrix(data=NA, nrow=5,

Hello R community members, I'm trying to install the 'rqpd' package which is developed by Roger Koenker and Stefan Bache. When I try to install the package using the command 'install.packages("rqpd",repos="http://R-Forge.R-project.org")' I'm getting the following two messages: i) package ?rqpd? is available as a source package but not as a binary

Dear Colleagues, I would like to provide a CITATION file for my package nat.nblast [1]. I have the correct citation in BibTeX format [2]. How can I convert this BibTeX to the format needed by R for a package CITATION file (I have a lot of other packages needing citations ...). I think what I need is the opposite of RefManageR::toBiblatex [3]. This seems like it should be a common need, so I

Hello, Does anyone know of an R package that can read Windows BMP format images? Many thanks, Greg Jefferis. -- Gregory Jefferis, PhD Division of Neurobiology MRC Laboratory of Molecular Biology, Hills Road, Cambridge, CB2 0QH, UK. http://www2.mrc-lmb.cam.ac.uk/group-leaders/h-to-m/g-jefferis

This seems to be a simple problem, and I feel that there ought to be a simple answer, but I can't seem to find it. I have a function that returns a number of values as a heterogeneous list - always the same length and same names(), but a number of different data types, including character. I want to apply it to many inputs, resulting in a list of lists. I would like to turn this list of

If I have a list I can set a sub-element as follows on the command line: people=list() people$tom$hair="brown" people But what if I have a string containing the name of the sub-element that I want to access? subel= "people$tom$hair" get(subel) # returns error assign(subel,"red") # silent but doesn't change list people The attempts above using assign/get

Hi, I am currently downloading CentOS-4.4-i386-bin1of4.iso and was wondering what is the minimum number of install disks I require to get a working install of Centos ? Do I need to download all four of the disks ? cheers, Jeffery -------------- next part -------------- A non-text attachment was scrubbed... Name: jeffery.vcf Type: text/x-vcard Size: 293 bytes Desc: not available URL: