similar to: tseries(arma) vs. stats(arima)

Hi, I'm trying to model an ARMA(1,[1,4]), i.e. I want only lags 1 and 4 of the Moving Average part. It's the '[1,4]' part that is giving me a problem. I've tried different arma's and arima's in different packages, namely: packages tseries, fArma, FinTS, timeSeries, TSA, Zelig, ds1, forecast For example, with package FinTS: > ( ARIMA(y, order=c(1,0,c(1,4))) )

I have weekly time series data with year, month, day, and price variables. The input data set for the weekly series takes the following form: Year month day price 1990 8 20 119.1 1990 8 27 124.5 1990 9 3 124.2 1990 9 10 125.2 1990 9 17 126.6 1990 9 24 127.2 1990 10 1 132.1 1990 10 8 133.3 1990 10 15 133.9 1990 10 22 134.5 1990 10 29 133.9 .. ... ... ... ... ... .... .... 2008 3 3 313.7 2008

Is anyone aware of an R procedure similar to STATA's "heckprob" procedure? "Heckprob" fits maximum likelihood probit models correcting for sample selection bias. Thanks, Richard Saba Department of Economics Auburn University Email: sabaric@auburn.edu [[alternative HTML version deleted]]

On Tue, Oct 13, 2009 at 5:06 PM, Rolf Turner <r.turner@auckland.ac.nz>wrote: > > Not clear to me what the OP really wants. Perhaps the seasonal > model is what's required; perhaps an arima(1,0,4) model with > theta_2 and theta_3 constrained to be 0. The latter can be > achieved with > > arima(x,order=c(1,0,4),fixed=c(NA,NA,0,0,NA,NA)) > > Or perhaps

Can someone direct me to a resource or resources that list the formulae used by R functions (i.e. predict.lm ) to calculate the statistic reported. I am not a programmer and studying the r code is extremely slow going. I have searched r-project.org and all the function help files without success. For example I have attempted to replicate by hand the se.fit calculation from a lm object

I am relatively new to R and object oriented programming. I have relied on SAS for most of my data analysis. I teach an introductory undergraduate forecasting course using the Diebold text and I am considering using R in addition to SAS and Eviews in the course. I work primarily with univariate or multivariate time series data. I am having a great deal of difficulty understanding and working with

I'd like to play with ARIMA models of stock prices, but I am a complete novice. Could some kind soul explain the relationship among packages "ts", "tseries", "dse", "dse2", and "fracdiff"? Are they 'competing' products or does one depend on another? Where would be the best place for a novice to begin? Thanks for any advice. PS. I

I have an ascii data set of monthly observation starting in Jan 1946 with a header. hstarts 57 65 95 103 103 97 94 . . . Which I read with the following code tab6.1<-ts(read.table(fname, header=TRUE),frequency=12,start=c(1946,1)) I would like to run a time series model with dummy variables for each month. If I had a variable which take values from 1 to 12 indicating the month I could use

Hello! I installed the fSeries package to get some farima time-series which i tried with farimaSim, but unfortunately i got always an error. I tried it this way: > farimaSim(n = 1000, model = list(ar = 0.5, d = 0.3, ma = 0.1), method="freq") Error in farimaSim(n = 1000, model = list(ar = 0.5, d = 0.3, ma = 0.1), : ... used in an incorrect context Some ideas? Regards, ___

Hi, I am trying to understand exactly what xreg does in arima. The documentation for xreg says:"xreg Optionally, a vector or matrix of external regressors, which must have the same number of rows as x." What does this mean with regard to the action of xreg in arima? Apparently somehow xreg made the following two arima fit equivalent in R: arima(x, order=c(1,1,1), xreg=1:length(x)) is

I would like to read the following hierarchical data set. There is a family record followed by one or more personal records. If col. 7 is "1" it is a family record. If it is "2" it is a personal record. The family record is formatted as follows: col. 1-5 family id col. 7 "1" col. 9 dwelling type code The personal record is formatted as follows: col.

There seems to be an error in the summary() function when applied to "ts" class objects. The results of a call to summary( ), on the R "ts" data set USAccDeaths , reports the wrong value for Max. The value reported by the summary function is 11320. The max( ) function returns the correct value 11317, the July 1993 value. Coercing the data to a data.frame and calling summary

Hello again! Is there a way to include a mean in the garch function in the library(tseries), please? I tried include.mean=T in the function statement but it didn't work thanks in advance! R Version 2.2.1 Windows Sincerely, Erin mailto: hodgess at gator.uhd.edu

Hello,Could somebody tell me what is the difference between  theese 3 calls of functionsarma(x,order=c(1,0)), arima(x,order=c(1,0,0)) ar(x,order=1)?I expected same residuals of theese three models,but unexpectably for the first two R requiredinitial value of something (what?)...Thanks in advance!  [[alternative HTML version deleted]]

Hi there, Does anyone know how to predict ARMA? It doesn?t have either predict or forecast methods. I found couple of packages called fbasic and fseries at http://www.itp.phys.ethz.ch/econophysics/R/, which has ?arma.predict? in it, but it doesn?t seem to be working. Any help in this regard would be appreciated. Thanks in advance. Regards Skanda Kallur "Prediction is very difficult,

R-guRus , ARMA function in tseries, seems to be calculating the AR coeff 's as coef <- lm(xx[,1]~xx[,lag$ar+1])$coef [*snipped* from around line 77,] I'd like to modify this model with another term somewhat in these lines lm(xx[,1] ~xx[,lag$ar+1]+mvgsignal)$coef where mvgsignal is a moving average signal based on some indicators, the question is could i simply hack into

Dear all R user, Please forgive me if my problem is too simple. Actually my problem is basically Statistical rather directly R related. Suppose I have return series ret with mean zero. And I want to fit a Garch(1,1) on this. my is r[t] = h[i]*z[t] h[t] = w + alpha*r[t-1]^2 + beta*h[t-1] I want to estimate the three parameters here; the R syntax is as follows: #

Hello, I am having some problems with fitting an ARMA model to my time series data (randomly generated numbers). The thing is I have tried many packages [tseries, fseries, FitARMA etc.] and all of them giving very different results. I would appreciate if someone could post here what the best package is for my purpose. Also, after having done the fitting, I would like to check for the model's

Hi, I'm running R 2.4.1 on Fedora Core 6 and am unable to install the tseries package. I've resolved a few problems getting to this point, by running a yum update, installing the gcc-gfortran dependency, but now I'm stuck. Could someone please point me in the right direction? ========R install.packages output ======= ==================================

In the summary of the output of arma, there's a number Pr(>|t|), however, I don't know what is its meaning - at least, it doesn't _seem_ to be a Student's t distribution. Reproducible test case: x <- c(0.5, sin(1:9)) reg <- arma(x, c(1,0)) summary(reg) <output> Call: arma(x = x, order = c(1, 0)) Model: ARMA(1,0) Residuals: Min 1Q Median 3Q