similar to: TIme Series AR to MA and (viceversa)

I would like to compute the next 15 observations for an ARFIMA(2,1,0) model along with confidence intervals. Can someone provide code? Many thanks. Jill ____________________________________________________________________________________ [[elided Yahoo spam]]

I am trying to read an ftable using read.ftable, but I get the following error message: > jobSatTable <- read.ftable("http://definetti.uark.edu/~gpetris/stat5333/jobSatisfaction.dat",skip=2) Error in seek(file, where = 0) : no applicable method for "seek" In addition: Warning messages: 1: no non-missing arguments to max; returning -Inf 2: no non-missing arguments to

Can this be re-implemented to run faster (without the loop) ? r <- list() n = nrow(prices) for (i in (w+1):n) { window <- prices[(i-w):(i-1),] if (prices[i,]$settle > max(window$high)) r <- append(r, 1) else if (prices[i,]$settle < min(window$low)) r <- append(r, -1)

Hello R Users! I have a question about the output of ARMAtoMA when used to calculate the variance of a model. I have a mixed model of the form ARMA(1,1). The actual model takes the form: X(t) = 0.75X(t-12) + a(t) - 0.4a(t-1) Given that gamma(0) takes the form [(1 + theta^2 - 2*theta*phi)/(1-phi^2)]*sigma(a), I would expect a process variance of 4.02*sigma(a) when I substitute 0.75 for phi and

Hi All, I have a problem trying to forecast using the dlm package, can anyone offer any advise? I setup my problem as follows, (following the manual as much as possible) data for example to run code CostUSD <- c(27.24031,32.97051, 38.72474, 22.78394, 28.58938, 49.85973, 42.93949, 35.92468) library(dlm) buildFun <- function(x) { dlmModPoly(1, dV = exp(x[1]), dW = exp(x[2])) } fit <-

Dear All, I have a database of 200 observations named myD. In the dataframe there are a column named code (with codes varying from 1 to 77), a column named "prevalence" with some quantitative measurements are given and an column named Pr_mean, with no values. I would like to set a cycle to compute the average of prevalence values for each different code and store the averages under the

Hi - I read through dse package manual a bit. I'm not quite certain how I can use it to estimate a time varying coefficient regression model? I might pick up an inappropriate package. Any suggestion would be greatly appreciated. Thank you. rh

How can I get the results of lm() into a list so I can loop through the results? e.g. myResults[1] <- lm(...) myResults[2] <- lm(...) myResults[3] <- lm(...) ... myResults[15] <- lm(...) myResults[16] <- lm(...) so far every attempt I've tried doesn't work throwing a "number of items to replace is not a multiple of replacement length" error or simply not

Hi all, ok, can someone tell me which header fields that the "Sender" function in sieve looks at ??? I've created a sieve rule to move a mail to a directory that looks like this : # rule:[File Message] if anyof (address :contains "Sender" "Dude") { fileinto "TestFolder"; } And the emails headers look like this : Received: (qmail 5811 invoked

Hi, I am studying Kalman Filter and it seems to be difficult for me to apply the filter on a simple ARMA. It is easy to construct the state-space model, for instance: dlmModARMA(ar=c(0.4,-0.2),ma=c(0.2,-0.1, sigma2=1) but applying the dlmFilter on it, it doesn't work... I don't know if my problem is clear but if anyone has already worked on Kalman filter, it could be great to advise me!

Dear R community, You may have the solution to how to construct a function using the function "source()" to build the function; i.e. myfunction <- function(...){ source('file1.r') source('file2.r') } After compiling and installing the "myfunction" in R, then calling the "myfunction" gives an error because the content of 'file1.r' and

On 5/13/2008 10:27 AM, Yukihiro Ishii wrote: > Hi Rusers! > > I am ashed of asking such a simple question. > > X<-matrix(rnorm(24), 4) > X0<-apply(X,2,mean) > > What I want is a matrix which consists of colums such as X[,1]--X0[1]. > > X-X0 doesn't work. > > Perhaps apply function? scale(X, scale=FALSE) ?scale > Thanks in advance. > >

Hi All, I am new to R, having used S-plus a number of years back. I would like to set up a time series object for forecasting, with data collected daily between 5th April 2001 and 16th September 2001. Any help would be most appreciated as I have been unable to find any suitable examples in the documentation. Thanks, Mark.

Hello to everyone! I have a question for you..I need to predict multivariate time series, for example sales of 2 products related one to the other, having the 2 prices like inputs.. Is there in R a function to do it? I saw dse package but I didn't find what a I'm looking for.. Could anyone help me? Thank you very much Giusy -- View this message in context:

I am using tapply to compute means by group, as in > tapply(y, z, mean) 1 2 3 21.00757 20.50031 NA where y is a vector and z a factor with levels 1:3. I would like to get a zero instead of the NA I get when a subset is empty. The problem is that FUN is not applied to empty cells of the ragged array. (Of course I am willing to define a function my.mean that returns

Hello, Recently updated Ubuntu to 10.04, I have installed the r-base and r-base-dev packages from CRAN with sudo apt-get install r-base But now when I want to start R I get the following error: gpetris at definetti:~$ R /usr/lib/R/bin/exec/R: symbol lookup error: /usr/local/lib/libreadline.so.6: undefined symbol: PC Any clue??? Has anybody seen that before? I have tried to reinstall the

Is there a way of `restarting' split.screen? This is what I am getting: > close.screen() [1] 10 11 12 13 > close.screen(all=TRUE) Error in par(args) : parameter "i" in "mfg" is out of range > graphics.off() > x11() > close.screen() [1] 10 11 12 13 > close.screen(all=TRUE) Error in par(args) : parameter "i" in "mfg" is out of range As

Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel. An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on

Hello! I am trying to put two math expressions in the title of a plot. As you can see below, I can place correctly one expression at a time, but not both. Ideally I would like to have them separated by a comma. Any suggestions? > k <- 1 > n.eff <- c(20, 30) > ### this works > plot(0,0, main = substitute(n == k, list(k = k))) > ### this works > plot(0,0, main =

We have tried Asterisk 1.0.9 on FC4 and have never been able to get CAPI (with Fritz card, fcpci) to work properly. Apart from that Asterisk works fine in switching internal calls. But it's useless if we can't make outgoing calls on our ISDN line. We are considering abandoning FC4 for Debian or SuSe. What is the general concensus on the best Linux to run Asterisk with CAPI? /Why Tea