similar to: Bootstrap Multivariate Times Series Forecast

Hi R-helpers, I'm new to time series modelling, but my requirement seems to fall just outside the capabilities of the arima function in R. I'd like to fit an ARMA model where the variance of the disturbances is a function of some exogenous variable. So something like: Y_t = a_0 + a_1 * Y_(t-1) +...+ a_p * Y_(t-p) + b_1 * e_(t-1) +...+ b_q * e_(t-q) + e_t, where e_t ~ N(0, sigma^2_t),

Dear useRs, I am trying to model a time series with a transfer function. I think it can be put into the ARMA framework, and estimated with the 'arima' function (and others have made similar comments on this list). I have tried to do that, but the results have so far been disappointing. Maybe I am trying to make 'arima' do something it can't... The data are time series of

Hi, I'm trying to plot the forecasts I generated using the Plot.Fore function of the BootPR package. But I got an error from zoo: My data: Time Series: Start = 1 End = 18 Frequency = 1 [1] 38731 38628 39117 92809 71984 31226 58613 72360 107956 92066 [11] 95208 99098 95848 120383 110717 105680 98469 101916 Script: y1<-ts(y1);

This posting has nothing to do with R (except maybe that I am using R very heavily in writing the paper to which the question pertains.) I simply wish to draw upon the impressive knowledge and wisdom of the R community. Since this question is way off topic, if anybody has the urge to reply, they should probably email me directly: rolf at math.unb.ca rather than via this list. My question

Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem is trying to figure out how the (model) identifier (e.g, EII, VII, VVI, etc.) relates to the covariance matrix. The parameterization of the covariance matrix makes use of the method of decomposition in Banfield and Rraftery (1993) and Fraley and Raftery (2002) where

I'm using the arima function to estimate coefficients and also using predict.Arima to forecast. This works nicely and I can see that the results are the same as using SAS's proc arima. I can also take the coefficent estimates for a simple model like ARIMA(2,1,0) and manually compute the forecast. The results agree to 5 or 6 decimal places. I can do this for models with and without

Dear All, I have six different campaign and 5 different agent have login on that campaign.*Same thing i have done using agi and database,i never use queue management on this scenario. Agent** can also shuffling one campaign to anther campaign. * Now i want to do some work with queue.I want to use single queue to managing this. Eg: campaign Agent Login A a_1,a_3

Hi R-helpers I would like to plot some planes which are perpendicular to the x-y plane, such as x=y. Is there a way to do this in wireframe? I realize that I am not plotting a function of x, y since there are infinite number of z's that satisfy the above relation.... Hmm... Somewhat related, I would also like to plot a line in 3 d space... Finally, if you are feeling really brave, I

Dear All, I would like to generate m positive real numbers c_i, I=1,...,m, such that (1) c_1 + c_2 + ... + c_m=m, (1) after being ordered into c_1 >= c_2 >= .... >=c_m>0, we have that c_m is of the same order of m^(-1/8), when m is sufficiently large. Thanks, -Chee [[alternative HTML version deleted]]

Dear all, I've only just joined this list, but I can't find any mention of this idea anywhere else, so I thought I'd just post here before getting too deep into programming and possibly reinventing the wheel. Here at Aber, we have around 30 unix and linux servers doing core services. Each one is maintaining its own logfiles and, for various reasons, we want to keep these on the

Hi all, I have n independent variables A_1, A_2, A_3,......,A_n, and each with known variances var(A_1), var(A_2),..., but unknown mean. How can I get the approximation of the variance of the product of the variables using numerical computation, i.e. var(A_1*A_2*A_3*.....*A_n)? Thanks. Sincerely, Yanwei Zhang Department of Actuarial Research and Modeling Munich Re America Tel: 609-275-2176

Dear All, I've been searching relevant topics about replacing values, none seemed to be applicable to me... I have a file with many many varieties, and want to replace some of them into different names. I tried various of ways, still don't know how to do that most efficiently.. Here is part of the example data: Gen Rep A_1 1 A_1 2 A_2 1 A_2 2 B_1 1 B_1

Hello, Apologies if this question has already arised, hope you can help me to the find the solution to this or point the place to look at. I have a multidimensional piece-wise regression linear problem, i.e. to find not only the regression coefficients for each "interval" but also the beginning and ends of the intervals. To simplify it to the one dimensional case and two intervals,

I've posted this before but have not been able to locate what I'm doing wrong. I cannot determine how the forecast is made using the estimated coefficients from a simple AR(2) model when there is an exogenous variable. Does anyone know what the problem is? The help file for arima doesn't show the model with any exogenous variables. I haven't been able to locate any documents

Hi! I want to write portions of my data (3573 columns at a time) to twenty folders I have available titled "A_1" to "A_20" such that the first 3573 columns will go to folder A_1, next 3573 to folder A_2 and so on. This code below ensures that the data is written into all 20 folders, but only the last iteration of the loop (last 3573 columns) is being written into ALL of the

Hi, > Both A and B are undef: > LHS = (undef & undef) | (undef & undef) = undef // Since ~undef = undef > RHS = undef > Thus transform is correct. LLVM documentation (http://llvm.org/docs/LangRef.html#undefined-values) suggests that it is unsafe to consider (a & undef = undef) and (a | undef = undef). "As such, it is unsafe to optimize or assume

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. The details of the model are: Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~ Beta(mu_{ij} * phi_{ij}; (1 - mu_{ij}) * phi_{ij}) y_{ij} = mu_{ij} + w_{ij} * *with* *logit(mu_{ij}) = Beta_{0i} + Beta_{1i} * x1_{ij} + b2 * x2_{ij}

*Dear R-help: My name is Rodrigo and I have a question with nlme package in R to fit a mixed beta regression model. I'm so sorry. In the last email, I forgot to say that W is also a unknown parameter in the mixed beta regression model. In any case, here I send you the correct formulation. ** Suppose that:* *j in {1, ..., J}* *(level 1)* *i in {1, ..., n_j}* *(level 2)* *y_{ij} ~

Dear R experts, I was wondering if there are any R functions that give the tail area of a sum of chisquare distributions of the type: a_1 X_1 + a_2 X_2 where a_1 and a_2 are constants and X_1 and X_2 are independent chi-square variables with different degrees of freedom. Thanks, Klaus -- "Feel free" - 5 GB Mailbox, 50 FreeSMS/Monat ...

Hi List l am interested in developing price model. I have found a research paper related to price model of corn in US market where it has taken demand & supply forces into consideration. Following are the equation: Supply equation: St= a0+a1Pt-1+a2Rt-1+a3St-1+a5D1+a6D2+a7D3+U1 -(1) Where D1,D2,D3=Quarterly Dummy Variables(Since quarterly data are considered) Here, Supply