similar to: weighted symmetric estimator

Dear R users, I am using the vars package to calculate the impulse response functions and the forecast error variance decomposition of a VAR model. Unfortunately I do not know whether these functions assume unit or one standard deviation shocks. I tried to look into the code of these functions, but in vain: neither irf, nor vars::irf, nor vars:::irf output the code of the functions. Does someone

Dear R users, I need perform structural analysis on a no intercept VAR model. Unfortunately the functions irf.VAR and dfev that come with the MSBVAR package only work with objects output by the reduced.form.var function, which seems to only evaluate VAR models with intercept. Is there a way to suppress the estimation of intercept term in reduced.form.var? Do I need to modify the code, and if I

Dear R users, I need to to the the following. Let a= 1 2 3 4 5 6 and b= -1 -2 -3 be (2x3) matrices. -4 -5 -6 I need to combine the two matrices into a new (2x6) matrix like this: ab = ( 1 -1 2 -2 3 -3 ) 4 -4 5 -5 6 -6 How can this be done in R? ----------------------------------------------------------------- ?????? ???

Dear R users, I need perform structural analysis on a no intercept VAR model. Unfortunately the functions irf.VAR and dfev that come with the MSBVAR package only work with objects output by the reduced.form.var function, which seems to only evaluate VAR models with intercept. Is there a way to suppress the estimation of intercept term in reduced.form.var? Do I need to modify the code, and if I

Hello! I would like to ask if there is in R a function that estimates the spectral density function of a stochastic series at frequency zero by the "plug-in method", advocated by Andrews in his paper "Heteroscedasticity and Autocorrelation Consistent Covariance Matrix Estimation", Econometrica, 59,817-858. I saw R has functions that employ Andrews' plug-in method using an

Dear nut users, I am running nut-2.6.0 on Slackware-Linux-13.37.0 with kernel 2.6.37.6. I am trying to get the software work for a repotec UPS with model name: RPT-1003AU. The UPS communicates with the computer via USB. I know that the model is not officially supported but I want to try out whether it will work with some of the repotec drivers. Here is the result with the genericups upstype=13

Thank you for your response. The point is that I have configured my system to shut down at the event "upsc Inform@localhost ups.status"!="OL", not at "LB". So there will be no /etc/killpower file created. At an ups.status != "OL" event, the following script named "notifyme" is run: echo `date` >> /root/upsLog echo " No power, shutting

R has a Dickey-Fuller Test implementation (adf.test) that tests for unit roots in an autoregressive process with a constant and linear trend. Is there a DF implementation that doesn't use the constant or trend? Thanks, Jake. -- View this message in context: http://www.nabble.com/Dickey-Fuller-Tests-with-no-constant-and-no-trend-tp23565210p23565210.html Sent from the R help mailing list

Hello, How can an interative unit root test be implemented in R? More specifically, given a time series, I wish to perform the Dickey Fuller Test on a daily basis for say the last 100 observations. It would be interative in the sense that this test would be repeated each day for the last 100 observations. Given the daily Dickey Fuller estimates of delta for the autoregressive process d(Y(t))

Hello! I need to perform weighted multidimensional scaling analysis(WMDS). I did rsitesearch, googled, but I could find no info on how to perform WMDS using R. In several places they say it is possible with the ALSCAL algorithm, but I could not find the relevant function to carry it out. ----------------------------------------------------------------- ???????? ?? ??????? ?? ??? ???????????!

Hello, I have a survivor curve that shows account cancellations during the past 3.5 months.  Fortunately for our business, but unfortunately for my analysis, the survivor curve doesn't yet pass through 50%.  Is there a safe way to extend the survivor curve and estimate at what time we'll hit 50%? We started a new program 3.5 months ago, and I believe that this set of accounts behaves

I have a question about constructing the likelihood function where there is censoring at level 1 in a two-level random effects sum. In a conventional solution, the likelihood function is constructed using the density for failures and the survivor function for (in this case, right) censored results. Within (for example) an R environment, this is easy to do and gives the same solution as survreg

Dear David, Dear R Users , Suppose that we want to regress for example a certain autoregressive model using SVM. We have our data and also some fixed kernels in libSVM behinde e1071 in front. The question: Where can we insert our certain autoregressive model ? During creating data frame ? Or perhaps we can make a relationship between our variables ended to desired autoregressive model ?

Hello! After I type "upsmon", in syslog the following message appears: Aug 16 14:11:09 marto upsmon[7731]: Master privileges unavailable on UPS [Inform@localhost] Aug 16 14:11:09 marto upsmon[7731]: Reason: Access denied Can you tell me where my error is, my configuration files are: ups.conf: [Inform] driver = powercom port = /dev/ttyS0 desc = "inform guard" linevoltage = 220

Hi everyone I've got a problem regarding the arima() and the ar() function for autoregressive series. I would simply like to combine them. To better understand my question, I first show you how I'm using these two functions individually (see file in the attachement). 1) apply(TSX,2, function(x) ar(na.omit(x),method="mle")$order # this function finds the optimal

There has been some recent discussion on this list about the value of using intervals with lme() to check for whether a model is ill-defined. My question is, what else can drive very large confidence intervals for the variance components (or cause the error message "Error in intervals.lme(Object) : Cannot get confidence intervals on var-cov components: Non-positive definite approximate

Hello everybody, I've a question about "autoregressive Regressionmodels". Let Y[1],.....,Y[n], be a time series. Given the model: Y[t] = phi[1]*Y[t-1] + phi[2]*Y[t-1] + ... + phi[p]*Y[t-p] + x_t^T*beta + u_t, where x_t=(x[1t],x[2t],....x[mt]) and beta=(beta[1],...,beta[m]) and u_t~(0,1) I want to estimate the coefficients phi and beta. Are in R any functions or packages for

Hi, I am trying to replicate Engle and Manganelli's paper Conditional Autoregressive Value at Risk (CAViaR) by Regression Quantiles. I have the Matlab code which I cannot get to work as I have never used Matlab before, does anyone know if there is the same code available to estimate the CAViaR models in R? Thanks, Shane -- View this message in context:

Hello, I am using R for fitting parameters of a time series model. The model is as below. Y(t) = mu + a*X(t) + YN(t) where YN(t) = b*YN(t-1) + innovation and Z(t) follows N(0,1). The main obstacle for me is the autoregressive error term, YN(t). I can't figure out how to estimate the parameters (mu, a, b) with usual 'arima' function in R. What I have tried is.... 1. Do the

Dear R users, I have a lot of experience with traditional R graphics, but I decided to turn to trellis as it was recommended for spatial graphs by the sp package. In traditional R graphics I always first set the size of the device region absolute units (e.g. mm) and then I firmly fix the inner margins with mai and the outer margins with oma also in absolute units. What is left from the device