similar to: Fwd: smart updates and rolling windows

Hi, Assuming I have a time series on which I will perform rolling-window MLE. In other words, if I stand at time t, I'm using points t-L+1 to t for my MLE estimate of parameters at time t (here L is my rolling window width). Next, at t+1, I'll do the same. My question is that is there anyway to avoid performing MLE each time like does the above. My impression is that rolling from point t

Hello R'ers, I have a large set of data which has many y samples for each unit x. The data might look like: Seconds Response_time ---------- ---------------- 0 0.150 0 0.202 0 0.065 1 0.110 1 0.280 2 0.230 2 0.156 3 0.070 3 0.185 3 0.255 3 0.311 3 0.120 4 .... and so on When I do a basic plot with type=l or the default of points it obviously plots every

Hello R'ers In previous versions of R (I now use 2.0.1) there was a package "bootstrap". I wrote some programs that depended heavily on this package but can unfortunately not find it on Cran today. Is there any way to find an older version of "bootstrap" and use it with the new version - 2.0.1? Sincerely Fredrik Lundgren Norrk??ping Sweden

Hello fellow R'ers, I have a simple calculation with a very large data set being generated (34.9 million values) on a somewhat unreliable XP box that will likely take ~ 74hrs. I wanted to know if there is a way to have my script automatically "save.image()" throughout the calculation in case of a crash. This could be on the basis of output generated or time elapsed. I checked the

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

On Mon, Mar 25, 2013 at 4:02 PM, Jakob Stoklund Olesen <stoklund at 2pi.dk>wrote: > > On Mar 25, 2013, at 2:51 PM, Akira Hatanaka <ahatanak at gmail.com> wrote: > > > Yes, it sounds like it will solve the problem. > > > > Using the following example where live ranges of accumulators $vreg_acc0 > and $vreg_acc1 conflict, > > > > MULT

> Quoting Stephane Lesage <stephane.lesage at ateis-international.com>: > > Is this a bug ? Is it possible to fix it ? > > (I use version speex 1.2beta2, because newer versions just > don't work > > on my > > platform) > > This is probable the cause. 1.2beta2 was the first release to > include the resampler and it had many bugs. I suggest trying

Dear R'ers, I use R version 2.10.1 (2009-12-14) Copyright (C) 2009 The R Foundation for Statistical Computing ISBN 3-900051-07-0 and Mac OS X 10.5.8 and use to work in different dierectories. When doing so it appears as if the .Rhistory file from my home directory always is loaded even when I want the file from the directory where I'm working. I have tried to fix it through preferences

On Mar 26, 2013, at 10:29 AM, Akira Hatanaka <ahatanak at gmail.com> wrote: > The size of general purpose integer registers for mips32 is 32-bit and accumulators are 64-bit registers consisting of 32-bit hi/lo register pairs. So you will need two instructions to copy two 32-bit GPR registers to a 64-bit accumulator register. If spilling to multiple registers is unsupported, perhaps I can

I just found out that this message got scrambled with other threads, so I trying to re-send... Dear R'ers, I am trying to build a composite plot (with several plots in one figure). I have tried, but I cannot use facetting, as I need to customize each plot using grid. Since all the plots are the same (with different data, but same layout and categories), I would like to have only one legend,

Dear R'ers, I would like to use barplot or a similar function to plot data demonstrating the distribution of the length of a kind of conservation in about 25000 DNA sequences. My data look like this: #Total sequences: 23873 0 19936 1 218 2 391 3 477 4 360 5 431 6 294 7 215 8 320 9 209 10 160 (.....) 99 0 100 1 101 0 102 0 103 1 104 0 105 0 106 0 107 0 108 0 109 1 Therefore, I would like

Dear R'ers, I've scanned available documentation and the web, but I can't seem to figure out where I've gone wrong in adding numbers to scale the colorkey in levelplot (Lattice package). For example, levelplot(z, contour=T, labels=T, cuts=10, region=T,... colorkey=list(T, space="bottom", at=seq(0,2000,length=20),

Hi R'ers- I'll believe this can be done and that I'm just not clever enough. I'm trying to do this in 2D, but i think we have the same problem in 1D. #Assume you have some 1D time series d<-rnorm(100) #Say you want to get the average over intervals of different lengths #like from time 10 to time 24 and time 51 to time 86 starts <- c(10,51) ends <- c(24,86) results

Dear R'ers, I am wondering what is the smallest geographicterritorial unit available for formatting in Canada. Provinces? I know that in the US it is the county so that I can color US counties any way I want, for example: ### Example for coloring US counties ### Creating an ARTIFICIAL criterion for coloring US counties: library(maps) allcounties<-data.frame(county=map('county',

Dear R'ers, I was not able to change the font size of axis label using cex.   Tried R archive solutions, though not solved yet. My code is: print(dotplot(ss ~ Response | bb*sr*tr, pp, pch = 21,               strip = FALSE, strip.left = TRUE,layout = c(3,8),               scales = list(y = list(relation = "free")), cex.axis=.2, cex.label=.2,                       ylab =

Dear R'ers, I have trouble installing tikzDevice in Ubuntu.  When I use install.packages("tikzDevice"), it gives error message: ERROR: dependency ‘filehash’ is not available for package ‘tikzDevice’ * removing ‘/usr/local/lib/R/site-library/tikzDevice’ Then I tried filehash installation, I get the message: "package ‘filehash’ is not available (for R version 2.13.1)"

Dear R'ers, I have a very large data frame (over 4000 rows and 2,500 columns). My task is very simple - I have to replace all NAs with a zero. My code works fine on smaller data frames - but I have to deal with a huge one and there are many NAs in each column. R runs out of memory on me ("Reached total allocation of 1535Mb: see help(memory.size)"). Is there any other, more efficient

On Mar 25, 2013, at 2:51 PM, Akira Hatanaka <ahatanak at gmail.com> wrote: > Yes, it sounds like it will solve the problem. > > Using the following example where live ranges of accumulators $vreg_acc0 and $vreg_acc1 conflict, > > MULT $vreg_acc0, $vreg_gpr0, $vreg_gpr1 > MULT $vreg_acc1, $vreg_gpr2, $vreg_gpr3 > > (consumer of $vreg_acc1) > (consumer of

It won't cause a negative effect, go for it! Dynamic_cast is realllly slow compared to dyn_cast, it is worth the memory. -Chris On Oct 3, 2012, at 9:39 PM, Sean Silva <silvas at purdue.edu> wrote: > Thanks for the feedback! > >>> Does anybody have anything else they think should go into TGContext or >>> any other responsibilities it should have? Or any

R'ers: I was hoping I could get some direction on this. I have a dataset of the form: Y1,Y2,...,YM = f(X1,X2,...,XN), where N is >>> M The response data (Y1,Y2,...,YM) is frequency data, such that the sum of all Yi = 1.0. Both Xj and Yi are continuous variables. I'm trying to figure out the best approach(es) to solving for the model f() -- any ideas? I could solve