similar to: R-beta: ar

I have just put up a new executable as a replacement for the one in rseptbeta.zip there have only been a few changes; mostly to the menu's. I am about to start on a major overhaul including getting survival to work and grabbing the 0.60 version once it's stable. Please let me know about other enhancements you want.... robert

I have just put up a new executable as a replacement for the one in rseptbeta.zip there have only been a few changes; mostly to the menu's. I am about to start on a major overhaul including getting survival to work and grabbing the 0.60 version once it's stable. Please let me know about other enhancements you want.... robert

I am trying to understand the function acf stats:::acf shows me the function I am having trouble understanding the usage "$acf" in the following acf <- array(.C(R_acf, as.double(x), as.integer(sampleT), as.integer(nser), as.integer(lag.max), as.integer(type == "correlation"), acf = double((lag.max + 1L) * nser * nser), NAOK =

--- On Mon, 1/17/11, Raymond Wong <raywong365@yahoo.ca> wrote: From: Raymond Wong <raywong365@yahoo.ca> Subject: Re: [R] help in calculating ar on ranked vector To: "Uwe Ligges" <ligges@statistik.tu-dortmund.de> Received: Monday, January 17, 2011, 11:56 AM Thanks Uwe:   Here is my code. the first set of print statements work, but not the second.   #

Hello! I have two questions. First of all, I have a problem dealing with acf (Autocovariance function) and need help. First I defined a time series, x, which is a vector created by x <- ts(rnorm(200)). So I plugged the series directly into the acf function, acf(x) and an error message popped up as: Error in .C("acf", as.double(x), as.integer(sampleT), as.integer(nser), :

Dear list, according to the documentation of acf{ts} "the partial correlation coefficient is estimated by fitting autoregressive models of successively higher orders up to lag.max. " However, R seems to return the Yule-Walker estimates of the PACF by default. You can check this using c(1:10) as the series: the YW estimates are 0.7000000 and -0.1527035 for lags 1 and 2 . If the PACF

I'm using R-2.9.1, so forgive me if this has already been resolved. One element of the object returned from the acf function is n.used, described in the man page as "The number of observations in the time series". However, I've noticed that this value is set to nrow(x) via the sampleT variable, i.e. the number of rows in the passed-in series. This forces an assumption that all

Hi there, we are beginners in R and we are trying to fit the following time series using ar(2): > x <- c(1.89, 2.46, 3.23, 3.95, 4.56, 5.07, 5.62, 6.16, 6.26, 6.56, 6.98, > 7.36, 7.53, 7.84, 8.09) The reason of choosing the present time series is that the we have previously calculated analitically the autoregressive coefficients using the direct inversion method as 1.1, 0.765, 0.1173.

Please pardon me if you see this message twice. The mail server has a bit problem. ***************************************************** Hello! I have two questions. First of all, I have a problem dealing with acf (Autocovariance function) and need help. First I defined a time series, x, which is a vector created by x <- ts(rnorm(200)). So I plugged the series directly into the acf

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

There is now a preliminary version of a time series package in the R-devel snapshots, and we would welcome feedback on it. It is based in part on the packages bats (Martyn Plummer) and tseries (Adrian Trapletti) and in part on code I had or have written. (Thanks for the contributions, Martyn and Adrian!) Some of the existing ts code has been changed, for example to plot multiple time series, so

Dear helpers, When I use ar to fit the data with length 180, I have the following error: ar(x,method="burg") Error in acf(x, type = "covariance",lag.max=order,plot=FALSE): lag.max must be at least 1 If I use ar(x), then I have Call (x=x) order selected 0 sigma^2 estimated as 5374 Obviously I missed some points for using ar. This is R 1.7.0 under Redhat Linux

I've been looking for an IACF() procedure in R for a long time (it's a very convenient function to check for overdifferencing time series), and eventually decided to write my own function. Here's what I came up with : 3 web-pages helped me estimate it : http://www.xycoon.com/inverse_autocorrelations.htm

Hi I'm getting the following error that do not make sense to me, what am Idoing wrong ? > acf(Recsim[1,], lag.max=1) Error in acf(Recsim[1, ], lag.max = 1) : 'lag.max' must be at least 1 Regards EJ

In some tests I am getting a difference in the 7th or 8th significant digit depending on whether I set DUP=T or DUP=F in the .Fortran call below. While this is not a huge error it is bigger than I would expect and may be a symptom of another problem. Any comments or suggestions? (R1.2.2 on Solaris.) Paul Gilbert _______ genD.ARMA <- function(model, data, d=0.01, eps=1e-4, r=6, warn=F){ n

Time Series functions in R ========================== I think a good basic S-like functionality for library(ts) in base R would include ts class, tsp, is.ts, as.ts plot methods start end window frequency cycle deltat lag diff aggregate filter spectrum, spec.pgram, spec.taper, cumulative periodogram, spec.ar? ar -- at least univariate by Yule-Walker arima -- sim, filter, mle, diag, forecast

Dear all, I have two different data frames, that have two common variables: date and sample. Here is a small extract of both of them > head(traffic) datet sessiont samplet buddleiat 1 07-08-2012 1 1 1 2 07-08-2012 1 1 1 3 07-08-2012 1 1 1 4 07-08-2012 1 2 3 5 07-08-2012 1 2

I tried rsept31.zip under Windows 3.11 yesterday and it didn't work at all ... the screen went black for a second and that was it. Any glue what I've done wrong? Should it work under 3.11 or only under 95/NT? .f =-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hi I am trying this > x <- read.table("/home/kenji/1245/GDS1_2grps_.cls", header = F, skip = 2) > x <- read.table("/home/kenji/1246/MYCset.cls", header = F, skip = 2) Warning message: In read.table("/home/kenji/1246/MYCset.cls", header = F, skip = 2) : incomplete final line found by readTableHeader on '/home/kenji/1246/MYCset.cls' Here are the