similar to: Lag() and lag()

Is it possible to fit a lagged regression, "y[t]=b0+b1*x[t-1]+e", using the function "lag"? If so, how? If not, of what use is the function "lag"? I get the same answer from y~x as y~lag(x), whether using lm or arima. I found it using y~c(NA, x[-length(x)])). Consider the following: > set.seed(1) > x <- rep(c(rep(0, 4), 9), len=9) > y <-

Dear R users, I have an unbalanced panel with an average of I=100 individuals and a total of T=1370 time intervals, i.e. T>>I. So far, I have been using the plm package. I wish to estimate a FE model like: res<-plm(x~c+v, data=pdata_frame, effect="twoways", model="within", na.action=na.omit) ?where c varies over i and t, and v represents an exogenous impact on x

Dear R-help-archive.. I am trying to figure out how to make arima prediction when I have a process involving multivariate time series input, and one output time series (output is to be predicted) .. (thus strictly speaking its an ARMAX process). I know that the arima function of R was not designed to handle multivariate analysis (there is dse but it doesnt handle arma multivariate analysis, only

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))) )

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

Hello, The "arma" function in the "tseries" package allows estimation of models with specific "ar" and "ma" lags with its "lag" argument. For example: y[t] = a[0] + a[1]y[t-3] +b[1]e[t-2] + e[t] can be estimated with the following specification : arma(y, lag=list(ar=3,ma=2)). Is this possible with the "arima" function in the

As I am reading ?arima, only NA entries in the argument fixed= imports. The following seems to indicate otherwise: x <- arima.sim(model=list(ar=0.8), n=100) + (1:100)/50 > t <- 1:100 > mod1 <- lm(x ~ t) > > init1 <- c(0, coef(mod1)[2]) > fixed1 <- c(as.numeric(NA), 0) > > arima(x, order=c(1,0,0), xreg=t, include.mean=FALSE, init=init1, fixed=fixed1)

When I run predict.Arima in my code, I get warnings like: Warning message: In cbind(intercept = rep(1, n), xreg) : number of rows of result is not a multiple of vector length (arg 1) I think this is because I'm not running predict.Arima in the same environment that I did the fit, so the data object used in the fit is no longer present. Looking at the predict.Arima source,

I am using arima to develop a time series regression model, I am using arima b/c I have autocorrelated errors. Several of my independent variables are categorical and I have coded them as factors . When I run ARIMA I don't get any warning or error message, but I do not seem to get estimates for all the levels of the factor. Can/how does ARIMA handle factors in xreg? here is some example

i am using sarima() function as below ___________________________________________________________________________________________ sarima=function(data,p,d,q,P=0,D=0,Q=0,S=-1,tol=.001){ n=length(data) constant=1:n xmean=matrix(1,n,1) if (d>0 & D>0) fitit=arima(data, order=c(p,d,q), seasonal=list(order=c(P,D,Q), period=S),

Hi, how to estimate a the following model in R: y(t)=beta0+beta1*x1(t)+beta2*x2(t)+...+beta5*x5(t)+beta6*y(t-1)+beta7*y(t-2)+beta8*y(t-3) 1) using "lm" : dates &lt;- as.Date(data.df[,1]) selection&lt;-which(dates&gt;=as.Date("1986-1-1") &amp; dates&lt;=as.Date("2007-12-31")) dep &lt;- ts(data.df[selection,c("dep")]) indep.ret1

Hello I want to fit an AR model were two of the coefficients are fixed to zero (the second and third ar-coefficients). I used the "arima" function with the "fixed" argument but the ar3 coefficient is not set to zero: ============================================== > arima(Y, order=c(4,0,0), xreg=1:23, fixed=c(NA,0,0,NA,NA,NA)) Call: arima(x = Y, order = c(4, 0, 0), xreg =

Hello, Does anybody know if the function arima with an external regressor (xreg) applies the auto correlation on the dependant variable or on the residuals. In comparison with SAS (proc autoreg), it seems that the auto correlation applies on the residuals but i'd like to have the confirmation. I want to estimate: Y[t] = a[1]*X[t] + a[2] + E[t] with E[t]=b[1]*E[t-1] Should I use : arima(Y,

Hello, Does anyone know how the parameter estimates are calculated for xreg variables when called as part of an arima() command, or know of any literature that provides this info? In particular, I was wondering if there is a quick way to compare different combinations of "xreg" variables in the arima() fit in the same way that you would in multiple regression (using AIC & R^2

Two time series questions: FITTING TRANSFER FUNCTIONS WITH LAGS: Consider the following toy example: > dates <- paste(11:21, "/01/2005", sep="") > Dates <- as.Date(dates, "%d/%m/%Y") > set.seed(1) > DF <- data.frame(date=Dates, y=rnorm(11), x=rnorm(11, 3)) > arima(DF$y, c(1,0,0), xreg=lag(DF$x, 1)) ar1 intercept lag(DF$x,

## The following corrects what I believe to have been an error in my previous post: #### Two time series questions: I. FITTING TRANSFER FUNCTIONS WITH LAGS: Consider the following toy example: > dates <- paste(11:21, "/01/2005", sep="") > Dates <- as.Date(dates, "%d/%m/%Y") > set.seed(1) > DF <- data.frame(date=Dates, y=rnorm(11),

When I do ARMA(2,2) using one lag of LCPIH data This is eview result > > *Dependent Variable: DLCPIH > **Method: Least Squares > **Date: 08/12/11 Time: 12:44 > **Sample (adjusted): 1970Q2 2010Q2 > **Included observations: 161 after adjustments > **Convergence achieved after 14 iterations > **MA Backcast: 1969Q4 1970Q1 > ** > **Variable Coefficient Std.

On 21 May 2015, at 12:49 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: >>>>>> peter dalgaard <pdalgd at gmail.com> >>>>>> on Thu, 21 May 2015 11:03:05 +0200 writes: > >> On 21 May 2015, at 10:35 , Martin Maechler <maechler at lynne.stat.math.ethz.ch> wrote: > >>>> >>>> I noticed that

I would like the text plotted with 'mtext' to be alighned like it is for printing on the console. Here is what I have: > print(emt) ME RMSE MAE MPE MAPE MASE original -1.034568e+07 1.097695e+08 2.433160e+07 -31.30554 37.47713 1.5100050 xreg 1.561235e+01 2.008599e+03 9.089473e+02 267.05490 280.66734

Hello all, I''m having fun again with the arima function. This time I read in: http://www.stat.pitt.edu/stoffer/tsa2/R_time_series_quick_fix.htm <<It has recently been suggested (by a reliable source) that using xreg in arima() does NOT fit an ARMAX model [insert slap head icon here]. This will be investigated as soon as time permits.>> (by R.H. Shumway & D.S. Stoffer)