similar to: warnings in ARMA with other regressor variables

Hi, I want to fit the following model to my data: Y_t= a+bY_(t-1)+cY_(t-2) + Z_t +Z_(t-1) + Z_(t-2) + X_t + M_t i.e. it is an ARMA(2,2) with some additional regressors X and M. [Z_t's are the white noise variables] How do I find the estimates of the coefficients in R? And also I would like to know what technique R employs to find the estimates? Any help is appreciated. Thanks,

Hello all, I tried to extract the p-value from adf.test in tseries; however, I got the error message such as > ht=adf.test(list.var$aa) > ht$p-value Error in ht$p - value : non-numeric argument to binary operator > ht Augmented Dickey-Fuller Test data: list.var$aa Dickey-Fuller = -2.3147, Lag order = 4, p-value = 0.4461 alternative hypothesis: stationary > ht$data [1]

Dear all, I have just sent an enquiry but probably I hadn’t expressed myself properly. Could anyone help me with the following? When I run the code on my data I get a boxplot with outliers identified by numbers 200 & 201. However, what I would like is to label these outliers with their corresponding “DATA$num” values of the data frame. In this example, the outliers should be labelled as:

Hello everybody, I have got a data set with about 400 companies. Each company has a score for its enviroment comportment between 0 and 100. These companies belong to about 15 different countries. I have e.g. 70 companies from UK and 5 from Luxembourg,- so the data set is pretty unbalanced and I want to do an ANOVA. Somthing like aov(enviromentscore~country). But the aov function is just for

Hi! I've been looking for a way to calculate std errors of marginal effects when I use interaction terms, but with no success. I pretty much have two cases: continuous variable * continuous variable, and continuous variable * binary variable. In both cases, I know how to calculate the marginal effects, even with simulation. But I still can't figure out how to calculate the std errors of

Hi all, Suppose I am fitting an arma(p,q) model to a time series y_t. So, my model should contain (q+1) white noise variables. As far as I know, each of them should have the same variance. How do I get the estimate of this variance by running the arma(y) function (or is there any other way)? Appreciate your help. Thanks, Preetam -- Preetam Pal (+91)-9432212774 M-Stat 2nd Year,

Hello all! I have a problem with a double exponential equation. This are my data's> structure(list(proc = c(1870.52067384719, 766.789388745793, 358.701545859122, 237.113777545511, 43.2726259059654, 148.985133316262, 92.6242882655781, 88.4521557193262, 56.6404686159112, 27.0374477259404, 34.3347291080268, 18.3226992991316, 15.2196612445747, 5.31600719692165, 16.7015717397302,

I have a data set to be analyzed using to binary logistic regression. The data set is iin grouped form. My question is: how I can compute Hosmer-Lemeshow test and measures like sensitivity and specificity? Any suggestion will be greatly appreciated. Thank you Endy [[alternative HTML version deleted]]

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,

Hi all, I have 4 matrices, each having 5 columns and 4 rows .....denoted by B1,B2,B3,B4. I have generated a vector of 7 indices, say (1,2,4,3,2,3,1} which refers to the index of the matrices to be chosen and then appended one on the top of the next: like, in this case, I wish to have the following mega matrix: B1over B2 over B4 over B3 over B2 over B3 over B1. 1> How can I achieve this?

Hello. I've got a problem with arma/xreg. I would like to get a better model-fit by implenting some external explanatory variable, so I thought I can implement it by expand the arima-function with an xreg-argument: I have two stationary data vectors y and x of length 201: y <-

Hi all, In my data analysis, I have created a random matrix M ( of order 500 X 7). I want to use the same matrix when I start a new session, or suppose I want to send this matrix to one of my friends (because this matrix is randomly generated, and I dont want to use any other 500X7 matrix randomly generated by R). How can I save and call this matrix in the later sessions as well? Appreciate

Hi all, I have a data set containing variables LOSS, GDP, HPI and UE. (I have attached it in case it is required). Having renamed the variables as l,g,h and u, I wish to run a Lasso Regression with l as the dependent variable and all the other 3 as the independent variables. data=read.table("data.txt", header=T) l=data$LOSS h=data$HPI u=data$UE g=data$GDP matrix=data.frame(l,g,h,u)

Hi all there I am enjoying R since 2 weeks and I come to my first deadlock, il am trying to use predict.Arima in the ts package. I get a "Error in cbind(...) : cannot create a matrix from these types" -- Start R session ----------------------------------------------------- > fitdiv <- arima(data, c(2, 0, 3), xreg = y ) ; print(fitdiv) Call: arima(x = data, order = c(2, 0, 3),

Hi all there I am enjoying R since 2 weeks and I come to my first deadlock, il am trying to use predict.Arima in the ts package. I get a "Error in cbind(...) : cannot create a matrix from these types" -- Start R session ----------------------------------------------------- > fitdiv <- arima(data, c(2, 0, 3), xreg = y ) ; print(fitdiv) Call: arima(x = data, order = c(2, 0, 3),

Hi all, 1>i have 3 vectors a,b and c, each of length 25....... i want to define a new data frame z such that z[1] = (a[1] b[1] c[1]), z[2] = (a[2] b[2] c[2]) and so on...how do i do it in R 2> Then i want to draw bootstrap samples from z. Kindly suggest how i can do this in R. Thanks, Preetam -- Preetam Pal (+91)-9432212774 M-Stat 2nd Year,

Hi guys, I am running glm(y~., data = history,family=binomial)-essentially, logistic regression for credit scoring (y = 0 or 1). The dataset 'history' has 14 variables, a few examples: history <- read.csv("history.csv". header = TRUE) 1> 'income = 100,200,300 (these are numbers in my dataset; however interpretation is that these are just tags or labels,for every

Hi all, I have run a ridge regression as follows: reg=lm.ridge(final$l~final$lag1+final$lag2+final$g+final$u, lambda=seq(0,10,0.01)) Then I enter : select(reg) and it returns: modified HKB estimator is 19.3409 modified L-W estimator is 36.18617 smallest value of GCV at 10 I think it means that it is advisable to

Hi all, I was running the adf test in R. CODE 1: adf.test(data$LOSS) Augmented Dickey-Fuller Test data: data$LOSS Dickey-Fuller = -1.9864, Lag order = 2, p-value = 0.5775 alternative hypothesis: stationary CODE 2: adf.test(diff(diff(data$LOSS))) Augmented Dickey-Fuller Test data: diff(diff(data$LOSS)) Dickey-Fuller = -6.9287, Lag order = 2, p-value = 0.01 alternative

Hi all, I have run a ridge regression as follows: reg=lm.ridge(final$l~final$lag1+final$lag2+final$g+final$g+final$u, lambda=seq(0,10,0.01)) Then I enter : select(reg) and it returns: modified HKB estimator is 19.3409 modified L-W estimator is 36.18617 smallest value of GCV at 10 I think it means that it is