similar to: CAPdiscrim error in BiodiversityR

Hello, Using the technique suggested here ( https://stat.ethz.ch/pipermail/r-devel/2007-September/047028.html), we can update a matrix in place (see example 1 below), with huge time and memory savings. Is it possible to tweak example 2 to update in place also? The only difference is that the matrix is held within a new environment. Thanks, Simon # example 1 - fast and doesn't

Hello, I am running an ARMA model to run forecast for changes in S&P 500 prices. My ARMA calculations look as follows armacal <- arma( spdata, order = c(0,4), lag = list(ma = c(1,2,4)) ) Output: Call: arma(x = spdata, order = c(0, 4), lag = list(ma = c(1, 2, 4)) ) Coefficient(s): ma1 ma2 ma4 intercept -0.073868 0.058020 -0.081292 0.007082 All's

Dear R users, I'm estimating a mixed effects model in which the spatial correlation is controlled for by the "corGaus" structure. I'm wondering if there is a document or paper that explains how the spatial correlation structure (such as "corExp" or "corGaus") works. Let me use the example and data posted on UCLA's R FAQ webpage to explain my problems.

I need to perform a QR-decomposition of a sparse matrix, so I've been trying to use the Matrix package. Unfortunately I don't seem to be getting exactly the same results as if I had used the qr() command from the base package. Here is an example of what I'm doing. > spdata <-rpois(50,1) > y <- rnorm(10,0,1) > S <-

Hello. I want to test a multi-site setup w/ replication to a central DC only. In a test scenario I have set up a central samba4 DC (sernet-samba 4.1.12-9). For provisioning I have used the --site=central, so no Default-First-Site-Name has been created on provisioning, just 'central'. In AD Site & Services I have unchecked the "Bridge all site links" in "Inter-Site

Hello R list, I have a list of data frames, in the form (dump output for an example list follows message): $site1 year data 1 2000 0.03685042 2 2001 0.02583885 6 2005 0.02480015 7 2006 0.03458745 $site2 year data 1 2002 4.071134e-03 2 2003 -4.513524e-08 3 2004 8.336272e-08 4 2005 2.270803e-08 5 2006 1.263597e-02 6 2007 NaN I would like to plot each

spData <- as.vector(getYahooData("SPY", start="19900101",end="20081231")[,"Close"]) > cr <- cReturns(spData, ttr = "sma", params=c(20)) Error in ind[t - k] <- pos[t - k + 1] - pos[t - k] : replacement has length zero I am getting the above error when running the cReturns function. Any ideas on this? How does one drill down into

Cada vez que intento cargar el paquete "agricolae" me da el siguiente mensaje: Error: package or namespace load failed for 'agricolae' in loadNamespace(i, c(lib.loc, .libPaths()), versionCheck = vI[[i]]): there is no package called 'spData' He intentado todo actualizar los paquetes, volverlos a instalar pero no hay forma. Si sabeis la solución ?? Un saludo. Juan